mshd.net
当前位置:首页 >> mAtlAB 灰色预测 GM(1,1) >>

mAtlAB 灰色预测 GM(1,1)

x0=[1.085 1.189 1.375 1.566 1.774 1.957 2.106 2.347]; t=1999:2006; %年份 GM10(x0,t) 请输入未来几年数T:3 ans = 1...

clc,clear x0=[692.5 649.6 669.8 596.0 607.9 617.9 635.2 652.2 667.8 682.8 675.1 695.7]; n=length(x0); lamda=x0(1:n-1)./x0(2:n) range=minmax(lamda) x1=cumsum(x0) for i=2:n z(i)=0.5*(x1(i)+x1(i-1)); end B=[-z(2:n)',ones(n-1,1)]; ...

function GM1_1(X0) %format long ; [m,n]=size(X0); X1=cumsum(X0); %累加 X2=[]; for i=1:n-1 X2(i,:)=X1(i)+X1(i+1); end B=-0.5.*X2 ; t=ones(n-1,1); B=[B,t] ; % 求B矩阵 YN=X0(2:end) ; P_t=YN./X1(1:(length(X0)-1)) %对原始数据序列X0...

%%%%%%%%%%%%%%%%把下面函数保存为gmcal.m文件%%%%%%%%%%% function gmcal=gm1(x) sizexd2 = size(x,2); %求数组长度 k=0; for y1=x k=k+1; if k>1 x1(k)=x1(k-1)+x(k); %累加生成 z1(k-1)=-0.5*(x1(k)+x1(k-1)); %z1维数减1,用于计算B yn1(k-1...

定义函数: function pre = GM11(x) %x,建模数据序列x = [x(1),x(2),...,x(n)] % pre,由GM11模型得到的预测值,即由x估计出的x(n+1) n = length(x); x1 = cumsum(x); z1 = (x1(1:n-1)+x1(2:n))/2; Y = x(2:length(x))'; B = [-z1',ones(n-1,1)];...

第一组数据,预测2010-2015年的值 2010 0.0495265145611443 2011 0.0500349983509087 2012 0.050548702693072 2013 0.0510676811865275 2014 0.0515919879804586 2015 0.0521216777799994 相对误差 0 2.00323906921688 2.89291599461445 2.69245...

结果:第 2013 年估值: 230109第 2014 年估值: 242828第 2015 年估值: 256250第 2016 年估值: 270414第 2017 年估值: 285361代码:% Purpose:% 解答网友 岑风吟 提出的问题:“求大神用matlab软件做GM(1.1)灰色预测”% Record of revision:...

y=[6134 6776 8151 9536.07 15911.96 18873.2 29286.8 33720.1 39131]; y1=y'; for i=1:9 j=1:i; y2(i)=sum(y1(j)); end for i=2:9 y3(i-1)=-0.5*(y2(i-1)+y2(i)); end b=ones(8,1); b=[y3',b]; yy=y1(2:end)'; a=(b'*b)^(-1)*b'*yy; yy2(1)=y1(...

第一个貌似是leslie矩阵 第二个才是你说的灰色模型 第二个的用法也很简单 把整个函数扔到M文件中 在主程序的command window中调用即可

把数据改一下应该就行了 clc clear x0=[1827 1729.3 1431 1231.7 1109.9 768.2 792.9 805.8]; %道路交通事故财产损失原始数据 for i=2:8 x1(1)=x0(1); x1(i)=x1(i-1)+x0(i); end x1 %对原始数据进行累加 for i=1:7 yn(i)=x0(i+1); end yn for i=...

网站首页 | 网站地图
All rights reserved Powered by www.mshd.net
copyright ©right 2010-2021。
内容来自网络,如有侵犯请联系客服。zhit325@qq.com