数字信号处理作业

成绩。作业与上机实验。

第三四章）班级： 13-电信

学号：姓名：

任课老师：李宏民

完成时间： 15.11.08

信息与通信工程学院。

2015—2016学年第 1学期。

第3-4章离散傅里叶变换及快速算法与应用。

一、实验内容及要求：

1、对信号在x(n)=，n=0，1，2...7,求其n=8点的dft。要求采用基于时间抽取算法编写fft实现程序，画出dft幅度谱与相位谱，并将计算结果与用matlab自带的fft函数计算结果进行比较。

2、一个由40hz和100 hz正弦信号构成的信号，受零均值随机噪声的干扰（噪声服从标准正态分布，由randn函数产生，n(t)= randn(m,n)返回一个m*n的随机矩阵），即。数据采样率为500hz，试用fft函数来分析其信号频率成分。

求其幅度频谱，从频谱图（横坐标以hz为单位，用plot函数画图）中能否观察出信号的2个频率分量？

提高采样点数，再求该信号的幅度频谱图，此时幅度频谱发生了什么变化？信号的2个模拟频率和数字频率各为多少？fft频谱分析结果与理论上是否一致？

3、研究高密度频谱与高分辨率频谱。频率分辨率是指所用的算法能将信号中两个靠得很近的谱峰分开的能力。信号末尾补零由于没有对原信号增加任何新的信息，因此不能提高频率分辨率，但可以减小栅栏效应，所得到的频谱称为高密度频谱。

在维持采样频率不变的情况下，为提高分辨率只能增加采样点数n，此时所得到的频谱称为高分辨率频谱。设有连续信号。

以采样频率对信号采样，分析下列几种情况的幅频特性。

1 采集数据长度n=16点，做n=16的fft；采集数据长度n=16点，补零到256点，做256点的fft；

2 采集数据长度n=64点，做n=64的fft；采集数据长度n=64点，补零到256点，做256点的fft；

3 采集数据长度n=256点，做n=16点的fft。

4 观察以上几幅不同的幅频特性曲线，分析和比较它们的特点，并说明形成的原因。

5 注：在matlab中，可用zeros函数实现填零运算，例如x1=[x(1:1:16),zeros(1,240)]。所有频谱图的横坐标以hz为单位，用plot函数画图。

6 一个连续信号含三个频率分量，经采样得以下序列。已知n =16，分别为1/16，1/64，观察其频谱；当n =64，128，不变，其结果有何不同，为什么？

4、语音信号频谱分析。

首先应用录音软件(见下图)，录音“我是xxx（学生姓名）”，得到语音数据“录音。w**”并将“录音。w**”文件复制到 matlab work 目录下。

其次在matlab环境下，x(n)=w**read('录音。w** ',2048)，得到2048点的语音数据x(n),对x(n)进行fft频谱分析，画出频谱图（横坐标以hz为单位，用plot函数画图）。

5、双音多频（dtmf）信号是将拨号盘上的0-f共16个数字，用音频范围的8个频率来表示的一种编码方式。8个频率分为高频群和低频群两组，分别作为列频和行频。每个字符的信号由来自列频和行频的两个频率的正弦信号叠加而成。

频率组合方式如表8-2所示。

双音多频（dtmf）信号频率组合方式。

找一部**机，在免提状态下拨号0～9中的任意7位号码，并录制双音多频拨号声音。采用快速傅里叶变换（fft）分析这7个号码拨号时的频谱，并与理论值进行比较。

三、波形图。

1．dft幅度谱与相位谱如图1可见。

xn=[1,2,3,4,5,6,7,8];

m=nextpow2(length(xn)),n=2^m,a=[xn,zeros(1,n-length(xn))]

g=0;for i=0:n-1;

if i t=a(i+1);a(i+1)=a(g+1);a(g+1)=t;

endk=n/2;

while g>=k;

g=g-k;k=k/2;

endg=g+k;

endwn=exp(-j*2*pi/n);

for l=1:m;

disp('级次运算'),disp(l);

b=2^(l-1);

for r=0:b-1;

p=2^(m-l)*r;

for k=r:2^l:n-2;

t=a(k+1)+a(k+b+1)*wn^p;

a(k+b+1)=a(k+1)-a(k+b+1)*wn^p;

a(k+1)=t;

endend

endxk=a,fftxn=fft(xn,n),k=0:7;wk=2*k/8;

subplot(2,2,1);stem(wk,abs(xk),'k.')

xlabel('hz');ylabel('幅度');

subplot(2,2,2);stem(wk,angle(xk),'k.')

xlabel('hz');ylabel('相位');

subplot(2,2,3);stem(wk,abs(fftxn),'k.')

xlabel('hz');ylabel('幅度');

subplot(2,2,4);stem(wk,angle(fftxn),'k.')

xlabel('hz');ylabel('相位');

答：由此可以看出信号经过fft函数和dft函数所得到的结果基本没有差异，说明fft函数能在复杂的运算中代替dft函数简化运算。

2.①幅度频谱如图2可见。

t=0:0.001:1;

s=sin(2*pi*100*t)+sin(2*pi*40*t)+randn(1,length(t));

figure(1)

plot(t,s);

y=fft(s,501);

f=500*(0:500)/501;

magy=abs(y);

figure(2);

plot(f,magy);

xlabel('hz');

ylabel('·幅度');

提高采样点数后的幅度频谱图如图3可见。

t=0:0.001:1;

s=sin(2*pi*100*t)+sin(2*pi*40*t)+randn(1,length(t));

figure(1)

plot(t,s);

y=fft(s,1001);

f=500*(0:1000)/1001;

magy=abs(y);

figure(2);

plot(f,magy);

xlabel('hz');

ylabel('幅度');

图2 频谱图。

答：提高采样点数之后，图像分辨率提高很多，更容易看出频率分量，模拟频率大概为0~3hz，数字频率大概为0~3hz,fft结果与理论上还是稍微有点区别。

3. ①16的fft和补零到256点的fft如图4可见。

magxa=abs(xa);

figure(1);

subplot(2,1,1);

plot(f,magxa);

xlabel('hz');

ylabel('幅度');

x=[xa(1:16),zeros(1,240)];

x=fft(x,256);

f=fs*(0:255)/256;

magx=abs(x);

subplot(2,1,2);

plot(f,magx);

xlabel('hz');

ylabel('幅度')

2) 64的fft和补零到256点的fft见图4

n=64;n=0:n-1;

fs=32000;

xa=cos(2*pi*6.5*10^3*n/fs)+cos(2*pi*7*10^3*n/fs)+cos(2*pi*9*10^3*n/fs);

xa=fft(xa,n);

f=fs*(0:63)/64;

magx=abs(xa);

figure(1);

subplot(2,1,1);

plot(f,magx);

xlabel('hz');

ylabel('幅度');

x=[xa(1:16),zeros(1,192)];

x=fft(x,256);

f=fs*(0:255)/256;

magx=abs(x);

subplot(2,1,2);

plot(f,magx);

xlabel('hz');

ylabel('幅度')

3) 16点的fft如图5可见。

n=256;

n=0:n-1;

fs=32000;

xa=cos(2*pi*6.5*10^3*n/fs)+cos(2*pi*7*10^3*n/fs)+cos(2*pi*9*10^3*n/fs);

xa=fft(xa,n);

f=fs*(0:255)/256;

magx=abs(xa);

figure(1);

plot(f,magx);

xlabel('hz');

ylabel('幅度');

采样数据为16的时候图像由于点数太少，很难反映出频谱的细节特征，只能分辨出两个频率分量，随着采样点的不段增加频率间隔缩小，连续曲线越来越平滑，并且由16点经过补零所得到的，频率分辨率并未提高，仍然只能看出两个频率分量。

6) n=16，为1/16时的源程序，如图6可见：

n=16;n=0:n-1;

df=1/16;

x=sin(2*pi*0.15*n)+cos(2*pi*(0.15+df)*n)+cos(2*pi*(0.15+2*df)*n);

数字信号处理作业

数字信号处理作业

数字信号处理作业

数字信号处理作业

其他用户还读了