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...

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...

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)]; ...

%%%%%%%%%%%%%%%%把下面函数保存为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)];...

结果:第 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(...

用bintprog求解,可通过help bintprog先熟悉其用法. 由于bintprog的解为向量形式,因此,若未知量x为矩阵形式,要先将x转换为向量形式进行求解,结果再转回矩阵形式. 为方便说明,下面的系数向量s由4*3的矩阵构造,你的原始问题只须将s换为20*...

把数据改一下应该就行了 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=...

t2 = 2014 x0_pre2 = 16518.7590276555 复制去Google翻译翻译结果

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