MATLAB图像处理-椭圆滤波算法

滤波器：

滤波器是一种二端口网络。它具有选择频率的特性，即可以让某些频率顺利通过，而对其它频率则加以阻拦，目前由于在雷达、微波、通讯等部门，多频率工作越来越普遍，对分隔频率的要求也相应提高；所以需用大量的滤波器。又如声卡中的语音合成输出，再如我们学校试验室中的合成信号发生器等，为了滤除谐波干扰，获得高精度的模拟信号，大多数要求衰减特性陡峭，而这种特性椭圆滤波器具有，因此椭圆滤波器作用很大，用途广泛。

 

椭圆滤波器的幅频特性

滤波器的比较：

 

椭圆滤波器程序：

1、信号采样

Fs=100; t=(1:100)/Fs; s1=sin(2*pi*t*5); s2=sin(2*pi*t*15); s3=sin(2*pi*t*30); s=s1+s2+s3; % 分别绘制四个函数 figure(1) plot(t,s,'r',t,s1,'m',t,s2,'g',t,s3,'b'); xlabel('时间(s)'); ylabel('幅值'); % 为绘制的曲线添加标签 legend({'s曲线','s1曲线','s2曲线','s3曲线'});

2、滤波器设计

[b,a]=ellip(5,0.1,40,5*2/Fs); % 求取离散系统频率响应函数 [H,w]=freqz(b,a,512); figure(2) % 绘制椭圆滤波器频率响应图 plot(w*Fs/(2*pi),abs(H)); xlabel('频率(Hz)'); ylabel('频率响应图');

3、信号进行分析和变换

S=fft(s,512); SF=fft(sf,512); w=(0:255)/256*(Fs/2); % 绘制实验结果图 figure(3) plot(w,abs(S(1:256)'),'r',w,abs( SF(1:256)'),'b'); xlabel('频率(Hz)');ylabel('傅立叶变换图'); grid; legend({'滤波前','滤波后'});