查看: 3323|回复: 8

谁有GM(1,1) matlab程序

[复制链接]

字体大小: 正常放大

yysclshi

5 主题	4 听众	218 积分

升级 59%

该用户从未签到

电梯直达

1^#

发表于 2009-8-26 13:15 |只看该作者 |倒序浏览

|招呼Ta 关注Ta

谁有GM(1,1) matlab程序，我非常的需要
哪位兄台有啊，发到我的邮箱yysclshi@qq.com，非常的感谢啊

zan

转播0 淘帖0 分享0 收藏0 支持0 反对0 微信

使用道具举报

longshui2848

12 主题	4 听众	133 积分

升级 16.5%

该用户从未签到

2^#

发表于 2009-8-26 13:31 |只看该作者 |招呼Ta 关注Ta

我有，加我QQ:841356436

我也说一句

发表

使用道具举报

longshui2848

12 主题	4 听众	133 积分

升级 16.5%

该用户从未签到

3^#

发表于 2009-8-26 13:32 |只看该作者 |招呼Ta 关注Ta

我需要你把你的题目发给我，咱俩交换~

我也说一句

发表

使用道具举报

yysclshi

5 主题	4 听众	218 积分

升级 59%

该用户从未签到

4^#

发表于 2009-8-26 14:00 |只看该作者 |招呼Ta 关注Ta

我没有题目，只是看了往年的例题，看到这个东西
想要，

我也说一句

发表

使用道具举报

mcm-dlu-edu

98 主题	230 听众	1万积分

mcm

升级 0%

TA的每日心情

	怒 2015-11-22 11:00

签到天数: 679 天

[LV.9]以坛为家II

网络挑战赛参赛者

群组: 2015美赛优秀论文解析

群组: 第六届国赛赛前冲刺培

群组: 2014年地区赛数学建模

群组: 2013年美赛优秀论文解

群组: 2014美国中学生数学建

5^#

发表于 2009-8-26 14:03 |只看该作者 |招呼Ta 关注Ta

本帖最后由 mcm-dlu-edu 于 2009-8-26 14:06 编辑

% GM(1,1) si short for gray model
% t  is short for time;
% E  is short for contrary error
% e  is short for absolute error
% [q,e,E]=gm(Q,t)

% by mcm-dlu-edu

function[q,e,E]=GM(Q,t)
Q1=cumsum(Q);                      % 一阶累加
B1=Q1';B1(1)=[];B2=Q1';B2(end)=[];
B=[-0.5*(B1+B2) ones(length(B1),1)];  % 构造矩阵B
Xn=Q';Xn(1)=[];                      % 构造矩阵Xn
c=(B'*B)^-1*B'*Xn;a=c(1);u=c(2)    % 计算a 和 u
q1=(Q(1)-u/a)*exp(-a.*t)+u/a;       % 计算出 q1
q2=diff(q1);q=[q1(1) q2];          % 还原成 q
e=Q-q(1:length(Q));                % 算绝对误差
E=e./Q;                            % 算相对误差