控制系统A实验报告.docVIP

《控制系统GAD》实验指导书 编者：刘辉 实验一 Mat lab使用方法和程序设计 一、实验目的 掌握Matlab软件使用的棊本方法； 熟悉Mat lab的基本运算和程序控制语句; 熟悉Matlab程序设计的基木方法。 二、实验内容 求多项式的根 — + + 3%2 + 5x + 4 程序：p=[1 23 54]; x=roots(p) 结果截图： ? p=[l 2 3 5 4]; x二roots (p) x = 0.2166 + 1.51991 0.2166 - 1.5199i -1.2166 + 0.4657i -1.2166 - 0.46571 2 己知 f = ab2-c2)+Z?4(c2-a2)+c4(a2 -的，试 使用符号运算的方法对其因式分解。 程序： Syms a b c; f=aA4*(bA2-cA2)+bA4*(cA2-aA2)+cA4*(aA2-bA2); R1 =factor(f) 结果截图： Syms a b c: ￡=aA4*(bA2-cA2)+bA4*(cA2-aA2)+cA4* (aA2_bA2): Rl = ￡actor (￡) R1 = (bc)* (b+c)* (ac)* (a+c)* (a—b)* (a+b) 结果?? 3.编与一*个函数，完成求和s=l+2+3+...+i+...+n。 程序：求1000个数相加的和 sum=0; for i=l:1000 sum=sum+i; end sum 结果截阁： sum=0; ? for i=l:1000 sum=sum+i: end sum sum 500500 5+2 4.已知一传递函数 4.已知一传递函数 5 + 5 5 + 4 ,试将其分解 为部分分式。 程序： num:[1 2]; den=[1 5 4]; [r，p，k]=residue(num，den) 结果截图： ruim=[l 2]: den=[l 5 4]: p, k]=resi due (num, den) 0.6667 0.3333 -4 -1 []结果: [] 结果: 实验二 一、 实验目的 掌握如何使用Matlab进行系统的时域分析。 掌握如何使用Matlab进行系统的频域分析。 掌握如何使用Matlab进行系统的根轨迹分析。 掌握如何使用Matlab进行系统的稳定性分析。 二、 实验内容 1.时域分析 根据下面传递函数模型：绘制其甲.位阶跃响应曲线并从图上读取最人 超调铽，绘制系统的单位脉冲响应。3 超调铽，绘制系统的单位脉冲响应。 3(52 + 55* + 6) 5 3 + 6 5 + 10 5 + 8 程序： num=[3 15 18]; den:[1 6 10 8]; G=tf (num, den); Grid;impulse(G) time=[0:0. 1:20]; step(G, time) 结果截阁： nvun=[3 15 18]; den=[l 6 10 8]; G=tf (num, den): time=[0：0. 1:20]; step (G, time) ? nujri=[3 15 K den=[l 6 10 8]： G=t￡ (rtuiTi, den): Grid: impulse (G: [Figure No. File Edit View Insert Tools Window Help a / | ^ ^ o a / | ^ ^ o Filo Edit View Insert Tools Window Help Q GB； BI ■ | k A ，Z J? o Fioure No.: I o 卜回4— Impulse Response Impulse Response 0 1 2 3 4 5 6 Time (sec) 结果分析： 超调量2.41 计算 频域分析 典堺二阶系统传递函数为:G(s)= 典堺二阶系统传递函数为: G(s)= + + vv/ 当C =0.7，con取6时的Bode Nichols Nyquist图的单位阶跃响应 (1) Bode图程序: num=36; den=[l 8.4 36]; g=tf(num,den); bode(g,{0.001, 100}); grid; num=36: den=[l 8.4 36]; g=t￡ (num, den): bode (g, {0.001, 100}); grid; i 结果截图： Eile Edit 乂iew Insert loots Window tlclp□ 穆 kA， Eile Edit 乂iew Insert loots Window tlclp □ 穆 kA，/ O Bod? Diagram Frequency (r^d/sec) 结果