椭圆带通滤波器设计程序

第五章仿真程序和仿真图5.1带通滤波器设计程序

Fs=100;

t=(1:100)/Fs;

s=sin(2*pi*t*5)+sin(2*pi*t*15)+sin(2*pi*t*30); figure(1);

subplot(111);

plot(t,s);

xlabel('时间(秒)');

ylabel('幅值');

[b,a]=ellip(6,0.1,40,[10 20]*2/Fs);

[H,w]=freqz(b,a,512);

figure(2);

subplot(111);

plot(w*Fs/(2*pi),abs(H));

xlabel('频率 (Hz)');ylabel('频率响应图'); grid;

sf=filter(b,a,s);

figure(3);

subplot(111);

plot(t,sf);

xlabel('时间 (s)');

ylabel('幅值');

axis([0 1 -1 1]);

S=fft(s,512);

SF=fft(sf,512);

w=(0:255)/256*(Fs/2);

figure(4);

subplot(111);

plot(w,abs([S(1:256)' SF(1:256)'])); xlabel('频率(Hz)');

ylabel('傅立叶变换图');

grid;

legend({'滤波前','滤波后'});

5.2信号的仿真图

1. 由Matlab软件可实现指定信号的输入，其波形如下图

图5.1连续信号波形图

2.这个信号包含三个分量，分别为5Hz，15Hz，30Hz的正弦波，其频谱如下图：

图5.2 连续信号的频谱图

3.通过Ellip函数设计出椭圆带通滤波器，它的频率响应图如下：

图5.3 椭圆带通滤波器频率响应

4.信号通过椭圆带通滤波器的仿真图，如下图所示

图5.4 信号通过滤波器后的频谱图

注：图中蓝色曲线代表滤波前的幅频曲线，绿线代表滤波后的幅频曲线。