4.1
num=[5];den=[1 1 44];[r,p,k]=residue(num,den); 结果:
r = -0.5000 - 0.2500i -0.5000 + 0.2500i 1.0000 + 0.0000i p =-0.0000 + 2.0000i -0.0000 - 2.0000i -1.0000 + 0.0000i k = [] 4.4
a)由图可知: H(s)=1/(s+2); num=[1];den=[1 2];
poles=roots(den) ;t=0:0.02:10; h=impulse(num,den,t); plot(t,h);
xlabel('t(s)');ylabel('h(t)');grid on; [H,w]=freqs(num,den);
figure(2);plot(w,abs(H));
xlabel('\\omega(rag/s)');ylabel('|H(j\\omega)|'); title('Magenitude Response');grid on; (b)由图可知:H(s)=1/((s+4j)(s-4j)); num=[1];den=conv([1 -4j],[1 4j]); poles=roots(den) t=0:0.02:10;
h=impulse(num,den,t); figure(1);plot(t,h);
xlabel('t(s)');ylabel('h(t)');grid on; [H,w]=freqs(num,den);
figure(2);plot(w,abs(H));
xlabel('\\omega(rag/s)');ylabel('|H(j\\omega)|');
title('Magenitude Response');grid on;