%--------------------------------------------------------------------------
% exa070903_remez_1.m,   for example 7.9.3 and 7.4.1;
% to test remez.m and to design Low-pass FIR filter;
%-------------------------------------------------------------------------
clear all;

f=[0 .6 .7 1];
% 給定頻率軸分點；
A=[1 1  0 0];
% 給定在這些頻率分點上理想的幅頻響應；
weigh=[1 10];	
% 給定在這些頻率分點上的加權；

b=remez(32,f,A,weigh);
% 設計出切比雪夫最佳一致逼近濾波器；
%
[h,w]=freqz(b,1,256,1);
h=abs(h);
h=20*log10(h);
figure(1)
stem(b,'.');grid;
figure(2)
plot(w,h);grid;