模糊RBF神经网络——程序

博大广阔

1. %模糊RBF神经网络
   
2. %参考刘金琨老师的智能控制
   
3. clc
   
4. clear all
   
5. ts=0.001;
   
6. y_1=0; u_1=0;
   
7. 
   
8. w=rands(25,1)/10; w_1=w; w_2=w;
   
9. 
   
10. c=[-4 -2 0 2 4;
    
11.     -4 -2 0 2 4].'; %五条规则，两个输入
    
12. B=10*ones(5,2); B_1=B;B_2=B;
    
13. c_1=c; c_2=c_1;
    
14. for K=1:1:10000
    
15. %xitong
    
16. u(K)=0.4*sin(4*pi*K*ts);
    
17. y(K)=u_1^3+y_1/(1+y_1^2);
    
18. %--------------------------------------%   
    
19. 
    
20. input0(K,:)=[u(K);y_1];
    
21. %%
    
22. input=input0(K);
    
23. for n=1:1:size(c,1)
    
24. f2(n,:)=exp(-(input.'-c(n,:)).^2./(B(n,:)).^2);   
    
25. end
    
26. f3=f2(:,1)*f2(:,2).';
    
27. f3=reshape(f3,1,25).';
    
28. output=w.'*f3;
    
29. %%
    
30. ym(K)=output;
    
31. Error=y(K)-ym(K);
    
32. w=w_1+0.01*Error*f3+0.01*(w_1-w_2); w_2=w_1; w_1=w;
    
33. 
    
34. deta2=-Error*w.'*f3;
    
35. for n=1:1:size(c,1)
    
36. d_b(n,:)=2*(input.'-c(n,:)).*((B(n,:)).^-3);   
    
37. end
    
38. B=B_1+0.01*d_b+0.01*(B_1-B_2); B_2=B_1; B_1=B;
    
39. for n=1:1:size(c,1)
    
40. d_c(n,:)=-2*(input.'-c(n,:)).*((B(n,:)).^-2);  
    
41. end
    
42. c=c_1+0.01*d_c+0.01*(c_1-c_2); c_2=c_1; c_1=c;  
    
43. time(K)=K*ts;
    
44. u_1=u(K);
    
45. y_1=y(K);
    
46. end
    
47. figure(1);plot(time,y,'b',time,ym,'r');grid on
    
48. %figure(2);plot(time,Error,'r');grid on

复制代码

不是钱林

liumingze

liubo1990

学习一下，顶一个