灰色预测MATLAB

下沙小僧

[p=130, null, left][size=130px]灰色系统预测程序

[p=145, null, left][size=130px]function yuce(x)

[p=145, null, left][size=130px]x=[ 36.22 37.66 39.09 40.53 41.76 42.99 43.90 44.94 45.69 46.59 50.00 ];

[p=145, null, left][size=130px]format long g

[p=145, null, left][size=130px]n=length(x);

[p=145, null, left][size=130px]x1=zeros(1,n);

[p=145, null, left][size=130px]x1(1)=x(1);

[p=145, null, left][size=130px]for i=2:n

* y  m9 k  W% \

& ~# O, \2 h; t& ~  u
[p=145, null, left][size=130px]x1(i)=x(i)+x1(i-1);

[p=145, null, left][size=130px]end

! e! O6 e5 i/ ^[p=145, null, left][size=130px]%

[p=145, null, left][size=130px]对序列

[p=145, null, left][size=130px]x

[p=145, null, left][size=130px]进行累加生成

[p=145, null, left][size=130px]z1=zeros(1,n-1);

[p=145, null, left][size=130px]for i=2:n

  R" {" Y- _- J8 J. n: ]  e1 c1 Y5 D[p=145, null, left][size=130px]z1(i-1)=(x1(i)+x1(i-1))/2;

[p=145, null, left][size=130px]end

[p=145, null, left][size=130px]%

[p=145, null, left][size=130px]对序列

[p=145, null, left][size=130px]x1

[p=145, null, left][size=130px]进行紧邻均值生成

, G' l& e( d0 O* t  n; p* `6 C[p=145, null, left][size=130px]z2=z1';

[p=145, null, left][size=130px]z3=z2.*(-1);

[p=145, null, left][size=130px]B=[z3,ones(n-1,1)];

[p=145, null, left][size=130px]y=zeros(n-1,1);

[p=145, null, left][size=130px]for i=1:n-1

3 _+ |' Y' j9 U/ I3 A8 H
/ A/ D7 L' Z- W7 n5 ]0 a
! w7 ^, w% P) H) o
% A# e  z% l) c4 }, {! p[p=145, null, left][size=130px]y(i)=x(i+1);

[p=145, null, left][size=130px]end

[p=145, null, left][size=130px]%

[p=145, null, left][size=130px]构造

[p=145, null, left][size=130px]b

[p=145, null, left][size=130px]矩阵和

[p=145, null, left][size=130px]y

[p=145, null, left][size=130px]矩阵

[p=145, null, left][size=130px]au=inv(B'*B)*B'*y; %

[p=145, null, left][size=130px]最小二乘法的参数估计

$ a, ~3 q0 C) R' b" C/ J  u% a[p=145, null, left][size=130px]x,x1,z1,B,y,au(1),au(2)%

[p=145, null, left][size=130px]输出原始序列、

[p=145, null, left][size=130px]一次累加生成序列、

[p=145, null, left][size=130px]紧邻均值生成、

[p=145, null, left][size=130px]B

[p=145, null, left][size=130px]矩阵、

[p=145, null, left][size=130px]y

[p=145, null, left][size=130px]矩阵、

[p=130, null, left][size=130px]参数估计

[p=145, null, left][size=130px]a=x1(1)-au(2)/au(1);

[p=145, null, left][size=130px]b=au(2)/au(1);

3 C7 m* u. F* H: X
1 u6 S$ y9 k' z, F3 ?[p=145, null, left][size=130px]X=zeros(1,n);

[p=145, null, left][size=130px]X1=zeros(1,n);

[p=145, null, left][size=130px]for i=1:n

, ^: y: g9 z* @. O

[p=145, null, left][size=130px]X1(i)=a*exp(-au(1)*(i-1))+b;

[p=145, null, left][size=130px]end

[p=145, null, left][size=130px]X(1)=X1(1);

[p=145, null, left][size=130px]for i=2:n

- s; {" Y0 T, W' d2 J
1 l; W) ?: `2 r0 O  o8 o2 Q
) v7 ]/ I# l( q, O$ I! `9 H
[p=145, null, left][size=130px]X(i)=X1(i)-X1(i-1);

[p=145, null, left][size=130px]end

[p=145, null, left][size=130px]X

[p=145, null, left][size=130px]%

[p=145, null, left][size=130px]计算模拟序列

[p=145, null, left][size=130px]X

[p=145, null, left][size=130px]dt=x-X;

[p=145, null, left][size=130px]v1=sum(x)/n;

[p=145, null, left][size=130px]v2=sum(dt)/n;

[p=145, null, left][size=130px]s1=0;

[p=145, null, left][size=130px]s2=0;

[p=145, null, left][size=130px]for i=1:n

7 H( K8 i& Z/ o/ @9 W; l! n
0 d! Z0 @; D# A9 T/ d" d

- G5 E- i4 q( f8 L% W[p=145, null, left][size=130px]s1=s1+(x(i)-v1)^2;

1 W# i, x% X1 \$ d- j  e8 ~, q: V

* v! z" r8 A* W, g% D; R
[p=145, null, left][size=130px]s2=s2+(dt(i)-v2)^2;

9 y9 W/ a- H+ E  F; `
6 e/ u2 G, e2 x. Z: d/ i* L' A9 M/ X+ o

& m2 \5 C. x, |2 J$ N3 A

6 P" l# G$ [- R

2014全国一级建造师资格考试备考资料真题集锦建筑工程经济 建筑工程项目管理 建筑工程法规 专业工程管理与实务
) F1 k  U; Y3 ?6 r* |

[p=145, null, left][size=130px]end

[p=145, null, left][size=130px]s1=s1/n;

[p=145, null, left][size=130px]s2=s2/(n-1);

[p=145, null, left][size=130px]C=(sqrt(s2))/(sqrt(s1));

[p=145, null, left][size=130px]C %

[p=145, null, left][size=130px]计算后验差比

/ b0 F& t) @2 O5 u% e[p=145, null, left][size=130px]m=0;

[p=145, null, left][size=130px]for i=1:n

7 W6 l  B$ U" C$ E) r
* d6 N6 f+ I" B/ V" s* N

: i  Z% p( f2 r# x. X7 l. \0 b[p=145, null, left][size=130px]if abs(dt(i)-v2)<0.6745*(sqrt(s1))

+ v8 a- x+ ?$ h9 c  `. I
/ s$ h& ^( {: {, e6 S, Y


0 q# o3 M8 P6 h* L
) d. {: I6 q1 O2 K6 m

0 z# r5 o6 Y# ^9 W; H[p=145, null, left][size=130px]m=m+1;

3 V4 s7 t. L3 D# N: |$ E. z& u% \+ s
& S* {) `1 p% }2 E6 G* ~9 O
[p=145, null, left][size=130px]end

[p=145, null, left][size=130px]end

[p=145, null, left][size=130px]p=m/n;

[p=145, null, left][size=130px]p%

[p=145, null, left][size=130px]计算小误差概率

4 x. \& a; J3 c1 T6 S# M( @0 E$ x[p=145, null, left][size=130px]if p>0.95&C<0.35

& G' f$ }: B$ U8 ~5 O

5 F; z0 y' A1 r: E( K. ?7 d6 ?[p=145, null, left][size=130px]disp('

[p=145, null, left][size=130px]预测精度好

[p=145, null, left][size=130px]');

[p=145, null, left][size=130px]else if p>0.8&C<0.5

1 q9 h* l( H, g7 m) Q3 L
" V. s# c3 v0 W6 t6 V* |( l# Q





[p=145, null, left][size=130px]disp('

[p=145, null, left][size=130px]预测合格

[p=145, null, left][size=130px]');

[p=145, null, left][size=130px]else if p>0.7&C<0.65

0 b8 ~: }8 O+ a- j$ g
# @- o7 L% z/ V  d" S: G- d
+ d  H. c' _7 J; i8 E/ T. R6 q) A
1 b1 F5 p) F; l% q) P
6 V! i2 t$ H2 i4 \( |; J9 d


( I: g" Z( I- o: Z

& K+ y. L  Q# h; n2 Z# X+ G
  Z  n. I0 z# x7 O" I3 R' w, L
; o+ B8 e7 W! W( i* p& `' N; d, M[p=145, null, left][size=130px]disp('

[p=145, null, left][size=130px]预测勉强合格

[p=145, null, left][size=130px]');

1 O  y8 T# f/ h# S. ]4 l9 `2 {5 l! ^
+ E0 @2 b5 _) K2 P

$ A1 y6 ^, v. t
+ f6 y* M7 ]" w) j

[p=145, null, left][size=130px]else disp('

[p=145, null, left][size=130px]预测不合格

[p=145, null, left][size=130px]');

7 g) N, l& y- F/ G7 i! U
# s% h  ]0 [7 L8 j5 v

/ _2 }: b5 V0 n+ u* \
2 X, k# j8 K6 I/ }* U
5 @" ~( Y  l) p# L2 }

% G2 w7 K. N6 b! B: a* K

: A. ^5 W+ D0 o' i8 V
5 M' d/ b7 u9 I; A[p=145, null, left][size=130px]end

4 ?' W( \: j% C( ]4 G

  j% g1 K5 q1 G3 A
% W) ^6 S# T0 c/ N9 b* P9 R[p=145, null, left][size=130px]end

[p=145, null, left][size=130px]end

[p=145, null, left][size=130px]for i=1:n+3

+ A  C. X" q2 c

/ z- P8 d+ @: F- _/ r/ Y
[p=145, null, left][size=130px]yc1(i)=a*exp(-au(1)*(i-1))+b;

[p=145, null, left][size=130px]end

[p=145, null, left][size=130px]yc(1)=yc1(1);

[p=145, null, left][size=130px]for i=2:n+3

[p=145, null, left][size=130px]yc(i)=yc1(i)-yc1(i-1);

[p=145, null, left][size=130px]end

[p=145, null, left][size=130px]yc

8 q! O4 {- h- ^' [' u7 R/ X& Y
[p=145, null, left][size=130px]%

[p=145, null, left][size=130px]预测三年数据

% r: z3 Q+ \* @5 S- B
8 j$ I8 w; U7 H3 o# t. x[p=130, null, left][size=130px]求一元回归方程和相关系数

[p=145, null, left][size=130px]y=[2112 2170 2250 2359 2778 3168 3367 3864 3800 4681 5032];

[p=145, null, left][size=130px]x=[

[p=145, null, left][size=130px]4984.05

[p=145, null, left][size=130px]6344.11

[p=145, null, left][size=130px]7790.92

[p=145, null, left][size=130px]10153.80

[p=145, null, left][size=130px]13158.30

[p=145, null, left][size=130px]15909.20

[p=145, null, left][size=130px]19422.90

[p=145, null, left][size=130px]25288.80

[p=145, null, left][size=130px]31203.20

[p=145, null, left][size=130px]36241.80 48267.00 ];

[p=145, null, left][size=130px]plot(x,y,'o')

[p=145, null, left][size=130px]xlabel('

[p=145, null, left][size=130px]金融政策（房产投资金额）

[p=145, null, left][size=130px]');

[p=145, null, left][size=130px]ylabel('

[p=145, null, left][size=130px]商品房平均销售价格

[p=145, null, left][size=130px]');

[p=145, null, left][size=130px]hold on

[p=145, null, left][size=130px]p=polyfit(x,y,1)

[p=145, null, left][size=130px]y1=polyval(p,x);

3 a1 ^8 `4 u* g( g9 A4 g+ }1 s


0 w( t! T" d6 p5 M1 q& W" ?



[p=145, null, left][size=130px]plot(x,y1)

[p=145, null, left][size=130px]a=corrcoef(x,y)

0 T+ |5 b# k& N
[p=130, null, left][size=130px]灰色关系程序

/ s" ~5 |1 h2 x  E[p=145, null, left][size=130px]Y=[2112 2170 2250 2359 2778 3168 3367 3864 3800 4681 5032];%

[p=145, null, left][size=130px]商品房平均销售价格

[p=145, null, left][size=130px]x1=[99214.6

[p=145, null, left][size=130px]109655.2

[p=145, null, left][size=130px]120332.7

[p=145, null, left][size=130px]135822.8

[p=145, null, left][size=130px]159878.3

[p=145, null, left][size=130px]184937.4

[p=145, null, left][size=130px]216314.4

[p=145, null, left][size=130px]265810.3

[p=145, null, left][size=130px]314045.4

[p=145, null, left][size=130px]340506.9 397983.0

% S3 ^7 u( B+ y4 i[p=145, null, left][size=130px]];%

[p=145, null, left][size=130px]国内生产总值

, K" z3 I9 k2 H3 X[p=145, null, left][size=130px]x2=[36.22 37.66 39.09 40.53 41.76 42.99 43.90 44.94 45.69 46.59 50.00 ];%

[p=145, null, left][size=130px]城市化率

[p=145, null, left][size=130px]x3=[542.44 599.65 589.0 622.0 606.4 519.2 696.8 820.9 970.9 1134.1 1349.9];%

[p=145, null, left][size=130px]经济使用房投资

[p=130, null, left][size=130px]额

[p=145, null, left][size=130px]x4=[101.50 98.60 98.20 99.70 105.10 103.10 101.90 103.00 109.50 101.10 102.30 ];%

[p=145, null, left][size=130px]物价水平

[p=130, null, left][size=130px]（建材价格）

[p=145, null, left][size=130px]x5=[101.00 98.60 107.70 112.40 111.60 110.30 106.00 113.70 111.30 106.40 127.60];%

[p=145, null, left][size=130px]居民用地

[p=130, null, left][size=130px]价格指数

  ~) S  @; V  L  [3 ^[p=145, null, left][size=130px]x6=[209.60 228.60 282.40 323.90 366.30 435.90 515.20 575.10 680.34 803.66 894.06 ];%

[p=145, null, left][size=130px]税收政

[p=130, null, left][size=130px]策（房地产税收）

+ H. P7 [7 B1 ~[p=145, null, left][size=130px]x7=[4984.05

[p=145, null, left][size=130px]6344.11

[p=145, null, left][size=130px]7790.92

[p=145, null, left][size=130px]10153.80

[p=145, null, left][size=130px]13158.30

[p=145, null, left][size=130px]15909.20

[p=145, null, left][size=130px]19422.90

[p=145, null, left][size=130px]25288.80

[p=145, null, left][size=130px]31203.20

[p=145, null, left][size=130px]36241.80 48267.00 ];%

[p=145, null, left][size=130px]金融政策（房地产投资额）

2 `( {1 d$ {/ x& ]1 F6 w[p=145, null, left][size=130px]Y_1=Y(:,1)./Y;

[p=145, null, left][size=130px]Y_2=Y./(sum(Y)/length(Y));

[p=145, null, left][size=130px]x1_1=x1(:,1)./x1;

[p=145, null, left][size=130px]x2_1=x2(:,1)./x2;

[p=145, null, left][size=130px]x3_1=x3(:,1)./x3;

[p=145, null, left][size=130px]x4_1=x4./(sum(x4)/11);

[p=145, null, left][size=130px]x5_1=x5./(sum(x5)/11);

[p=145, null, left][size=130px]x6_1=x6(:,1)./x6;

[p=145, null, left][size=130px]x7_1=x7(:,1)./x7;

[p=145, null, left][size=130px]X_1=[Y_1' x1_1' x2_1'

% j1 f+ K! t/ V# g: X) H" W6 G[p=145, null, left][size=130px]x3_1' x4_1'

! v. i% n7 ]& D[p=145, null, left][size=130px]x5_1' x6_1' x7_1' ];%

[p=145, null, left][size=130px]初值化矩阵

- i1 Q9 l# `, P. q+ y4 j; g! A[p=145, null, left][size=130px]X_2=[ x1_1' x2_1'

[p=145, null, left][size=130px]x3_1' x4_1'

  D7 B6 q0 Z- ^[p=145, null, left][size=130px]x5_1' x6_1' x7_1' ];

[p=145, null, left][size=130px]X_3=[Y_1' Y_1' Y_1' Y_2' Y_2' Y_1' Y_1' ];

[p=145, null, left][size=130px]X_4=abs(X_2-X_3);%

[p=145, null, left][size=130px]绝对差矩阵

[p=145, null, left][size=130px]p=0.7;%

[p=145, null, left][size=130px]分辨系数

7 i( G+ R1 `2 G# A9 k[p=145, null, left][size=130px]mn=min(min(X_4));

[p=145, null, left][size=130px]mx=max(max(X_4));

[p=145, null, left][size=130px]X_5=(mn+p*mx)./(X_4+p*mx);%

[p=145, null, left][size=130px]各影响因素关联系数

[p=145, null, left][size=130px][a,b]=size(X_5);

[p=145, null, left][size=130px]for i=1:b

* S. s6 Q. N+ B# ?3 o* b9 [


[p=145, null, left][size=130px]X_6(i)=sum(X_5(:,i))./a;%

[p=145, null, left][size=130px]各影响因素关联度

) C- w; j# _% ?$ {9 y" b[p=145, null, left][size=130px]end

[p=145, null, left][size=130px]X_6

8 e  L5 y0 J1 o! F$ X
) _* R, Q1 _! Z  z% l8 d( W
6 W! z  t/ y7 |* s2 C: M% n


+ [* F- U5 h! h3 R% U0 A9 }2 X
3 V. l7 ?5 F3 O1 E2 [3 A0 o8 B

6 V& _" h; F; s: b, m* ]
[p=130, null, left][size=130px]多元回归方程和残差分析

[p=145, null, left][size=130px]x1=[3161.0 3710.5

[p=145, null, left][size=130px]5023.8 6060.3 6886.3 9846.8 11115.0 12153.0 14113.6 ];%

[p=145, null, left][size=130px]国内生产总值

[p=145, null, left][size=130px]x2=[77.54 78.06

  E# P+ r9 O( S% x, g[p=145, null, left][size=130px]79.05 79.22 83.62

: X4 h, I5 ]: Q8 R[p=145, null, left][size=130px]84.50 84.90 85.00 85.96];%

[p=145, null, left][size=130px]城市化率

+ N1 N1 F8 n* T' k' s' k$ z7 p/ T[p=145, null, left][size=130px]x3=[34.19 58.20

[p=145, null, left][size=130px]69.07 72.57 44.76

9 S6 q. o% W% |; B4 p[p=145, null, left][size=130px]28.30 35.95 67.76 48.98];%

[p=145, null, left][size=130px]经济使用房投资额

[p=145, null, left][size=130px]x4=[99.5 97.6

( w$ Q$ C: m% W' A4 k[p=145, null, left][size=130px]105.8 101.8 99.5

9 J  j9 k+ L8 ?2 J) I9 a[p=145, null, left][size=130px]103.5 115.2 99.4 102.7 ];%

[p=145, null, left][size=130px]物价水平（建材价格）

[p=145, null, left][size=130px]x5=[100.0 100.0

[p=145, null, left][size=130px]101.0 102.0 103.7 105.3 104.6 104.5 107.8 ];%

[p=145, null, left][size=130px]居民用地价格指数

[p=145, null, left][size=130px]x6=[19.01 21.43

+ z4 m' f2 ^0 p7 C" ^9 k6 q& p& Y2 \[p=145, null, left][size=130px]30.79 31.97 35.21

[p=145, null, left][size=130px]51.75 63.84 73.98 83.83];%

[p=145, null, left][size=130px]税收政策（房地产税收）

2 Q: o3 M, Q% h7 {; ^# e1 I[p=145, null, left][size=130px]x7=[4919 5062 4764 5053 6788 11553 12418 13799 17782 ];%

[p=145, null, left][size=130px]金融政策（房地产投资额）

[p=145, null, left][size=130px]X=[ones(9,1),x1',x2',x3',x4',x5',x6',x7'];

[p=145, null, left][size=130px]Y=[2112 2170 2359 2778 3168

[p=145, null, left][size=130px]3864 3800 4681 5032];%

[p=145, null, left][size=130px]商品房平均销售价格

4 X/ _. U- q- ?) Y[p=145, null, left][size=130px]alpha=0.05;

[p=145, null, left][size=130px][b,bint,r,rint,stats]=regress(Y',X,alpha)

[p=145, null, left][size=130px]y1=X*b;

[p=145, null, left][size=130px]e=Y'-y1;

[p=145, null, left][size=130px]plot(y1,e,'*');%

[p=145, null, left][size=130px]画出以

[p=145, null, left][size=130px]y^i

[p=145, null, left][size=130px]为横坐标的残差图

( }$ V# _% \: y# q[p=145, null, left][size=130px]figure;

[p=145, null, left][size=130px]rcoplot(r,rint);%

[p=145, null, left][size=130px]画出以观测值序号为横坐标的残差图

1 ~1 p& h  t! x: r' p! _+ b- U



1 H9 t! g% \1 E: E+ Y0 F! J


/ n, Y- ]* N) `2 v% r+ d/ {

至子星

表示深深看不懂

至子星

表示深深看不懂

戴上耳机