语音信号的FIR滤波器处理华南农业大学

（农业畜牧行业）语音信号的FIR滤波器处理(华南农业大学)摘要DSP(DigitalSignalProcessing)也就是我们常说的数字信号处理,它是利用计感谢阅读算机或专用处理设备,以数字形式对信号进行采集变换滤波估值增强压谢谢阅读缩识别等处理,以得到符合人们需要的信号形式DSP芯片在信号处理通信、感谢阅读雷达等许多领域得到广泛的应用。Matlab矩阵计算科学数据可视化等诸精品文档放心下载多强大功能集成在一个易于使用的视窗环境中，为科学研究工程设计以及必须谢谢阅读进行有效数值计算的众多科学领域提供了一种全面的解决方案。感谢阅读CSS集成开发环境使用CCS内置的软件仿真simulator对程序进行编译，谢谢阅读,并能通过仿真器与目谢谢阅读标板连接，在目标板上实时观察效果。在本次设计中，我们选择的课题是基于DSP的语言信号的FIR滤波处理。精品文档放心下载首先利用MATLAB进行了仿真，得到滤波前后的时域波形和频谱。然后通过调谢谢阅读用MATLAB的分析工具FDATOOL，根据仿真结果导出了滤波器的相关参数，谢谢阅读将原始信号数据和滤波器参数输入CCS进行DSP编程最后在DSP中实现了谢谢阅读FIR带通滤波，并通过CCS的频谱分析功能查看了最终DSP的滤波效果。精品文档放心下载关键词:语音信号DSPFIR滤波MATLABCCS感谢阅读目录1.前言42.设计原理42.1数字信号处理器42.2滤波器简介52.3fir滤波器原理62.4窗函数简介73.FIR滤波器的MATLAB设计83.1总体方案的设计83.2语音信号的采集93.3MATLAP的具体处理103.5FIR滤波器的系统参数144.FIR滤波器的CCS设计与仿真164.1CCS程序流程图174.2CCS仿真结果205.问题及解决办法226.设计感想23参考文献237.附录241前言精品文档放心下载技术领域之一它在通信语音图像自动控制雷达军事航空航天医精品文档放心下载疗和家用电器等众多领域得到了广泛的应用在数字信号处理的基本方法中，通感谢阅读常会涉及到变换滤波频谱分析调制解调和编码解码等处理其中滤波是应谢谢阅读用非常广泛的一个环节，数字滤波器的理论和相关设计也一直都是人们研究的重感谢阅读点之一。FIR滤波器的是非递归的，稳定性好，精度高；更重要的是，FIR滤波感谢阅读器在满足幅频响应要求的同时，可以获得严格的线性相位特征因此，它在高保感谢阅读真的信号处理，如数字音频图像处理数据传输和生物医学等领域得到广泛应谢谢阅读用。在数字信号处理中，滤波占有极其重要的地位。数字滤波是语音信号处理、精品文档放心下载图像处理、模式识别、频谱分析等应用的基本处理算法。用DSP芯片实现数字谢谢阅读滤波除了具有稳定性好精确度高不受环境影响等优点外，还具有灵活性好等谢谢阅读特点。本文介绍了数字滤波器的设计基础及用窗函数法设计FIR滤波器的方法，感谢阅读运用MATHLAB语言实现了低通滤波器的设计并用CCS3.3进行观察效果。谢谢阅读2设计原理2.1数字信号处理器DSP（DigitalSignalProcessor）是一种独特的微处理器，是以数字信号来感谢阅读处理大量信息的器件。其工作原理是接收模拟信号，转换为0或1的数字信号。谢谢阅读再对数字信号进行修改删除强化，并在其他系统芯片中把数字数据解译回模感谢阅读拟数据或实际环境格式它不仅具有可编程性，而且其实时运行速度可达每秒数谢谢阅读以千万条复杂指令程序，远远超过通用微处理器，是数字化电子世界中日益重要精品文档放心下载的电脑芯片。它的强大数据处理能力和高运行速度，是最值得称道的两大特色。感谢阅读由Ti公司提供专业的开发工具CCSDSP/BIOS谢谢阅读编写适合DSP开发工程及文件，满足DSP程序设计要求。谢谢阅读由MathWorks公司和TI公司联合开发的精品文档放心下载DSPMATLABLinkforCCSDevelopmentTools（简称CCSLink）是MATLAB6.5谢谢阅读Release13MATLABCCS和DSP精品文档放心下载MATLAB变量一样来操作DSP器件的存感谢阅读MATLAB环境下完成对DSP感谢阅读提高DSP应用系统的开发进程。本设计主要是使用DSP对语音噪声信号进行处理，使用CCS开发工具编谢谢阅读写程序，同时利用MATLAB操作DSP器件的存储器和寄存器。谢谢阅读2.2滤波器简介“波”是一个非常广泛的物理概谢谢阅读念，在电子技术领域，“波”被狭义地局限于特指描述各种物理量的取值随时间起精品文档放心下载伏变化的过程谢谢阅读因为自变量时间‘是连续取值的，精品文档放心下载所以称之为连续时间信号，又习惯地称之为模拟信号(AnalogSignal)。随着数字精品文档放心下载式电子计算机(一般简称计算机)技术的产生和飞速发展，为了便于计算机对信号精品文档放心下载进行处理，产生了在抽样定理指导下将连续时间信号变换成离散时间信号的完整谢谢阅读的理论和方法也就是说，可以只用原模拟信号在一系列离散时间坐标点上的样精品文档放心下载本值表达原始信号而不丢失任何信息，波波形信号这些概念既然表达的是客感谢阅读观世界中各种物理量的变化，自然就是现代社会赖以生存的各种信息的载体信谢谢阅读息需要传播，靠的就是波形信号的传递信号在它的产生转换传输的每一个感谢阅读环节都可能由于环境和干扰的存在而畸变，有时，甚至是在相当多的情况下，这精品文档放心下载种畸变还很严重，以致信号及其所携带的信息被深深地埋在噪声当中了。谢谢阅读感谢阅读的过程。滤波器特性可以用其频率响应来描述，按其特性的不同，可以分为低通滤波器，精品文档放心下载高通滤波器，带通滤波器和带阻滤波器等。带通滤波器（英语：Band-passfilter）是指能通过某一频率范围内的频率分量、精品文档放心下载但将其他范围的频率分量衰减到极低水平的滤波器，与带阻滤波器的概念相对。感谢阅读一个模拟带通滤波器的例子是电阻-电感-电容电路(RLCcircuit)。这些滤波器也可精品文档放心下载以用低通滤波器同高通滤波器组合来产生。2.3FIR滤波器原理FIR滤波器(有限长单位冲激响应滤波器)是在数字信号处理（DSP）中经常精品文档放心下载IIR滤波器IIR滤波器是无限冲激响应精品文档放心下载滤波器。h(n)感谢阅读能在时域中实现频域中的滤波。在频域，当其输入信号为X(ejω)时，如滤波器的频率响应为H(ejω)，则其精品文档放心下载输出信号为Y(ejω)=X(ejω)H(ejω)。精品文档放心下载在时域，设滤波器的单位冲激响应h(n)为一N点序列，即0≤n≤N-1时h(n)

x（n应y（n）表示为x（n）与h(n)的卷积和，即：即求出时域的谢谢阅读h(n)后，便可通过卷积来实现频域的滤波。卷积和运算主要有以下几个步骤：精品文档放心下载(1)h(n)序列N个点数值的存储由于h(n)是根据滤波性能要求已经设计好的有限长单位冲激响应，故其N精品文档放心下载个点的数值是已知的，因此可以存放在ROM或RAM当中，且对应着N个不同谢谢阅读的地址，便于寻址。(2)输入序列x（n）的移位寄存输入序列x（n）是不断变化的，因此只能对其进行移位寄存，寄存器的个数为谢谢阅读NN个寄存器中分别存放着x（nx（n-1）……x（n-N+1n感谢阅读的变化而变化。(3)乘法器用以完成两个数值的乘法，即h(m)x（n-mN地址所对应精品文档放心下载的N个固定数值h(m)分别与N个移位寄存器中的不断变化的N个变化数值x谢谢阅读（n-m）相乘。(4)累加器用以实现N个乘积的累加，即将当前x（n）所对应的N个乘积进行累加，感谢阅读所得到的和就是y（n当滤波器的下一个输入值即x（n+1精品文档放心下载零，并重新将下一组x（n+1）所对应的N个乘积进行累加，所得到的和就是y感谢阅读（n+1。2.4窗函数简介数字信号处理的主要数学工具是傅里叶变换。而傅里叶变换是研究整个时间精品文档放心下载域和频率域的关系不过，当运用计算机实现工程测试信号处理时，不可能对无谢谢阅读限长的信号进行测量和运算，而是取其有限的时间片段进行分析做法是从信号谢谢阅读中截取一个时间片段，然后用观察的信号时间片段进行周期延拓处理，得到虚拟谢谢阅读的无限长的信号，然后就可以对信号进行傅里叶变换相关分析等数学处理无谢谢阅读线长的信号被截断以后，其频谱发生了畸变，原来集中在f(0)处的能量被分散到精品文档放心下载两个较宽的频带中去了（这种现象称之为频谱能量泄漏）。精品文档放心下载为了减少频谱能量泄漏，可采用不同的截取函数对信号进行截短，截断函数精品文档放心下载称为窗函数，简称为窗信号截短以后产生的能量泄漏现象是必然的，因为窗函谢谢阅读数w(t)是一个频带无限的函数，所以即使原信号x(t)是有限带宽信号，而在截短谢谢阅读以后也必然成为无限带宽的函数，即信号在频域的能量与分布被扩展了又从采精品文档放心下载样定理可知，无论采样频率多高，只要信号一经截短，就不可避免地引起混叠，谢谢阅读因此信号截短必然导致一些误差。泄漏与窗函数频谱的两侧旁瓣有关，如果两侧瓣的高度趋于零，而使能量相精品文档放心下载对集中在主瓣，就可以较为接近于真实的频谱，为此，在时间域中可采用不同的谢谢阅读窗函数来截短信号。3FIR滤波器的Matlab设计3.1.总体方案设计本次课程设计所采用的主要软件为Matlab。设计中用到的函数主要有谢谢阅读wavread,sound,fft,subplot,plot,fir1,freqz,filter。谢谢阅读课程设计利用窗函数法设计FIR数字滤波器，并利用MATLAB作为辅助工具仿精品文档放心下载真滤波前后的时域波形和频谱对于DSPMatlab给语音信谢谢阅读号添加高频噪音，即将原始信号调制到某一高频上（本设计用10000hz感谢阅读matlab谢谢阅读并存入后缀为dat的文件中。之后根据仿真结果使用MATLAB的设计分析工具感谢阅读FDATOOL产生滤波器参数。将混叠后的信号数据和滤波器参数输入CCS进行

DSP编程，在DSP中实现FIR带通滤波，并使用CCS的频谱分析功能，查看感谢阅读最终DSP的滤波效果。3.2语音信号的采集使用PC机上的声卡和WINDOWS操作系统可以进行数字信号的采集将话感谢阅读筒插入计算机的语音输入插口上,启动MOOo录音专家按下录音按钮,接着对话谢谢阅读筒录音,录音结束后屏幕左侧将显示所录声音的长度。点击放音按钮,可以实现所谢谢阅读录音的重现如图3.1“voice”保存入matlab当前工作路径中可精品文档放心下载以看到,文件存储器的后缀默认为.wav,这是WINDOWS操作系统规定的声音文感谢阅读件存的标准。3.3matlab具体处理3.3.1原始音频抽样Fs=10000;bits=16;谢谢阅读T=1;ts=1/Fs;N=38;[x,Fs,bits]=wavread('');%x=x(:,1);m=length(x);%sound(x,Fs,bits);谢谢阅读y=fft(x,m);%对录制的信号进行傅立叶变换f=(Fs/m)*[1:m];t=[1:m]/Fs;figure(1);精品文档放心下载subplot(211);plot(t,x);title('原始信号波形');xlabel('time(s)');subplot(212);plot(f,abs(y));title('原始信号频谱');xlabel('frequency(hz)');noise4=0.01*sin(2*pi*5000*t)+0.01*sin(2*pi*200*t);带通噪声%精品文档放心下载3.3.2加噪声并带通滤波fpl4=500;fph4=2000;fsl4=700;fsh4=1800;wp4=[2*fpl4/Fs2*fph4/Fs];谢谢阅读ws4=[2*fsl4/Fs2*fsh4/Fs];X4=x+noise4';%加了带通噪声的信号%sound(X4,Fs,bits);%回放噪声信号精品文档放心下载Y41=fft(X4,m);%对加了带通噪声的信号进行傅立叶变换figure(2);精品文档放心下载subplot(211);plot(t,X4);title('加了带通噪声的信号时域图');xlabel('time(s)');

subplot(212);plot(f,abs(Y41));title('加了带通噪声的信号频谱图

');xlabel('frequency(hz)');谢谢阅读figure(3);wc4=(wp4+ws4)/2;感谢阅读b4=fir1(N,wc4,kaiser(N+1));%调用kaiser计算带通FIRDF的b4（n）[hn4wp4]=freqz(b4,1,512);谢谢阅读subplot(211);plot(wp4/pi,20*log10(abs(hn4)));title('带通滤波器幅频');xlabel('w');

subplot(212);plot(wp4/pi,unwrap(angle(hn4)));title('带通滤波器相频');xlabel('w');

figure(4);x4get=filter(b4,1,X4);Y4get=fft(x4get,m);%sound(x4get,Fs,bits);谢谢阅读subplot(211);plot(t,x4get);title('加了带通噪声的信号滤波后时域图');xlabel('time(s)');谢谢阅读subplot(212);plot(f,abs(Y4get));title('加了带通噪声的信号滤波后频谱图');xlabel('frequency(hz)');感谢阅读3.3.3导出混合信号供DSP处理xout=X4/max(X4);xto_ccs=round(32767*xout);感谢阅读%由于TMS32054X是16FIR滤波器系精品文档放心下载数需要采用Q15格式。myfile=fopen('fir1.dat','w');%新建fprintf(myfile,'16511e5112c\n');感谢阅读fprintf(myfile,'%d\n',xto_ccs);%把调制信号+原始信号输出文件fclose(myfile);感谢阅读3.5FIR滤波器的系统参数在本设计中，使用MATLAB模拟产生合成信号，然后利用CCS进行滤波。谢谢阅读因此设计一个FIR带通滤波器，其参数为：采样频率：Fs=44100Hz通带1：Fp=700Hz阻带1：Fc=500Hz通带2：Fp=1800Hz阻带2：Fp=2000Hz通带最大衰减：1dB阻带最少衰减：100dB滤波器级数：N=38（由MATLAB根据前述参数求得）感谢阅读其中，采样频率由[y,fs,bits]=wavread()求出。其中fs为采样频率。执行程谢谢阅读序可得fs为22050HZ滤波器级数由%据ωｐ１＝２*ｆｐ１／ｆｓ；ωc１＝２*精品文档放心下载ｆcωｐ１－ωc精品文档放心下载38。滤波器系数使用matlab滤波器的设计/分析工具FDATOOL设计FIR滤波谢谢阅读器。步骤如下图所示：图1Matlab声音信号调制示意图图2Matlab滤波器设计”DesignFilter”按钮完成滤波器设感谢阅读计。设计成功后的结果如图下图所示。在fdatool中，选择精品文档放心下载Targets->CodeComposerStudio™IDE在出现的对话框中选择输出文件类型为

C.headerfile,输出系数类型为signed16-bitinteger11—13点击Generate感谢阅读按钮，选择路径，即可输出前一步设计出的FIR滤波器系数表。谢谢阅读Error!Nobookmarknamegiven.4FIR滤波器的CCS设计与仿真精品文档放心下载CCS是一种针对TMS320系列DSP的集成开发环境，在windows操作系精品文档放心下载统下，采用图形接口界面，提供环境配置源文件编辑程序调试跟踪和分析谢谢阅读等工具，提供了配置建立调试跟踪和分析程序的工具，从而完成编辑编谢谢阅读译、链接、调试和数据分析等工作。在CCSFIR滤波的功能本课程设精品文档放心下载计使用CCS开发应用程序的一般步骤为：（1）创建一个工程项目文件fir2.h，用汇编语言编写处理主程序fir2.asm。另外谢谢阅读根据板上的存储器配置方式,编写存储器配置文件文件,将MATLAB生成的fir2.h精品文档放心下载和文件拷贝到工程文件夹下。（2。使用CCS提供的集成编辑环境，对头文件、链接命令文件感谢阅读和源程序进行编辑。（3）对工程项目进行编译，生成可执行文件。如出现语法错误，将在构建窗口感谢阅读中显示错误信息。用户可以根据显示信息定位错误信息，更改错误。谢谢阅读（4）下载程序、输入数据，执行程序，对结果和数据进行分析和算法评估。利精品文档放心下载用CCS提供的探测点、图形显示等工具，对运行结果、输出数据进行分析，评感谢阅读估算法性能。4.1CCS程序流程图4.1.1程序分析FIR滤波器的DSP实现主要有线性缓冲区法和循环缓冲区法两种。感谢阅读线性缓冲区法又称延迟线法。其特点：（1）对于N级的FIR滤波器，在数据存储器中开辟一个N单元的缓冲区（滑精品文档放心下载N个输入样本；（2）从最老样本开始取数，每取一个样本后，将此样本向下移位；谢谢阅读（3）读完最后一个样本后，输入最新样本存入缓冲区的顶部。感谢阅读循环缓冲区法的特点如下：（1N级FIR滤波器，在数据存储器中开辟一个N谢谢阅读用来存放最新的N个输入样本；（2）从最新样本开始取数；（3）读完最后一个样本（最老样本）后，输入最新样本来代替最老样本，而其谢谢阅读他数据位置不变；（4BK感谢阅读区地址首尾相邻。本次设计的FIR滤波器所采用的就是循环缓冲区法。谢谢阅读对DSP进行初始化，并设置好存储单元。为数字滤波作准备，将预先设计精品文档放心下载好的有N个抽头的FIR数字滤波器的冲激响应序列h(n)中的N个数值放入存储精品文档放心下载单元bnindataxn精品文档放心下载中。主要循环程序如下：STM#-1,AR0RPTZA,#Ｎ;将ACC累加器清零，循环Ｎ次MAC*AR2+0%,*AR3+0%,ASTHA,*(yn);保存输出MVKD*(yn),*AR5+将累加器清零，将存储单元bn与xn值对应相乘并累加到ACC中，循环Ｎ谢谢阅读次将ACC中的值保存输出。再从indata数据存储段中读入一个新的数据，代替谢谢阅读原来的最老样本MVDD*AR4+,*AR2+0%,重复循环直到数据全部处理完。精品文档放心下载最后将处理完的数据输入到outdata数据存储段中。精品文档放心下载4.1.2数据输入MATLAB辅助DSP实现FIR滤波器滤波,其总体过程是为DSP中编写处理程序精品文档放心下载提供滤波器系数在MATLAB中利用滤波器设计分析工具(FDATOOL),根据指谢谢阅读定的滤波器性能快速设计一个FIR,然后把滤波器系数以头文件形式导入CCS中,谢谢阅读头文件中包含MATLAB辅助DSP实现FIR数字滤波器的阶数和系数数组,在感谢阅读MATLAB中调试、运行DSP程序并显示、分析处理后的数据。精品文档放心下载打开CCSIDE中file的data加载初始化dataAddress中输入数据存放感谢阅读的起始地址indata，并设置数据栈长度。4.1.3运行程序仿真运行程序，查看输入输出波形，修改相应参数进行调试。在View的Graph中单击Time/frequency出现graphpropertydialog框。谢谢阅读将显示类型，图形名称，起始地址，抽样点数，数据类型等分别进行设置，谢谢阅读输出各种波形输出滤波前的信号时域波形首先，设置如下图所示，输出滤波谢谢阅读前的信号时域波形。图3滤波前语音信号的时域GraphPropertyDialog设置参数图感谢阅读将DisplayType改为FFTMagnitude就可以看到滤波前信号的频谱，将谢谢阅读StartAddress改为outdata就可以看到滤波后信号的时域波形和频谱。4.2CCS仿真结果感谢阅读图4滤波前语音信号的时域波形图图5滤波前信号频谱波形图4.2.2FIR滤波器滤波后的信号感谢阅读图6滤波后信号时域波形图图7滤波后信号频谱图图17波前的信号，由高频信号和原信号叠加而成的波形较为杂乱，加噪后谢谢阅读的信号幅值较大。图19波后的信号，从图上可以看到信号幅度较低，波形较为谢谢阅读清晰。对比图18图20可以看到加躁后信号幅值较大，信号的能量集中在高频段，谢谢阅读滤波后，高频成分大部分被滤去，信号的能量主要集中在低频段。谢谢阅读5问题及解决方法1.为什么存入DAT文件的信号数据要乘上32768？感谢阅读由于TMS32054X是16位定点处理器，在进行汇编程序设计时，FIR滤波谢谢阅读器系数需要采用Q15格式，即必须将上述系数据转化位Q15定点式。.

2.matlab编译错误出现如下错误？谢谢阅读Thecarrierfrequencymustbelessthanhalfthesamplingfrequency.故降低载波频率感谢阅读3.在View的Graph中单击Time/frequency出现graphpropertydialog框，感谢阅读未修改抽样点数，显示的图形出现差异.4.在CCS编译时，发现fir.h文件老是出错，最后发现滤波系数要以整数体感谢阅读现，而且不能加分号。6设计感想通过本次课程设计，我们发现在课本上学到的知识实在是不能使我们有能力精品文档放心下载驾驭这个课程设计。一开始决定做这个课题，我们就去查找了大量资料，通过查找资料我们都清感谢阅读楚了大概要做些什么。首先要进行MATLAB仿真，并通过MATLAB设计的FIR感谢阅读滤波器输出滤波器系数接着便是通过CCSDSP上运行我们觉得感谢阅读思路是很清晰的，但是实际操作起来却困难重重。FIR感谢阅读滤波器的相关介绍，重点看了凯塞窗。然后查资料找出与本次设计有关的谢谢阅读MATLAB知识谢谢阅读转了几转。动工，这时才深刻体会到“纸上学来终觉浅，用到实际才叫真！”开始精品文档放心下载设计从噪声选取，到低通，高通，带阻，带通滤波器的设计，其中是错漏百出。精品文档放心下载当然，作为老手自然不会再去犯什么程序设计上的低级错误，主要还是对滤波器感谢阅读的认知不够深刻有的只是照抄书本，最后是连自己都为自己所犯的错误惊愕不感谢阅读已其中许多错误差点让自己从头来过，一度让自己有放弃的打算最后还是坚感谢阅读持了下来，自己大叹不容易。。精品文档放心下载拥有扎实的理论知识，认真细致的处事态度，随机应变的实际应用方法，不拘一精品文档放心下载格的思路都将使得自己能有好的未来。参考文献[1]邹彦，DSP2012．

[2]王宏.MATLAB6.5及其在信号处理中的应用.[M].北京:清华大学出版社,2004

[3]高西全,丁玉美,阔永红.数字信号处理原理实现及应用.[M].北京:电子工业出谢谢阅读版社,2006[4]胡广书编著，数字信号处理理论、算法与实现.2005；感谢阅读附录１.Matlab源程序%%%%%%%%%%%原始音频抽样%%%%%%%%%%%%%Fs=10000;bits=16;精品文档放心下载T=1;ts=1/Fs;N=38;[x,Fs,bits]=wavread('');%利用Windows下的录音机，录制一段自己的语音信号，精品文档放心下载时间为1sx=x(:,1);m=length(x);%sound(x,Fs,bits);精品文档放心下载y=fft(x,m);%对录制的信号进行傅立叶变换f=(Fs/m)*[1:m];t=[1:m]/Fs;figure(1);精品文档放心下载subplot(211);plot(t,x);title('原始信号波形');xlabel('time(s)');subplot(212);plot(f,abs(y));title('原始信号频谱');xlabel('frequency(hz)');noise4=0.01*sin(2*pi*5000*t)+0.01*sin(2*pi*200*t);%带通噪声%%%%%%%%%%%%%带通滤波%%%%%%%%%%%%%fpl4=500;fph4=2000;fsl4=700;fsh4=1800;wp4=[2*fpl4/Fs2*fph4/Fs];ws4=[2*fsl4/Fs2*fsh4/Fs];X4=x+noise4';%加了带通噪声的信号

%sound(X4,Fs,bits);%回放噪声信号谢谢阅读Y41=fft(X4,m);%对加了带通噪声的信号进行傅立叶变换figure(2);感谢阅读subplot(211);plot(t,X4);title('加了带通噪声的信号时域图');xlabel('time(s)');

subplot(212);plot(f,abs(Y41));title('加了带通噪声的信号频谱图

');xlabel('frequency(hz)');感谢阅读figure(3);wc4=(wp4+ws4)/2;感谢阅读b4=fir1(N,wc4,kaiser(N+1));%调用kaiser计算带通FIRDF的b4（n）[hn4wp4]=freqz(b4,1,512);精品文档放心下载subplot(211);plot(wp4/pi,20*log10(abs(hn4)));title('带通滤波器幅频');xlabel('w');

subplot(212);plot(wp4/pi,unwrap(angle(hn4)));title('带通滤波器相频');xlabel('w');

figure(4);x4get=filter(b4,1,X4);感谢阅读Y4get=fft(x4get,m);%sound(x4get,Fs,bits);精品文档放心下载subplot(211);plot(t,x4get);title('加了带通噪声的信号滤波后时域图');xlabel('time(s)');精品文档放心下载subplot(212);plot(f,abs(Y4get));title('加了带通噪声的信号滤波后频谱图');xlabel('frequency(hz)');精品文档放心下载%%%%%%%%%