用matlab设计低通带通高通和带阻fir数字滤波器

《用matlab设计低通带通高通和带阻fir数字滤波器》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。

1、抽样频率为f为1KHZ的数字带通滤波器，性能要求为：通带范围从200HZ到250HZ，在此两频率处衰减不大于3dB，在100HZ和400HZ频率处频率衰减不得小于20dB，采用巴特沃斯滤波器4阶Butterworth滤波器源代码n=2;Wn=[200250]/500;[b,a]=butter(n,Wn);freqz(b,a,512,1000);这个滤波器100Hz、400Hz处衰减在-30db以上。可以通过增大n增加衰减。用MATLAB设计低通带通高通和带阻FIR数字滤波器（1）低通滤波器的技术指标

2、：H(ejw)=1,0

3、验目的1．了解数字信号处理系统的一般构成；2．掌握奈奎斯特抽样定理。1.2实验仪器1．YBLD智能综合信号源测试仪1台2．双踪示波器1台3．MCOM－TG305数字信号处理与现代通信技术实验箱1台4．PC机（装有MATLAB、MCOM－TG305配套实验软件）1台1.3实验原理一个典型的DSP系统除了数字信号处理部分外，还包括A/D和D/A两部分。这是因为自然界的信号，如声音、图像等大多是模拟信号，因此需要将其数字化后进行数字信号处理，模拟信号的数字化即称为A/D转换。数字信号处理后的数据可能需还原

4、为模拟信号，这就需要进行D/A转换。一个仅包括A/D和D/A两部分的简化数字信号处理系统功能如图1所示。A/D转换包括三个紧密相关的过程，即抽样、量化和编码。A/D转换中需解决的以下几个重要问题：抽样后输出信号中还有没有原始信号的信息？如果有能不能把它取出来？抽样频率应该如何选择？奈奎斯特抽样定理（即低通信号的均匀抽样定理）告诉我们，一个频带限制在0至fx以内的低通信号x(t)，如果以fs≥2fx的抽样速率进行均匀抽样，则x(t)可以由抽样后的信号xs(t)完全地确定，即xs(t)包含有x(t)的成

5、分，可以通过适当的低通滤波器不失真地恢复出x(t)。最小抽样速率fs=2fx称为奈奎斯特速率。低通译码编码量化抽样输入信号样点输出滤波输出A/D（模数转换）D/A（数模转换）图1低通采样定理演示为方便实现，实验中更换了一种表现形式，即抽样频率固定（10KHz），通过改变输入模拟信号的频率来展示低通抽样定理。我们可以通过研究抽样频率和模拟信号最高频率分量的频率之间的关系，来验证低通抽样定理。1.4实验内容1．软件仿真实验：编写并调试MATLAB程序，分析有关参数，记录有关波形。2．硬件实验：输入不同频

6、率的正弦信号，观察采样时钟波形、输入信号波形、样点输出波形和滤波输出波形。1.5MATLAB参考程序和仿真内容%*******************************************************************%%f—余弦信号的频率%M—基2FFT幂次数N=2^M为采样点数，这样取值是为了便于作基2的FFT分析%2.采样频率Fs%*******************************************************************%

7、functionsamples(f,Fs,M)N=2^M;%fft点数=取样总点数Ts=1/Fs;%取样时间间隔T=N*Ts;%取样总时间=取样总点数*取样时间间隔n=0:N-1;t=n*Ts;Xn=cos(2*f*pi*t);subplot(2,1,1);stem(t,Xn);axis([0T1.1*min(Xn)1.1*max(Xn)]);xlabel('t-->');ylabel('Xn');Xk=abs(fft(Xn,N));subplot(2,1,2);stem(n,Xk);axis([0

8、N1.1*min(Xk)1.1*max(Xk)]);xlabel('frequency-->');ylabel('!Xk!');%*******************************************************************%假如有一个1Hz的余弦信号y=cos(2*π*t)，对其用4Hz的采样频率进行采样，共采样32点，只需执行samples(1,4,5)，即可得到仿真结果。软件仿真实验内容如下表所示：仿真参数fFsWo(计算