磁盘驱动读取系统的分析设计.doc

磁盘驱动读取系统的分析设计 一、闭环系统的性能分析 （1）确定使闭环系统稳定的Ka的取值范围 G1=tf([5000],[1,1000]); G2=tf([1],conv([1,0],[1,20])); G=series(G1,G2) Transfer function: 5000 ------------------------ s^3 + 1020 s^2 + 20000 s g3= 一一开环传函 G3= 一一闭环传函 syms K den den=[1 1020 20000 5000*K]; K=den(2)*den(3)/den(1)/5000 K = 4080 有劳斯判据可得k的范围是0 K 4080 （2）在上述取值范围内取较小和较大的两个Ka值，仿真闭环系统的阶跃响应，并进行分析 K=100时 g=100*G g1=feedback(g,1) C=dcgain(g1) Transfer function: 500000 ------------------------ s^3 + 1020 s^2 + 20000 s Transfer function: 500000 --------------------------------- s^3 + 1020 s^2 + 20000 s + 500000 C = 1 [c,t]=step(g1); [y,k]=max(c); percentovershoot=100*(y-C)/C percentovershoot = 21.6918 t=setllingtime(g1) t = 0.3697 K=1000时 g=1000*G Transfer function: 5e006 -------------------------------- s^3 + 1020 s^2 + 20000 s g2=feedback(g,1) Transfer function: 5e006 ------------------------------------------ s^3 + 1020 s^2 + 20000 s + 5e006 [c,t]=step(g2); C=dcgain(g2) C = 1 [y,k]=max(c) y = 1.7109 k = 11 percentovershoot=100*(y-C)/C percentovershoot = 71.0891 t=setllingtime(g2) t = 0.4989 超调量 调节时间0.3697(s) K=1000时 超调量 调节时间0.4989 (s) （3）考察扰动信号为单位阶跃时，上述两个Ka取值情况下，系统的抗干扰能力，并进行分析 g2=tf([1],conv([1 0],[1 20])) Transfer function: 1 ---------- s^2 + 20 s g1=tf([5000],[1 1000]) Transfer function: 5000 ----------- s + 1000 syms k g3=feedback(g2,g1,1) Transfer function: s + 1000 ------------------------------- s^3 + 1020 s^2 + 20000 s - 5000 g3=feedback(g2,-g1,1) Transfer function: s + 1000 ------------------------------- s^3 + 1020 s^2 + 20000 s + 5000 g4=-g3 Transfer function: -s - 1000 ------------------------------- s^3 + 1020 s^2 + 20000 s + 5000 --------------------扰动输入的传递函数 当K=100时 g=tf([-1 -1000],[1 1020 20000 500000]) Transfer function: -s - 1000 --------------------------------- s^3