数学建模社区-数学中国

标题: 求帮助…… [打印本页]

作者: 苏简若的爱 时间: 2012-8-4 15:00
标题: 求帮助……
%灰色系统做预测的matlab程序
clear;clc
syms a b;
c=[a b]';
A=[174 179 183 189 207 234 220.5 256 270 385];%数据行，竞赛时直接替换数据就行
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+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=1995:2004;
t2=1995:2014;
G;a,b
plot(t1,A,'o',t2,G)

??? Error using ==> mtimes
Inner matrix dimensions must agree.

Error in ==> huisexitong at 15
c=inv(E*E')*E*D;
要求矩阵内部维数一样，但是不怎么理解

作者: 苏简若的爱 时间: 2012-8-5 11:11
求帮助呀……非常感谢。

作者: darker50 时间: 2012-8-5 12:16
那个代码用编辑器的上面的“<>”这个符号就能使用让代码不出现表情的了。

作者: 苏简若的爱 时间: 2012-8-5 14:28

darker50 发表于 2012-8-5 12:16
那个代码用编辑器的上面的“”这个符号就能使用让代码不出现表情的了。

恩恩，谢谢！但是现在不是你说的那个问题，我本身编辑的时候上面没有错，只是最后出现
??? Error using ==> mtimes
Inner matrix dimensions must agree.

Error in ==> huisexitong at 15
c=inv(E*E')*E*D;

我不知道究竟出现在哪里
麻烦给解答一下。

作者: 我能 时间: 2012-8-12 15:53
有很多的问题要自己去看比如你的算除法时直接用 E\d 你那有写问题就是线代中的矩阵除法你一步步试下
我发了个我做深圳夏令营的题的程序你看下

附件1 深圳历年人口数据.xls

19.5 KB, 下载次数: 0, 下载积分: 体力 -2 点

HUSE(chang).M

1.01 KB, 下载次数: 0, 下载积分: 体力 -2 点

售价: 2 点体力 [记录] [购买]

欢迎光临数学建模社区-数学中国 (http://www.madio.net/)