DFT变换.docVIP

地球物理信息处理基础 DFT--实验作业报告 一、 实验目的： 1、 信号：f(t)=2+3cos(2∏*50t-30)+1.5cos(2∏*75t+90)。对该信号波形有直观的认识。 2、 利用DFT将时间信号变换至频率域， 并绘制振幅谱、相位谱。 3、 将50Hz频率成分滤去，绘制滤波后的振幅谱、相位谱。 4、 将频谱反变换回时间域，对比前后波形。 5、 分析如何获得高分辨率谱与减弱截断效应。 二、实验方法： 1、利用matlab编写程序。 2、根据原理，不直接使用函数，实现dft变换，dft反变换以及滤波过程。 3、将这部分利用欧拉公式分成实虚部，然后分别做dft。那么变换到频率域的振幅谱就可以直接利用实部的平方加虚部的平方再开根号来求，相位谱也可以直接利用虚部除以实部再求反正切得到。这里需要注意时间域的离散造成频率域的周期化，出现对称现象。 4、滤波过程的关键在于求取50Hz对应的K值，以及对称的K1值，将滤波器充零。 5、反变换得到的频谱既有实部也有虚部，我们只取实部。 实验过程： 具体见如下程序 %地球物理信息处理11月15日上机作业 clc clf clear all %绘制波形图 N=500; %500个采样点 ddt=0.002; %采样间隔是0.002 dt=0:ddt:N*ddt-ddt; df=1/(N*ddt); t=N.*dt; f1=