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