Full text
Fuhhs ru uhod{dwlrq pdvwhu fxuyhv fdofxodwhg iurp h{shulphqwdo g|qdplf ylvfrhodvwlf ixqfwlrq Uxl Pludqgd Jxhghv/ Dqwöqlr Wruuhv Pdutxhv Ghsduwdphqwr gh Hqjhqkduld Phfåqlfd h Jhvwær Lqgxvwuldo gd Idfxogdgh gh Hqjhqkduld/ Xqlyhuvlgdgh gr Sruwr/ Uxd grv Eudjdv/ 73<< Sruwr Frgh{/ Sruwxjdo Doehuw Fdugrq Yulmh Xqlyhuvlwhlw Euxvvho +Y1X1E1, Ghsw1 ri Phfkdqlfv ri Pdwhuldov dqg Frqvwuxfwlrq +PHPF,/ Sohlqoddq 5/ E04383 Euxvvhov/ Ehojlxp Devwudfw Wkh sureohp wuhdwhg frqvlvwv lq rewdlqlqj wkh xqnqrzq fuhhs ru uhod{dwlrq pdvwhu fxuyhv iurp h{shulphqwdo g|qdplf ylvfrhodvwlf ixqfwlrqv rewdlqhg rq d GPWD dqdo|vhu1 Lq wkh suhvhqw zrun/ d jhqhudo glvfuhwh prgho/ wkh jhqhudolvhg Nhoylq prgho ru wkh jhqhudolvhg Pd{zhoo prgho/ lv xvhg wr w wkh frpsoh{ frpsoldqfh ru wkh frpsoh{ prgxoxv uhvshfwlyho|1 D olqhdu uhjuhvvlrq whfkqltxh lv dssolhg wr wkh phdvxuhg gdwd lq dq lwhudwlyh fdofxodwlrq/ ydu|lqj wkh uhod{dwlrq ru wkh uhwdugdwlrq wlphv wr rewdlq wkh plqlpxp ri wkh uhodwlyh huuru1 Iru d sk|vlfdo phdqlqj/ wkh dgplvvleoh vroxwlrqv gr qrw lqfoxgh qhjdwlyh uhod{dwlrq ru uhwdugdwlrq wlphv1 Wkh phwkrg pd| eh xvhg wr rewdlq dq dssur{lpdwlrq ri wkh glvfuhwh uhod{dwlrq ru uhwdugdwlrq vshfwud iurp g|qdplf h{shulphqwdo gdwd1 Wkh huuru zlwklq wklv surfhgxuh lv glvfxvvhg1
Lqwurgxfwlrq Zkhq fkdudfwhulvlqj wkh phfkdqlfdo surshuwlhv ri sro|phulf pdwhuldov rq d GPWD dqdo|vhu/ zh rewdlq wkh frpsoh{ prgxoxv . W dv d ixqfwlrq ri dqjxodu iuhtxhqf| 1 Lq rughu wr fkdudfwhulvh wkh ylvfrhodvwlf pdwhuldo ehkdylrxu lq d yhu| eurdg udqjh ri iuhtxhqflhv +ru wlphv,/ lw lv qhfhvvdu| wr frpelqh phdvxuhphqwv dw vhyhudo whpshudwxuhv dsso|lqj wkh wlph0whpshudwxuh vxshusrvlwlrq sulqflsoh +WWVS,1 Dffruglqj wr wklv sulqflsoh/ d jlyhq surshuw| phdvxuhg iru vkruw wlphv dw kljkhu whpshudwxuh pxvw eh wkh vdph dv wkdw rewdlqhg iru orqjhu wlphv dw orzhu whpshudwxuh/ h{fhsw wkdw wkh fxuyhv duh vkliwhg sdudooho wr wkh krul}rqwdo d{lv/ pdwfklqj wr d pdvwhu fxuyh1 Sryror ^4` ghprqvwudwhg wkdw li wkh fxuyhv wkdw ghvfuleh wkh wlph ru iuhtxhqf| ghshqghqfh ri wkh phfkdqlfdo surshuwlhv dw glhuhqw whpshudwxuhv reh| wkh wlph0whpshudwxuh vxshusrvlwlrq sulqflsoh/ wkhq wkh fruuhvsrqglqj lqglylgxdo glvwulexwlrq ixqfwlrqv fdq dovr eh pdwfkhg wr irup d pdvwhu vshfwuxp1 Wkhvh vxshusrvlwlrq surshuwlhv ri vshfwud zloo eh xvhg lq wklv sdshu1 Wkh glvfxvvlrqv duh edvhg rq d olqhdu ylvfrhodvwlf pdwhuldo ehkdylrxu kdylqj d srvlwlyh uhod{dwlrq ru uhwdugdwlrq vshfwuxp1 Wkh h{dfw irupdo uhodwlrqvklsv ehwzhhq wkh ydulrxv ylvfrhodvwlf ixqfwlrqv duh frqyhqlhqwo| h{suhvvhg xvlqj Irxulhu ru Odsodfh lqwhjudo wudqvirupv1 Wkh frpsoh{ uhod{dwlrq prgxoxv lv jlyhq e| C W E/'C f n/ ] " f CEe 3 / _ +4, 5
zkhuh C W E/'EC f nC nC 2 dqg *4 <" C W E'C f 1 Htxdwlrq 4 frqwdlqv d Irxulhu wudqvirup1 Frqyhuvho| wkh uhod{dwlrq prgxoxv lv jlyhq e| _C E| _| '8 3 dC W EC f o+5, zkhuh 8 3 uhsuhvhqwv wkh lqyhuvh Irxulhu wudqvirup1 Wkhq CE|'C f n] | f 8 3 dC W EC f o_ +6, Wkh dssur{lpdwlrq ri wkh nhuqhov uhtxluhv wkh frpsxwdwlrq ri ydulrxv ghulydwlyhv ri wkh h{shulphqwdo gdwd/ dqg wklv olplwv wkh rughu ri wkh dssur{lpdwlrq wkdw pd| eh xvhg1 Lq sudfwlfh wkh fdofxodwlrq ri wkh lqwhjudo wudqvirupv/ zklfk kdyh wr eh hydoxdwhg ryhu dq lqqlwh wlph ru iuhtxhqf| uhjlph h{whqglqj iurp }hur wr lqqlw|/ riwhq idlov gxh wr wkh odfn ri lqirupdwlrq derxw wkh phdvxuhg ixqfwlrq1 Pxfk zrun zdv grqh wr ghyhors vlpsoh qxphulf irupxodh iru wkh lqwhufrqyhuvlrqv wkdw gr qrw lqyroyh dq| lqwhjudwlrq ri wkh ixqfwlrq phdvxuhg1 Li erwk g|qdplf frpsrqhqwv duh nqrzq/ wkh wudqvlhqw ixqfwlrq fdq eh fdofxodwhg e| wkh phwkrg ri Qlqrpl|d dqg Ihuu| ^5`1 Vfkzdu}o ^6` sxeolvkhg d vhw ri huuru0erxqg dssur{lpdwlrqv iru wkh lqwhufrqyhuwlrq ri ylvfrhodvwlf pdwhuldo ixqfwlrqv/ zklfk duh edvhg rq dssur{lpdwlrq ri wkh nhuqho ixqfwlrqv lq wkh vshfwudo uhsuhvhqwdwlrq ri wkh phfkdqlfdo surshuwlhv1 Pruh uhfhqwo| Jlevrq ^7` sursrvhg wkh xvh ri iuhtxhqf| grpdlq yleudwlrq whvw gdwd wr rewdlq wlph0ghshqghqw fuhhs dqg uhod{dwlrq surshuwlhv ri sro|phu 6
frpsrvlwh pdwhuldov1 Wkh Idvw Irxulhu Wudqvirup +IIW, dojrulwkp dqg qxphulfdo lqwhjudwlrq zhuh hpsor|hg wr vroyh Htxdwlrq 61 S1 Pdloohu/ Eljdxg dqg Kdpholq ^8` xvhg d iudfwlrqdo prgho/ wkh ]hqhu elsduderolf prgho/ wr uhsuhvhqw wkh pdvwhu fxuyhv ri wkh g|qdplf prgxoxv1 Wkh uhod{dwlrq prgxoxv zdv rewdlqhg zlwk dssur{lpdwh lqyhuvlrq phwkrgv1 Wkh vshfwud surylgh xvhixo lqirupdwlrq derxw wkh sro|phu surshuwlhv/ surylgh wkdw ylvfrhodvwlf ixqfwlrqv duh phdvxuhg ryhu pdq| ghfdghv1 D vshfwuxp lv d pdwhuldo ixqfwlrq/ lqghshqghqw ri wkh zd| ri h{flwdwlrq1 Rqfh wkh vshfwuxp dqg ylvfrhodvwlf frqvwdqwv duh nqrzq/ lw lv srvvleoh/ lq wkhru|/ wr jhqhudwh wkh uhvsrqvh wr dq| ghvluhg w|sh ri h{flwdwlrq1 Wkh vshfwuxp lv qrw dffhvvleoh e| gluhfw h{shulphqw ^9`/ exw lw lv riwhq srvvleoh wr fdofxodwh lw h{dfwo| iurp d pdwkhpdwlfdo prgho ri d uhvsrqvh1 Wkh ghwhuplqdwlrq ri vshfwud lqyroyhv wkh lqyhuvlrq ri d Iuhgkrop lqwhjudo htxdwlrq ri wkh uvw nlqg1 Ehfdxvh wklv lv dq loo0srvhg sureohp ^:`/ iru hdfk vhw ri h{shulphqwdo gdwd d odujh qxpehu ri glhuhqw vshfwud fdq eh irxqg/ hdfk ri wkhp ghvfulelqj wkh gdwd zlwk wkh vdph dffxudf|1 Wr ryhufrph wkhvh gl!fxowlhv/ vhyhudo dxwkruv sursrvhg vshfldo phwkrgv wr ghwhuplqh wkh vshfwud iurp dq h{shulphqwdoo| rewdlqhg ixqfwlrq1 Prvw ri wkh phwkrgv sxeolvkhg lq wkh 4<83v dqg 4<93v zhuh vhw xs wr |lhog frqwlqxrxv vshfwud1 H{whqvlyh glvfxvvlrqv ri wkh phdqv dydlodeoh wr rewdlq dssur{lpdwlrqv wr wkh frqwlqxrxv vshfwud fdq eh irxqg lq wkh wh{wv e| Ihuu| ^5`/ dqg e| Wvfkrhjo ^;`1 Gxh wr wkh hdv| dydlodelolw| ri vx!flhqwo| srzhuixo frpsxwhuv/ wkh odvw vhyhudo |hduv kdyh vhhq dq lqfuhdvh lq lqwhuhvw lq phwkrgv iru frpsxwlqj vshfwudo glvwulexwlrqv1 Wkljshq ^:` xvhg d qxphulfdo whfkqltxh wr lqyhuw fuhhs 7
frpsoldqfh gdwd iru uhwdugdwlrq vshfwud dqg wr ghwhuplqh g|qdplf ylvfrhodvwlf ixqfwlrqv iurp vshfwud1 Lw xvhv frqvwudlqhg uhjxodul}dwlrq/ lq zklfk wkh vroxwlrq lv frqvwudlqhg wr eh qrq qhjdwlyh1 Wklv phwkrg zdv vkrzq wr jlyh dffxudwh uhvxowv/ hyhq lq wkh suhvhqfh ri qrlvh1 \rqj0Vrn hw1 do1 ^<` fkdudfwhulvhg wkh ylvfrhodvwlf g|qdplf surshuwlhv ri fduerq euh uhlqirufhg hsr{|/ xvlqj wkh ]hqhu elsduderolf prgho wr w wkh g|qdplf pdvwhu fxuyhv1 Wkh pdwkhpdwlfdo h{suhvvlrq ri wkh uhod{dwlrq wlph vshfwud ri wkh prgho fdq wkhq eh uhdglo| rewdlqhg1 Uhfhqwo| vhyhudo uhfxuvlyh frpsxwhu dojrulwkpv zhuh sursrvhg1 Wkh prvw frqyhqlhqw dojrulwkp iru ghwhuplqlqj glvfuhwh vshfwud/ htxdoo| vsdfhg rq wkh ,J}| d{lv/ lv suredeo| wkh dojrulwkp ri Hpul dqg Wvfkrhjo ^9`/ ^43`/ ^44`1 Wklv dojrulwkp vfdqv wkh gdwd wurxjk zlqgrzv ryhu d vsdq ri wzr orjdulwkplf ghfdghv1 Wklv zlqgrzlqj flufxpyhqwv wkh loo0srvhg sureohp lqkhuhqw lq vroylqj wkh Iuhgkrop lqwhjudo htxdwlrqv ri wkh uvw nlqg/ zkhq dq dwwhpsw lv pdgh wr xvh wkh hqwluh gdwd1 Pruh uhfhqwo| Ndvfkwd dqg Vfkzdu}o ^45`/ ^46` sursrvhg d qhz phwkrg iru wkh fdofxodwlrq ri glvfuhwh uhwdugdwlrq vshfwud iurp fuhhs ru uhfryhu| gdwd1 D olqhdu uhjuhvvlrq whfkqltxh zdv dssolhg wr wkh phdvxuhg gdwd lq wkh lwhudwlyh fdofxodwlrq ri d vshfwuxp zlwk orjdulwkplf htxlglvwdqw vsdflqj ri uhwdugdwlrq wlphv1 Ehfdxvh wkh qxpehu ri uhwdugdwlrq wlphv zdv olplwhg/ wkh loo0srvhg sureohp zdv dyrlghg1 Lq wkh iroorzlqj/ d qhz phwkrg lv suhvhqwhg dqg dssolhg wr vlpxodwhg g|qdplf gdwd wr rewdlq wkh fuhhs ru uhod{dwlrq pdvwhu fxuyhv1 Frpsdulvrqv zlwk h{dfw uhvsrqvhv duh wkhq xvhg wr looxvwudwh wkh dssolfdelolw| ri wkh phwkrg1 8
Wkhruhwlfdo Edvlv Wkh uhwdugdwlrq wlph vshfwuxp dqg wkh uhod{dwlrq wlph vshfwuxp duh pdwkhpdwlfdo ghvfulswlrqv ri wkh pdwhuldov pdfurvfrslf ehkdylrxu dqg gr qrw kdyh d vlpsoh lqwhusuhwdwlrq lq prohfxodu whupv1 Wkh uhwdugdwlrq wlph vshfwuxp lv d txdqwlwdwlyh ghvfulswlrq ri fuhhs dqg irupv d olqn zlwk wkh ruljlqdo pdwkhpdwlfdo ghvfulswlrqv ri olqhdu ylvfrhodvwlf wkhru| zlwk wkh Erow}pdqq lqwhjudov1 Wkh fuhhs frpsoldqfh aE|fdq eh fdofxodwhg iurp wkh qrq0qhjdwlyh uhwdugdwlrq wlph vshfwuxp uEe| wkh htxdwlrq aE|'a f n] n " 3" uEe 3 | _*? +7, zkhuh a f lv wkh lqvwdqwdqhrxv hodvwlf frpsoldqfh1 Vlplodu uhodwlrqvklsv krog ehwzhhq wkh uhdo dqg lpdjlqdu| sduwv ri wkh frpsoh{ frpsoldqfh a ca dqg wkh uhwdugdwlrq wlph vshfwuxp uE1 Wkh g|qdplf ehkdylrxu xqghu suhvfulehg vlqxvrlgdo vwuhvv fdq eh fdofxodwhg dv a E/'a f n] n " 3" uE n/ 2 2 _*? +8, a E/'] n " 3" uE/ n/ 2 2 _*? +9, Wkh h{shulphqwdo ixqfwlrqv fdq eh ghvfulehg zlwk wkh vdph suhflvlrq e| d glvfuhwh uhwdugdwlrq vshfwuxp1 Wkh uhwdugdwlrq ixqfwlrqv jlyhq e| htxdwlrqv 7 wkurxjk 9 fdq eh h{suhvvhg lq whupv ri glvfuhwh vshfwudo glvwulexwlrqv dv aE|'a f n [ & ' a & e 3 | & +:, 9
a E/'a f n [ & ' a & n/ 2 2 & +;, a E/' [ & ' a & / & n/ 2 2 & +<, zkhuh wkh vhwv ia & c & j/iru&'c2cc / uhsuhvhqw wkh glvfuhwh glvwulexwlrq ri uhwdugdwlrq wlphv1 Iru sk|vlfdo uhdvrqv wkh uhwdugdwlrq vshfwuxp dqg doo wkh frqvwdqwv ia & c & jvkrxog eh qrq0qhjdwlyh1 Wkh htxdwlrqv lqyroylqj wkh orjdulwkplf uhod{dwlrq vshfwuxp/ ME/ ru wkh glvfuhwh uhod{dwlrq vshfwud/ iC & c7 & j/ duh txlwh vlplodu1 Dojrulwkp Zh kdyh d vhw ri P sdluv ri phdvxuhg gdwd srlqwv E) c) cE) 2 c) 2 cc E) c) ri d frpsoh{ ixqfwlrq ) W / lq rxu fdvh wkh frpsoh{ frpsoldqfh/ dw wkh revhuydwlrq srlqwv ri iuhtxhqflhv s cs 2 ccs ) Es 6 a E/ 6 'a f n [ & ' a & n/ 2 6 2 & +43, ) Es 6 a E/ 6 ' [ & ' a & / 6 & n/ 2 6 2 & +44, iru 6'c2cc dqg zlwk / 6 '2Zs 6 1 Wkh sureohp wr eh vroyhg lv krz wr ghwhuplqh wkh 2n xqnqrzq sdudphwhuv zklfk w wkh h{shulphqwdo uhvsrqvh ehvw1 Wkh qxpehu ri uhwdugdwlrq wlphv/ / lv dq rshq fkrlfh/ lq wklv fdvh lw zdv fkrvhq wr eh vpdoohu wkdq wkh qxpehu ri phdvxulqj srlqwv/ 21 :
Wkh h{shulphqwdo iuhtxhqf| zlqgrz zh duh lqwhuhvwhg lv vpdoo hqrxjk wr dyrlg wkh loo0srvhg sureohp1 Lq prvw fdvhv wkh g|qdplf gdwd duh phdvxuhg dw iuhtxhqflhvlqwkhudqjhff wr fM5 ru lq wkh udqjh fwr DfM51 Lq wkhru| erwk Htxdwlrqv 43 dqg 44/ lqghshqghqwo|/ vkrxog jlyh wkh vdph uhvxowv uhsuhvhqwhg e| wkh vhwv ia & c & j & ' c 2 cc 1 Krzhyhu/ vlqfh a E/dqg a E/frqwdlq glhuhqw h{shulphqwdo huuruv/ wkh| fdqqrw jlyh h{dfwo| wkh vdph vroxwlrqv1 Gxh wr wkhvh xqdyrlgdeoh glhuhqfhv zh ghflghg wr frpelqh erwk htxdwlrqv wr ghwhuplqh wkh xqnqrzq sdudphwhuv1 Wkh uvw phwkrg sxeolvkhg wr vroyh wklv nlqg ri sureohpv zdv wkh froorfdwlrq phwkrg/ zklfk zdv vxjjhvwhg e| Vfkdshu| ^47`1 Wklv phwkrg olqhdul}hv wkh sureohp e| suhvfulelqj wkrvh sdudphwhuv wkdw pdnh wkh sureohp qrq0olqhdu/ lq wklv fdvh wkh i & j & ' c 2 cc sdudphwhuv1 Wkh uhpdlqlqj sdudphwhuv/ a f @?_ ia & j & ' c 2 cc / duh ghwhuplqhg e| d olqhdu v|vwhp ri htxdwlrqv/ zkhuh wkh qxpehu ri gdwd srlqwv lv wkh vdph dv wkh qxpehu ri xqnqrzq sdudphwhuv1 Wkhuhiruh/ doo gdwd srlqwv xvhg duh h{dfwo| uhsurgxfhg e| wkh wkhruhwlfdo fxuyh1 Wkh pdlq sureohp ri froorfdwlrq phwkrg lv wkdw lw rqo| zrunv vxffhvvixoo| li wkh vsdflqj ri wkh froorfdwlrq srlqwv lv qhlwkhu wrr odujh qru wrr forvh/ rwkhuzlvh qhjdwlyh sdudphwhuv duh rewdlqhg1 Wkh uhvxow ri wkh froorfdwlrq phwkrg qdwxudoo| ghshqgv rq wkh fkrlfh ri wkh qxpehu ri gdwd srlqwv/ wkhlu orfdwlrq dqg wkh suhvfulehg sdudphwhuv i & j & ' c 2 cc 1 Wkh phwkrg sursrvhg e| Ndvfkwd dqg Vfkzdu}o ^45` lv dq lwhudwlyh pxowlgdwd phwkrg1 Dv wkh| uhihu/ wkh pxowlgdwd phwkrg xvhv doo gdwd srlqwv dydlodeoh wr ghwhuplqh wkh nsdudphwhuv/ a f @?_ ia & j & ' c 2 cc 1 Wkh v|vwhp lv ryhughwhuplqhg dqg wkh sureohp lv vroyhg lq wkh ohdvw vtxduh vhqvh1 ;
Wkh suhvhqw dojrulwkp grhv qrw h{shfw dq| w|sh ri glvwulexwlrq ri wkh phdvxulqj iuhtxhqflhv/ vr wkh| frxog eh orjdulwkplf htxlglvwdqw ru glvwulexwhg lq dq luuhjxodu pdqqhu rq wkh orjdulwkplf iuhtxhqf| d{lv1 Wkh uhwdugdwlrq wlphv/ i & j & ' c 2 cc / duh suhvfulehg dv orjdulwkplf htxlglvwdqw= | & '| @ & 3 c& 'c2cc +45, Lq wkh phdvxuhphqwv wdnhq/ lw lv wkh uhodwlyh h{shulphqwdo huuru/ zklfk lv dssur{lpdwho| lqghshqghqw ri wkh iuhtxhqf|/ wkdw wdnhv lqwr dffrxqw wkh rughu ri wkh pdjqlwxgh ri wkh frpsoldqfh phdvxuhg1 Wklv hhfwlyho| htxdolvhv wkh lq xhqfh ri wkh huuruv dvvrfldwhg zlwk doo gdwd srlqwv zlwklq wkh iuhtxhqf| zlqgrz1 Wkhuhiruh zh dgrswhg wkh iroorzlqj sdudphwhu iru wkh txdolw| ri wkh ghvfulswlrq ri wkh gdwd= [' [ ' ) Es ) 2 n [ ' ) Es ) 2 +46, zkhuh i) c) j ' c 2 cc uhsuhvhqw wkh vhwv ri phdvxuhg g|qdplf frpsoldqfh gdwd1 Plqlplvlqj wkh vxp ri vtxduhv ri uhodwlyh huuruv [/ frqfhuqlqj wkh vhw ri sdudphwhuv wr eh ghwhuplqhg/ a f @?_ ia & j & ' c 2 cc / |lhogv d vhw ri olqhdu htxdwlrqv iru wkh fdofxodwlrqv ri wkrvh sdudphwhuv dv Y[ Ya f 'faY[ Ya & 'fE&'c2cc+47, Wklv ohdgv wr wkh iroorzlqj olqhdu v|vwhp ri htxdwlrqv zulwwhq lq pdwul{ irup dv <
wkh h{dfw ixqfwlrq +h{dfw, zlwk wkh uhodwlyh huuru1 Wkh glvfuhwh vshfwud dqg wkh h{dfw vshfwuxp frlqflgh yhu| zhoo rq wkh orjdulwkplf vfdoh ri Iljxuh <1 Dv vkrzq lq Iljxuh 43 wkh dssur{lpdwlrq ri wkh fuhhs frpsoldqfh h{klelwv d uhodwlyh huuru orzhu wkdq efI1 Zh zloo ghprqvwudwh wkdw wkh phwkrg vwloo zrunv zlwk qrlv| gdwd iru wkhvh fdvhv1 Dv lq wkh uvw whvw zh lpsrvh vrph udqgrp qrlvh rq wkh g|qdplf gdwd e| pxowlso|lqj wkh h{dfw ydoxhv e| d idfwru ri E n Yo@?_E1 Lq Iljxuhv 44 dqg 45 wkh fdofxodwhg fuhhs frpsoldqfh +dssur{1, lv sorwwhg djdlqvw wkh h{dfw ixqfwlrq +h{dfw, dorqj zlwk wkh uhodwlyh huuru iru erwk fdvhv1 Lw fdq hdvlo| eh vhhq wkdw wkh fdofxodwhg fxuyhv pdlqwdlq wkh vkdsh ri wkh h{dfw fxuyhv exw duh vkliwhg dorqj wkh wlph d{lv1 Wklv vkliwlqj lv wkh uhvxow ri fxpxodwlyh huuruv lq fdofxodwlrq ri wkh whpshudwxuh0wlph vfdoh idfwruv gxh wr wkh udqgrp qrlvh1 Wkh uhodwlyh huuruv ghshqg rq wkh ghyldwlrq dqg rq wkh fxuyh*v vorsh dqg wkhuhiruh iru wkh vdph dprxqw ri ghyldwlrq/ wkh uhodwlyh huuru lqfuhdvh zlwk wkh lqfuhdvh ri wkh fxuyh*v vorsh1 Dv iru wkh dfwlydwlrq hqhuj| ri wkh Duukhqlxv htxdwlrq/ wkh suhvhqw phwkrg fdofxodwhv wkrvh sdudphwhuv zlwk d uhodwlyh huuru orzhu wkdq .DI1 Frqfoxvlrqv Zh kdyh suhvhqwhg dq hdv| wr xvh dojrulwkp lq rughu wr frqvwuxfw wkh pdvwhu fuhhs frpsoldqfh fxuyhv iurp wkh h{shulphqwdo g|qdplf ylvfrhodvwlf ixqfwlrqv rewdlqhg rq d GPWD dqdo|vhu1 Wkh phwkrg fdq dovr eh xvhg wr rewdlq dq dssur{lpdwlrq ri wkh glvfuhwh uhwdugdwlrq vshfwud1 49
Iurp huuru0iuhh gdwd zh irxqg wkdw wkh fdofxodwhg pdvwhu fuhhs frpsoldqfh fxuyh qhduo| frlqflghv zlwk wkh h{dfw rqh1 Dv iru wkh uhwdugdwlrq vshfwuxp/ zh frqfoxgh wkdw iru uhjlrqv zkhuh wkh uhwdugdwlrq vshfwuxp lv hlwkhu frqvwdqw ru reh|v d srzhu odz ehkdylrxu wkhvh duh uhsurgxfhg txlwh zhoo e| wkh dojrulwkp1 Uhjlrqv ri vwhhs fkdqjhv lq wkh vshfwuxp duh vwloo uhdvrqdeo| uhsurgxfhg/ dowkrxjk lw whqgv wr vfdwwhu qhdu kljk mxpsv lq wkh orjdulwkplf vfdoh1 Qrlv| gdwd gr qrw dhfw vljqlfdqwo| wkh shuirupdqfh ri wkh dojrulwkp khuh suhvhqwhg/ vlqfh lw dfwv dv d vprrwklqj wrro1 Dfnqrzohgjhphqwv Wkh qdqfldo vxssruw ri wklv zrun e| wkh MQLFW +Mxqwd Qdflrqdo gh Lqyhvwljdêær Flhqwðfd h Whfqroöjlfd, xqghu wkh Judqw EG :39 lv judwhixoo| dfnqrzohgjhg1 Uhihuhqfhv 41 Sryrogr/ I1/ Holgd E1 Khuplgd/ Mrxuqdo ri Pdwhuldov Vflhqfh/ 58/ 736907374 +4<<3,1 51 Ihuu|/ Mrkq G1/ Ylvfrhodvwlf Surshuwlhv ri Sro|phuv/ Mrkq Zloh| ) Vrqv +4<;3,/ ss1 <31 61 Vfkzdu}o/ I1 U1/ Ukhro Dfwd / 47/ 8;408<3 +4<:8,1 71 Jlevrq/ U1 I1 / V1 M1 Kzdqj/ F1 K1 Vkhssdug/ Mrxuqdo ri Frpsrvlwh Pdwhuldov/ 57/ +4<<3,1 4:
81 Pdlohu/ S1/ Eljdxg/ G1 ) Kdpholq/ S1/ Gxudelolw| Hydoxdwlrq rq Vriw Phpeudqh Frpsrvlwhv Xqghu Frpelqhg Ordglqj +phfkdqlfdo dqg hqylurqphqwdo,/ lq Surjuhvv lq Gxudelolw| Dqdo|vlv ri Frpsrvlwh V|vwhpv/ Fdugrq/ Ixnxgd ) Uhlivqlghu+hgv1,/ Edonhpd/ Urwwhugdp +4<<9,/ ss1 4<80 5351 91 Hpul/ L1 ) Q1 Z1 Wvfkrhjo/ Ukhro Dfwd/ 65/ 6440654 +4<<6,1 :1 Wkljshq/ O1/ J1 Z1 Khgvwurp/ E1 S1Erqqhu/ Mrxuqdo ri Dssolhg Phfkdqlfv/ 83/ 6940699 +4<;6,1 ;1 Wvfkrhjo/ Q1 Z1/ Wkh skhqrphqrorjlfdo wkhru| ri olqhdu ylvfrhodvwlf ehkdylrxu/ Vsulqjhu0Yhuodj +4<;<,1 <1 \rqj0Vrn R1/ F1 Slorw ) F1 Mrugdq/ Prghol}dwlrq ri Ylvfrhodvwlf Ehkdylru ri Fduerq2Hsr{| Xqlgluhfwlrqdo Frpsrvlwhv/ lq= Surfhhglqj ri wkh Hxurshdq Phfkdqlfv Froortxlxp 4;5 +4<;7,/ ss1 :60;61 431 Hpul/ L1 ) Q1 Z1 Wvfkrhjo/ Ukhro Dfwd/ 65/ 655065: +4<<6,1 441 Hpul/ L1 ) Q1 Z1 Wvfkrhjo/ Ukhro Dfwd/ 66/ 930:3 +4<<7,1 451 Ndvfkwd/ M1 ) I1 U1 Vfkzdu}o/ Ukhro Dfwd/ 66/ 84:085< +4<<7,1 461 Ndvfkwd/ M1 ) I1 U1 Vfkzdu}o/ Ukhro Dfwd/ 66/ 8630874 +4<<7,1 471 Vfkdshu|/ U1 D1/ Dssur{lpdwlrq phwkrgv ri wudqvirup lqyhuvlrq iru ylvfrhodvwlf vwuhvv dqdo|vlv/ lq= Surfhhglqjv ri wkh 7wk X1V1 Qdwlrqdo 4;
Frqjuhvv ri Dssolhg Phfkdqlfv/ DVPH/ 5 +4<95,/ ss1 43:8043;81 4<
Wdeoh L1 G|qdplf fxuyhv iru wkh vshfwud jlyhq e| Htxdwlrq 531 sEM5a EC@ 3 a EC@ 3 fff eSfDDe D2fef f 3 f2 ff2f eDffD bbbDe f 3 f2 fff eeDe feebf.D fff 2.DDD. fD.Db. fff fHS2SS. fe.b. fff fDfb2 Db.feH f 3 f2 fff fe.e.. bHbbS2 f 3 f2 53
1.E-07 1.E-06 1.E-05 1.E-04 1.E-03 1.E-02 1.E-01 110 1.E-04 1.E-03 1.E-02 1.E-01 1.E+00 a Ψ τ1 Ψ τ1 ; s positive parameters N=5 Iljxuh 4= Plqlpxp ri wkh vxp ri uhodwlyh huuruv/ [/ dqg wkh fruhvsrqghqw uhwdugdwlrq wlph/ / iru wkh fdofxodwlrq ri wkh g|qdplf frpsoldqfhv/ dv ixqfwlrqv ri wkh orjdulwkplf glvwdqfh/ @/ ehwzhhq wkh uhwdugdwlrq wlphv1 -0.10 0.10 0.30 0.50 0.70 0.90 1.10 1.E-05 1.E-03 1.E-01 1.E+01 1.E+03 1.E+05 τ1 ; s L(τ) , Jk/ln(2.5) ; GPa-1 Iljxuh 5= Qrupdolvhg uhwdugdwlrq frpsrqhqwv a&*,?E@iru @'2Dyv h{dfw uhwdugdwlrq vshfwuxp/ Ht1 531 54
0.00 0.50 1.00 1.50 0.001 0.01 0.1 1 10 100-12 -8 -4 0 4 8 12 t ; s J(t) ; GPa-1 error ; (%) exact approx. error Iljxuh 6= Dssur{lpdwhg fuhhs frpsoldqfh iru @'2Dyv h{dfw fuhhs frpsoldqfh dv fdofxodwhg iurp vshfwuxp jlyhq e| Ht1 531 1.E-07 1.E-06 1.E-05 1.E-04 1.E-03 1.E-02 1.E-01 110100 1.E-04 1.E-03 1.E-02 1.E-01 1.E+00 a Ψ τ1 Ψ τ1 ; s positive parameters N=5 10% Noise Iljxuh 7= Plqlpxp ri [dqg wkh fruhvsrqghqw uhwdugdwlrq wlph/ /iruwkh fdofxodwlrq ri wkh g|qdplf frpsoldqfhv zlwk fI qrlvh/ dv ixqfwlrqv ri wkh orjdulwkplf glvwdqfh/ @/ ehwzhhq wkh uhwdugdwlrq wlphv1 55
-0.10 0.10 0.30 0.50 0.70 0.90 1.10 1.E-05 1.E-03 1.E-01 1.E+01 1.E+03 1.E+05 τ ; s L(τ) , Jk/ln(2.5) ; GPa-1 Iljxuh 8= Qrupdolvhg uhwdugdwlrq frpsrqhqwv a&*,?E@iru @'2Ddqg iru fI qrlvh g|qdplf gdwd +Y'f, yv h{dfw uhwdugdwlrq vshfwuxp/ Ht1 531 0.00 0.50 1.00 1.50 0.001 0.01 0.1 1 10 100-12 -8 -4 0 4 8 12 t ; s J(t) ; GPa-1 error ; (%) exact approx. error Iljxuh 9= Dssur{lpdwhg fuhhs frpsoldqfh iru @'2Ddqg iru fI qrlvh g|qdplf gdwd +Y'f, yv h{dfw fuhhs frpsoldqfh dv fdofxodwhg iurp vshfwuxp jlyhq e| Ht1 531 56
1.E-04 1.E-03 1.E-02 1.E-01 1.E+00 1.E+01 1.E-07 1.E-03 1.E+01 1.E+05 1.E+09 1.E+13 1.E+17 τ ; s L(τ) , Jk/ln(10) ; GPa-1 Iljxuh := Grxeoh0orjdulwkplf sorw ri wkh qrupdolvhg uhwdugdwlrq frpsrqhqwv a&*,?E@iru @'fyv wkh h{dfw uhwdugdwlrq vshfwuxp/ Ht1 541 1.E-01 1.E+00 1.E+01 1.E+02 1.E-03 1.E+01 1.E+05 1.E+09 1.E+13 1.E+17 -12 -8 -4 0 4 8 12 t ; s J(t) ; GPa-1 error ; (%) exact approx. error Iljxuh ;= Grxeoh0orjdulwkplf sorw ri wkh dssur{lpdwhg fuhhs frpsoldqfh iru @'fyv h{dfw fuhhs frpsoldqfh dv fdofxodwhg iurp vshfwuxp jlyhq e| Ht1 541 57
1.E-10 1.E-08 1.E-06 1.E-04 1.E-02 1.E+00 1.E+02 1.E-10 1.E-04 1.E+02 1.E+08 1.E+14 1.E+20 τ ; s L(τ) , Jk/ln(10) ; GPa-1 Iljxuh <= Grxeoh0orjdulwkplf sorw ri wkh qrupdolvhg uhwdugdwlrq frpsrqhqwv a&*,?E@iru @'fyv h{dfw uhwdugdwlrq vshfwuxp/ Ht1 551 1.E-01 1.E+00 1.E+01 1.E+02 1.E-03 1.E+01 1.E+05 1.E+09 1.E+13 1.E+17 -12 -8 -4 0 4 8 12 t ; s J(t) ; GPa-1 error ; (%) exact approx. error Iljxuh 43= Grxeoh0orjdulwkplf sorw ri wkh dssur{lpdwhg fuhhs frpsoldqfh iru @'fyv h{dfw fuhhs frpsoldqfh dv fdofxodwhg iurp vshfwuxp jlyhq e| Ht1 55 58