模糊控制

jiantao310

模糊控制，附带说明，希望对大家有用哈！

语言变量X，Y，Z的隶属度函数.

论域      语   言值	-4	-3	-2	-1	0	1	2	3	4
NB	1	0.5	0	0	0	0	0	0	0
NS	0	0.5	1	0.5	0	0	0	0	0
ZO	0	0	0	0.5	1	0.5	0	0	0
PS	0	0	0	0	0	0.5	1	0.5	0
PB	0	0	0	0	0	0	0	0.5	1

模糊控制规则表

Z		X
		NB	NS	ZO	PS	PB
Y	NB	PB	PB	PS	PS	ZO
	NS	PB	PS	PS	ZO	ZO
	ZO	PS	PS	ZO	ZO	NS
	PS	PS	ZO	ZO	NS	NS
	PB	ZO	ZO	NS	NS	NB

设计带有纯延迟的一阶惯性环节（假设T=6，τ=0.02）：



的模糊控制器，观察仿真结果。

编程如下：

1. 
   
2. %被控系统建模
   
3. num=1;
   
4. den=[6,1];
   
5. [a1,b,c,d]=tf2ss(num,den);%传递函数转换到状态空间
   
6. x=[0];
   
7. %系统参数
   
8. T=0.01;h=T;td=0.02;N=1000;
   
9. nd=td/T;%系统纯延迟
   
10. R=ones(1,N);%参考输入
    
11. %定义输入和输出变量及隶属度函数
    
12. a=newfis('Simple');
    
13. a=addvar(a,'input','e',[-4 4]);
    
14. a=addmf(a,'input',1,'NB','trimf',[-4,-4,-2]);
    
15. a=addmf(a,'input',1,'NS','trimf',[-4,-2,0]);
    
16. a=addmf(a,'input',1,'ZO','trimf',[-2,0,2]);
    
17. a=addmf(a,'input',1,'PS','trimf',[0,2,4]);
    
18. a=addmf(a,'input',1,'PB','trimf',[2,4,4]);
    
19. a=addvar(a,'input','de',[-4 4]);
    
20. a=addmf(a,'input',2,'NB','trimf',[-4,-4,-2]);
    
21. a=addmf(a,'input',2,'NS','trimf',[-4,-2,0]);
    
22. a=addmf(a,'input',2,'ZO','trimf',[-2,0,2]);
    
23. a=addmf(a,'input',2,'PS','trimf',[0,2,4]);
    
24. a=addmf(a,'input',2,'PB','trimf',[2,4,4]);
    
25. a=addvar(a,'output','u',[-4 4]);
    
26. a=addmf(a,'output',1,'NB','trimf',[-4,-4,-2]);
    
27. a=addmf(a,'output',1,'NS','trimf',[-4,-2,0]);
    
28. a=addmf(a,'output',1,'ZO','trimf',[-2,0,2]);
    
29. a=addmf(a,'output',1,'PS','trimf',[0,2,4]);
    
30. a=addmf(a,'output',1,'PB','trimf',[2,4,4]);
    
31. %模糊规则矩阵
    
32. rr=[5 5 4 4 3
    
33. 5 4 4 3 3
    
34. 4 4 3 3 2
    
35. 4 3 3 2 2
    
36. 3 3 2 2 1];
    
37. r1=zeros(prod(size(rr)),3);k=1;
    
38. for i=1:size(rr,1)
    
39. for j=1:size(rr,2)
    
40. r1(k,:)=[i,j,rr(i,j)];
    
41. k=k+1;
    
42. end
    
43. end
    
44. [r,s]=size(r1);
    
45. r2=ones(r,2);
    
46. rulelist=[r1,r2];
    
47. a=addrule(a,rulelist);
    
48. %采用模糊控制器的二阶系统仿真
    
49. e=0;de=0;
    
50. ke=30;kd=5;ku=1;
    
51. for k=1:N
    
52. %输入变量变换至论域
    
53. e1=ke*e;
    
54. de1=kd*de;
    
55. if e1>=4
    
56. e1=4;
    
57. elseif e1<=-4
    
58. e1=-4;
    
59. end
    
60. if de1>=4
    
61. de1=4;
    
62. elseif de1<=-4
    
63. de1=-4;
    
64. end
    
65. %模糊推理，计算出被控对象的控制输入
    
66. in=[e1 de1];
    
67. u=ku*evalfis(in,a);
    
68. uu(1,k)=u;
    
69. if k<=nd
    
70. u=0;
    
71. else
    
72. u=uu(1,k-nd);
    
73. end
    
74. 
    
75. %控制作用于被控系统，计算系统输出；利用龙格—库塔法进行系统仿真
    
76. k0=a1*x+b*u;
    
77. k1=a1*(x+h*k0/2)+b*u;
    
78. k2=a1*(x+h*k1/2)+b*u;
    
79. k3=a1*(x+h*k2)+b*u;
    
80. x=x+(k0+2*k1+2*k2+k3)*h/6;
    
81. y=c*x+d*u;
    
82. yy(1,k)=y;
    
83. %计算系统输出误差及误差导数
    
84. e1=e;e=y-R(1,k);
    
85. de=(e-e1)/T;
    
86. end
    
87. %典型二阶环节的模糊控制输出曲线
    
88. kk=[1:N]*T;
    
89. figure(1);
    
90. plot(kk,R,'k',kk,yy,'r');grid on
    
91. xlabel('时间（秒）');ylabel('输出');
    

复制代码

fengziche

模糊控制的相关资料很多，你这个是基础内容还是有什么特别之处

csqi060106

谢啦 哈哈哈

07047224

很好。就是不知道这个规则表如何建立。能给详细回答最好