基于matlab的数字音效处理器——数字信号处理课设报告.doc

实验题目 专 业： 通信工程 班 级： 114班 姓 名： XXXXX 学 号： XXXXXXXX 指导教师： XXXXX 一、实验类型 综合研究性实验 二、实验目的 1．掌握采样定理及FFT谱分析的基本原理及其利用Matlab的实现方法； 2. 掌握数字滤波器的设计原理和方法； 3. 学习用MATLAB编程实现语音数字滤波系统。 三、实验要求 1. 利用Windows下的录音机设备采集语音信号； 2. 对语音信号进行采样并混进加性噪声，作频谱分析； 3. 通过频谱分析选择合适的滤波器性能指标，设计合适的数字滤波器，并对含噪音的语音信号进行数字滤波； 4. 设计处理系统的用户界面（GUI），在所设计的系统界面上可以选择滤波器的参数，显示滤波器的频率响应，选择信号等。 四、数字滤波器的设计原理 数字滤波器可以理解为是一个计算程序或算法，将代表输入信号的数字时间序列转化为代表输出信号的数字时间序列，并在转化过程中，使信号按预定的形式变化。数字滤波器有多种分类，根据数字滤波器冲激响应的时域特征，可将数字滤波器分为两种，即无限长冲激响应(IIR)滤波器和有限长冲激响应(FIR)滤波器。IIR数字滤波器具有无限宽的冲激响应，与模拟滤波器相匹配，所以IIR滤波器的设计可以采取在模拟滤波器设计的基础上进一步变换的方法。其设计方法主要有经典设计法、直接设计法和最大平滑滤波器设计法。FIR数字滤波器的单位脉冲响应是有限长序列。它的设计问题实质上是确定能满足所要求的转移序列或脉冲响应的常数问题设计方法主要有窗函数法、频率采样法和等波纹最佳逼近法等。在对滤波器实际设计时，整个过程的运算量是很大的。设计阶数较高的IIR滤波器时，计算量更大，设计过程中改变参数或滤波器类型时都要重新计算。设计完成后对已设计的滤波器的频率响应要进行校核。要得到幅频、相频响应特性，运算量也是很大的。平时所要设计的数字滤波器，阶数和类型并不一定是完全给定的，很多时候要根据设计要求和滤波效果不断地调整，以达到设计的最优化。在这种情况下，滤波器设计就要进行大量复杂的运算，单纯的靠公式计算和编制简单的程序很难在短时间内完成。利用MATLAB强大的计算功能进行计算机辅助设计，可以快速有效地设计数字滤波器，大大地简化了计算量。 fp2=4800;fs2=5000; wp2=2*fp2/Fs; ws2=2*fs2/Fs; rp=1; as=100; [N2,wp2]=ellipord(wp2,ws2,rp,as); [B2,A2]=ellip(N2,rp,as,wp2,high); y2=filter(B2,A2,x); figure; freqz(B2,A2); figure; y2t=y_2(t); subplot(2,1,1); t=n*T; plot(t,y2); xlabel(t/s);ylabel(y2t); axis([0,t(end),min(y2),1.2*max(y2)]) subplot(2,1,2); plot(f,abs(fft(y2))); 5、 带通滤波器的设计 fpl=1200; fpu=3000; fsl=1000; fsu=3200; wp3=[2*fpl/Fs,2*fpu/Fs]; ws3=[2*fsl/Fs,2*fsu/Fs]; rp=1; as=100; [N3,wp3]=ellipord(wp3,ws3,rp,as); [B3,A3]=ellip(N3,rp,as,wp3); y3=filter(B3,A3,x); figure; freqz(B3,A3); y3t=y_3(t); figure; subplot(2,1,1); t=n*T; plot(t,y3); xlabel(t/s);ylabel(y3t); axis([0,t(end),min(y3),1.2*max(y3)]) subplot(2,1,2); plot(f,abs(fft(y3))); 6、 音频放大 figure o=10 for xa=1:16384; p(xa,1)=x(xa,1)*xa*o; end subplot(2,1,1); plot(x); xlabel(原始波形) subplot(2,1,2); plot(n