东南大学信号与系统MATLAB实践第二次作业
(18页)
本资源提供全文预览,点击全文预览即可全文预览,如果喜欢文档就下载吧,查找使用更方便哦!
19.9 积分
. . . .练习二实验六一.用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-1p = -1k = 12. >>Ns=[10]Ds=[1,-5,0]sys2=tf(Ns,Ds)实验结果:Ns = 10Ds = 1 -5 0sys2 = 10 --------- s^2 - 5 s>>[z,p,k]=tf2zp([10],[1,-5,0])z = Empty matrix: 0-by-1p = 0 5k =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 1Ds = 1 0 -1sys2 = 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 0Ds = 1 6 8sys3 = 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);endmag = 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);endmag = 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);endmag = 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);endmag = 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);zplane(z,p)>> num=[1];den=[1,0];[z,p,k]=tf2zp(num,den)zplane(z,p)实验结果:z = Empty matrix: 0-by-1p = 0k = 13. >> num=[1,0,1];den=[1,2,5];[z,p,k]=tf2zp(num,den)zplane(z,p)实验结果:z = 0 + 1.0000i 0 - 1.0000ip = -1.0000 + 2.0000i -1.0000 - 2.0000ik = 14. >> num=[1.8,1.2,1.2,3];den=[1,3,2,1];[z,p,k]=tf2zp(num,den)zplane(z,p)实验结果:z = -1.2284 0.2809 + 1.1304i 0.2809 - 1.1304ip = -2.3247 -0.3376 + 0.5623i -0.3 关 键 词: 东南大学 信号 系统 MATLAB 实践 第二次 作业
天天文库所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。