资源描述:
《智能控制--专家PID控制仿真程序.doc》由会员上传分享,免费在线阅读,更多相关内容在应用文档-天天文库。
1、智能控制作业第2章2-3:专家PID控制仿真程序:%ExpertPIDControllerclearall;%清理数据库中所有数据closeall;%关闭所有界面图形ts=0.001;%对象采样时间,1mssys=tf(133,[1,25,0]);%受控对象的传递函数,transferfunctiondsys=c2d(sys,ts,'z');%连续系统转化为离散系统,continoustodiscrete[num,den]=tfdata(dsys,'v');%离散化后参数,得到num和den值u_1=0
2、;u_2=0;%设定初值,u_1是第(k-1)步控制器输出量y_1=0;y_2=0;%设定初值,y_1是第(k-1)步系统对象输出量x=[0,0,0]';%设定误差x1、误差导数x2、误差积分x3三个变量的初值,x2_1=0;%设定误差导数x2_1的初值error_1=0;%设定误差error_1的初值kp=0.6;%设定比例环节系数kd=0.01;%设定微分环节系数ki=0.03;%设定积分环节系数fork=1:1:5000%for循环开始,k从1变化到500,每步的增量为1time(k)=k*ts;
3、%仿真时长[0.0010.5]sr(k)=1.0;%系统输入信号,TracingStepSignalu(k)=kp*x(1)+kd*x(2)+ki*x(3);%PID控制器,PIDController%开始设计专家PID控制器及规则,Expertcontrolrule%规则1:开环控制,Rule1:Unclosedcontrolruleifabs(x(1))>0.8%if循环开始,产生式规则,if...then...;误差的绝对值大于..u(k)=0.45;%控制器输出量等于..elseifabs(x(
4、1))>0.40u(k)=0.40;elseifabs(x(1))>0.20u(k)=0.12;elseifabs(x(1))>0.01u(k)=0.10;end%if循环结束%规则2,Rule2ifx(1)*x(2)>0
5、(x(2)==0)%if循环开始,如果误差增大或不变ifabs(x(1))>=0.05%内嵌if循环开始,如果误差绝对值大于..u(k)=u_1+2*kp*x(1);%控制器输出量施加较强控制else%否则u(k)=u_1+0.4*kp*x(1);%控制器输出量施加一般控制end%内
6、嵌if循环结束end%if循环结束%规则3,Rule3if(x(1)*x(2)<0&x(2)*x2_1>0)
7、(x(1)==0)%if循环开始,如果误差减小或消除u(k)=u(k);%控制器输出量不变end%if循环结束%规则4,Rule4ifx(1)*x(2)<0&x(2)*x2_1<0%if循环开始,如果误差处于极值状态ifabs(x(1))>=0.05%内嵌if循环开始,如果误差绝对值大于..u(k)=u_1+2*kp*error_1;%控制器输出量施加较强控制else%否则u(k)=u_1+0.
8、6*kp*error_1;%控制器输出量施加一般控制end%内嵌if循环结束end%if循环结束%规则5,运用PI控制来消除误差Rule5:IntegrationseparationPIcontrolifabs(x(1))<=0.001%if循环开始,如果误差绝对值小于..(很小)u(k)=0.5*x(1)+0.010*x(3);%控制器输出量用比例和积分输出end%if循环结束%对控制输出设限,Restrictingtheoutputofcontrollerifu(k)>=10u(k)=10;%设控制
9、器输出量上限值endifu(k)<=-10u(k)=-10;%设控制器输出量下限值end%z变化后系统的线性模型,Linearmodely(k)=-den(2)*y_1-den(3)*y_2+num(1)*u(k)+num(2)*u_1+num(3)*u_2;error(k)=r(k)-y(k);%系统误差error的表达式,等于系统输入减去输出%----------Returnofparameters------------%每步计算时的参数更新u_2=u_1;u_1=u(k);%u(k)代替u_1,
10、..y_2=y_1;y_1=y(k);%y(k)代替y_1,.x(1)=error(k);%赋误差error值于x1,CalculatingPx2_1=x(2);%赋当前步计算时的误差导数x2的值于x2_1x(2)=(error(k)-error_1)/ts;%求误差导数x2,用于下一步的计算,CalculatingDx(3)=x(3)+error(k)*ts;%求误差积分x3;CalculatingIerror_1=error(k);%%