计算机控制系统课程设计论文
计算机课控制系统程设计论文
班级:0311410
姓名:王超
学号:031141010
1.设计背景 (3)
2.设计过程 (3)
2.1.设计目的 (3)
2.2.设计过程 (3)
3.系统仿真程序及相关结论 (4)
3.1.系统仿真流程图 (4)
3.2.校正前离散系统伯德图绘制程序 (5)
3.3.校正前离散系统伯德图 (6)
3.4.校正后离散系统伯德图绘制程序 (6)
3.5.校正前离散系统伯德图 (7)
3.6.校正前后离散系统伯德图绘制程序 (8)
3.7.校正前后离散系统伯德图 (9)
3.8.校正后系统时域响应流程图 (9)
3.9.校正后系统时域响应程序 (10)
3.10.校正后系统时域响应曲线 (11)
4.增大采样周期和开环增益后系统实现（T=2,K=2） (11)
4.1.改变后伯德图绘制程序 (11)
4.2.系统校正前后伯德图对比 (13)
5.心得体会.................................................................................................. 错误！未定义书签。
6.参考文献 (13)
1.设计背景
由于滞后-超前校正适用于对校正后系统的动态和静态性能有更多更高要求的场合。施加滞后-超前校正环节，主要是利用其超前部分增大系统的相位裕度，以改善系统的动态性能；利用其滞后部分改善系统的静态性能。
所以我们在此利用超前校正的目的是改善系统的动态性能，实现在系统静态性能不受损的前提下，提高系统的动态性能。通过加入超前校正环节，利用其相位超前特性来增大系统的相位裕度，改变系统的开环频率特性。一般使校正环节的最大相位超前角出现在系统新的穿越频率点。同时利用滞后校正通过加入滞后校正环节，使系统的开环增益有较大幅度增加，同时又使校正后的系统动态指标保持原系统的良好状态。它利用滞后校正环节的低通滤波特性，在不影响校正后系统低频特性的情况下，使校正后系统中高频段增益降低，从而使其穿越频率前移，达到增加系统相位裕度的目的。 2.设计过程 2.1.设计目的
采用W 变换的Bode 图法设计离散控制系统的数字控制器，使得系统满足一定的性能指标，并且采用仿真实验对校正后的系统进行校验。 2.2.设计过程
()
D z 0()
G s ZOH
+
-
()e t *()r t ()e t ()u t *()
y t ()R z ()E z ()
Y z ()
U z T T ()
G z
图1 离散控制系统结构图
某控制系统如图1所示，已知被控对象的传递函数为
0()(1)
k
G s s s =
+
（1）设采样周期1T s =，要求系统的相位裕度为45，且为了满足准确度要求所需要的v k 值而要求1k =，试用W 变换的Bode 图法确定数字控制器()D z ； （2）使用MATLAB 编程序绘制校正前后离散系统Bode 图；
（3）使用MATLAB 编程序绘制校正后系统时域响应，并求取响应的各项性能指标； （4）增大采样周期（如2T s =）和开环增益（如2k =），重复上述校正过程及仿真校验，比较分析采样周期和开环增益对系统校正效果的影响；
（5）根据校正后的系统时域响应，确定其超调量，若使得超调量减小30%，重新确定数字控制器()D z 。 3.系统仿真程序及相关结论 3.1.系统仿真流程图
V=0.01:0.01:100
遍取?
输入传递函数
G （W ）
结束
开始
W=jv 将G （W ）变为
频域
L(v)=20lg|G(jv)|Phi(v)=Phi(G(jv))
Semilogx(v,L(v))Semilogx(v,Phi(v))
V=0.01:0.01:100
Yes No
图2 系统仿真流程图
3.2.校正前离散系统伯德图绘制程序
w=0.01:0.01:100
L=20.*log10(((0.082.*w).^2+1).^0.5)+20.*log10(((-0.5.*w).^2+1).^0.5)-20.*log10(w)-20.*log10(((1.082.*w).^2+1).^0.5)
Phi=atan(0.082.*w).*180/pi+atan(-0.5.*w).*180/pi-90-atan(1.082.*w).*1 80/pi
figure(1)
subplot(2,1,1)
semilogx(w,L)
grid on
title('幅频曲线')
xlabel('w/rad/s')
ylabel('Magnitude(dB)')
hold on
subplot(2,1,2)
semilogx(w,Phi)
grid on
title('相频曲线')
xlabel('w/rad/s')
ylabel('Phase(deg)')
3.3.校正前离散系统伯德图
图3 校正前系统伯德图
由伯德图可知，相位裕度<45°，不满足设计要求，需要进行校正。
3.4.校正后离散系统伯德图绘制程序
v=0.01:0.01:100;
w=j.*v;
D1w=(1.368.*w+1)./(0.555.*w+1);
D2w=(0.972.*w+1)./(0.648.*w+1);
Gw=(0.082.*w+1).*(-0.5.*w+1)./(w.*(1.082.*w+1));
G1=Gw.*D1w.*D2w;
L1=20.*log10(abs(G1));
Phi1=atan(0.082.*v)*180/pi+atan(-0.5.*v))*180/pi+atan(1.368.*v))*180/ pi+atan(0.972.*v))*180/pi-90-atan(1.082.*v))*180/pi-atan(0.555.*v)*18 0/pi+atan(0.648.*v))*180/pi;
figure(1);
subplot(2,1,1);
semilogx(v,L1);
grid on;
title('幅频曲线');
xlabel('v(rad/s)');
ylabel('Magnitude(dB)');
subplot(2,1,2);
semilogx(v,Phi1);
grid on;
title('相频曲线');
xlabel('v(rad/s)');
ylabel('Phase(deg)');
3.5.校正前离散系统伯德图
图4 校正后系统伯德图
由伯德图可知，经过校正，相位裕度>45°，系统的各方面性能满足要求。
3.6.校正前后离散系统伯德图绘制程序
clear
x=0.01:0.1:0.01:100;
y1befor=20*(log10(sqrt(1+(0.5^2).*(x.^2)))+log10(sqrt((0.082^2).*(x.^ 2)+1))-log10(sqrt(x.^2))-log10(sqrt((1.082^2).*(x.^2)+1)));
y1after=20*(log10(0.66)+log10(sqrt(1+(1.758^2).*(x.^2)))+log10(sqrt(1 +(0.5^2).*(x.^2)))+log10(sqrt((0.082^2).*(x.^2)+1))-log10(sqrt(1+(0.8 62^2).*(x.^2)))-log10(sqrt(x.^2))-log10(sqrt((1.082^2).*(x.^2)+1))); subplot(211)
semilogx(x,y1befor,'-',x,y1after,':')
legend('befor','after')
hold on
grid
y2befor=-90-atan(1.082.*x)*180/pi-atan(0.5.*x)*180/pi+atan(0.082.*x)* 180/pi;
y2after=-90-atan(1.082.*x)*180/pi-atan(0.5.*x)*180/pi+atan(0.082.*x)* 180/pi+atan(1.758.*x)*180/pi-atan(0.862.*x)*180/pi;
subplot(212)
semilogx(x,y2befor,'-',x,y2after,':')
legend('befor','after')
grid
3.7.校正前后离散系统伯德图
图5 校正后系统伯德图由伯德图可知，校正后系统性能得到提高。
3.8.校正后系统时域响应流程图
开始
初始化变量h y_0 Time u_1
及向量a b c d x
定义循环嵌套，一个外循环
和内循环
绘制校正后系统时域响应
曲线
输出超调量
结束
图6 校正后系统时域响应流程图
3.9.校正后系统时域响应程序
lc,clear,close all
num=1;den=[1 1 0] ; %G0传递函数
[A,B,C,D]=tf2ss(num,den); %传递函数变换为状态方程
r=1;x=[0;0]; %单位阶跃输入,状态变量
T=1;h=0.1;Tk=T/h;Time=60; %采样周期，连续部分仿真步长，仿真时间
y_0=0;ym=1; %输出初始化
e_1=0;e_2=0;u_1=0;u_2=0; %误差，控制器输出初始化
tp=0;mp=0;
td=0;md=0.5;
ts=0;ms=0.02;
tr=0;
for t=1:Time %控制器输出仿真
e(t)=r-y_0;
u(t)=2.1187*e(t)-0.9544*e_1-0.014*e_2-0.1615*u_1+0.0111*u_2;
for k=1:Tk
k1=A*x+B*u(t);
k2=A*(x+h/2*k1)+B*u(t);
k3=A*(x+h/2*k2)+B*u(t);
k4=A*(x+h*k3)+B*u(t);
x=x+(k1+2*k2+2*k3+k4)*h/6;
y((t-1)*Tk+k)=C*x+D*u(t);
y_1=y((t-1)*Tk+k); %G0输出
y_0=y_1;
end
u_2=u_1;u_1=u(t);e_2=e_1;e_1=e(t);
end
plot((1:Time*Tk)*h,y)
3.10.系统时域响应曲线
图7 校正后系统时域响应曲线
由时域响应曲线可知，系统较快的进入稳态，但超调量偏大。
4.大采样周期和开环增益后系统实现（T=2,K=2）
4.1.改变后伯德图绘制程序
clear
x=0.01:0.1:0.01:100:1000;
y1befor=20*(log10(2.03)+log10(sqrt(1+(1.00^2).*(x.^2)))+log10(sqrt((0 .32^2).*(x.^2)+1))-log10(sqrt(x.^2))-log10(sqrt((1.31^2).*(x.^2)+1))) ;
y1after=20*(log10(2.03)+log10(sqrt(1+(0.66^2).*(x.^2)))+log10(sqrt(1+
(1.00^2).*(x.^2)))+log10(sqrt((0.32^2).*(x.^2)+1))-log10(sqrt(1+(0.05 ^2).*(x.^2)))-log10(sqrt(x.^2))-log10(sqrt((1.31^2).*(x.^2)+1))); subplot(211)
semilogx(x,y1befor,'-',x,y1after,':')
legend('befor','after')
hold on
grid
y2befor=-90-atan(1.31.*x)*180/pi-atan(1.*x)*180/pi+atan(0.32.*x)*180/ pi;
y2after=-90-atan(1.31.*x)*180/pi-atan(1.*x)*180/pi+atan(0.32.*x)*180/ pi+atan(0.68.*x)*180/pi-atan(0.05.*x)*180/pi;
subplot(212)
semilogx(x,y2befor,'-',x,y2after,':')
legend('befor','after')
grid
4.2.系统校正前后伯德图对比
图8 改变后系统校正前后离散系统伯德图
由图像可知，当开环增益，采样周期增大时，改变了系统的动态性能，稳态性能。
6.参考文献
1 姜学军，刘新国，李晓静. 计算机控制技术.北京：清华大学出版社，2009
2 胡寿松. 自动控制原理.北京：科学出版社，2007
3 薛定宇.控制系统计算机辅助设计――MATLAB语言与应用.北京：清华大学出版社，2006