matlab ellipticf,[Matlab]椭圆滤波器设计：低通、高通、带通和带阻

椭圆滤波器(Elliptic filter)又称考尔滤波器(Cauer filter)：

这是在通带和阻带等波纹的一种滤波器。

椭圆滤波器相比其他类型的滤波器，在阶数相同的条件下有着最小的通带和阻带波动。

它在通带和阻带的波动相同，这一点区别于在通带和阻带都平坦的巴特沃斯滤波器，以及通带平坦、阻带等波纹或是阻带平坦、通带等波纹的切比雪夫滤波器。

测试代码：

% EllipticFilter.m

% 椭圆滤波器的设计

%

clear;

close all;

clc;

fs = 1000; %Hz 采样频率

Ts = 1/fs;

N = 1000; %序列长度

t = (0:N-1)*Ts;

delta_f = 1*fs/N;

f1 = 50;

f2 = 100;

f3 = 200;

f4 = 400;

x1 = 2*0.5*sin(2*pi*f1*t);

x2 = 2*0.5*sin(2*pi*f2*t);

x3 = 2*0.5*sin(2*pi*f3*t);

x4 = 2*0.5*sin(2*pi*f4*t);

x = x1 + x2 + x3 + x4; %待处理信号由四个分量组成

X = fftshift(abs(fft(x)))/N;

X_angle = fftshift(angle(fft(x)));

f = (-N/2:N/2-1)*delta_f;

figure(1);

subplot(3,1,1);

plot(t,x);

title('原信号');

subplot(3,1,2);

plot(f,X);

grid on;

title('原信号频谱幅度特性');

subplot(3,1,3);

plot(f,X_angle);

title('原信号频谱相位特性');

grid on;

%设计一个切比雪夫低通滤波器,要求把50Hz的频率分量保留,其他分量滤掉

wp = 55/(fs