计算机控制系统专业课程设计东北电力大学

计算机控制系统专业课程设计东北电力大学LOGO集团标准化工作小组[Q8QX9QT-X8QQB8Q8-NQ8QJ8-M8QMN]1.设计背景由于滞后-超前校正适用于对校正后系统的动态和静态性能有更多更高要求的场合。施加滞后-超前校正环节，主要是利用其超前部分增大系统的相位裕度，以改善系统的动态性能；利用其滞后部分改善系统的静态性能。所以我们在此利用超前校正的目的是改善系统的动态性能，实现在系统静态性能不受损的前提下，提高系统的动态性能。通过加入超前校正环节，利用其相位超前特性来增大系统的相位裕度，改变系统的开环频率特性。一般使校正环节的最大相位超前角出现在系统新的穿越频率点。同时利用滞后校正通过加入滞后校正环节，使系统的开环增益有较大幅度增加，同时又使校正后的系统动态指标保持原系统的良好状态。它利用滞后校正环节的低通滤波特性，在不影响校正后系统低频特性的情况下，使校正后系统中高频段增益降低，从而使其穿越频率前移，达到增加系统相位裕度的目的。2.设计题目介绍.设计目的采用W变换的Bode图法设计离散控制系统的数字控制器，使得系统满足一定的性能指标，并且采用仿真实验对校正后的系统进行校验。.设计过程图离散控制系统结构图某控制系统如图1所示，已知被控对象的传递函数为(1)要求用W变换的Bode图法设计数字控制器Z)(z),使校正后的系统满足以下性能指标：幅值裕度为"10四，相位裕度7245。，^,>4,采样周期T=0.2s；(2)使用MATLAB编程序绘制校正前后离散系统Bode图；(3)使用MATLAB编程序绘制校正后系统单位阶跃响应，并求取各项性能指标；(4)增大采样周期(如T=2s),重复上述校正过程及仿真校验，比较分析采样周期对系统校正效果的影响.4系统仿真程序与仿真曲线T二校正前后离散系统bode图clearw=::100yl=20.*logl0(4)+20.*logl0((.*w)2+1).+20.*logl0((.*w).八2+1)*logl0(w)-20.*logl0((.*w)2+1)."y2=20.*logl0(4)+20.*logl0((.*w)/2+1).+20.*logl0((.*w)/2+1).+20.*logl0((.*w),2+1).八一20.*logl0(w)-20.*logl0((.*w).C2+1),八-20.*logl0(((68.*w).-2+1).八;subplot(2,1,1)semilogx(w,yl,'-',w,y2,':')legendbefore','after，)title('幅频曲线')xlabel('w(rad/s)’)ylabelCMagnitude(dB)')holdongridonphil=atan.*w).*180/pi+atan.*w).*180/pi-90-atan.*w).*180/pi;phi2=atan.*w).*180/pi+atan.*w).*180/pi+atan.*w).*180/pi-90-atan.*w).*180/pi-atan(68.*w).*180/pi;subplot(2,1,2)semilogx(w,phil,'-'，w,phi2,':')legendbefore','after，)title('幅频曲线')xlabel('w(rad/s)’)ylabel(?Magnitude(dB)')gridon图T=离散系统校正前后bode图T=2s校正前后离散系统bode图clearw=::100yl=20.*logl0(4)+20.*logl0(((1.*w)."2+1).+20.*logl0((.*w).2+1)20.*logl0(w)-20.*logl0((.*w).-2+1)."y2=20.*logl0(4)+20.*logl0(((1.*w)2+1).\+20.*logl0((.*w)2+1).>20.*logl0(((25.*w).八2+1).-20.*logl0(w)-20.*logl0((.*w).^2+1)20.*log!0(((240.*w)「2+l)」;subplot(2,1,1)semilogx(w,yl,'-，w,y2,':')legendCbefore','after')

title('幅频曲线')xlabel('w(rad/s)’)ylabel(JMagnitude(dB)')holdongridonphil=atan(-1.*w).*180/pi+atan.*w).*180/pi-90-atan.*w).*180/pi;phi2=atan(-1.*w).*180/pi+atan.*w),*180/pi+atan(25.*w).*180/pi-90-atan.*w).*180/pi-atan(240.*w).*180/pi;subplot(2,1,2)semilogx(w,phil,',，w,phi2,':')legendCbefore','after5)title('幅频曲线')xlabel('w(rad/s)J)ylabelCMagnitude(dB)')gridon图T=2s离散系统校正前后bode图T二校正后离散系统单位阶跃响应clc,clear,closeallnum=4;den=[l10];[A,B,C,D]=tf2ss(num,den);r=l;x=[0;0];clc,clear,closeallnum=4;den=[l10];[A,B,C,D]=tf2ss(num,den);r=l;x=[0;0];T=l;h=;Tk=T/h;Time=100;y_0=0;ym=l;e_l=0;e_2=0;u_l=0;u_2=0;fort=l:Timee(t)=r-y0;u(t)=*e(t)*e_l+*u_l;clc,clear,closeallnum=4;den=[l10];[A,B,C,D]=tf2ss(num,den);r=l;x=[0;0];T=l;h=;Tk=T/h;Time=100;y_0=0;ym=l;clc,clear,closeallnum=4;den=[l10];[A,B,C,D]=tf2ss(num,den);r=l;x=[0;0];T=l;h=;Tk=T/h;Time=100;y_0=0;ym=l;e_l=0;e_2=0;u_l=0;u_2=0;fort=l:Timee(t)=r-y0;u(t)=*e(t)*e_l+*u_l;%G0传递函数%传递函数变换为状态方程%单位阶跃输入,状态变量%采样周期，连续部分仿真步长，仿真时间%输出初始化%误差，控制器输出初始化%控制器输出仿真fork=l:Tkfork=l:Tkkl=A*x+B*u(t);fork=l:Tkk2=A*(x+h/2*kl)+B*u(t);k3=A*(x+h/2*k2)+B*u(t);k4=A*(x+h*k3)+B*u(t);x=x+(kl+2*k2+2*k3+k4)*h/6;y((t-1)*Tk+k)=C*x+D*u(t);y_l=y((t-l)*Tk+k);%G0输出y_O=y_l;endu_2=u_l;u_l=u(t);e_2=e_l;e_l=e(t);endplot((1:Time*Tk)*h,y)Function[pos,tr,ts,tp]=stepchar(gO,delta)[y,t]=step(gO);_mp,ind]=max(y);dimt=length(t);yss=y(dimt);pos=100*(mp-yss)/yss;tp=t(ind);fori=l:dimtify(i)>=ltr=t⑴；break;endend;fori=l:length(y)ify(i)<=(l-delta)*yss|y(i)>=(1+delta)*yssts=t(i);endend校正后系统阶跃响应5结论参考文献[1]姜学军，刘新国，李晓静