利用MATLAB实现信号DFT的计算

1、07 级电信（ 2）班刘坤洋 24试验一利用 matlab实现信号dft的运算一、试验目的：1、熟识利用 matlab运算信号 dft的方法2、把握利用 matlab实现由 dft运算线性卷积的方法二、试验设备： 电脑、 matlab 软件三、试验内容：1、练习用 matlab 中供应的内部函数用于运算dft（1） fft（ x），fft（ x， n），ifft（x），ifft（ x， n）的含义及用法（2） 在进行 dft时选取合适的时域样本点数n请举例，并编程实现题目：利用 dft 分析xt cos4f1tcos4f2t 的频谱f1100hz, f 2120hz, fs600hz源程序：

2、>> n=30;%数据的长度>>l=512;%dft的点数>>f1=100; f2=120;>>fs=600;%抽样频率>>t=1/fs;%抽样间隔>>ws=2*pi*fs;>>t=0:n-1*t;>>f=cos4*pi*f1*t+cos4*pi*f2*t;>>f=fftshiftfftf,l;>>w=-ws/2+0:l-1*ws/l/2*pi;>>hd=plotw,absf;>>ylabel'幅度谱 '>> xlabel&

3、#39;频率/hz'>> title'my picture'结果图：（3） 在对信号进行 dft时挑选 hamming窗增加频率辨论率请举例，并编程实现题目：利用 dft 分析 xt cos4f1t0.15cos4f2 t的频谱f1100 hz , f 2150 hz , f s600hz源程序： >> n=50;%数据的长度>>l=512;%dft的点数>>f1=100;f2=150;>>fs=600;%抽样频率>>t=1/fs;%抽样间隔>>ws=2*pi*fs;>>t=

4、0:n-1*t;>>f=cos4*pi*f1*t+0.15*cos4*pi*f2*t;>>wh=hammingn'>>f=f.*wh;>>f=fftshiftfftf,l;>>w=-ws/2+0:l-1*ws/l/2*pi;>>plotw,absf;>>ylabel'幅度谱 '>> xlabel'频率/hz'>> title'my picture'>> legend'n=50'结果图：2、增加 dft点数

5、 m以显示更多频谱细节请举例，并编程实现题目：x kcos4rk/ n ,n32,r4利用 matlab运算 16 点序列 x k 的 512 点 dft；源程序： >> n = 32;>> k = 0:n-1;>> l = 0:511;>> x = cos4*pi*k*4./n;>> x = fftx;plotk/n,absx,'o'>> hold on>> xe = fftx,512;plotl/512,absxe ;>> legend'n=32'结果图：3、利用

6、 matlab实现由 dft 运算线性卷积请举例，并编程实现题目：利用ma tlab由 dft 运算 x k* hk；xk=1, 2, 1, 1,h k=2, 2, 1, 1源程序 ： %利用 dft 运算线性卷积>> x = 1 2 1 1;>>h = 2 2 1 1;>> %确定卷积结果的长度>>l = lengthx+lengthh-1;%运算序列的 l 点 dft>>xe = fftx,l;>.he = ffth,l;>> %由 idft运算卷积结果>>y1 = ifftxe.*he;>> % 绘制卷积结果及误差图>> k=0:l-1;>> subplot1,2,1;>> stemk,realy1;axis0 6 0 7;>> title'result of linear convolution'>> xlabel'time index k'ylabel'amplitude'>> y2=convx,h;error=y1-y2;>> subplot1,2,2;stemk,abserr