这是我的程序
麻烦你帮我看一下吧~~
p =[-60.5 -63.5 -43.67 -44.94 -52.76 -94.71
-43.4 -0.34 14.93 40.33 58 50.69
29.34 23.24 5.83 13.9 46.44 58.48
55.79 51.16 49.31 62.97 70.94 71.26
80.55 84.9 49.52 24.09 6.33 -9.33
-8.67 -9.5 -15 18.17 -10.77 29.31
-20.63 -15.83 -55.67 5.44 9.33 21
27.94 36.45 52.33 71.39 87.29 82.93
82.33 79.33 100.37 129.18 134.48 112.5
151.54 108.53 97.65 103.81 100.43 85.28
58.52 47.67 64.94 99.29 54.33 62
41.12 17.88 -6.13 -15.86 -43.86 -10.17
-24.73 -38.09 -68.82 -71.14 -82.5 -63.14
-76.39 -45.06 -60.29 -66.92 -65.83 -14.39
0.71 43.57 116.91 185.73 143.24 120.27
100.1 104.03 114.12 142.62 127.9 127.19
117.88 112.15 130.58 98 117.69 131.81
118.29 106.86 90.57 52.24 44.42 56.55
114.37 122.43 108.33 98.18 70.86 50.08
-0.83 -48.48 -29.06 -75.52 -76.57 -81.8
-54.84 -81.06 -99.87 -93.53 -20.12 12.53
0.24 -5 -13 -15 -58.85 -36.47
74.18 51.76 76.69 22.84 -47.56 -57.12
-30.31 -47.43 -0.63 -15.57 -67.94 -78.63
-46.76 -90.56 -71.69 -86.71 19.14 -6.13
-0.3 14.76 -6.94 10.86 35.29 48.97
131.64 101.85 65.94 45.94 73.94 65.43
87.09 65.06 113.41 112.61 101.73 101.65
46.67 26.29 -28 -16.64 11.25 -16.29
-20.16 0.84 -199.42 17.82 -12.19 -81.33
-58.06 7.56 -11.64 13.81 -31.63 -71.31
-46.61 -22.41 12.88 83.97 68.26 69.28
143 97.48 97 35.94 194 163.79
171.56 173.68 158 145.45 112.24 93.33
84.45 99.72 110.55 103.47 83.97 62.06];
t=[47.9 -28.67 9.91 62.17 44.37 8.67];
[pn,minp,maxp,tn,mint,maxt]=premnmx(p,t);
net=newff([minp,maxp],[8,1],{'tansig','purelin'},'traingd');
net.trainParam.epochs=5000;%设置最大训练循环次数
net.trainParam.lr=0.001; %设置学习系数
net.trainParam.goal=0.00001;
net.trainParam.show=50;
net=train(net,pn,tn);
an=sim(net,pn);
a=postmnmx(an,mint,maxt);
net.IW{1};
net.LW{2};
pnew=[71.26 80.55 84.9 49.52 24.09 6.33
-9.33 -8.67 -9.5 -15 18.17 -10.77
29.31 -20.63 -15.83 -55.67 5.44 9.33
21 27.94 36.45 52.33 71.39 87.29
82.93 82.33 79.33 100.37 129.18 134.48
112.5 151.54 108.53 97.65 103.81 100.43
85.28 58.52 47.67 64.94 99.29 54.33
62 41.12 17.88 -6.13 -15.86 -43.86
-10.17 -24.73 -38.09 -68.82 -71.14 -82.5
-63.14 -76.39 -45.06 -60.29 -66.92 -65.83
-14.39 0.71 43.57 116.91 185.73 143.24
120.27 100.1 104.03 114.12 142.62 127.9
127.19 117.88 112.15 130.58 98 117.69
131.81 118.29 106.86 90.57 52.24 44.42
56.55 114.37 122.43 108.33 98.18 70.86
50.08 -0.83 -48.48 -29.06 -75.52 -76.57
-81.8 -54.84 -81.06 -99.87 -93.53 -20.12
12.53 0.24 -5 -13 -15 -58.85
-36.47 74.18 51.76 76.69 22.84 -47.56
-57.12 -30.31 -47.43 -0.63 -15.57 -67.94
-78.63 -46.76 -90.56 -71.69 -86.71 19.14
-6.13 -0.3 14.76 -6.94 10.86 35.29
48.97 131.64 101.85 65.94 45.94 73.94
65.43 87.09 65.06 113.41 112.61 101.73
101.65 46.67 26.29 -28 -16.64 11.25
-16.29 -20.16 0.84 -199.42 17.82 -12.19
-81.33 -58.06 7.56 -11.64 13.81 -31.63
-71.31 -46.61 -22.41 12.88 83.97 68.26
69.28 143 97.48 97 35.94 194
163.79 171.56 173.68 158 145.45 112.24
93.33 84.45 99.72 110.55 103.47 83.97
62.06 47.9 -28.67 9.91 62.17 44.37
8.67 -1 -24.53 -45.91 -50.81 -78.68
-27.25 -116.59 -132.18 -144.61 -163.63 -147.39
-104.5 -51.91 -124.5 -48.69 -41.94 6.08];
z=sim(net,pnew);
pnewn=tramnmx(pnew,minp,maxp);
anewn=sim(net,pnew);
anew=postmnmx(anewn,mint,maxt);
[本帖最后由 Jaredcau 于 2010-4-15 15:40 编辑]
2010-4-15 09:52 上传
点击文件名下载附件
1.41 KB, 下载次数: 2636