ex0824.asv

数字信号处理（机械工业出版社）的源码

% Chapter 08: Example 8.24 Design a filter
n=135;
f = [0 0.2 0.4 0.7 0.9 1]; 
a = [0 0.5 0 1 0];
up=[0.005 0.51 0.03 1.02 0.05]
lo=
b=firls(n,f,a,'hilbert');
[h,w]=freqz(b,1,512,2);
plot(w,abs(h),f,a,':'); grid on
xlabel('Normalized Frequency'); ylabel('magnitude(dB)');
legend('Ideal','firls Design')