Question

一个MATLAB基础问题 为什么他会显示 错误使用 plot 矢量长度必须相同。

yms a b;
c=[a b]';
A=[11.501 13.595 13.957 10,721 8.571 10.813 9.467];% 原始序列
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+1)
F(i)=(A(1)-b/a)/exp(a*(i-1))+b/a;
end
G=[];G(1)=A(1);
for i=2:(n+1)
G(i)=F(i)-F(i-1);
end
t1=2010:2016;
t2=2010:2019;
G
plot(t1,A,'o',t2,G)

Answer 1

问题关键的错误是 i=2:(n+1)中(n+1)，应该为(n+3)。因为你想预测2016到2019的数据，2019-2016=3，而不是1。

修改后运行可以得到如下图形。

追问

已知2010-2016年数据分别为：[11.501 13.595 13.957 10,721 8.571 10.813 9.467]
需要求得灰参数a和u的值；预测模型；累加序列预测值和原始序列预测值；c和p值；外推2016-2019年的预测值。还有残差修正以后进行的模型运算。总之就是灰色预测模型所需要的一切，能不能帮忙整一下，谢谢啊！