freqz_m.m

来自「FIR设计滤波！！绝对可以运行！非常的棒」· M 代码 · 共 7 行

function[db,mag,pha,w]=freqz_m(b,a);
[H,w]=freqz(b,a,59000,'whole');
H=(H(1:1:(59000/2+1)))';
w=(w(1:1:(59000/2+1)))';
mag=abs(H);
db=20*log10((mag+eps)/max(mag));
pha=angle(H);