资源描述：
《DSP软件实验报告_让你最牛》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、北京邮电大学DSP软件课程实验报告姓名学号：班级：院系：电子工程学院实验一：数字信号的FFT分析1、实验内容及要求(1)离散信号的频谱分析：设信号此信号的0.3pi和0.302pi两根谱线相距很近，谱线0.45pi的幅度很小，请选择合适的序列长度N和窗函数，用DFT分析其频谱，要求得到清楚的三根谱线。(2)DTMF信号频谱分析用计算机声卡采用一段通信系统中电话双音多频（DTMF）拨号数字0～9的数据，采用快速傅立叶变换（FFT）分析这10个号码DTMF拨号时的频谱。2、实验目的通过本次实验，应该掌握：(a)用傅立叶变换进行信号分析时基本参数的选择。
2、(b)经过离散时间傅立叶变换（DTFT）和有限长度离散傅立叶变换（DFT）后信号频谱上的区别，前者DTFT时间域是离散信号，频率域还是连续的，而DFT在两个域中都是离散的。(c)离散傅立叶变换的基本原理、特性，以及经典的快速算法（基2时间抽选法），体会快速算法的效率。(d)获得一个高密度频谱和高分辨率频谱的概念和方法，建立频率分辨率和时间分辨率的概念，为将来进一步进行时频分析（例如小波）的学习和研究打下基础。(e)建立DFT从整体上可看成是由窄带相邻滤波器组成的滤波器组的概念，此概念的一个典型应用是数字音频压缩中的分析滤波器，例如DVDAC3和MP
3、EGAudio。3、实验步骤及结果(1)离散信号频谱分析程序如下：closeall;%关闭文件clear;%清除内存N=1000;n=0:N-1;%n=0,1,.....999xn=0.001*cos(0.45*n*pi)+sin(0.3*n*pi)-cos(0.302*n*pi-pi/4);%生成含有1000个元素的信号序列yn=fft(xn,N);%快速傅里叶变换stem(2*n/N,log(abs(yn)+1),'b');%abs(yn)表示求模，作脉冲图axis([0.25,0.5,-1,10])%限制坐标范围title('FFT频谱分析'
4、);%标题ylabel('对X(k)的模取对数');%Y轴标签text(0.35,8,'电子4班，许金良，10210929')text(0.3,2,'注意：幅度用log函数进行了对数处理')运行结果：图一（左边两个频率分量相距很近，右边一个频率分量很小）(2)DTMF信号频谱分析程序如下：closeall;clear;tm=[1,2,3,65;4,5,6,66;7,8,9,67;42,0,35,68];%DTMF信号代表的16个数N=205;f1=[697,770,852,941];%行频率向量f2=[1209,1336,1477,1633];TN
5、=input('键入1位0-9的数字=');forp=1:4forq=1:4iftm(p,q)==abs(TN),break,end%检测码相符的列号qendiftm(p,q)==abs(TN),break,endendn=0:1023;x=sin(2*pi*n*f1(p)/8000)+sin(2*pi*n*f2(q)/8000);%构成双频信号X=fft(x',1024);X=abs(X);y=sort(X);fori=1:1024ifX(i)<=y(1020)X(i)=0;end;endx=[1:512]./512*4000;X=X(1:512
6、);fo=[];fori=1:512ifX(i)~=0fo=[fo,i];endendfo=round(fo*8000/1024);stem(x,X)wucha=(fo-[f1(p),f2(q)])./[f1(p),f2(q)]*100;wucha=strcat(['频率误差：','[',num2str(wucha),']','%']);text(2000,350,['所拨打的号码：',int2str(TN)])text(2000,300,['检测到的频率：',int2str(fo)]);text(2000,250,['查表所得频率：',int2s
7、tr([f1(p),f2(q)])]);text(2000,200,wucha,'color','r');title('DTMF频谱分析与误差','color','r')disp(['检测到的频率为：',num2str(fo)])输入0:图二输入1:图三输入2:图四输入3:图五输入4:图六输入5：图七输入:6：图八输入7：图九输入8：图十输入9：图十一实验二：DTMF信号的编码1、实验内容及要求（1）把您的联系电话号码通过DTMF编码生成为一个.wav文件。技术指标：根据ITUQ.23建议，DTMF信号的技术指标是：传送/接收率为每秒10个号码，或
8、每个号码100ms。每个号码传送过程中，信号存在时间至少45ms，且不多于55ms，100ms的其余时间是静音。在每个频率
