高通滤波器(matlab编写)切比雪夫

设计一个高通滤波器，并检验它的性能

采样率为10kHZ

阻带边缘为1.5Khz,衰减为40bB

通带边缘为2kHz，波纹为3Db

采用切比雪夫2型滤波器

Fs=1e4 ; fs=1.5e3; fp=2e3; As=40; Rp=3 ;wp=2*fp/Fs; ws=2*fs/Fs;

用MA TLAB工具：

设计的滤波器性能

>>Fs=1e4;

>>fs=1.5e3;

>>fp=2e3;

>>As=40;

>>Rp=3;

>>wp=2*fp/Fs;

>>ws=2*fs/Fs;

>>[N,wn]=cheb2ord(wp,ws,Rp,As);

>>[b,a]=cheby2(N,As,wn,'high');

>>[db,mag,pha,grd,w]=freqz_m(b,a);

>> subplot(2,2,1);plot(w/pi,mag);

>> axis([0,1,0,1]);

>> setX([0 0.3 0.4 1]);

>>setY([0.01 0.7279 1])

>> title('Magnitude Response');

>> subplot(2,2,2);plot(w/pi,db);

>> axis([0 1 -70 0])

>> setX([0 0.3 0.4 1])

>> setY([-40 -2.7589])

>> title('Magnitude Response in dB');

输入信号x=cos(0.6*pi*n); 取200个抽样值，为了便于观看，在画图是裁剪掉100个值

不直观，直接做频谱图

为了说明该滤波器的性能，做一组对比。换上另一种频率的信号x=cos(0.2*pi*n);

由上图可以看到该信号完全被滤去了

可以看出在20和180点的信号几乎淹没了。

注意在此时选取点数时一定要是周期的整数倍，此题N=2*pi/0.6pi=10/3,周期是10 如果选取201点，可以看出频率已经泄露

>>n=0:200;

>>x=cos(0.6*pi*n);

>>y=filter(b,a,x);

>>subplot(2,2,3);plot(n,x);

>>subplot(2,2,4);plot(n,y);

>>x1=fft(x,201);

>>x11=abs(x1);

>>subplot(2,2,1);stem(n,x11);

>>y1=fft(y,201);

>>y11=abs(y1);

>>subplot(2,2,2);stem(n,y11);

>>setX([0 60 140 201])

>>title('FFT of y');

>>subplot(2,2,1);stem(n,x11);

>>setX([0 60 140 201])

>>title('FFT of x');

>>g=x11-y11;

>>subplot(2,2,3);stem(n,g);