zfirls.m

来自「用MATLAB实现数字信号处理的程序」· M 代码 · 共 14 行

%Firls函数：通过最小二乘算法设计线性相位的FIR滤波器
N=35;
F=[0 0.4 0.5 0.6 0.8 0.9];
A=[0 1 0 0 0.5 0.5];
b=firls(N,F,A,'hilbert');
for j=1:1:5
    plot([F(j) F(j+1)],[A(j) A(j+1)],'r*');
    hold on
end
[H,W]=freqz(b,1,512,2);
plot(W,abs(H),'g-');
grid on;
xlabel('频率/F'),ylabel('幅度/A');