灰色预测MATLAB
灰色系统预测程序function yuce(x) x=[ 36.22 37.66 39.09 40.53 41.76 42.99 43.90 44.94 45.69 46.59 50.00 ]; format long g n=length(x); x1=zeros(1,n); x1(1)=x(1); for i=2:n
x1(i)=x(i)+x1(i-1); end
%对序列x进行累加生成
z1=zeros(1,n-1); for i=2:n
z1(i-1)=(x1(i)+x1(i-1))/2; end
%对序列x1进行紧邻均值生成
z2=z1'; z3=z2.*(-1); B=; y=zeros(n-1,1); for i=1:n-1
y(i)=x(i+1); end
%构造b矩阵和y矩阵
au=inv(B'*B)*B'*y; %最小二乘法的参数估计
x,x1,z1,B,y,au(1),au(2)%输出原始序列、一次累加生成序列、紧邻均值生成、B矩阵、y矩阵、参数估计
a=x1(1)-au(2)/au(1); b=au(2)/au(1);
X=zeros(1,n); X1=zeros(1,n); for i=1:n
X1(i)=a*exp(-au(1)*(i-1))+b; end X(1)=X1(1); for i=2:n
X(i)=X1(i)-X1(i-1); end X
%计算模拟序列X dt=x-X; v1=sum(x)/n; v2=sum(dt)/n; s1=0; s2=0; for i=1:n
s1=s1+(x(i)-v1)^2;
s2=s2+(dt(i)-v2)^2;
2014全国一级建造师资格考试备考资料真题集锦建筑工程经济 建筑工程项目管理 建筑工程法规 专业工程管理与实务
end s1=s1/n; s2=s2/(n-1); C=(sqrt(s2))/(sqrt(s1)); C %计算后验差比
m=0; for i=1:n
if abs(dt(i)-v2)<0.6745*(sqrt(s1))
m=m+1;
end end p=m/n; p%计算小误差概率
if p>0.95&C<0.35
disp('预测精度好'); else if p>0.8&C<0.5
disp('预测合格');
else if p>0.7&C<0.65
disp('预测勉强合格');
else disp('预测不合格');
end
end end for i=1:n+3
yc1(i)=a*exp(-au(1)*(i-1))+b; end yc(1)=yc1(1); for i=2:n+3
yc(i)=yc1(i)-yc1(i-1); end yc
%预测三年数据
求一元回归方程和相关系数
y=; x=[ 4984.05 6344.11 7790.92 10153.80 13158.30 15909.20 19422.90 25288.80 31203.20 36241.80 48267.00 ]; plot(x,y,'o') xlabel('金融政策（房产投资金额）'); ylabel('商品房平均销售价格'); hold on p=polyfit(x,y,1) y1=polyval(p,x);
plot(x,y1) a=corrcoef(x,y)
灰色关系程序
Y=;%商品房平均销售价格
x1=109655.2 120332.7 135822.8 159878.3 184937.4 216314.4 265810.3 314045.4 340506.9 397983.0
];%国内生产总值
x2=;%城市化率
x3=;%经济使用房投资额
x4=;%物价水平（建材价格）
x5=;%居民用地价格指数
x6=;%税收政策（房地产税收）
x7=6344.11 7790.92 10153.80 13158.30 15909.20 19422.90 25288.80 31203.20 36241.80 48267.00 ];%金融政策（房地产投资额）
Y_1=Y(:,1)./Y; Y_2=Y./(sum(Y)/length(Y)); x1_1=x1(:,1)./x1; x2_1=x2(:,1)./x2; x3_1=x3(:,1)./x3; x4_1=x4./(sum(x4)/11); x5_1=x5./(sum(x5)/11); x6_1=x6(:,1)./x6; x7_1=x7(:,1)./x7; X_1=
x3_1' x4_1'
x5_1' x6_1' x7_1' ];%初值化矩阵
X_2=[ x1_1' x2_1'
x3_1' x4_1'
x5_1' x6_1' x7_1' ]; X_3=; X_4=abs(X_2-X_3);%绝对差矩阵
p=0.7;%分辨系数
mn=min(min(X_4)); mx=max(max(X_4)); X_5=(mn+p*mx)./(X_4+p*mx);%各影响因素关联系数
=size(X_5); for i=1:b
X_6(i)=sum(X_5(:,i))./a;%各影响因素关联度
end X_6
多元回归方程和残差分析
x1=
5023.8 6060.3 6886.3 9846.8 11115.0 12153.0 14113.6 ];%国内生产总值
x2=
79.05 79.22 83.62
84.50 84.90 85.00 85.96];%城市化率
x3=
69.07 72.57 44.76
28.30 35.95 67.76 48.98];%经济使用房投资额
x4=
105.8 101.8 99.5
103.5 115.2 99.4 102.7 ];%物价水平（建材价格）
x5=
101.0 102.0 103.7 105.3 104.6 104.5 107.8 ];%居民用地价格指数
x6=
30.79 31.97 35.21
51.75 63.84 73.98 83.83];%税收政策（房地产税收）
x7=;%金融政策（房地产投资额）
X=; Y=
3864 3800 4681 5032];%商品房平均销售价格
alpha=0.05; =regress(Y',X,alpha) y1=X*b; e=Y'-y1; plot(y1,e,'*');%画出以y^i为横坐标的残差图
figure; rcoplot(r,rint);%画出以观测值序号为横坐标的残差图
{:3_41:}{:3_41:}{:3_41:}{:3_41:}{:3_41:}{:3_41:}{:3_41:} 表示深深看不懂 表示深深看不懂
页:
[1]