没有合适的资源?快使用搜索试试~ 我知道了~
是运行过的MATLAB灰色模型,十分有用
资源详情
资源评论
资源推荐
function aa
X0=[2068664 2188917 2200947 2298368 2523943 2820911 2857872 3103836
1851000 3155000];
AU=aa(X0);a=AU(1);u=AU(2);m2=length(X0);
for k=1:1:m2-1
xx1(k+1)=(X0(1)-u/a)*exp(-a*k)+u/a;
end
s=0;xx0(1)=X0(1);
for jj=2:1:m2;
xx0(jj)=xx1(jj)-xx1(jj-1);
end
disp('GM(1,1)¶ÔÊýÁнøÐÐÔ¤²â½á¹û');xx0
disp('ÊýÁÐ1Ôʼ¹Û²âÊý¾Ý');X0
disp('a');AU(1)
disp('u');AU(2)
function au=aa(X0)
m=length(X0);
s1=0;
for jj=1:1:m;X1(jj)=s1+X0(jj); s1=X1(jj); end
for ii=1:1:m-1; B(ii)=-(X1(ii)+X1(ii+1))/2; end
B=[B(:),ones(m-1,1)];
y=X0([2:m])';
au=inv((B'*B))*B'*y;
%%灰色预测
function [feval,au,ec,ex,r,C,P]=GM1_1(x,myear)
%%
%%
%%
%%
if nargin<2
myear=0;
end
[mx,nx]=size(x);
if mx==1
x=x';
end
SUISUI61
- 粉丝: 0
- 资源: 1
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功
评论0