MATLAB课程设计报告--信号的频域分析及MATLAB实现.docx

MATLAB电子信息应用课程设计

设计五

信号的频域分析及MATLAB实现

学院：信息学院专业班级：电信1201

学号：201

姓名：

一、设计目的

通过该设计，理解傅里叶变换的定义及含义，掌握对信号进行频域分析的方法。

二、设计内容相关的信号知识：

信号的离散傅里叶变换从连续到离散：

连续时间信号x(t)以及对应的连续傅里叶变换x(w)都是连续函数。将x和×都离散化，并且建立对应的傅里叶变换。

假设x(t)时限于[0,L],再通过时域采样将x(t)离散化，就可以得到有限长离散信号，记为xdscrete(t)。设采样周期为T,则时域采样点数N=L/T。

它的傅里叶变换为

这就是x(t)在时域采样后的连续傅里叶变换，也就是离散时间傅里叶变换，它在频域依然是连续的。

下面将频域信号转化为有限长离散信号。与对时域信号的处理类似，假设频域信号是带限的，再经过离散化，即可得到有限长离散信号。依据采样定理，时域采样若要能完全重建原信号，频域信号x(の)应当带限

于(0,1/T)。由于时域信号时限于[0,L],由采样定理以及时频对偶的关系，频域的采样间隔应为1/L。故，频域采样点数为：

即频域采样的点数和时域采样同为N,频域采样点为

{wk=2πk/NT}o≤kN

令T=1,将其归一化，就得到离散傅里叶变换。因此，DFT就是先将信

号在时域离散化，求其连续傅里叶变换后，再在频域离散化的结果。离散傅里叶变换：

1、有限长序列的离散傅里叶变换公式

MATLAB函数：fft功能是实现快速傅里叶变换，fft函数的格式为：y=fft(x),返回向量x的不连续fourier变化。ifft功能是实现快速反傅里叶变换，ifft函数的格式为：y=ifft(x),返回向量x的不连续inverse

fourier变化。

是一个N=12的有限序列，利用MATLAB计算x(n)它的傅里叶变换X(k)并画出图形，然后再对X(k)进行离散傅里叶反变

换，并求出x(n)画出其波形。

2、频率分辨率与DFT参数的选择

在DFT问题中，频率分辨率是指在频率轴上所能得到的最小频率间隔

即最小频率间隔反比于数据的长度N。若在x(n)中的两个频率分别为f和f?的信号，对x(n)用矩形窗截断，要分辨出这两个频率，N必须满足

通过下面实验，验证上面的结论：设一序列中含有两种频率成分，

f?=2Hz,f?=2.05Hz,采样频率取为f,=10Hz,表示为

x(n)=sin(2nf?n/f)=sin(2nf?n/fg)

根据上面结论，要区分出着两种频率成分，必须满足N400。

1)取x(n)(0≤n≤128)时，计算x(n)的DFTX(k);

2)取x(n)(0≤n≤512)时，计算x(n)的DFTX(k)。

三、总体方案设计

1.信号的离散傅里叶变换

有限长序列的离散傅里叶变换公式

MATLAB函数：fft功能是实现快速傅里叶变换，fft函数的格式为：y=fft(x),返回向量x的不连续fourier变化。ifft功能是实现快速反傅里叶变换，ifft函数的格式为：y=ifft(x),返回向量x的不连续inversefourier变化。

)是一个N=12的有限序列，利用MATLAB计算x(n)它

的傅里叶变换X(k)并画出图形(见仿真结果中图1和图2),然后再对

X(k)进行离散傅里叶反变换，并求出x(n)画出其波形(见仿真结果中图3和图4)。

2.频率分辨率与DFT参数的选择

在DFT问题中，频率分辨率是指在频率轴上所能得到的最小频率间隔,即最小频率间隔反比于数据的长度N。若在x(n)中的两个频率

分别为f和f?的信号，对x(n)用矩形窗截断，要分辨出这两个频率，N

必须满足

通过下面实验，验证上面的结论：设一序列中含有两种频率成分，

f?=2Hz,f?=2.05Hz,采样频率取为f=10Hz,表示为

x(n)=sin(2πf?n/fg)=sin(2πf?n/f)

根据上面结论，要区分出着两种频率成分，必须满足N400。

1)取x(n)(0≤n≤128)时，计算x(n)的DFTX(k),并画出x(n)和x(n)的DFTX(k)的图形(见仿真结果中图5、图6)

2)取x(n)(