东南大学信号与系统MATLAB实践第二次 作业.docVIP

练习二 实验六 一．用MATLAB语言描述下列系统，并求出极零点、 1. Ns=[1]; Ds=[1,1]; sys1=tf(Ns,Ds) 实验结果： sys1 = 1 ----- s + 1 [z,p,k]=tf2zp([1],[1,1]) z = Empty matrix: 0-by-1 p = -1 k = 1 2. Ns=[10] Ds=[1,-5,0] sys2=tf(Ns,Ds) 实验结果： Ns = 10 Ds = 1 -5 0 sys2 = 10 --------- s^2 - 5 s [z,p,k]=tf2zp([10],[1,-5,0]) z = Empty matrix: 0-by-1 p = 0 5 k = 10 二．已知系统的系统函数如下，用MATLAB描述下列系统。 1． z=[0]; p=[-1,-4]; k=1; sys1=zpk(z,p,k) 实验结果： sys1 = s ----------- (s+1) (s+4) Continuous-time zero/pole/gain model. 2. Ns=[1,1] Ds=[1,0,-1] sys2=tf(Ns,Ds) 实验结果： Ns = 1 1 Ds = 1 0 -1 sys2 = s + 1 ------- s^2 - 1 Continuous-time transfer function. 3． Ns=[1,6,6,0]; Ds=[1,6,8]; sys3=tf(Ns,Ds) 实验结果： Ns = 1 6 6 0 Ds = 1 6 8 sys3 = s^3 + 6 s^2 + 6 s ----------------- s^2 + 6 s + 8 Continuous-time transfer function. 六．已知下列H(s)或H(z)，请分别画出其直角坐标系下的频率特性曲线。 1. clear; for n = 1:400 w(n) = (n-1)*0.05; H(n) = (1j*w(n))/(1j*w(n)+1); end mag = abs(H); phase = angle(H); subplot(2,1,1) plot(w,mag);title(幅频特性) subplot(2,1,2) plot(w,phase);title(相频特性) 实验结果： 2. clear; for n = 1:400 w(n) = (n-1)*0.05; H(n) = (2*j*w(n))/((1j*w(n))^2+sqrt(2)*j*w(n)+1); end mag = abs(H); phase = angle(H); subplot(2,1,1) plot(w,mag);title(幅频特性) subplot(2,1,2) plot(w,phase);title(相频特性) 实验结果： 3. clear; for n = 1:400 w(n) = (n-1)*0.05; H(n) = (1j*w(n)+1)^2/((1j*w(n))^2+0.61); end mag = abs(H); phase = angle(H); subplot(2,1,1) plot(w,mag);title(幅频特性) subplot(2,1,2) plot(w,phase);title(相频特性) 实验结果： 4. clear; for n = 1:400 w(n) = (n-1)*0.05; H(n) =3*(1j*w(n)-1)*(1j*w(n)-2)/(1j*w(n)+1)*(1j*w(n)+2); end mag = abs(H); phase = angle(H); subplot(2,1,1) plot(w,mag);title(幅频特性) subplot(2,1,2) plot(w,phase);title(相频特性) 实验结果： 实验七 三．已知下列传递函数H(s)或H(z)，求其极零点，并画出极零图。 1. z=[1,2]; p=[-1,-2]; zplane(z,p) 实验结果： 2. z=[1,2]; p=[-1,-2]; zplane(z,p) num=[1]; den=[1,0]; [z,p,k]=tf2zp(num,den)