?? greyforecasting.m
字號:
clear
close all
%灰色預測模型
%X0原始數據
%計算原始數據的累加生成值
X0=[2.655 2.819 3.006 3.093 3.277 3.515 3.77 4.107 4.496 4.973 5.451 5.847 6.213 6.775 7.542 8.384 9.138 9.942 10.65 11.198 11.431 12.176 13.51];
n=length(X0);
X1=zeros(1,n);
X1(1)=X0(1);
for i=2:n
X1(i)=X0(i)+X1(i-1);
end
%計算數據矩陣B和數據向量Yn
%采用GM(1,1)模型所對應的數據矩陣為
for i=1:n-1
B(i,1)=-(X1(i)+X1(i+1))/2;
B(i,2)=1;
end
for i=1:n-1
Yn(i,1)=X0(i+1);
end
%計算GM(1,1)微分方程的參數a和u
A=inv(B'*B)*B'*Yn;
%計算預測結果
Xp1=zeros(n,1);
for i=0:n-1+5
Xp1(i+1)=(X0(1)-A(2,1)/A(1,1))*exp(-A(1,1)*i)+A(2,1)/A(1,1);
end
Xp0(1,1)=Xp1(1);
for i=1:n-1+5
Xp0(i+1,1)=Xp1(i+1)-Xp1(i);
end
?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -