matlab程序之——滤波器(带通-带阻)教学内容Word文档下载推荐.docx

matlab程序之——滤波器(带通-带阻)教学内容Word文档下载推荐.docx

《matlab程序之——滤波器(带通-带阻)教学内容Word文档下载推荐.docx》由会员分享，可在线阅读，更多相关《matlab程序之——滤波器(带通-带阻)教学内容Word文档下载推荐.docx（10页珍藏版）》请在冰豆网上搜索。

%f3：

通带右边界

%fs1：

衰减截止左边界

%fsh：

衰变截止右边界

%rp：

边带区衰减DB数设置

%rs：

截止区衰减DB数设置

%FS：

序列x的采样频率

%f1=300;

f3=500;

%通带截止频率上下限

%fsl=200;

fsh=600;

%阻带截止频率上下限

%rp=0.1;

rs=30;

%通带边衰减DB值和阻带边衰减DB值

%Fs=2000;

%采样率

%

wp1=2*pi*f1/Fs;

wp3=2*pi*f3/Fs;

wsl=2*pi*fsl/Fs;

wsh=2*pi*fsh/Fs;

wp=[wp1wp3];

ws=[wslwsh];

%设计切比雪夫滤波器；

[n,wn]=cheb1ord（ws/pi,wp/pi,rp,rs）;

[bz1,az1]=cheby1（n,rp,wp/pi）;

%查看设计滤波器的曲线[h,w]=freqz（bz1,az1,256,Fs）;

h=20*log10（abs（h））;

figure;

plot（w,h）;

title（'

所设计滤波器的通带曲线'

）;

gridon;

y=filter（bz1,az1,x）;

end

带通滤波器使用例子

%--

%带通滤波器测试程序fs=2000;

t=（fs;

ff1=100;

ff2=400;

ff3=700;

x=sin（2*pi*ff1*t）+sin（2*pi*ff2*t）+sin（2*pi*ff3*t）;

figure;

subplot（211）;

plot（t,x）;

subplot（212）;

hua_fft（x,fs,1）;

%y=filter（bz1,az1,x）;

y=bandp（x,300,500,200,600,0.1,30,fs）;

subplot（211）;

plot（t,y）;

hua_fft（y,fs,1）;

%调用到的hua_fft（）函数代码如下

functionhua_fft（y,fs,style,varargin）

%当style=1,画幅值谱；

当style=2,画功率谱;

当style=其他的，那么花幅值谱和功率谱

%当style=1时，还可以多输入2个可选参数

%可选输入参数是用来控制需要查看的频率段的

%第一个是需要查看的频率段起点

%第二个是需要查看的频率段的终点

%其他style不具备可选输入参数，如果输入发生位置错误

nfft=2^nextpow2（length（y））;

%找出大于y的个数的最大的2的指数值（自动进算最佳FFT步长nfft）

%nfft=1024;

%人为设置FFT的步长nffty=y-mean（y）;

%去除直流分量

y_ft=fft（y,nfft）;

%对y信号进行DFT，得到频率的幅值分布

y_p=y_ft.*conj（y_ft）/nfft;

%conj（）函数是求y函数的共轭复数，实数的共轭复数是他本身。

y_f=fs*（2-1）/nfft;

�T变换后对应的频率的序列

%y_p=y_ft.*conj（y_ft）/nfft;

ifstyle==1

ifnargin==3

plot（y_f,2*abs（y_ft

（2））/length（y））;

%matlab 的帮助里画FFT的方法

%ylabel（'

幅值'

xlabel（'

频率'

信号幅值谱'

%plot（y_f,abs（y_ft

（2）））;

% 论坛上画FFT的方法

elsef1=varargin{1};

fn=varargin{2};

ni=round（f1*nfft/fs+1）;

na=round（fn*nfft/fs+1）;

plot（y_f（ni:

na）,abs（y_ft（ni:

na）*2/nfft））;

end

elseifstyle==2plot（y_f,y_p

（2））;

功率谱密度'

信号功率谱'

elsesubplot（211）;

ylabel（'

plot（y_f,y_p

（2））;

ylabel（'

运行结果如下图，第一幅是滤波前测试信号的时频图，第二幅是滤波器的滤波曲线图，第三幅是经滤波后的测试信号时频图。

2.带阻滤波器

functiony=bands（x,f1,f3,fsl,fsh,rp,rs,Fs）

%带阻滤波

[bz1,az1]=cheby1（n,rp,wp/pi,'

stop'

使用例子

%带阻滤波器测试fs=1000;

y=sin（2*pi*100*t）+sin（2*pi*150*t）+sin（2*pi*200*t）;

z=bands（y,110,190,140,160,0.1,30,fs）;

hua_fft（z,fs,1）;

运行结果如下图，第一幅是滤波前测试信号的频谱图，第二幅是滤波器的滤波曲线图，第三幅是经滤波后的测试信号频谱图。