谁有GM(1,1) matlab程序
谁有GM(1,1) matlab程序，我非常的需要哪位兄台有啊，发到我的邮箱yysclshi@qq.com，非常的感谢啊 我有，加我QQ:841356436 我需要你把你的题目发给我，咱俩交换~ 我没有题目，只是看了往年的例题，看到这个东西
想要， 本帖最后由 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
% =gm(Q,t)
% by mcm-dlu-edu
function=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=; % 还原成 q
e=Q-q(1:length(Q)); % 算绝对误差
E=e./Q; % 算相对误差 我有,jiawo dddddddddd ddddddd
页:
[1]