设计性实验参考程序.docVIP

１.１ k=2; b=[1]*k; a=[1 0]; sys=tf(s); %根据零极点构造系统函数H(s) sys=tf(b,a) p=roots(a); z=roots(b); subplot(1,2,1); plot(real(p),imag(p), *);hold on; plot(real(z),imag(z),o); xlabel(t);title(零极点图); grid on %系统的冲激响应 T=0:0.01:3; y=impulse(b,a,T); subplot(1,2,2) plot(T,y); xlabel(t);title(冲激响应h(t)); grid on 1.2 k=2; b=[1]*k; a=conv([1 7j],[1 -7j]); sys=tf(s); %根据零极点构造系统函数H(s) sys=tf(b,a) p=roots(a); z=roots(b); subplot(1,2,1); plot(real(p),imag(p),*);hold on; plot(real(z),imag(z),o); xlabel(t);title(零极点图); grid on %系统的冲激响应 T=0:0.01:3; y=impulse(b,a,T); subplot(1,2,2); plot(T,y); xlabel(t);ti