[问题求助] 请问各位高手，下面的这个代码哪些地方错了？怎么改啊

[复制链接]

字体大小: 正常放大

Mr.Q

7 主题	4 听众	950 积分

升级 87.5%

TA的每日心情

	开心 2018-4-19 21:34

签到天数: 111 天

[LV.6]常住居民II

2012挑战赛参赛者

群组: 第二届数模基础实训

群组: 数学建模培训课堂1

群组: 2011年第一期数学建模

群组: 数学建模培训课堂2

电梯直达

1^#

发表于 2012-4-14 10:51 |只看该作者 |倒序浏览

|招呼Ta 关注Ta |邮箱已经成功绑定

本帖最后由 Mr.Q 于 2012-4-14 10:52 编辑

clear
syms a b;
c=[a b]';
A=[26.2
23.2
23.4
23.3
23.2
23.2
24.0
25.2
24.6
25.2
26.3
25.9
21.9
22.4
20.1
20.1
22.5
20.4
20.7
19.9
19.4
19.0
19.0
20.4
20.9
21.1
20.4
20.8
18.7
20.2
];
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+30)
F(i)=(A(1)-b/a)/exp(a*(i-1))+b/a ;
end
G=[];G(1)=A(1);
for i=2: (n+30)
G(i)=F(i)-F(i-1); %得到预测出来的数据
end
t1=1:30;
t2=1:61;
G, a, b % 输出预测值，发展系数和灰色作用量
plot(t1,A,'o',t2,G)  %原始数据与预测数据的比较

zan

数据