Errorusing==>mtimesInnermatrixdimensionsmustagree.clearsymsau;c=[au]';%灰参数cA=[489.9518.6507.9519.9527.4520.3521.82515.9514.1501.2];%原始数据Ago=cumsum(A);%原始数据一次累加n=length(A);%原始数据个数fori=1:(

Question

Error using ==> mtimes Inner matrix dimensions must agree.
clear
syms a u;
c=[a u]'; %灰参数c
A=[489.9
518.6
507.9
519.9
527.4
520.3
521.82
515.9
514.1
501.2
]; %原始数据
Ago=cumsum(A); %原始数据一次累加
n=length(A); %原始数据个数
for i=1:(n-1)
C(i)=(Ago(i)+Ago(i+1))/2; %生成累加矩阵
end % 计算待定参数的值
Yn=A;
Yn(1)=[];
Yn=Yn';
E=[-C;ones(1,n-1)];
c=inv(E*E')*E*Yn;
c=c';
a=c(1);u=c(2); % 预测后续数据
F=[];F(1)=A(1);
for i=2:(n+10)
F(i)=(A(1)-u/a)/exp(a*(i-1))+u/a ;
end
G=[];
G(1)=A(1);
for i=2:(n+10)
G(i)=F(i)-F(i-1); %得到预测出来的数据
end
t1=2002:2011;
t2=2002:2018;
G;a;u; %输出预测值、发展系数和灰色作用量
plot(t1,A,'o',t2,G) %原始数据与预测数据的比较

Answer 1

c=inv(E*E')*E*Yn;
这一步 E是2*9的矩阵,Yn是1*9的向量.矩阵运算维数不匹配,请检查E和Yn是否设错.