DSP软件实验报告--matlab北邮

实验一数字信号的FFT分析1实验目的通过本次试验，应该掌握（A）用傅里叶变换进行信号分析时基本参数的选择（B）经过离散时间傅里叶变换和有限长度离散傅里叶变换后信号频谱上的区别，前者DTFT时间域是离散信号，频率域还是连续的，而DFT在两个域中都是离散的。（C）离散傅里叶变化的基本原理、特性，以及经典的快速算法（基2时间抽选法），体会快速算法的效率。（D）获得一个高密度频谱和高分辨率频谱的概念和方法，建立频率分辨率和时间分辨率的概念，为将来进一步进行时频分析（例如小波）的学习和研究打下基础。（E）建立DFT从整体上可看成是由窄带相邻滤波器组成的滤波器组的概念，此概念的一个典型应用时数字音频压缩中的分析滤波器，例如DVDAC3和MPEGAUDIO。2实验内容、要求及结果。（1）离散信号的频谱分析设信号XN0001COS045NSIN03NCOS0302N4此信号的03谱线相距很近，谱线045的幅度很小，请选和0302两根择合适的序列长度N和窗函数，用DFT分析其频谱，要求得到清楚的三根谱线。【实验代码】CLEARCLOSEALLN5000N11NX0001COS045NPISIN03NPICOS0302NPIPI/4YFFTX,NAABSY11N/21K01N/2W2PI/NKSTEMW/PI,AAXIS029,046,0,10【实验结果图】（2）DTMF信号频谱分析用计算机声卡采集一段通信系统中电话双音多频（DTMF）拨号数字09的数据，采用快速傅里叶变换（FFT）分析这10个号码DTMF拨号时的频谱。【实验代码】CLEARCLOSEALLCOLUMN1209,1336,1477,1633LINE697,770,852,941FS10000N1024TS1/FSN0N1F0FS/NFS/NN1KEYZEROS16,NKEY1,COS2PICOLUMN1NTSCOS2PILINE1NTSKEY2,COS2PICOLUMN2NTSCOS2PILINE1NTSKEY3,COS2PICOLUMN3NTSCOS2PILINE1NTSKEY4,COS2PICOLUMN1NTSCOS2PILINE2NTSKEY5,COS2PICOLUMN2NTSCOS2PILINE2NTSKEY6,COS2PICOLUMN3NTSCOS2PILINE2NTSKEY7,COS2PICOLUMN1NTSCOS2PILINE3NTSKEY8,COS2PICOLUMN2NTSCOS2PILINE3NTSKEY9,COS2PICOLUMN3NTSCOS2PILINE3NTSKEY10,COS2PICOLUMN2NTSCOS2PILINE4NTSFIGUREFORI110SUBPLOT4,4,IPLOTF,ABSFFTKEYI,GRIDEND【实验结果图】实验二DTMF信号的编码1实验目的（A）复习和巩固IIR数字滤波器的基本概念；（B）掌握IIR数字滤波器的设计方法；（C）掌握IIR数字滤波器的实现结构；（D）能够由滤波器的实现结构分析滤波器的性能（字长效应）；（E）了解通信系统电话DTMF拨号的基本原理和IIR数字滤波器的实现方法。2实验内容、要求及结果1）把你的联系电话号码通过DTMF编码生成一个WAV文件。技术指标根据ITUQ23建议，DTMF信号的技术指标是传送/接收率为每秒10个号码，或每个号码100MS。每个号码传送过程中，信号存在时间至少45MS，且不多于55MS，100MS的其余时间是静音。在每个频率点上允许有不超过15的频率误差。任何超过给定频率35的信号，均被认为是无效的，拒绝接收。（其中关键是不同频率的正弦波的产生。可以使用查表方式模拟产生两个不同频率的正弦波。正弦表的制定要保证合成信号的频率误差在15以内，同时是取样点数尽量少）【实验代码】DINPUT请键入电话号码,SSUMLENGTHDTOTAL_XSUM_XSUM_XSUM_X,ZEROS1,800FORA1SUMSYMBOLABSDATM49,50,51,6552,53,54,6655,56,57,6742,48,35,68FORP14FORQ14IFTMP,QABSDABREAK,ENDENDIFTMP,QABSDABREAK,ENDENDF1697,770,852,941F21209,1336,1477,1633N1400XSIN2PINF1P/8000SIN2PINF2Q/8000XX,ZEROS1,400SUM_XSUM_XXTOTAL_XTOTAL_X,XENDWAVWRITETOTAL_X,SOUNDWAVESOUNDTOTAL_XT18800/8000SUBPLOT2,1,1PLOTT,TOTAL_XAXIS0,12,2,2XLABEL时间/STITLEDTMF信号时域波形XKFFTXMXKABSXKSUBPLOT2,1,2K1800SUM8000/800PLOTK,MXKXLABEL频率TITLEDTMF信号频谱DISP双频信号已生成并发出【实验结果图】2）对所生成的DTMF文件进行解码。DTMF信号解码可以采用FFT计算N点频率处的频谱值，然后估计出所拨号码。但FFT计算了许多不需要的值，计算量太大，而且为保证频率分辨率，FFT的点数较大，不利于实时实现。因此，FFT不适合于DTMF信号解码的应用。由于只需要知道8个特定点的频谱值，因此采用一种称为GOERTZEL算法的IIR滤波器可以有效的提高计算效率。其传递函数为HKZ12112212【实验代码】K1820222431343842N205DISP接收端检测到的号码为FORA1SUMM800A1XGOERTZELTOTAL_XM1MN,K1VALABSXXK2VAL2XK2LIMIT80FORS58IFVALSLIMIT,BREAK,ENDFORR14IFVALRLIMIT,BREAK,ENDENDDISPSETSTRTMR,S4END【实验结果图】实验三FIR数字滤波器的设计和实现1实验目的通过本次试验，掌握一下知识FIR数字滤波器窗口设计法的原理和设计步骤；GIBBS效应发生的原因和影响；不同类型的窗函数对滤波效果的影响，以及窗函数和长度N的选择。（效果，耳机听前后声音，或者看前后的频谱图）2实验内容、要求及结果录制一段自己的声音，长度为十几秒，取样频率32KHZ，然后叠加一个高斯白噪声，（知道噪声分布，知道噪声功率，只要知道输入信号功率），使得信噪比为20DB。请采用窗口法，设计一个FIR带通滤波器，滤除噪声提高质量。提示滤波器指标参考通带边缘频率为4KHZ，阻带边缘频率为45KHZ，阻带衰减大于50DB；MATLAB函数YAWGNX,SNR,MEASURED,首先测量输入信号X的功率，然后对其叠加高斯白噪声。【实验代码】CLC,CLEARFS32000X,FS,BITSWAVREADWOWAVSNR20X2AWGNX,SNR,MEASURED,DBWAVWRITEX2,FS,16,WO1WAVT01/FSSIZEX21/FSWP8000PI/32000WS9000PI/32000WDELTAWSWPNCEIL11PI/WDELTA取整WNWSWP/2BFIR1N,WN/PI,BLACKMANN1选择窗函数，并归一化截止频率FIGURE1FREQZB,1,512F2FILTERB,1,X2TITLE滤波器幅频、相频特性FIGURE2SUBPLOT2,1,1PLOTT,X2TITLE滤波前时域波形SUBPLOT2,1,2PLOTT,F2TITLE滤波后时域波形F0FFTF2,1024FFS0511/1024FIGURE3Y2FFTX2,1024SUBPLOT2,1,1PLOTF,ABSY21512TITLE滤波前频谱XLABELHZYLABEL幅度SUBPLOT2,1,2F2PLOTF,ABSF01512TITLE滤波后频谱XLABELHZYLABEL幅度WAVWRITEF2,FS,16,WO2WAV【实验结果图】实验总结这次的MATLAB实验让学到了许多，以前是学过MATLAB的，但是已经忘的差不多了，正好有这次机会让我重新复习相关的知识。把MATLAB与数字信号处理的相关知识很好的联系到一起。在数字信号处理的理论课上，对于傅里叶变换、数字滤波器等知识的了解并不是很深刻，但是通过这次试验，我更好的理解了用傅里叶变换进行信号分析时基本参数的选择、离散时间傅里叶变换和有限长度离散傅里叶变换后信号频谱上的区别、离散傅里叶变化的基本原理、特性，以及经典的快速算法的效率，也更好的掌握了IIR数字滤波器的设计方法、实现结构及其性能，还有不同类型的窗函数对滤波效果的影响，以及窗函数和