[推荐]一个经典的BP网络原程序

lxq

对于初学者是挺有用的
边看边学

1. %
   
2. clf; figure(1); setfsize(500,300); echo on; clc; pause % 键入任意键继续
   
3. clc; P = [-3.0 +2.0]; T = [+0.4 +0.8]; pause
   
4. clc; wv = -4:0.4:7.6; bv = -4:0.4:7.6; es = errsurf(P,T,wv,bv,'logsig');
   
5. plotes(wv,bv,es,[60 30]); pause
   
6. clc; net=newff(minmax(P),[1],{'logsig'},'traingd','learngd','sse');
   
7. net.iw{1,1}; net.b{1}; pause; echo off; k = pickic(1);
   
8. if k == 2, net.iw{1,1} = -2.1617; net.b{1} = -1.7862;
   
9. elseif k == 3, subplot(1,2,2); h=text(2,2,sprintf('Click on me'))
   
10. [net.iw{1,1},net.b{1}] = ginput(1); delete(h)
    
11. end
    
12. echo on; clc
    
13. 
    
14. net.trainParam.show = 5; % 学习过程显示频率
    
15. net.trainparam.goal = 0.01; % 误差指标
    
16. net.trainParam.lr = 2;% 学习率
    
17. me=100; A=sim(net,P); sse=sumsqr(T-A);
    
18. h=plotep(net.iw{1,1},net.b{1},sse)
    
19. for i=1:me
    
20. if sse<net.trainParam.goal,i=i-1;break,end
    
21. net.trainParam.epochs=1;
    
22. figure(2); setfsize(50,50); [net,tr] = train(net,P,T);
    
23. trp(i)=tr.perf(1,2); A=sim(net,P); sse=sumsqr(T-A);
    
24. if rem(i,net.trainParam.show)==0
    
25. figure(1); h=plotep(net.iw{1,1},net.b{1},sse,h);
    
26. end
    
27. end
    
28. tr.p=trp(1:i);
    
29. message=sprintf('Train:%%g/%g epochs,sse=%%g.\n',me);
    
30. fprintf(message,i,sse)
    
31. figure(1); plotep(net.iw{1,1},net.b{1},sse,h); pause
    
32. subplot(1,2,2); plot(1:i,tr.p); hold on
    
33. plot(1:i,net.trainParam.goal,'r--'); hold off
    
34. title('Error Signals'); xlabel('Epochs'); ylabel('Error'); pause
    
35. p = -1.2; a = sim(net,p); echo off

复制代码

[ 本帖最后由 ChaChing 于 2009-12-28 14:49 编辑 ]

rolo

??? Undefined command/function 'setfsize'.<BR>能把这个m文件给出来吗？

yu1983619

setfsize是什么函数<BR>还有pickic提示是作废的函数

hustbill

却是，没见过这个函数<BR>

wxr1000

好多函数都没有见过阿，行不行啊，不过 仍然谢谢

dreambreez

Undefined function or variable 'setfsize'.

.m
On line 4  ==> setfsize(500,300);

wateropt

我也试了试,有好多函数每见过,程序通不过

suffer

1. function h=setfsize(wx,wy)
   
2. % h=setfsize (wx,wy)
   
3. %
   
4. % Set figure position.       
   
5. %
   
6. % h = setfsize(X,Y)
   
7. %   X - Horizontal length.
   
8. %   Y - Vertical height.
   
9. % h is the figure handle to the figure just handled
   
10. % Sets current figure lower left position to X by Y pixels (measured from
    
11. % lower left corner of screen. On image, max pixels are about
    
12. % 1280, 1024, you can use  get (0,'ScreenSize') to get the screen size
    
13. %
    
14. %  adapted from function setfsize.m
    
15. %
    
16. %                                Ziad Saad Sept 27 97
    
17. 
    
18. 
    
19. z = get(gcf,'position');
    
20. set(gcf,'position',[z(1) z(2) wx wy]);
    
21. h = gcf;
    
22. return;

复制代码

suffer

把setfsize这个函数加上就可正常运行了

不过里边存在若干作废函数

lxq

呵呵  suffer真强!

suffer

原帖由 lxq 于 2006-11-7 21:20 发表
呵呵  suffer真强!

这个在网络上就能够找到