发布时间 : 星期二 文章DSP软件实验分析报告--matlab-北邮更新完毕开始阅读34cf225c7a3e0912a21614791711cc7930b7784f
DSP软件实验报告--matlab-北邮
———————————————————————————————— 作者： ———————————————————————————————— 日期：
2
实验一：数字信号的FFT分析 1.实验目的
通过本次试验，应该掌握：
（a） 用傅里叶变换进行信号分析时基本参数的选择
（b） 经过离散时间傅里叶变换和有限长度离散傅里叶变换后信号频谱上的区别，前者DTFT
时间域是离散信号，频率域还是连续的，而DFT在两个域中都是离散的。 （c） 离散傅里叶变化的基本原理、特性，以及经典的快速算法（基2时间抽选法），体会快
速算法的效率。
（d） 获得一个高密度频谱和高分辨率频谱的概念和方法，建立频率分辨率和时间分辨率的
概念，为将来进一步进行时频分析（例如小波）的学习和研究打下基础。
（e） 建立DFT从整体上可看成是由窄带相邻滤波器组成的滤波器组的概念，此概念的一个
典型应用时数字音频压缩中的分析滤波器，例如DVD AC3和MPEG Audio。
2.实验内容、要求及结果。 （1）离散信号的频谱分析：
设信号x(n)=0.001*cos(0.45nπ)+sin(0.3nπ)-cos(0.302nπ-) 4π
此信号的0.3π和0.302π两根谱线相距很近，谱线0.45π的幅度很小，请选
择合适的序列长度N和窗函数，用DFT分析其频谱，要求得到清楚的三根谱线。
【实验代码】：
clear;close all; N=5000; n=1: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:1:N/2; w=2*pi/N*k; stem(w/pi,a);
axis([0.29,0.46,0,10]);
【实验结果图】：
3
（2）DTMF信号频谱分析
用计算机声卡采集一段通信系统中电话双音多频（DTMF）拨号数字0～9的数据，采用快速傅里叶变换（FFT）分析这10个号码DTMF拨号时的频谱。
【实验代码】：
clear;
close all;
column=[1209,1336,1477,1633]; line=[697,770,852,941] fs=10000; N=1024; ts=1/fs; n=0:N-1;
f=0:fs/N:fs/N*(N-1);
key=zeros(16,N);
key(1,:)=cos(2*pi*column(1)*n*ts)+cos(2*pi*line(1)*n*ts); key(2,:)=cos(2*pi*column(2)*n*ts)+cos(2*pi*line(1)*n*ts); key(3,:)=cos(2*pi*column(3)*n*ts)+cos(2*pi*line(1)*n*ts);
4