北邮dsp软件实验报告

Matlab仿真实验

实验报告

学院：电子工程学院

专业：电子信息科学与技术

班级：

学号：

姓名：

时间：2015年12月23日

实验一：数字信号的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分析其频谱，要求得到清楚的三根谱线。【实验代码】：

k=2000;

n=[1:1:k];

x=0.001*cos(0.45*n*pi)+sin(0.3*n*pi)-cos(0.302*n*pi-pi/4);

subplot(2,1,1);

stem(n,x,'.');

title(‘时域序列');

xlabel('n');

ylabel('x(n)');

xk=fft(x,k);

w=2*pi/k*[0:1:k-1];

subplot(2,1,2);

stem(w/pi,abs(xk));

axis([0 0.5 0 2]);

title('1000点DFT');

xlabel('数字频率');

ylabel('|xk(k)|');

【实验结果图】：

（2）DTMF信号频谱分析

用计算机声卡采集一段通信系统中电话双音多频（DTMF）拨号数字0～9的数据，采用快速傅里叶变换（FFT）分析这10个号码DTMF拨号时的频谱。

【实验代码】及【实验结果图】见实验二（2）；

实验二：DTMF信号的编码

1.实验目的：

（a）复习和巩固IIR数字滤波器的基本概念；

（b）掌握IIR数字滤波器的设计方法；

（c）掌握IIR数字滤波器的实现结构；

（d）能够由滤波器的实现结构分析滤波器的性能（字长效应）；

（e）了解通信系统电话DTMF拨号的基本原理和IIR数字滤波器的实现方法。2.实验内容、要求及结果：

1）把你的联系电话号码通过DTMF编码生成一个.wav文件。

✧技术指标：

➢根据ITU Q.23建议，DTMF信号的技术指标是：传送/接收率为每秒10个号码，或每个号码100ms。

➢每个号码传送过程中，信号存在时间至少45ms，且不多于55ms，100ms 的其余时间是静音。

➢在每个频率点上允许有不超过±1.5%的频率误差。任何超过给定频率±3.5%的信号，均被认为是无效的，拒绝接收。

（其中关键是不同频率的正弦波的产生。可以使用查表方式模拟产生两个不同频率的正弦波。正弦表的制定要保证合成信号的频率误差在±1.5%以内，同时是取样点数尽量少）

【实验代码】：

N=800;

total_x=[];

tm=[49,50,51,65;52,53,54,66;55,56,57,67;42,48,35,68];

n=1:N;

fs=8000;

lf=[697 770 852];

hf=[1209 1336 1477];

x1=sin(2*pi*lf(1)*n/fs)+sin(2*pi*hf(1)*n/fs);%1

x2=sin(2*pi*lf(1)*n/fs)+sin(2*pi*hf(3)*n/fs);%3

x3=sin(2*pi*lf(1)*n/fs)+sin(2*pi*hf(2)*n/fs);%2

x4=sin(2*pi*lf(2)*n/fs)+sin(2*pi*hf(3)*n/fs);%6

x5=sin(2*pi*lf(3)*n/fs)+sin(2*pi*hf(3)*n/fs);%9

x6=sin(2*pi*lf(2)*n/fs)+sin(2*pi*hf(3)*n/fs);%6

x7=sin(2*pi*lf(1)*n/fs)+sin(2*pi*hf(2)*n/fs);%2

x8=sin(2*pi*lf(3)*n/fs)+sin(2*pi*hf(2)*n/fs);%8

x9=sin(2*pi*lf(3)*n/fs)+sin(2*pi*hf(3)*n/fs);%9

xa=sin(2*pi*lf(3)*n/fs)+sin(2*pi*hf(3)*n/fs);%9

xb=sin(2*pi*lf(2)*n/fs)+sin(2*pi*hf(2)*n/fs);%5

total_x=[x1,x2,x3,x4,x5,x6,x7,x8,x9,xa,xb];

x11=[x1,zeros(1,800)];

x22=[x2,zeros(1,800)];

x33=[x3,zeros(1,800)];

x44=[x4,zeros(1,800)];

x55=[x5,zeros(1,800)];

x66=[x6,zeros(1,800)];

x77=[x7,zeros(1,800)];

x88=[x8,zeros(1,800)];

x99=[x9,zeros(1,800)];

xaa=[xa,zeros(1,800)];

xbb=[xb,zeros(1,800)];

x=[x11,x22,x33,x44,x55,x66,x77,x88,x99,xaa,xbb];

x=x/max(abs(x));