资源描述：
《dsp软件实验报告--matlab北邮》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、Matlab仿真实验实验报告学院：电子工程学院专业：电子信息科学与技术班级：2012211201学号：2012210815J1•Aa名：李昊轩时间：2024年22月19日实验一：数字信号的FFT分析1.实验目的通过本次试验，应该掌握：(a)用傅里叶变换进行信号分析时基本参数的选择(b)经过离散时间傅里叶变换和有限长度离散傅里叶变换后信号频谱上的区别，前者DTFT时间域是离散信号，频率域还是连续的，而DFT在两个域中都是离散的。(0)离散傅里叶变化的基本原理、特性，以及经典的快速算法(基2吋间抽选法)，体会快速算法的效率。(d)获
2、得一个高密度频谱和高分辨率频谱的概念和方法，建立频率分辨率和时间分辨率的概念，为将来进一步进行时频分析(例如小波)的学习和研究打下基础。(e)建立DFT从整体上可看成是由窄带相邻滤波器组成的滤波器组的概念，此概念的一个典型应用时数字咅频压缩中的分析滤波器，例如DVDAC3和MPEGAudio。2.实验内容、要求及结果。(1)离散信号的频谱分析：设信号x(n)=0.001*cos(0.45nn)+sin(0.3n兀)-cos(0.302n兀-寸)此信号的0.3k和0.302Ji两根谱线相距很近,谱线0.45n的幅度很小，请选择合适
3、的序列长度N和窗函数，用DFT分析其频谱，要求得到清楚的三根谱线。【实验代码】：clear;closeall;N=5000;n=l:1:N;x=0.001*cos(0.45*n*pi)+sin(0.3*n*pi)-cos(0.302*n*pi-pi/4);y=fft(x,N);a=abs(y(1:1:N/2+1));k=0:l:N/2;w=2*pi/N*k;stem(w/pi,a);axis([0.29,0.46,0,10]);【实验结果图】：'余工作区'余工作区tststrtr;closer；*cos(0:,N);r/2；N*
4、k;i.29,0.46,a)close(2)DTMF信号频谱分析用计算机声卡采集一段通信系统中电话双音多频(DTMF)拨号数字0〜9的数据,采用快速傅里叶变换(FFT)分析这10个号码DTMF拨号时的频谱。【实验代码】：num=[1,2,3,65;4Z5,6,66;7,8,9,67;42,0z35,68];fl=[697,770,852,941];f2=[1209z1336,1477,1633];in=input(*键入丄位0-9的数字=');forp=l:4forq=l:4ifnum(pzq)==inzbreak,endend
5、ifnum(p,q)==infbreaksendendn=0:1023;x=sin(2*pi*n*fl(p)/8000)+sin(2*pi*n*f2(q)/8000);X=fft(x1,1024);X=abs(X);y=sort(X);fori=l:1024ifX(i)<=y(1020)X(i)=0;end;endx=[l:512]./512*4000;X=X(1:512);f0=[];fori=l:512ifX(i)~=0f0=[f0zi];endendf0=round(fO*8000/1024);stem(x,X);text
6、(2000,350,『所拨打的号码1,int2str(in)]);text(2000,300,『检测到的频率1,int2str(fO)]);text(2000,250,『查表所得频率1/int2str([f1(p),f2(q)])]);title(fDTMF信号频谱分析1);曹二：囂Simulmk.S代讶S1MUUNK>►TOM►Documents►MATLABable*52.941•.1477I■ableuble站令行IB口n=u:】UQ：x=sin(2*pi*n*f1(p)/8000)+sm(2*pi*n*f2(q)X=ff
7、t(xJ,1024);X二abs(X);尸sort(X);fori=l:1024ifX(i)<=y(1020)X(i)=0;end;endx=[l:512]./512*4000.X=X(1:512);f0=()：fori=l:512ifX(i)-=0fO=(fO,i];endendf0=round(f0«8000/1024):stea(x,X):text(2000,350,「所拨打的号玛:int2str(in)]);text(2000,300,['检测到的频苹：r,mt2str(f0)]):tert(2000,250,V直夷所
8、谓频软int2str([f1(p),£2(q)])]):titleCDIMF信号频谱分析')；瀝入1位0-9的敎豪I»MFigure11►Documents►MATLAB命令行帝口n=u:iuzo:x=sin(2*pi*n*f1(p)/8000)+sin(2*p
