请教 灰色理论GM(1,1) 在MATLAB下的编程

simily

请教完整的灰色理论GM(1,1)在MATLAB下的程序 多谢！

small_bridge

我以前做过
这个其实很简单，我还把源代码发到simwe上了

aiwa

这是我以前做过的一个，可以参考一下

更多原创程序，请访问GreenSim团队的主页
http://blog.sina.com.cn/greensim

1. clc
   
2. clear
   
3. %原始数据输入；
   
4. x0=[55.1;54.7;52.2;54.1;54.1;54.5;55.0];
   
5. n=7;
   
6. %由数列x0生成数列x1;
   
7. s=0;
   
8. for i=1:n
   
9.     s=s+x0(i);
   
10.     x1(i)=s;
    
11. end
    
12. %由数列x1生成矩阵C和A；
    
13. for j=1:n-1
    
14.     C(j,1)=x1(j+1);
    
15.     A(j,1)=x1(j);
    
16.     A(j,2)=1;
    
17. end
    
18. %计算出估计值a1,a和u;
    
19. beta=inv(A'*A)*A'*C;
    
20. beta1=beta(1);
    
21. beta2=beta(2);
    
22. %求a,u值；
    
23. a=-log(beta1);
    
24. u=beta2/(1-beta1)*a;
    
25. %求出数列x1的预测值数列x2;
    
26. m=input('m=');
    
27. for k=0:(m-1)
    
28.     x2(k+1)=(x0(1)-u/a)*exp(-a*k)+u/a;
    
29. end
    
30. %求出原始数列x0的预测值数列x3;
    
31. x3(1)=x0(1);
    
32. for k=1:(m-1)
    
33.     x3(k+1)=(1-exp(a))*(x0(1)-u/a)*exp(-a*k);
    
34. end
    
35. %残差数据输入；
    
36. d0=[-0.3437   -0.2226   -0.3409   -0.2997   -0.0000];
    
37. k0=3;
    
38. %由数列d0生成数列d1;
    
39. sum=0;
    
40. for i=1:(n-k0+1)
    
41.     sum=sum+d0(i);
    
42.     d1(i)=sum;
    
43. end
    
44. %由数列d1生成矩阵C1和A1；
    
45. for j=1:(n-k0)
    
46.     C1(j,1)=d1(j+1);
    
47.     A1(j,1)=d1(j);
    
48.     A1(j,2)=1;
    
49. end
    
50. %计算出估计值b,a1和u1;
    
51. b=inv(A1'*A1)*A1'*C1;
    
52. b1=b(1);
    
53. b2=b(2);
    
54. a1=-log(b1);
    
55. u1=b2/(1-b1)*a1;
    
56. %求出数列d1的预测值数列d2;
    
57. for k=0:(n-k0)
    
58.     d2(k+1)=(d0(1)-u1/a1)*exp(-a1*k)+u1/a1;
    
59. end
    
60. %求出原始数列d0的预测值数列d3;
    
61. d3(1)=d0(1);
    
62. for k=1:(m-1)
    
63.     d3(k+1)=(1-exp(a1))*(d0(1)-u1/a1)*exp(-a1*k);
    
64. end
    
65. %求出原始数列x0的修正预测值数列x4;
    
66. x4(1)=x0(1);
    
67. for k=1:(m-1)
    
68.     if   k>=k0
    
69.          x4(k+1)=x2(k+1)-a1.*d3(k-k0+1)+u1;
    
70.     else x4(k+1)=x2(k+1);
    
71.     end
    
72. end
    
73. %求出原始数列x0的预测值数列x3;
    
74. x3(1)=x0(1);
    
75. for k=1:m-1
    
76.     x3(k+1)=(1-exp(a))*(x0(1)-u/a)*exp(-a*k);
    
77. end
    
78. x3
    
79. %残差检验；
    
80. for k=1:n
    
81. d(k)=x1(k)-x4(k);
    
82. end
    
83. d
    
84. %相对误差；
    
85. for k=1:n
    
86.     e(k)=abs(d(k)/x0(k))*100;
    
87. end
    
88. e
    
89. s5=0;
    
90. for i=1:n
    
91.     s5=s5+d(i);   
    
92. end
    
93. p=1/n*s5
    
94. %后验检验；
    
95. s1=0;s3=0;
    
96. for k=1:n
    
97. s1=s1+x0(k);
    
98. end
    
99. for k=1:n
    
100.     s3=s3+d(k);
     
101. end
     
102. mean1=1/n*s1;
     
103. mean2=1/n*s3;
     
104. s2=0;s4=0;
     
105. for k=1:n
     
106. s2=s2+(x0(k)-mean1)^2;
     
107. end
     
108. for k=1:n
     
109. s4=s4+(d(k)-mean2)^2;
     
110. end
     
111. sigma1=sqrt(1/n*s2)
     
112. sigma2=sqrt(1/n*s4)
     
113. c=sigma2/sigma1
     
114. p1=0.6745*sigma1;
     
115. sum=0;
     
116. for k=1:n
     
117.     if abs(d(k)-mean2)<p1
     
118.         sum=sum+1;
     
119.     else sum;
     
120.     end
     
121. end
     
122. p2=sum/n

复制代码

【GreenSim团队简介】

——更多原创程序，请访问GreenSim团队的主页→http://blog.sina.com.cn/greensim

    GreenSim团队长期从事建模仿真、算法设计、代写程序等业务，为您提供实验、课题、论文、项目等方面的仿真服务。
    GreenSim团队擅长的学科领域主要有
【智能算法】遗传算法、模拟退火、神经网络、蚁群算法、免疫优化、支持向量机、拟物拟人
【运筹优化】数学建模、数值优化算法、大规模/非线性问题、组合优化、NP完全问题
【信息与通信】无线通信、信号处理、计算机网络
【其它学科】任何问题，只要能归结为数学问题或者算法仿真，我们都有实力和信心解决
    GreenSim团队的联系方式：
QQ: 330264704，535115369  Email：greensim@163.com（来 信 必 复）

[ 本帖最后由 bainhome 于 2007-3-18 20:06 编辑 ]

snow_man_0

对线性方程组的欠定和超定问题是否要考虑？

liangbaojian

怎么没有人顶呀,我顶一下

franklings

[code]%由数列d0生成数列d1;
sum=0;
for i=1: (n-k0+1)
    sum=sum+d0(i);
    d1(i)=sum;
end[/code]
这一段直接用d1=cumsum(d0)就得到d1序列了！

azure0

我是菜鸟 弱弱问一句 这个是存为M文件再用 还是直接输入？

lichanghui

都是高手啊 菜鸟学习中