急求GM 1,N 的matlab程式

時間 2021-09-01 19:33:35

1樓:匿名使用者

% gm(1,1)模型計算及檢驗、作圖。檔名fungry1.m

function gm1=fungry1(x0) %輸入原始資料x0

t=input('t=')

x1=zeros(1,length(x0));

b=zeros(length(x0)-1,2);

yn=zeros(length(x0)-1,1);

hatx0=zeros(1,length(x0)-1,2);

hatx00=zeros(1,length(x0));

hatx1=zeros(1,length(x0)+t);

epsilon=zeros(length(x0),1);

omega=zeros(length(x0),1);

for i=1:length(x0)

for j=1:i

x1(i)=x1(i)+x0(j);

endend

for i=1:length(x0)-1

b(i,1)=(-1/2)*(x1(i)+x1(i+1));

b(i,2)=1;

yn(i)=x0(i+1);

endhata=(inv(b'*b))*b'*yn;% gm(1,1)模型引數估計

a=hata(1)

b=hata(2)

for k=1:length(x0)+t

hatx1(k)=(x0(1)-b/a)*exp(-a*(k-1))+b/a

endhatx0(1)=hatx1(1);

for k=2:length(x0)+t

hatx0(k)=hatx1(k)-hatx1(k-1) % 累減還原得到歷史資料的模擬值

endfor i=1:length(x0) % 開始模型檢驗

epsilon(i)=x0(i)-hatx0(i);

omega(i)=(epsilon(i)/x0(i))*100;

endc=std(epsilon)/std(x0);p=0;

for i=1:length(x0)

if abs(epsilon(i)-mean(epsilon))<0.6745*std(x0)

p=p+1;

endp=p/length(x0)

if p>0.95 & c<0.35

disp('the model is good,and the forecast is:');

disp(hatx0(length(x0)+t))

elseif p>0.85 & c<0.5

disp('the model is eligibility,and the forecast is:');

disp(hatx0(length(x0)+t))

elseif p>0.70 & c<0.65

disp('the model is not good,and the forecast is:');

disp(hatx0(length(x0)+t))

else p<=0.70 & c>0.65

disp('the model is bad and try again')

endfor i=1:length(x0)

hatx00(i)=hatx0(i);

endz=1:length(x0);

plot(z,x0,'-',z,hatx00,':') %將原始資料和模擬值畫在乙個圖上幫助觀察end

2樓:落地生根發芽

難度太大了 我不會。

an 2n 1 n為奇an 3n 1 n為偶求sn急

當n為奇數時,由已知可得an為等差為4,首項為1,n n 2的等差數列,當n為偶數時,由已知可得an為等差為6,首項為7,n n 2的等比數列,把兩種情況的sn加起來,就為sn 這裡有送你們的一首歌 喜歡請收聽哦 解當n為偶數時 sn a1 a2 a3 a4 an a1 a3 a5 a n 1 a2...

用VB編寫程式,輸入正整數n,求所有大於1且小於n,同時與n互質的數,存入陣列並每行的格式輸出

同時與n互質的數 這個不太明白?能說一下嗎?下面1000以內質數演算法。private sub command1 click text1.text for i 1 to 1000 k 1 for j 2 to i 2 if i mod j 0 then k 0next j if k 1 then p...

求古代或穿越的1女n男小說,求1女N男穿越小說

俊男坊 果子文筆也很好的 這個情節超不俗 一女七男 美男十二宮 同下面兩個都是狼大 逍遙紅塵 的文,這個不到最後看不出結局,設計的太好了 有一點點玄幻色彩 對了,這個不是穿越 但一點不影響精彩 同下面的兩篇 都是女尊 只是文筆極好 不捨得舍掉啦 如果樓主不是特別不能接受女尊 可以看一下 尤其是這個 ...