请教一论文中matlab编程问题

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

小弟最近在看《GM（1,1）模型灰色作用量的优化》这篇论文，对其中的计算用matlab，代码如下，(鄙人水平极其有限，对该软件只会简单操作）求得a,b,c与文中结果相同，但最终运行结果差别巨大，很久也没发现错在哪，今特意注册账号恳请各位前辈指点纠正，感激不尽：
clear;clc;
y1=[5425.1 5854 6280 6859.6 7702.8 8472.2];
n=6;
for i=1:n
x1(i)=sum(y1(1:i));
end
y=[y1(2:n)]';
for i=2:n
    x(i)=0.5*x1(i-1)+0.5*x1(i);
end
x=x(2:n);
x=[-x;ones(1,n-1);2:n]';
A=x\y;
a=A(1);
b=A(2);
c=A(3);
for i=1:n+2
    xy(i+1)=(1-exp(a))*(y1(1)-b/a-c/a+c/(a^2))*exp(-i*a)+c/a;%问题应该在这，这里我加上了c/a，论文中没有，核算过一次，好像文中式子漏掉了这个，也可能我算错了，但不管加与否，都对不上
end
xy(1)=y1(1);
xy

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：MATLAB编程 MATLAB matla atlab Lab matlab 论文 模型 软件

好几天了，有木有人能给点有用的信息呀。。。万分期待中

代码没仔细看，不过我这里有段灰色预测的现成代码，供你参考。
clear
syms a b;
c=[a b]';
A=[89677,99215,109655,120333,135823,159878,182321,209407,246619,300670];
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=1999:2008;
t2=1999:2018;
G
plot(t1,A,'o',t2,G)  %原始数据与预测数据的比较
xlabel('年份')
ylabel('利润')

你这是传统灰色模型的代码。。。不过还是谢谢你的回复