matlab中灰色模型代码,灰色模型matlab代码

灰色预测模型Matlab代码

sr(1) = 11.985;

sr(2) = 12.1121;

sr(3) = 12.2389;

sr(4) = 12.3626;

sr(5) = 12.481;

x1=[0 0 0 0 0];

for j = 1 : 5

for i = 1 : j

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

end

end

for k = 2 : 5

z1(k) = 0.5 * (x1(k) + x1(k - 1));

end

B=[-z1(2) 1;-z1(3) 1;-z1(4) 1;-z1(5) 1];

Yn=[sr(2);sr(3);sr(4);sr(5)];

ajg=inv(B'*B)*B'*Yn;

a=ajg(1,1);

b=ajg(2,1);

for i = 0 : 100

sc1(i+1) = (sr(1) - b / a) * exp(-(a*i)) + b / a; end

for j = 1 : 100

sc0(j + 1) = sc1(j + 1) - sc1(j);

end

sc0

sc1;

等维递补灰色预测模型Matlab代码

format compact

sr(1) = 11.985;

sr(2) = 12.1121;

sr(3) = 12.2389;

sr(4) = 12.3626;

sr(5) = 12.481;

p=0;

while p<=100

x1=[0 0 0 0 0];

for j = 1 : 5

for i = 1 : j

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

end

end

for k = 2 : 5