模糊控制

模糊控制

基本内容

利用模糊数学的基本思想和理论的控制方法。在传统的控制领域里，控制系统动态模式的精确与否是影响控制优劣的最主要关键。系统动态的信息越详细，则越能达到精确控制的目的。
然而，对于复杂的系统，由于变量太多，往往难以正确的描述系统的动态，于是工程师便利用各种方法来简化系统动态，以达成控制的目的，但却不尽理想。换言之，传统的控制理论对于明确系统有强而有力的控制能力，但对于过于复杂或难以精确描述的系统，则显得无能为力了。因此便尝试着以模糊数学来处理这些控制问题。
使用模糊控制器处理
模糊控制器包括四部分：
(1)模糊化。主要作用是选定模糊控制器的输入量，并将其转换为系统可识别的模糊量，具体包含以下三步：
第一，对输入量进行满足模糊控制需求的处理;
第二，对输入量进行尺度变换;
第三，确定各输入量的模糊语言取值和相应的隶属度函数。
(2)规则库。根据人类专家的经验建立模糊规则库。模糊规则库包含众多控制规则，是从实际控制经验过渡到模糊控制器的关键步骤。
(3)模糊推理。主要实现基于知识的推理决策。
(4)解模糊。主要作用是将推理得到的控制量转化为控制输出。
代码实现

%模糊控制器设计
a=newfis('fuzzf');                   %创建新的模糊推理系统

%输入1
f1=1; 
a=addvar(a,'input','温度',[-3*f1,3*f1]);                   
 %添加 e 的模糊语言变量
a=addmf(a,'input',1,'0-20℃','zmf',[-3*f1,-1*f1]);          
 %添加 e 的模糊语言变量的隶属度函数（z型）
a=addmf(a,'input',1,'20-40℃','trimf',[-3*f1,-2*f1,0]);      
  %隶属度函数为三角形
a=addmf(a,'input',1,'40-60℃','trimf',[-3*f1,-1*f1,1*f1]); 
a=addmf(a,'input',1,'60-80℃','trimf',[-2*f1,0,2*f1]); 
a=addmf(a,'input',1,'80-100℃','smf',[1*f1,3*f1]); 

%输入2
f2=1;
a=addvar(a,'input','湿度',[-3*f2,3*f2]);                   
 %添加 ec 的模糊语言变量
a=addmf(a,'input',2,'0-0.1','zmf',[-3*f2,-1*f2]); 
a=addmf(a,'input',2,'0.1-0.2','trimf',[-3*f2,-2*f2,0]);
a=addmf(a,'input',2,'0.2-0.3','trimf',[-3*f2,-1*f2,1*f2]);
a=addmf(a,'input',2,'0.3-0.4','trimf',[-2*f2,0,2*f2]);
a=addmf(a,'input',2,'0.4-0.5','smf',[1*f2,3*f2]); 

%输出
f3=1.5;
a=addvar(a,'output','乳酸菌存活情况',[-3*f3,3*f3]);                 
   %添加 u 的模糊语言变量
a=addmf(a,'output',1,'0-0.2','zmf',[-3*f3,-1*f3]); 
a=addmf(a,'output',1,'0.2-0.4','trimf',[-3*f3,-2*f3,0]);
a=addmf(a,'output',1,'0.4-0.6','trimf',[-3*f3,-1*f3,1*f3]);
a=addmf(a,'output',1,'0.6-0.8','trimf',[-2*f3,0,2*f3]);
a=addmf(a,'output',1,'0.8-1及以上','smf',[1*f3,3*f3]);

%规则库
rulelist=[     1 1 1 1 1;             %编辑模糊规则，后俩个数分别是规则权重和AND OR选项
               1 2 1 1 1;
               1 3 1 1 1;
               1 4 2 1 1;
               1 5 2 1 1;

          
               2 1 1 1 1;
               2 2 2 1 1;
               2 3 4 1 1;
               2 4 5 1 1;
               2 5 2 1 1;
               
               3 1 2 1 1;
               3 2 4 1 1;
               3 3 3 1 1;
               3 4 2 1 1;
               3 5 2 1 1;
               
               4 1 2 1 1;
               4 2 3 1 1;
               4 3 3 1 1;
               4 4 2 1 1;
               4 5 2 1 1;
               
               5 1 1 1 1;
               5 2 1 1 1;
               5 3 1 1 1;
               5 4 1 1 1;
               5 5 1 1 1;
               ];
           
a=addrule(a,rulelist);                %添加模糊规则函数
showrule(a)                             %显示模糊规则函数
a1=setfis(a,'DefuzzMethod','centroid');                  %设置解模糊方法
writefis(a1,'fuzzf');                       %保存模糊系统
a2=readfis('fuzzf');   %从磁盘读出保存的模糊系统
disp('fuzzy Controller table:e=[-3,+3],ec=[-3,+3]');%显示矩阵和数组内容

%推理
Ulist=zeros(5,5);                                   %全零矩阵
for i=1:5
       for j=1:5
           e(i)=-4+i;
           ec(j)=-4+j;
           Ulist(i,j)=evalfis([e(i),ec(j)],a2);    %完成模糊推理计算
       end
   end
%   Ulist=ceil(Ulist)                               %朝正无穷方向取整
   Ulist                               %朝正无穷方向取整
   
%画出模糊系统
figure(1); plotfis(a2);  
figure(2);plotmf(a,'input',1);
figure(3);plotmf(a,'input',2);
figure(4);plotmf(a,'output',1);





           
           

结果展示

'1. If (温度 is 0-20℃) and (湿度 is 0-0.1) then (乳酸菌存活情况 is 0-0.2) (1) ’
'2. If (温度 is 0-20℃) and (湿度 is 0.1-0.2) then (乳酸菌存活情况 is 0-0.2) (1) ’
'3. If (温度 is 0-20℃) and (湿度 is 0.2-0.3) then (乳酸菌存活情况 is 0-0.2) (1) ’
'4. If (温度 is 0-20℃) and (湿度 is 0.3-0.4) then (乳酸菌存活情况 is 0.2-0.4) (1) ’
'5. If (温度 is 0-20℃) and (湿度 is 0.4-0.5) then (乳酸菌存活情况 is 0.2-0.4) (1) ’
'6. If (温度 is 20-40℃) and (湿度 is 0-0.1) then (乳酸菌存活情况 is 0-0.2) (1) ’
'7. If (温度 is 20-40℃) and (湿度 is 0.1-0.2) then (乳酸菌存活情况 is 0.2-0.4) (1) ’
'8. If (温度 is 20-40℃) and (湿度 is 0.2-0.3) then (乳酸菌存活情况 is 0.6-0.8) (1) ’
‘9. If (温度 is 20-40℃) and (湿度 is 0.3-0.4) then (乳酸菌存活情况 is 0.8-1及以上) (1)’
‘10. If (温度 is 20-40℃) and (湿度 is 0.4-0.5) then (乳酸菌存活情况 is 0.2-0.4) (1)’
'11. If (温度 is 40-60℃) and (湿度 is 0-0.1) then (乳酸菌存活情况 is 0.2-0.4) (1) ’
‘12. If (温度 is 40-60℃) and (湿度 is 0.1-0.2) then (乳酸菌存活情况 is 0.6-0.8) (1)’
‘13. If (温度 is 40-60℃) and (湿度 is 0.2-0.3) then (乳酸菌存活情况 is 0.4-0.6) (1)’
‘14. If (温度 is 40-60℃) and (湿度 is 0.3-0.4) then (乳酸菌存活情况 is 0.2-0.4) (1)’
‘15. If (温度 is 40-60℃) and (湿度 is 0.4-0.5) then (乳酸菌存活情况 is 0.2-0.4) (1)’
'16. If (温度 is 60-80℃) and (湿度 is 0-0.1) then (乳酸菌存活情况 is 0.2-0.4) (1) ’
‘17. If (温度 is 60-80℃) and (湿度 is 0.1-0.2) then (乳酸菌存活情况 is 0.4-0.6) (1)’
‘18. If (温度 is 60-80℃) and (湿度 is 0.2-0.3) then (乳酸菌存活情况 is 0.4-0.6) (1)’
‘19. If (温度 is 60-80℃) and (湿度 is 0.3-0.4) then (乳酸菌存活情况 is 0.2-0.4) (1)’
‘20. If (温度 is 60-80℃) and (湿度 is 0.4-0.5) then (乳酸菌存活情况 is 0.2-0.4) (1)’
'21. If (温度 is 80-100℃) and (湿度 is 0-0.1) then (乳酸菌存活情况 is 0-0.2) (1) ’
'22. If (温度 is 80-100℃) and (湿度 is 0.1-0.2) then (乳酸菌存活情况 is 0-0.2) (1) ’
'23. If (温度 is 80-100℃) and (湿度 is 0.2-0.3) then (乳酸菌存活情况 is 0-0.2) (1) ’
'24. If (温度 is 80-100℃) and (湿度 is 0.3-0.4) then (乳酸菌存活情况 is 0-0.2) (1) ’
'25. If (温度 is 80-100℃) and (湿度 is 0.4-0.5) then (乳酸菌存活情况 is 0-0.2) (1) ’

fuzzy Controller table:e=[-3,+3],ec=[-3,+3]

Ulist =

-3.6509 -3.6509 -2.9123 -2.6521 -2.4164
-2.9123 -1.3632 -0.1099 0.3138 -0.1866
-2.6521 -0.9285 -0.1609 -0.2580 -0.1866
-2.4998 -1.0353 -1.0353 -1.8743 -2.4164
-2.4164 -1.6722 -1.6722 -1.6722 -2.4164

结论

本次我对乳酸菌在不同温度和湿度下的生存情况进行模糊控制
其中有两输入（e，ec）分别代表湿度和温度，一输出（u）代表乳酸菌的存活情况
两输入及一输出均有五个模糊集，总的会有25种模糊规则。