全部版块 我的主页
论坛 计量经济学与统计论坛 五区 计量经济学与统计软件
5952 0
2017-10-10
本人在做一个灰度预测的模型,但是做出来值可以出来就是图出不来,不知道哪里出错了,怎么解决,求助各位大神
clc,clear;
syms a b;
c=[a b]';
A=[1.0123,1.2512,1.4454,1.6205,1.7690,1.8022,1.9600];
B=cumsum(A);  %原始数据累加
n=length(A);
for i=1:(n-1)
    C(i)=(B(i)+B(i+1))/2; %生成累加矩阵
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+6)  %只推测后6个数据,
    F(i)=(A(1)-b/a)/exp(a*(i-1))+b/a;
end
G=[];G(1)=A(1);
for i=2:(n+6)  %只推测后6个数据,
    G(i)=F(i)-F(i-1);  %得到预测出来的数据
end
t1=2010:2016;
t2=2017:2022;  %多6组数据
G
h=plot(t1,A,'o',t2,G,'-'); %原始数据与预测数据的比较
set(h,'LineWidth',1.5);  

二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群