用MATLAB设计低通带通高通.docx
- 文档编号:2482170
- 上传时间:2022-10-30
- 格式:DOCX
- 页数:18
- 大小:26.05KB
用MATLAB设计低通带通高通.docx
《用MATLAB设计低通带通高通.docx》由会员分享,可在线阅读,更多相关《用MATLAB设计低通带通高通.docx(18页珍藏版)》请在冰豆网上搜索。
用MATLAB设计低通带通高通
用MATLAB设计低通,带通,高通和带阻FIR数字滤波器急!
!
!
!
悬赏分:
30-解决时间:
2006-11-120:
56
(1)低通滤波器的技术指标:
H(ejw)=1,0 (2)带通滤波器的技术指标: H(ejw=0,0 (3)高通滤波器的技术指标: H(ejw)=0,0 (4)带阻滤波器的技术指标: H(ejw)=1,0 提问者: wwce520-魔法学徒一级 最佳答案 低通采样定理实验 1.1实验目的 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转换。 数字信号处理后的数据可能需还原为模拟信号,这就需要进行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)的成分,可以通过适当的低通滤波器不失真地恢复出x(t)。 最小抽样速率fs=2fx称为奈奎斯特速率。 低通 译码 编码 量化 抽样 输入信号样点输出滤波输出 A/D(模数转换)D/A(数模转换) 图1低通采样定理演示 为方便实现,实验中更换了一种表现形式,即抽样频率固定(10KHz),通过改变输入模拟信号的频率来展示低通抽样定理。 我们可以通过研究抽样频率和模拟信号最高频率分量的频率之间的关系,来验证低通抽样定理。 1.4实验内容 1.软件仿真实验: 编写并调试MATLAB程序,分析有关参数,记录有关波形。 2.硬件实验: 输入不同频率的正弦信号,观察采样时钟波形、输入信号波形、样点输出波形和滤波输出波形。 1.5MATLAB参考程序和仿真内容 %*******************************************************************% %f—余弦信号的频率 %M—基2FFT幂次数N=2^M为采样点数,这样取值是为了便于作基2的FFT分析 %2.采样频率Fs %*******************************************************************% 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([0N1.1*min(Xk)1.1*max(Xk)]); xlabel('frequency-->'); ylabel('! Xk! '); %*******************************************************************% 假如有一个1Hz的余弦信号y=cos(2*π*t),对其用4Hz的采样频率进行采样,共采样32点,只需执行samples(1,4,5),即可得到仿真结果。 软件仿真实验内容如下表所示: 仿真参数 f Fs Wo(计算) Xn(图形) Xk(图形) (1,4,5) 另外记录图形,并标图号 (1,8,5) (2,8,6) 自选 1.6硬件实验步骤 本实验箱采样频率fs固定为10KHz,低通滤波器的截止频率约为4.5KHz。 1、用低频信号源产生正弦信号,正弦信号源频率f自定,并将其接至2TP2(模拟输入)端,将示波器通道一探头接至2TP6(采样时钟)端观察采样时钟波形,示波器通道二探头接至2TP2观察并记录输入信号波形。 2、将示波器通道二探头接至2TP3观察并记录样点输出波形。 3、将示波器通道二探头接至2TP4观察并记录滤波输出波形。 4、根据采样定理,分f=fs/8、f=fs/4、f=fs/2等3种情况更改正弦信号频率,重复步骤2至步骤3。 5、用低频信号源产生方波信号,重复步骤1至步骤4。 1.7思考题 1、讨论在仿真实验中所计算的数字域频率Wo和Xk的图形中非零谱线位置之间的对应关系。 2、讨论在仿真实验中自选参数的意义。 3、将在2TP2端加方波信号后的恢复波形,与相同频率的正弦信号的恢复波形相比,能够得出哪些结论? 2FFT频谱分析实验 2.1实验目的 1.通过实验加深对快速傅立叶变换(FFT)基本原理的理解。 2.了解FFT点数与频谱分辨率的关系,以及两种加长序列FFT与原序列FFT的关系。 2.2实验仪器 1.YBLD智能综合信号源测试仪1台 2.双踪示波器1台 3.MCOM-TG305数字信号处理与现代通信技术实验箱1台 4.PC机(装有MATLAB、MCOM-TG305配套实验软件)1台 2.3实验原理 离散傅里叶变换(DFT)和卷积是信号处理中两个最基本也是最常用的运算,它们涉及到信号与系统的分析与综合这一广泛的信号处理领域。 实际上卷积与DFT之间有着互通的联系: 卷积可化为DFT来实现,其它的许多算法,如相关、滤波和谱估计等都可化为DFT来实现,DFT也可化为卷积来实现。 对N点序列x(n),其DFT变换对定义为: 在DFT运算中包含大量的重复运算。 FFT算法利用了蝶形因子WN的周期性和对称性,从而加快了运算的速度。 FFT算法将长序列的DFT分解为短序列的DFT。 N点的DFT先分解为2个N/2点的DFT,每个N/2点的DFT又分解为2个N/4点的DFT。 按照此规律,最小变换的点数即所谓的“基数(radix)。 ”因此,基数为2的FFT算法的最小变换(或称蝶形)是2点DFT。 一般地,对N点FFT,对应于N个输入样值,有N个频域样值与之对应。 一般而言,FFT算法可以分为时间抽取(DIT)FFT和频率抽取(DIF)两大类。 在实际计算中,可以采用在原来序列后面补0的加长方法来提高FFT的分辨率;可以采用在原来序列后面重复的加长方法来增加FFT的幅度。 2.4实验内容 1.软件仿真实验: 分别观察并记录正弦序列、方波序列及改变FFT的点数后的频谱;分别观察并记录正弦序列、方波序列及2种加长序列等信号的频谱。 2.硬件实验: 分别观察并记录正弦信号、方波信号及改变FFT的点数后的频谱。 2.5MATLAB参考程序和仿真内容 %*******************************************************************% function[x]=ffts(mode,M) Nfft=2^M; x=zeros(1,Nfft);%定义一个长度为Nfft的一维全0数组 ifmode==1forn=0: Nfft-1x(n+1)=sin(2*pi*n/Nfft);end end%定义一个长度为Nfft的单周期正弦序列 ifmode==2forn=0: Nfft-1x(n+1)=sin(4*pi*n/Nfft);end end%定义一个长度为Nfft的双周期正弦序列 ifmode==3forn=0: Nfft/2-1x(n+1)=sin(4*pi*n/Nfft);end end%定义一个长度为Nfft/2的正弦序列,后面一半为0序列。 ifmode==4forn=0: Nfft-1x(n+1)=square(2*pi*n/Nfft);end end ifmode==5forn=0: Nfft-1x(n+1)=square(2*pi*n/Nfft);end end ifmode==6forn=0: Nfft/2-1x(n+1)=square(4*pi*n/Nfft);end end n=0: Nfft-1; subplot(2,1,1); stem(n,x); axis([0Nfft-11.1*min(x)1.1*max(x)]); xlabel('Points-->'); ylabel('x(n)'); X=abs(fft(x,Nfft)); subplot(2,1,2); stem(n,X); axis([0Nfft-11.1*min(X)1.1*max(X)]); xlabel('frequency-->'); ylabel('! X(k)! '); %*******************************************************************% 假设需观察方波信号的频谱,对一个周期的方波信号作32点的FFT,则只需在MATLAB的命令窗口下键入: [x]=ffts(21,5),程序进行模拟,并且输出FFT的结果。 关于软件仿真实验内容,建议在完成大量仿真例子的基础上,选择能够体现实验要求的4个以上的例子进行记录。 例如要观察后面补0的加长方法来提高FFT的分辨率的现象,可以仿真ffts(4,5)和ffts(6,6)两个例子。 2.6硬件实验步骤 1.将低频信号源输出加到实验箱模拟通道1输入端,将示波器探头接至模拟通道1输出端。 2.在保证实验箱正确加电且串口电
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- MATLAB 设计 低通带通高通
![提示](https://static.bdocx.com/images/bang_tan.gif)