samp7_14.m

来自「数字信号处理的Matlab实现的随书光盘」· M 代码 · 共 21 行

%Samp7_14
clf;n=30;                 %滤波器阶数
f=[0 0.9];a=[0 0.9];          %理想滤波器响应
b=firls(n,f,a,'d');          %采用 firls 设计滤波器
[h,w]=freqz(b);           %计算滤波器频率响应
bb=remez(n,f,a,'d');       % 采用 remez 设计滤波器
[hh,w]=freqz(bb);           %计算滤波器脉冲响应
figure(1)
plot(w/pi,abs(h),'b',w/pi, abs(hh),'r--',f,a,'m:'); %绘幅频响应
xlabel('归一化频率');ylabel('振幅');
legend('firls', 'remez','理想特性',4);  %绘制图例，4表示图例位置在右下角
grid on ;
f1=5;                     %输入信号频率
t=0:1/1000:1;               %时间序列
x=sin(2* pi*f1* t);         %输入正弦信号
y=fftfilt(bb,x);          %滤波
figure(2)
subplot(2,1,1), plot(t,x),title ('输入信号')% 绘输入信号
subplot(2,1,2), plot(t,y),            %绘输出信号
ylim([-0.01 0.01]),xlabel('时间/s'),title('输出信号')