%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');