ex0822.asv

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

% Chapter 08: Example 8.22 Design a filter
n=28;
f = [0 0.3 0.4 0.6 0.7 0.9]; 
a = [0 1 0 0 1 1];
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')