matlab运用报告yx.docVIP

08电信（1）班 袁旭 20081305043 MATLAB在信号与系统中的运用 例：求x[n]=(-0.9)^n,-10=n=10的离散时间傅里叶变换。 n=-5:5; x=(0.9).^n; k=-200:200; w=(pi/100)*k; X=x*(exp(-j*pi/100)).^(n*k); magX=abs(X); angX=angle(X); subplot(2,1,1); plot(w/pi,magX); grid axis([-2,2,0,15]) xlabel(frequence in pi units); ylabel(abs(X)); title(Magnitude Part) subplot(2,1,2) ; plot(w/pi,angX/pi); grid axis([-2,2,-1,1]) xlabel(frequence in pi units); ylabel(Radians/pi); title(Angle Part) 例题的幅频和相频曲线 MATLAB在自动控制原理中的运用 例：设系统的传递函数为：H(s)=50/((s-1.5)*(s+1)*(s+6)) 这是一个开环不稳定的系统。画出其Nyquist曲线，判别其闭环稳定性，并用MATLAB的其他函数加以检验。在此系统上加一个零点(s+0.8)后，再做同样的工作，把两种情况进行比较讨论。 s1=zpk([],[-6,-1,1.5],50); figure(1) subplot(2,2,1),nyquist(s1),grid sb1=feedback(s1,1) 程序运行后，s1闭环后构成的系统sb1零极增益模型为： Zero/pole/gain: 50 -------------------------------- (s+6.985) (s^2 - 1.485s + 5.87) subplot(2,2,2),impulse(s1),grid subplot(2,2,3),margin(s1),grid subplot(2,2,4),impulse(sb1),grid s2=zpk([-0.8],[-6,-1,1.2],50); figure(2) subplot(2,2,1),nyquist(s2),grid sb2=feedback(s2,1) s2闭环后构成的系统sb2零极增益模型为： Zero/pole/gain: 50 (s+0.8) ---------------------------------- (s+0.7485) (s^2 + 5.051s + 43.82) subplot(2,2,2),impulse(s2),grid subplot(2,2,3),margin(s2),grid subplot(2,2,4),impulse(sb2),grid 开环不稳定系统输出图形 加一个零点后系统输出图形