disp('可以用GM(1,1)建模') end x1 = cumsum(x0) %累加运算 B=[-0.5*(x1(1:n-1)+x1(2:n)),ones(n-1,1)] Y=x0(2:n); u=inv(B'*B)*B'*Y %拟合参数 %上述语句可替换为 u=B\Y 【有点不太理解】 % k=10; forecast1 = (x1(1)-u(2)./u(1)).*exp(-u(1).*([0:n-...
p=1;c=0.0148;预测等级为:好 从运行结果看,对于线性的数据使用GM(1,1)预测,其拟合效果还是不错。 【2】GM(1,n)模型及Matlab实现 1.GM(1,n)模型的预测原理:与GM(1,1)类似,不同在于输入数据变量是n个。 2. GM(1,n)模型的建模过程: 设系统有特征数据序列: 相关因素序列: (1) 令 的1-AGO序列为...
(1)原始值-模型值-残差-相对误差0级比偏差表格 (2)根据表格作出说明 一个实例: 北方某城市 1986～1992 年道路交通噪声平均声级数据如下: 求解分析: 求解MATLAB程序如下: %此程序原著pdf上程序是有bug的,以下已经调通 1clc,clear2x0=[71.172.472.472.171.472.071.6]';%注意这里为列向量3n=length(x0);4lamda...
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进行准光滑性检验，序列x0的光滑比P(t)=X0(t)/X1(t-1)A=inv(B.'*B)*B.'*YN.' ;a=A(1)u=A...
灰色模型预测GM（1,1）MATLAB程序代码 灰⾊模型预测GM（1,1）MATLAB程序代码 版权所有引⽤请注明出处 function gmcal=gm1(x)%% ⼆次拟合预测GM(1,1)模型 %x = [5999,5903,5848,5700,7884];sizexd2 = size(x,2);%求数组长度 k=0;for y1=x k=k+1;if k>1 x1(k)=x1(k-1)+x(k);%累加...
灰色预测;GM ( 1,1) 模型;M atl ab;算法;精度检验中图分类号:N 941.5文献标识码:A关键词:灰色系统理论是20世纪80年代由我国邓聚龙教授提出的一种数学方法,可用来解决小样本的信息不完备系统的复杂问题,在岩土工程中得到了广泛关注¨ q] 。如果有一个程序能够使工程人员只需输入数据,而不需编写计算公式,就...
灰色预测[GM(1 ,1 ) ]MATLAB 程序 % 本程序主要用来计算根据灰色理论建立的模型的预测值。 % 应用的数学模型是 GM(1 ,1 )。 % 原始数据的处理方法是一次累加法。 y=input('请输入数据 ');%输入数据请用如例所示形式: [48.7 57.1 7 68.76 92.1 5] n=length(y); yy=ones(n,1 ); yy(1 )=y(...
三、matlab程序 // 灰色预测g(1,1)模型,负荷预测clcclose allclear all%% 原始数据%年份。。。负荷S=[2000 4149 1149 2003.07 4851 2396 208.152001 4186 1273 2175.68 5221 2500 222.282002 4222 1360 2450.48 5829 2651 246.572003 4254 1447 2807.41 6624 2739 299.532004 4284 1524 3456.7 8097 3353 389.22005 ...
1.6837 2.0485 2.4435 2.8581 3.3117 (3)由于GM(1,1)的灰微分方程为X (t)-+aZ (k)=b 设 为待估参数, ,利用最小二乘法得到 , matlab程序为 >>B=[-z1(2:12)',ones(11,1)]; >> Y=(x0(2:12))'; >> alpha=inv(B'*B)*B'*Y; >> alpha alpha = 1.0e+04 * -0.0000 1.4000©...
2014年4月JiangsuScience&TechnologyInformationApril,2014改进的灰色预测GM(1,1)模型的MATLAB实现杨旭(郑州大学水利与环境学院,河南郑州450001)摘要:灰色系统理论中的灰色预测理论已得到了广泛的应用,文章简单介绍了改进的灰色预测GM(1,1)模型,使用MATLAB语言给出了建立模型的算法程序,为高效地利用MATLAB强大的科学计算功能...