课程设计任务书学生姓名：________ 专业班级： _______________指导教师：_______ 工作单位： ____________题目：飞行器控制系统设计初始条件：飞行器控制系统的开环传递函数为：G（s） -^500^s（s 361.2）控制系统性能指标为调节时间0.01s，单位斜坡输入的稳态误差0.000521，相角裕度大于84度。
要求完成的主要任务:（包括课程设计工作量及其技术要求，以及说明书撰写等具体要求）（1）设计一个控制器，使系统满足上述性能指标；（2）画出系统在校正前后的奈奎斯特曲线和波特图；（3）用Matlab画出上述每种情况的阶跃响应曲线，并根据曲线分析系统的动态性能指标；（4）对上述任务写出完整的课程设计说明书，说明书中必须写清楚分析计算的过程，给出响应曲线，并包含Matlab源程序或Simulink仿真模型，说明书的格式按照教务处标准书写时间安排：指导教师签名：系主任（或责任教师）签名:目录1串联滞后—超前校正的原理............ 错误! 未定义书签。
2 飞行器控制系统的设计过程. ................. 错误! 未定义书签。
2.1 飞行器控制系统的性能指标............... 错误! 未定义书签。
2.2 系统校正前的稳定情况................. 错误! 未定义书签。
2.2.1 校正前系统的波特图............. 错误! 未定义书签。
2.2.2 校正前系统的奈奎斯特曲线 (2)2.2.3 校正前系统的单位阶跃响应曲线......... 错误! 未定义书签。
2.3 飞行器控制系统的串联滞后—超前校正 (4)2.3.1 确定校正网络的相关参数 (4)2.3.2 验证已校正系统的性能指标 (6)2.4 系统校正前后的性能比较 (8)2.4.1 校正前后的波特图 (8)2.4.2 校正前后的奈奎斯特曲线 (9)2.4.3 校正前后的单位阶跃响应曲线 (11)3 设计总结与心得体会 (12)参考文献 (13)摘要根据被控对象及给定的技术指标要求，设计自动控制系统，既要保证所设计的系统有良好的性能，满足给定技术指标的要求，还有考虑方案的可靠性和经济性。
本说明书介绍了在给定的技术指标下，对飞行器控制系统的设计。
为了达到给定要求，主要采用了串联滞后—超前校正。
在对系统进行校正的时，采用了基于波特图的串联滞后—超前校正，对系统校正前后的性能作了分析和比较，并用MATLA进行了绘图和仿真。
对已校正系统的高频特性有要求时，采用频域法校正较其它方法更为方便。
关键词：飞行器控制系统校正MATLAB飞行器控制系统设计1 串联滞后—超前校正的原理如果系统设计要求满足的性能指标属频域特征量，则通常采用频域校正方法。
在开环系统对数频率特性基础上，以满足稳态误差、开环系统截止频率和相角裕度等要求为出发点，进行串联校正的方法。
在伯德图上虽然不能严格定量地给出系统的动态性能，但却能方便地根据频域指标确定校正装置的形式和参数，特别是对已校正系统的高频特性有要求时，采用频域法校正较其它方法更方便。
串联滞后—超前校正兼有滞后校正和超前校正的优点，当待校正系统不稳定，且要求校正后系统的响应速度、相角裕度和稳态精度要求较高时，应采用串联滞后—超前校正。
其基本原理是利用滞后—超前网络的超前部分来增大系统的相角裕度，同时利用滞后部分来改善系统的稳态性能。
串联滞后—超前校正的设计步骤如下：（1）根据稳态性能要求确定开环增益K。
（2）绘制待校正系统的对数幅频特性曲线，求出待校正系统的截止频率 c 相角裕度及幅值裕度h（dB）。
（3）在待校正系统对数幅频特性上，选择斜率从-20dB/dec 变为-40dB/dec 的交接频率作为校正网络超前部分的交接频率b。
b的这种选法，可以降低已校正系统的阶次，且可保证中频区斜率为期望的-20dB/dec ，并占据较宽的频带。
11（4）根据响应速度要求，选择系统的截止频率c和校正网络衰减因子1/a 。
要I I保证已校正系统的截止频率为所选的c，下列等式成立20lga L'（c''） 20lg T b c'' 0 （1）（5）根据相角裕度要求，估算校正网络滞后部分的交接频率a（6）校验已校正系统的各项性能指标。
根据滞后—超前校正的原理和步骤，可以在纯超前校正及纯滞后校正都不宜采用时，对系统进行串联滞后—超前校正。
奈奎斯特曲线，如图2-2所示2飞行器控制系统的设计过程2.1飞行器控制系统的性能指标飞行器控制系统的开环传递函数为G(s)4500K s(s 361.2)控制系统性能指标为调节时间t s0.01s，单位斜坡输入的稳态误差0.000521,相角裕度大于84度。
根据单位斜坡输入的稳态误差e ss0.000521，可以得出K v lim sG(s)s 04500K 1(3)361.2e ssK 1552.2系统校正前的稳定情况2.2.1校正前系统的波特图根据原有的飞行器控制系统的开环传递函数，在 波特图，如图2-1所示。
绘制校正前波特图的MATLA 源程序如下：num=697500; den=[1,361.2,0]; %校正前系统参数 bode( num,de n);%绘制波特图grid;2.2.2校正前系统的奈奎斯特曲线MATLAB 中绘制出校正前的根据原有的飞行器控制系统的开环传递函数，在MATLAB 中绘制出校正前的绘制校正前波特图的MATLA 源程序如下:num=697500; den=[1,361.2,0]; nyq uist (nu m,de n)%校正前系统参数 %绘制奈奎斯特曲线(2)图2-2系统校正前的奈奎斯特曲RealAxis(gp 〕ID10 10Frequency (radfcec)图2-1系统校正前的波特图册盪 +」沖匚一E 4Q E-Bode DsagrarriNyquist Diagram-4-1 04(4)223校正前系统的单位阶跃响应曲线校正刖系统的单位反馈闭环传递函数为C(s)697500~2R(s) s 361.2s 697500用MATLAB 绘制系统校正前的的单位阶跃响应曲线如图 2-3所示。
MATLAB 源程序如下所示:num=697500;图2-3系统校正前的单位阶跃响应曲线2.3飞行器控制系统的串联滞后一超前校正2.3.1确定校正网络的相关参数通过编写MATLAB 源程序求系统校正前的稳定裕度，程序如下:num= 697500;den=[1,361.2,697500]; %校正前系统参数 step( nu m,de n)%绘制阶跃响应曲线□ Step Response0.005 0.01 0015 0.02 0.025 0.03 0 035lime (sscjaa bden=[1,361.2,0];%系统校正前的参数[mag,phase,w]=bode （ nu m,de n ） [gm,pm,wcg,wcp]=margi n （ mag,phase,w ）% 求系统校正前的稳定裕度运行后，得出相角裕度 pm=24.4° ,截止频率wcp=797rad/s 。
由此可得，若 采用超前校正，需补偿超前角m 为m ”850 24.40 50 65.6° 60°（ 5）显然一级串联超前网络不能达到要求。
又由于要求校正后系统的响应速度、 相角裕度要求较高，所以采用串联滞后一超前校正。
从图2-1及校正前系统的开环传递函数可以看出，斜率从 -20dB/dec 变为 -40dB/dec 的交接频率为361.2rad/s,即校正网络超前部分的交接频率b361.2 rad/s 。
根据飞行器控制系统对调节时间的要求，可以得出校正以后系统的截止频率 为c630在校正前的波特图2-1可以读出待校正系统对数幅频特性的-20dB/dec 延长 线在:处的数值为24,则有20lg a a 24 015.9（9）设校正网络的传递函数为“ s s 、(1 -)(1 )G c (s)a b（10）由高阶系统频域指标与时域的关系, 有如下的公式t sK 。
(6) K 0 2 1.5(M r1) 2.5(M r 1)2 (7) M r1 sin(8)as s(1 —)(1 —)a a b(11)根据飞行器控制系统的相角裕度要求，得出递函数为2.3.2验证已校正系统的性能指标根据校正后系统的开环传递函数，验证校正后系统的相角裕度。
编写MATLAB 远程序如下：num=[73421,697500]; den=[0.105,606.88,361.2,0]; [mag,phase,w]=bode( nu m,de n) [gm,pm,wcg,wcp]=margi n( mag,phase,w)运行后得出校正后系统的相角裕度pm=84.6°，符合给定的相位裕度要求源程序如下:num=[73421,697500]; den=[0.105,606.88,361.2,0]; bode( nu m,de n) grid编写MATLA 程序，绘制已校正系统的奈奎斯特曲线，如图 2-5所示。
相应的MATLAB!程序如下：num=[73421,697500]; den=[0.105,606.88,361.2,0]; Nyquist (nu m,de n) %校正后系统参数%绘制校正后的奈奎斯特曲线编写MATLAB?序，绘制已校正系统的波特图,如图2-4所示。
相应的MATLAB考虑到180° arctg —90°aa carctgarctg361.2 rad/s,求得 〜9.5 rad/so 求得已校正系统的开环传697500(1G c (s)G °(s)s 361.2s(s 361.2)(1 些炎)(19.5s 15.9 361.2(12)G c (s)G 。
(s)73421s 6975000.105s 3 606.88s 2 361.2s(13)%校正后系统参数%求系统校正后的稳定裕度%校正后系统参数%绘制校正后的波特图图2-5系统校正后的奈奎斯特曲线Bode Diagram(Bap 】10 10 10 10 10 10Frequency (radfsec)图2-4系统校正后的波特图Nyquist Diagram10 12 3 i l a s也空卜」嗎匚-503|1|-■200(1 -1500 Real Axis^IDCO -500. 0编写MATLAB?序，绘制已校正系统的单位阶跃响应曲线，如图2-6所示相应的MATLA源程序如下：num=[73421,697500];den=[0.105,606.88,737822.2,697500];step( nun,den)grid图2-6校正后的单位阶跃响应曲线2.4系统校正前后的性能比较2.4.1校正前后的波特图确定了校正网络的各种参数，经过验证已校正系统的技术指标，基本达到标准后，可以将校正前后的性能指标进行对比。