|
加个循环行不行？for N=1:6[z,p,k]=buttap(N); [b,a]=zp2tf(z,p,k); [H,w]=freqs(b,a); magH2=abs(H).^2; hold on; figure(1) plot(w,magH2); xlabel('w/wc'); ylabel('|H(jw)|^2'); title('Butterword滤波器幅频特性'); grid on; endfigure(2) hold on; [z,p,k]=butter(N,0.2); zplane(z,p); xlabel('实部');ylabel('虚部');grid on; |
