ellipb.m

来自「使用Matlab实现了FIR与IIR滤波器」· M 代码 · 共 18 行

function [ output_args ] = Untitled8( input_args )
%UNTITLED8 Summary of this function goes here
%  Detailed explanation goes here
ft=22050;
Wp1 =2*1200/ft;
Wp2=2*3000/ft;
Ws1 = 2*1000/ft;
Ws2=2*3200/ft;
Wp=[Wp1,Wp2];
Ws=[Ws1,Ws2];
Rp = 1;
Rs = 100;
[N,Wn] = ellipord(Wp,Ws,Rp,Rs);
[b,a] = ellip(N,Rp,Rs,Wn);
[h,w]= freqz(b,a);
plot (w/pi,20*log10(abs(h)));grid;
xlabel('w/pi'); ylabel('Gain, dB');
title('Type Ellip Bandpass Filter');