DSP软件实验报告.docVIP

PAGE 2 DSP软件实验报告 班级： 姓名： 学号： 班内序号： 实验一：数字信号的 FFT 分析 实验内容及要求 离散信号的频谱分析 设信号 此信号的0.3pi和0.302pi两根谱线相距很近，谱线0.45pi的幅度很小，请选择合适的序列长度 N 和窗函数，用DFT分析其频谱，要求得到清楚的三根谱线。 DTMF 信号频谱分析 用计算机声卡采集一段通信系统中电话双音多频（DTMF）拨号数字0～9的数据，采用快速傅立叶变换（FFT）分析这10个号码DTMF拨号时的频谱。 实验目的 通过本次实验，应该掌握： 用傅立叶变换进行信号分析时基本参数的选择。 经过离散时间傅立叶变换（DTFT）和有限长度离散傅立叶变换（DFT） 后信号频谱上的区别，前者 DTFT 时间域是离散信号，频率域还是连续的，而 DFT 在两个域中都是离散的。 离散傅立叶变换的基本原理、特性，以及经典的快速算法（基2时间抽选法），体会快速算法的效率。 获得一个高密度频谱和高分辨率频谱的概念和方法，建立频率分辨率和时间分辨率的概念，为将来进一步进行时频分析（例如小波）的学习和研究打下基础。 建立 DFT 从整体上可看成是由窄带相邻滤波器组成的滤波器组的概念，此概念的一个典型应用是数字音频压缩中的分析滤波器，例如 DVD AC3 和MPEG Audio。 程序源代码 离散信号的频谱分析 % w/(2*pi) = k/N % 0.3pi/(2*pi) = 150/1000 0.302pi/(2*pi) = 151/1000 0.45pi/(2*pi) = 225/1000 N=1000; n=[0:1:N-1]; x=0.001*cos(0.45*n*pi)+sin(0.3*n*pi)-cos(0.302*n*pi-pi/4); y=fft(x,N); stem(n,abs(y),.); title(信号频谱); xlabel(k); ylabel(|X(k)|) DTMF 信号频谱分析 clear;close all; %% DTMF编码 Fs=8000; %采样频率8kHz N=800; %每个号码100ms,800个点 n=0:N/2-1; %信号时间50ms z=zeros(1,N/2); %静音时间50ms y=zeros(10,800); %编码信号初始化 freq=[941,1336; 697,1209; 697,1336; 697,1477; 770,1209; 770,1336; 770,1477; 852,1209; 852,1336; 852,1477]; %0~9对应频率 for i=0:9 %10次循环编码 f=freq(i+1,:); %按电话号码读取对应频率 x=sin(2*pi*f(1)/Fs*n)+sin(2*pi*f(2)/Fs*n); %产生信号 y(i+1,:)=[x,z]; %输出信号后接静音 end; %% DTMF解码 xn=y; %把y重构成800行10列的矩阵，每一列为100ms信号 xk=abs(fft(xn,Fs)); %DFT变换 for i=0:9 %打印号码0~9DTMF编码后信号频谱 subplot(2,5,i+1); plot(xk(:,i+1)); title([号码,num2str(i),频谱]); axis([0,Fs,0,200]); end 程序运行结果 离散信号的频谱分析 0.3pi和0.302pi两根谱线 0.45pi谱线 DTMF 信号频谱分析 实验二： DTMF 信号的编码 实验内容及要求 把您的联系电话号码通过DTMF编码生成为一个.wav文件。 技术指标：（时域） 根据 ITU Q.23 建议，DTMF 信号的技术指标是：传送/接收率为每秒 10 个号码，每个号码 100ms。 每个号码传送过程中，信号存在时间至少 45ms，且不多于 55ms，100ms 的其余时间是静音。 在每个频率点上允许有不超过 ±1.5% 的频率误差。任何超过给定频率 ±3.5% 的信号，均被认为是无效的，拒绝接收。 （其中关键是不同频率的正弦波的产生。可以使用查表方式模拟产生两个不同频率的正弦波。正弦表的制定要保证合成信号的频率误差在±1.5%以内，同时使取样点数尽量少） 对所生成的DTMF文件进行解码。 DTMF 信号解码可以采用 FFT 计算 N 点频率处的频谱值，然后估计出所拨号码。但 FFT计算了许多不需要的值，计算量太大，而且为保证频率分辨率，FFT的点数较大，不利于实