资源描述：
《倒立摆离散型模糊控制器设计与仿真.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、倒立摆离散型模糊控制器设计与仿真1.数学模型图1倒立摆示意图如上图1倒立摆示意图所示，为杆与垂线的夹角，为作用力，杆的质量，杆和小车的总质量，半杆长，重力加速度，采样周期。其数学模型为：2.离散模糊控制器的设计（1）确定输入、输出变量，选择维数本文为倒立摆控制系统设计了一个双输入、单输出结构的模糊控制器。模糊控制器系统的输入变量为系统设定与实际测得之间的偏差及其偏差变化率，输出语言变量为控制量。（2）确定基本论域和论域（比例因子），，都是实数域上的连续变量，在该例中它们的变化范围分别为：的范围为，的范围为，的范围为，将它们全部变换到离散论
2、域，得到离散论域上的输入、输出变量，，。（1）定义模糊集合及隶属函数表对，，分别定义7个模糊集合，,，，，,。对，，定义相同的隶属度函数，如表1所示。表1，，隶属度函数表（4）设计模糊控制规则总结工作经验得出52条控制规则，见表2。表2控制规则表（5）求模糊关系矩阵R和模糊控制表在Matlab中新建文件，编写M文件计算模糊控制表，M文件如下。Input1_Terms=[1,2,3,4,5,6,7];Input2_Terms=[1,2,3,4,5,6,7];Output_Terms=[1,2,3,4,5,6,7];Input1_Terms_
3、Membership=[0,0,0,0,0,0,0,0,0,0,0.2,0.7,1;0,0,0,0,0,0,0,0,0.2,0.8,1,0.8,0.2;0,0,0,0,0,0,0,0.8,1,0.8,0.2,0,0;0,0,0,0,0,0.5,1,0.5,0,0,0,0,0;0,0,0.2,0.8,1,0.8,0,0,0,0,0,0,0;0.2,0.8,1,0.8,0.2,0,0,0,0,0,0,0,0;1,0.7,0.2,0,0,0,0,0,0,0,0,0,0];Input2_Terms_Membership=Input1_Terms_
4、Membership;Output_Terms_Membership=Input1_Terms_Membership;Rule=[0,0,7,7,7,7,7;7,7,6,6,6,6,6;7,5,5,5,3,2,2;6,5,5,4,3,2,2;6,6,5,3,3,3,1;2,2,2,2,2,3,1;1,1,1,1,1,0,0];fori=1:169forj=1:13R(i,j)=0;endendforInput1_Terms_Index=1:7forInput2_Terms_Index=1:7if(Rule(Input1_Terms_Ind
5、ex,Input2_Terms_Index)~=0)Output_Terms_Index=Rule(Input1_Terms_Index,Input2_Terms_Index);A=Input1_Terms_Membership(Input1_Terms_Index,:);B=Input2_Terms_Membership(Input2_Terms_Index,:);C=Output_Terms_Membership(Output_Terms_Index,:);fori=1:13forj=1:13R1(i,j)=min(A(i),B(j)
6、);endendR2=[];fork=1:13R2=[R2;R1(k,:)'];endfori=1:169forj=1:13R3(i,j)=min(R2(i),C(j));endendR=max(R,R3);endendendforInput1_value_index=1:13forInput2_value_index=1:13Input1_value_membership=Input1_Terms_Membership(:,Input1_value_index);[Max_Input1_value,Max_Input1_index]=m
7、ax(Input1_value_membership);Ad=Input1_Terms_Membership(Max_Input1_index,:);Input2_value_membership=Input2_Terms_Membership(:,Input2_value_index);[Max_Input2_value,Max_Input2_index]=max(Input2_value_membership);Bd=Input2_Terms_Membership(Max_Input2_index,:);fori=1:13forj=1
8、:13Rd1(i,j)=min(Ad(i),Bd(j));endendRd2=[];fork=1:13Rd2=[Rd2;Rd1(k,:)'];endforj=1:13Cd(j)=max(min
