myfreqz.m

来自「用矩形窗设计线性相位FIR低通滤波器,给了源代码」· M 代码 · 共 7 行

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