设计性实验参考程序.docVIP

１.１ k=2; b=[1]*k; a=[1 0]; sys=tf(s); %根据零极点构造系统函数H(s) sys=tf(b,a) ; tf是传递函数，在s域中，你要输入G（s）=1/（s^2+2s+1），在matlab中输入G=tf（[1],[1 2 1]）就了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) ; 求取脉冲响应的函数 xlabel(t);title(冲激响应h(t)); grid on 关于impulse的使用： 先构建系统，再直接调用，如： a = [-0.5572 -0.7814;0.7814 0]; b = [1 -1;0 2]; c = [1.9691 6.4493]; sys = ss(a,b,c,0); impulse(sys)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