| | shǒu xiān shì yòu huá zhōng kē jì dà xué de dèng jù lóng jiào shòu tí chū de lǐ lùn。
tōng guò shǎo liàng de、 bù wán quán de xìn xī， jiàn lì huī sè wēi fēn yù cè mó xíng， duì shì wù fā zhǎn guī lǜ zuò chū mó hú xìng de cháng qī miáo shù（ mó hú yù cè lǐng yù zhōng lǐ lùn、 fāng fǎ jiào wéi wán shàn de yù cè xué fēn zhī）。
huī sè lǐ lùn rèn wéi xì tǒng de xíng wéi xiàn xiàng jìn guǎn shì méng lóng de， shù jù shì fù zá de， dàn tā bì jìng shì yòu xù de， shì yòu zhěng tǐ gōng néng de。 huī shù de shēng chéng， jiù shì cóng zá luàn zhōng xún zhǎo chū guī lǜ。 tóng shí， huī sè lǐ lùn jiàn lì de shì shēng chéng shù jù mó xíng， bù shì yuán shǐ shù jù mó xíng， yīn cǐ， huī sè yù cè de shù jù shì tōng guò shēng chéng shù jù de gm(1,1) mó xíng suǒ dé dào de yù cè zhí de nì chǔlǐ jiēguǒ。
<1>、 guān lián dù
tí chū xì tǒng de guān lián dù fēn xī fāng fǎ， shì duì xì tǒng fā zhǎn tài shì de liàng huà bǐ jiào fēn xī。 guān lián dù de yī bān biǎo dá shì wéi：
n
ri=1/n ∑ xi(k)
i=1
ri shì qū xiàn xi duì cān kǎo qū xiàn x0 de guān lián dù。
<2>、 shēng chéng shù
tōng guò duì yuán shǐ shù jù de zhěng lǐ xún zhǎo shù de guī lǜ， fēn wéi sān lèi：
a、 lěi jiā shēng chéng： tōng guò shù liè jiān gè shí kè shù jù de yǐ gè lěi jiā dé dào xīn de shù jù yǔ shù liè。 lěi jiā qián shù liè wéi yuán shǐ shù liè， lěi jiā hòu wéi shēng chéng shù liè。
jī běn guān xì shì：
jì x(0) wéi yuán shǐ shù liè
x(0)=(x(0)(k)xk=1,2, … ,n)=(x(0)(1),x(0)(2),… ,x(0)(n))
jì x(1) wéi shēng chéng shù liè
x(1)=(x(1)(k)xk=1,2,… ,n)=(x(1)(1),x(1)(2),… ,x(1)(n))
rú guǒ x(0) yǔ x(1) zhī jiān mǎn zú xià liè guān xì， jí
k
x(1)(k)=∑ x(0)(i)
i=a
chēng wéi yī cì lěi jiā shēng chéng。
b、 lěi jiǎn shēng chéng： qián hòu liǎng gè shù jù zhī chā， lěi jiā shēng chéng de nì yùn suàn。 lěi jiǎn shēng chéng kě jiāng lěi jiā shēng chéng hái yuán chéng fēi shēng chéng shù liè。
c、 yìng shè shēng chéng： lěi jiā、 lěi jiǎn yǐ wài de shēng chéng fāng shì。
<3>、 jiàn lì mó xíng
a、 jiàn mó jī lǐ
b、 bǎ yuán shǐ shù jù jiā gōng chéng shēng chéng shù；
c、 duì cán chā ( mó xíng jì suàn zhí yǔ shí jì zhí zhī chā ) xiū dìng hòu， jiàn lì chā fēn wēi fēn fāng chéng mó xíng；
d、 jī yú guān lián dù shōu liǎn de fēn xī；
e、 gm mó xíng suǒ dé shù jù xū jīng guò nì shēng chéng hái yuán hòu cái néng yòng。
f、 cǎi yòng “ wǔ bù jiàn mó（ xì tǒng dìng xìng fēn xī、 yīn sù fēn xī、 chū bù liàng huà、 dòng tài liàng huà、 yōu huà） ” fǎ， jiàn lì yī zhǒng chā fēn wēi fēn fāng chéng mó xíng gm(1,1) yù cè mó xíng。
jī běn suàn shì wéi：
lìng x(0)=(x(0)(1),x(0)(2),… ,x(0)(n))
zuò yī cì lěi jiā shēng chéng， k
x(1)(k)=∑ x(0)(m)
m=1
yòu x(1)=(x(1)(1),x(1)(2),… ,x(1)(n))
=(x(0)(1),x(1)(1)+x(0)(2),… ,x(1)(n-1)+x(0)(n))
x(1) kě jiàn lì bái huà fāng chéng： dx(1)/dt+ax(1)=u jí gm(1,1).
gāi fāng chéng de jiě wéi :x(1)(k+1)=(x(1)(1)-u/a)e-ak+u/a
<4>、 yù cè fāng fǎ：
a、 shù liè yù cè
b、 zāi biàn yù cè
c、 jì jié zāi biàn yù cè
d、 tuò pū yù cè
e、 xì tǒng zōng hé yù cè
f、 mó hú yù cè
duì yú yī gè mó hú xì tǒng lái shuō， chuán tǒng de yù cè fāng fǎ jiù huì shī qù zuò yòng。 chǔlǐ mó hú yù cè wèn tí de shù xué fāng fǎ shì mó hú shù xué。 mó hú shù xué de jī chǔ shì mó hú jí hé lùn， ér mó hú jí hé shì pǔ tōng jí hé de kuò zhǎn。 měi guó xué zhě l.a.zadeh jiào shòu jiàn lì de mó hú jí hé lùn， wéi mó hú yù cè lǐ lùn yǔ fāng fǎ de yán jiū diàn dìng liǎo lǐ lùn jī chǔ。 tā yòng jiǎn jié yòu lì de fāng fǎ chǔlǐ fù zá xì tǒng， zài mǒu zhǒng chéng dù shàng mí bǔ liǎo jīng diǎn shù xué yǔ tǒng jì shù xué de bù zú。
zài yù cè yìng yòng shàng， rú qì xiàng yù bào、 dì zhèn yù bào、 bìng chóng hài yù bào děng， guó nèi xué zhě zuò chū liǎo xǔ duō yòu yì de yán jiū。
4、 jī yú hùn dùn lǐ lùn de fēn xī yù cè
hùn dùn lǐ lùn shì jìn nián lái cháng zú fā zhǎn de yī mén xué kē。 hùn dùn xiàng shì jiè guī lǜ yùn dòng de jiǎ dìng xìng tí chū tiǎo zhàn。 yī fāng miàn， tā gào sù wǒ men， yǔ zhòu yuǎn bǐ wǒ men xiǎng dé dào de yào guài yì， tā shǐ xǔ duō chuán tǒng de kē xué fāng fǎ shòu dào huái yí。 lìng yī fāng miàn， hùn dùn rèn wéi xǔ duō wú guī zé de shì wù shí jì shàng kě néng shì jiǎn dān guī lǜ de jiēguǒ。 hùn dùn zhǎn xiàn gěi wǒ men de shì yī xiē xīn de guī lǜ。 tí chū： zūn cóng jiǎn dān guī lǜ de xì tǒng huì yǐ lìng rén jīng yà de fù zá fāng shì biǎo xiàn qí xíng wéi。 hùn dùn shì yǐn mì xíng shì de zhì xù。
hùn dùn xì tǒng shì zhǐ mǐn gǎn dì yǐ lài yú chū shǐ tiáo jiàn de nèi zài biàn huà xì tǒng， duì wài lái biàn huà de mǐn gǎn xìng běn shēn bìng bù yì wèi zhe hùn dùn。 hùn dùn lǐ lùn zuì lìng rén xīng fèn de shì： yī gè fēi cháng jiǎn dān de jué dìng lùn xì tǒng néng gòu chǎn shēng yì cháng fù zá de shū chū jiēguǒ。 gěi dìng yī gè jiǎn dān guī zé hé chū shǐ tiáo jiàn， xì tǒng jiāng chǎn shēng fù zá lián xù xì liè， zhè yī diǎn lèi sì“ wú zhōng shēng yòu”。
měi guó kē xué jiā pà kǎ dé hé tā de tóng shì jī yú hùn dùn hé shēng wù jìn huà lǐ lùn， jiè zhù jì suàn jī， zhì lì yú yòng tú xíng lái miáo shù jīn róng shì chǎng de hùn dùn xiàn xiàng。 pà kǎ dé rèn wéi， shì jiè shàng yòu dà liàng bù tóng de suí jī xiàn xiàng， tā suǒ yán jiū de shì dà tǐ zhǐ xū jǐ gè biàn liàng jiù néng miáo shù xì tǒng xíng wéi de yī zhǒng hùn dùn xiàn xiàng。 tā shì tú jiàn lì yī zhǒng xué xí suàn fǎ， duì jìn huà mó xíng jìn xíng chǔlǐ。 ér duì yú zhòng duō de mó xíng， pà kǎ dé cǎi yòng yī zhǒng chēng wéi yí chuán suàn fǎ de fāng fǎ chǔlǐ shù jù。 tā yòng lèi sì shēng wù fán zhí zhōng tū biàn hé zá jiāo xiàn xiàng de fāng fǎ lái gǎi biàn mó xíng。 zhè zhǒng fāng fǎ de hé xīn shì， jì suàn jī bù duàn shè dìng xīn de jiǎ shè huán jìng， cóng 'ér shǐ xué xí suàn fǎ gèng jù yòu shì yìng xìng。 rèn wéi yī gè hǎo de xué xí suàn fǎ bù jǐn néng jiàn lì shì yìng mó xíng， tā hái néng shí kè guān cè shù jù de biàn huà。 suǒ wèi“ xué xí suàn fǎ” shì yī zhǒng tè bié de chéng xù， tā shàn cháng duì dà liàng de、 gè zhǒng gè yàng de mó xíng jìn xíng bǐ jiào yán jiū， zhǎo chū nǎ gè mó xíng zuì shì yòng yú fēn xī mù qián hé wèi lái de shù jù。 | | shǒu xiān shì yòu huá zhōng kē jì dà xué de dèng jù lóng jiào shòu tí chū de lǐ lùn。
tōng guò shǎo liàng de、 bù wán quán de xìn xī， jiàn lì huī sè wēi fēn yù cè mó xíng， duì shì wù fā zhǎn guī lǜ zuò chū mó hú xìng de cháng qī miáo shù（ mó hú yù cè lǐng yù zhōng lǐ lùn、 fāng fǎ jiào wéi wán shàn de yù cè xué fēn zhī）。
huī sè lǐ lùn rèn wéi xì tǒng de xíng wéi xiàn xiàng jìn guǎn shì méng lóng de， shù jù shì fù zá de， dàn tā bì jìng shì yòu xù de， shì yòu zhěng tǐ gōng néng de。 huī shù de shēng chéng， jiù shì cóng zá luàn zhōng xún zhǎo chū guī lǜ。 tóng shí， huī sè lǐ lùn jiàn lì de shì shēng chéng shù jù mó xíng， bù shì yuán shǐ shù jù mó xíng， yīn cǐ， huī sè yù cè de shù jù shì tōng guò shēng chéng shù jù de gm(1,1) mó xíng suǒ dé dào de yù cè zhí de nì chǔlǐ jiēguǒ。
<1>、 guān lián dù
tí chū xì tǒng de guān lián dù fēn xī fāng fǎ， shì duì xì tǒng fā zhǎn tài shì de liàng huà bǐ jiào fēn xī。 guān lián dù de yī bān biǎo dá shì wéi：
n
ri=1/n∑ xi(k)
i=1
ri shì qū xiàn xi duì cān kǎo qū xiàn x0 de guān lián dù。
<2>、 shēng chéng shù
tōng guò duì yuán shǐ shù jù de zhěng lǐ xún zhǎo shù de guī lǜ， fēn wéi sān lèi：
a、 lěi jiā shēng chéng： tōng guò shù liè jiān gè shí kè shù jù de yǐ gè lěi jiā dé dào xīn de shù jù yǔ shù liè。 lěi jiā qián shù liè wéi yuán shǐ shù liè， lěi jiā hòu wéi shēng chéng shù liè。
jī běn guān xì shì：
jì x(0) wéi yuán shǐ shù liè
x(0)=(x(0)(k)xk=1,2,… ,n)=(x(0)(1),x(0)(2),… ,x(0)(n))
jì x(1) wéi shēng chéng shù liè
x(1)=(x(1)(k)xk=1,2,… ,n)=(x(1)(1),x(1)(2),… ,x(1)(n))
rú guǒ x(0) yǔ x(1) zhī jiān mǎn zú xià liè guān xì， jí
k
x(1)(k)=∑ x(0)(i)
i=a
chēng wéi yī cì lěi jiā shēng chéng。
b、 lěi jiǎn shēng chéng： qián hòu liǎng gè shù jù zhī chā， lěi jiā shēng chéng de nì yùn suàn。 lěi jiǎn shēng chéng kě jiāng lěi jiā shēng chéng hái yuán chéng fēi shēng chéng shù liè。
c、 yìng shè shēng chéng： lěi jiā、 lěi jiǎn yǐ wài de shēng chéng fāng shì。
<3>、 jiàn lì mó xíng
a、 jiàn mó jī lǐ
b、 bǎ yuán shǐ shù jù jiā gōng chéng shēng chéng shù；
c、 duì cán chā ( mó xíng jì suàn zhí yǔ shí jì zhí zhī chā ) xiū dìng hòu， jiàn lì chā fēn wēi fēn fāng chéng mó xíng；
d、 jī yú guān lián dù shōu liǎn de fēn xī；
e、 gm mó xíng suǒ dé shù jù xū jīng guò nì shēng chéng hái yuán hòu cái néng yòng。
f、 cǎi yòng“ wǔ bù jiàn mó（ xì tǒng dìng xìng fēn xī、 yīn sù fēn xī、 chū bù liàng huà、 dòng tài liàng huà、 yōu huà）” fǎ， jiàn lì yī zhǒng chā fēn wēi fēn fāng chéng mó xíng gm(1,1) yù cè mó xíng。
jī běn suàn shì wéi：
lìng x(0)=(x(0)(1),x(0)(2),… ,x(0)(n))
zuò yī cì lěi jiā shēng chéng， k
x(1)(k)=∑ x(0)(m)
m=1
yòu x(1)=(x(1)(1),x(1)(2),… ,x(1)(n))
=(x(0)(1),x(1)(1)+x(0)(2),… ,x(1)(n-1)+x(0)(n))
x(1) kě jiàn lì bái huà fāng chéng： dx(1)/dt+ax(1)=u jí gm(1,1).
gāi fāng chéng de jiě wéi :x(1)(k+1)=(x(1)(1)-u/a)e-ak+u/a
<4>、 yù cè fāng fǎ：
a、 shù liè yù cè
b、 zāi biàn yù cè
c、 jì jié zāi biàn yù cè
d、 tuò pū yù cè
e、 xì tǒng zōng hé yù cè
f、 mó hú yù cè
duì yú yī gè mó hú xì tǒng lái shuō， chuán tǒng de yù cè fāng fǎ jiù huì shī qù zuò yòng。 chǔlǐ mó hú yù cè wèn tí de shù xué fāng fǎ shì mó hú shù xué。 mó hú shù xué de jī chǔ shì mó hú jí hé lùn， ér mó hú jí hé shì pǔ tōng jí hé de kuò zhǎn。 měi guó xué zhě l.a.zadeh jiào shòu jiàn lì de mó hú jí hé lùn， wéi mó hú yù cè lǐ lùn yǔ fāng fǎ de yán jiū diàn dìng liǎo lǐ lùn jī chǔ。 tā yòng jiǎn jié yòu lì de fāng fǎ chǔlǐ fù zá xì tǒng， zài mǒu zhǒng chéng dù shàng mí bǔ liǎo jīng diǎn shù xué yǔ tǒng jì shù xué de bù zú。
zài yù cè yìng yòng shàng， rú qì xiàng yù bào、 dì zhèn yù bào、 bìng chóng hài yù bào děng， guó nèi xué zhě zuò chū liǎo xǔ duō yòu yì de yán jiū。
4、 jī yú hùn dùn lǐ lùn de fēn xī yù cè
hùn dùn lǐ lùn shì jìn nián lái cháng zú fā zhǎn de yī mén xué kē。 hùn dùn xiàng shì jiè guī lǜ yùn dòng de jiǎ dìng xìng tí chū tiǎo zhàn。 yī fāng miàn， tā gào sù wǒ men， yǔ zhòu yuǎn bǐ wǒ men xiǎng dé dào de yào guài yì， tā shǐ xǔ duō chuán tǒng de kē xué fāng fǎ shòu dào huái yí。 lìng yī fāng miàn， hùn dùn rèn wéi xǔ duō wú guī zé de shì wù shí jì shàng kě néng shì jiǎn dān guī lǜ de jiēguǒ。 hùn dùn zhǎn xiàn gěi wǒ men de shì yī xiē xīn de guī lǜ。 tí chū： zūn cóng jiǎn dān guī lǜ de xì tǒng huì yǐ lìng rén jīng yà de fù zá fāng shì biǎo xiàn qí xíng wéi。 hùn dùn shì yǐn mì xíng shì de zhì xù。
hùn dùn xì tǒng shì zhǐ mǐn gǎn dì yǐ lài yú chū shǐ tiáo jiàn de nèi zài biàn huà xì tǒng， duì wài lái biàn huà de mǐn gǎn xìng běn shēn bìng bù yì wèi zhe hùn dùn。 hùn dùn lǐ lùn zuì lìng rén xīng fèn de shì： yī gè fēi cháng jiǎn dān de jué dìng lùn xì tǒng néng gòu chǎn shēng yì cháng fù zá de shū chū jiēguǒ。 gěi dìng yī gè jiǎn dān guī zé hé chū shǐ tiáo jiàn， xì tǒng jiāng chǎn shēng fù zá lián xù xì liè， zhè yī diǎn lèi sì“ wú zhōng shēng yòu”。
měi guó kē xué jiā pà kǎ dé hé tā de tóng shì jī yú hùn dùn hé shēng wù jìn huà lǐ lùn， jiè zhù jì suàn jī， zhì lì yú yòng tú xíng lái miáo shù jīn róng shì chǎng de hùn dùn xiàn xiàng。 pà kǎ dé rèn wéi， shì jiè shàng yòu dà liàng bù tóng de suí jī xiàn xiàng， tā suǒ yán jiū de shì dà tǐ zhǐ xū jǐ gè biàn liàng jiù néng miáo shù xì tǒng xíng wéi de yī zhǒng hùn dùn xiàn xiàng。 tā shì tú jiàn lì yī zhǒng xué xí suàn fǎ， duì jìn huà mó xíng jìn xíng chǔlǐ。 ér duì yú zhòng duō de mó xíng， pà kǎ dé cǎi yòng yī zhǒng chēng wéi yí chuán suàn fǎ de fāng fǎ chǔlǐ shù jù。 tā yòng lèi sì shēng wù fán zhí zhōng tū biàn hé zá jiāo xiàn xiàng de fāng fǎ lái gǎi biàn mó xíng。 zhè zhǒng fāng fǎ de hé xīn shì， jì suàn jī bù duàn shè dìng xīn de jiǎ shè huán jìng， cóng 'ér shǐ xué xí suàn fǎ gèng jù yòu shì yìng xìng。 rèn wéi yī gè hǎo de xué xí suàn fǎ bù jǐn néng jiàn lì shì yìng mó xíng， tā hái néng shí kè guān cè shù jù de biàn huà。 suǒ wèi“ xué xí suàn fǎ” shì yī zhǒng tè bié de chéng xù， tā shàn cháng duì dà liàng de、 gè zhǒng gè yàng de mó xíng jìn xíng bǐ jiào yán jiū， zhǎo chū nǎ gè mó xíng zuì shì yòng yú fēn xī mù qián hé wèi lái de shù jù。
matlab chéng xù：
% xià miàn chéng xù shì huī sè mó xíng GM(1,1) chéng xù 'èr cì nǐ hé hé děng wéi xīn chén dài xiè gǎi jìn yù cè chéng xù ,matlab6.5, shǐ yòng běn chéng xù qǐng zhù míng， chéng xù cún chǔ wéi gm1.m
%x=[5999,5903,5848,5700,7884];gm1(x); cè shì shù jù
% èr cì nǐ hé yù cè GM(1,1) mó xíng
functiongmcal=gm1(x)
ynyz=[
124113154.1666667240.0909091224.3333333
180.9155.0909091223.4166667274.8333333327.25
33.3333333392.557.5382.5101.3333333
6069.7562.66666667136.6666667202.75
48.6666666760.3333333349.7597.25139.6666667
111.498.18181818118.75168.8243.5454545
280.25296.125316.7777778354458.8888889
62.6666666761.3333333395.7539.6666666779
223.1205235.3636364303.8355
148.25162161.1666667243361.6
650522.8629.7142857915.66666671270.4
780.3333333726.66666679241483.6666671656.666667
140104.75168.2322337.3333333
292.6666667344.3333333454.3333333416.3333333943.5
168.3333333161.5339415.6666667262.8
290.3333333270.25317548494
313.25543565.6666667696.75713
308328.4285714379432.1428571521.8333333
1723.25206021073643.6666672416.6
285.75508364.25526.251142
386.3333333511.2558.1666667549.5674.5
348.2857143412.3461538502.96694.88747.56
4091.9142864137.0555564823.6857146175.7222227706.305556
520619547.57861142
643.2352941489.8095238699.66666671133.251298.533333
964.259721220.8095241408.5652172286.909091
210.6234.8304428.5544.8571429
79.75110.6666667107.2283415
265.7333333286.5714286400.8125415425.6875
1873.28571416231854.2222222042.92665.5
219.175311.7406.8928571480.75642.0769231
88.2171233178.8438
33.33333333183223269733
105.4122.8125130.5625265.5235.9444444
134.4104.3333333247.25293385.1666667
134.8147216.8571429186.6666667246.25
33.3333333344.25116.759654.6
124.2147253544530.75
114.5109.5150.5219.3333333300.8571429
57.6666666751.689.6166.6666667253
959.751120.2521883439.54430
2461.8571433007.7142864046.57142947836518.5
1871.21785.8333333962.2524897272
2396.8333332322.7142862950.14285734084164.25
5154.2550928304.3333331957313888.33333
27372675.42857132273299.756860.4
1500.751813.8571432180.16666727223583.857143
14171205.81414.51211.7777781276.9
255278.6363636305.1666667386.4615385596.0769231
2033.33333333100132217.5
894614.3333333545425.3333333563
285.8308.4615385443.75511.2710.1
299.6316443893.33333331023.666667
240278.5714286361.6666667704848.6
5061.33333333100273841
119.2192.25251445715
159.6142.4225.5396.3333333357.5
33.3333333363.6666666771177.5137
256.6666667257375.6666667664.5677
215.1666667194329.6666667361.1428571385
296.6666667265244.25249240.6666667
292.5355.75415.75407.2574.4
317.5306.2408.5389.2543.6
315409.4504.8472.6666667711.4
84.33333333116.5153128253
6321013855.33333331089.3333331361.333333
227.25275.75348361520
225.3333333260.4264.5305426.5
237.5351.8378375.6666667550.8333333
77106.6104.3333333166.8346.8
148443189.3333333379.5604.5
474.51182125010191263.5
];
x=ynyz(72,:);
sizexd2=size(x,2);
% qiú shù zǔ cháng dù
k=0;
fory1=x
k=k+1;
ifk>1
x1(k)=x1(k-1)+x(k);
% lěi jiā shēng chéng
z1(k-1)=-0.5*(x1(k)+x1(k-1));
%z1 wéi shù jiǎn 1， yòng yú jì suàn B
yn1(k-1)=x(k);
else
x1(k)=x(k);
end
end
%x1,z1,k,yn1
sizez1=size(z1,2);
%size(yn1);
z2=z1';
z3=ones(1,sizez1)';
YN=yn1';% zhuǎn zhì
%YN
B=[z2z3];
au0=inv(B'*B)*B'*YN;
au=au0';
%B,au0,au
afor=au(1);
ufor=au(2);
ua=au(2)./au(1);
%afor,ufor,ua
% shū chū yù cè de au hé u/a de zhí
constant1=x(1)-ua;
afor1=-afor;
x1t1='x1(t+1)';
estr='exp';
tstr='t';
leftbra='(';
rightbra=')';
%constant1,afor1,x1t1,estr,tstr,leftbra,rightbra
strcat(x1t1,'=',num2str(constant1),estr,leftbra,num2str(afor1),tstr,rightbra,'+',leftbra,num2str(ua),rightbra)
% shū chū shí jiān xiǎng yìng fāng chéng
%******************************************************
% èr cì nǐ hé
k2=0;
fory2=x1
k2=k2+1;
ifk2>k
else
ze1(k2)=exp(-(k2-1)*afor);
end
end
%ze1
sizeze1=size(ze1,2);
z4=ones(1,sizeze1)';
G=[ze1'z4];
X1=x1';
au20=inv(G'*G)*G'*X1;
au2=au20';
%z4,X1,G,au20
Aval=au2(1);
Bval=au2(2);
%Aval,Bval
% shū chū yù cè de A,B de zhí
strcat(x1t1,'=',num2str(Aval),estr,leftbra,num2str(afor1),tstr,rightbra,'+',leftbra,num2str(Bval),rightbra)
% shū chū shí jiān xiǎng yìng fāng chéng
nfinal=sizexd2-1+1;
% jué dìng yù cè de bù zhòu shù 5 zhè gè bù zhòu kě yǐ tōng guò hán shù chuán rù
%nfinal=sizexd2-1+1;
% yù cè de bù zhòu shù 1
fork3=1:nfinal
x3fcast(k3)=constant1*exp(afor1*k3)+ua;
end
%x3fcast
% yī cì nǐ hé lěi jiā zhí
fork31=nfinal:-1:0
ifk31>1
x31fcast(k31+1)=x3fcast(k31)-x3fcast(k31-1);
else
ifk31>0
x31fcast(k31+1)=x3fcast(k31)-x(1);
else
x31fcast(k31+1)=x(1);
end
end
end
x31fcast
% yī cì nǐ hé yù cè zhí
fork4=1:nfinal
x4fcast(k4)=Aval*exp(afor1*k4)+Bval;
end
%x4fcast
fork41=nfinal:-1:0
ifk41>1
x41fcast(k41+1)=x4fcast(k41)-x4fcast(k41-1);
else
ifk41>0
x41fcast(k41+1)=x4fcast(k41)-x(1);
else
x41fcast(k41+1)=x(1);
end
end
end
x41fcast,x
% èr cì nǐ hé yù cè zhí
%*** jīng dù jiǎn yàn pC************//////////////////////////////////
k5=0;
fory5=x
k5=k5+1;
ifk5>sizexd2
else
err1(k5)=x(k5)-x41fcast(k5);
end
end
%err1
% jué duì wù chā
xavg=mean(x);
%xavg
%x píng jūn zhí
err1avg=mean(err1);
%err1avg
%err1 píng jūn zhí
k5=0;
s1total=0;
fory5=x
k5=k5+1;
ifk5>sizexd2
else
s1total=s1total+(x(k5)-xavg)^2;
end
end
s1suqare=s1total./sizexd2;
s1sqrt=sqrt(s1suqare);
%s1suqare,s1sqrt
%s1suqare cán chā shù liè x de fāng chā s1sqrt wéi x fāng chā de píng fāng gēn S1
k5=0;
s2total=0;
fory5=x
k5=k5+1;
ifk5>sizexd2
else
s2total=s2total+(err1(k5)-err1avg)^2;
end
end
s2suqare=s2total./sizexd2;
%s2suqare cán chā shù liè err1 de fāng chā S2
Cval=sqrt(s2suqare./s1suqare);
Cval
%nnn=0.6745*s1sqrt
%CvalC jiǎn yàn zhí
k5=0;
pnum=0;
fory5=x
k5=k5+1;
ifabs(err1(k5)-err1avg)<0.6745*s1sqrt
pnum=pnum+1;
%ppp=abs(err1(k5)-err1avg)
else
end
end
pval=pnum./sizexd2;
pval
%p jiǎn yàn zhí
%arr1=x41fcast(1:6)
% yù cè jiēguǒ wéi qū jiān fàn wéi yù cè bù cháng hé shù jù cháng dù kě diào zhěng chéng xù cān shù jìn xíng gǎi jìn | | | huī sè yù cè fǎ | shì gù huī sè yù cè mó xíng | |
|
|
|