%灰色预测公司的利润
clc,clear
syms a b;
c=[a,b]';
A=[89677,99215,109655,120333,135823,159878,182321,209407,246619,300670];
B=cumsum(A); %原始数据累加
n=length(A);
for i=1:(n-1)
C(i)=0.5*(B(i)+B(i+1)); %生成累加矩阵
end
%计算待定参数的值
D=A;D(1)=[];
D=D';
E=[-C;ones(1,n-1)];
c=inv(E*E')*E*D;
c=c';
a=c(1);b=c(2);
%预测后续数据
F=[];F(1)=A(1);
for i=2:(n+10)
F(i)=(A(1)-b/a)/exp(a*(i-1))+b/a;
end
G=[];G(1)=A(1);
for i=2:(n+10)
G(i)=F(i)-F(i-1); %预测所得数据
end
t1=1999:2008;
t2=1999:2018;
G
plot(t1,A,'o',t2,G);
xlabel('年份');ylabel('利润');
数模---灰色预测模型代码
最新推荐文章于 2024-04-17 23:21:46 发布
![](https://img-home.csdnimg.cn/images/20240711042549.png)