离散信号X(n)变换的matlab实现.docVIP

实验四离散信号的频域分析

实验目的

掌握序列的傅里叶变换、离散傅里叶级数、离散傅里叶变换、快速傅里叶变换的Matlab实现；

学习用FFT对连续信号和离散信号进行谱分析的方法，了解可能出现的分析误差及其原因，以便在实际中正确应用FFT。



实验内容及步骤

1.计算序列的DTFT和DFT，观察栅栏效应

设，要求用MATLAB实现：

〔1〕计算的傅里叶变换，并绘出其幅度谱；

〔2〕分别计算的4点DFT和8点DFT，绘出其幅度谱。并说明它们和的关系。

〔提示：DFT变换可用MATLAB提供的函数fft实现，也可以自己用C语言或matlab编写〕

2.计算序列的FFT，观察频谱泄漏

周期为16的信号。

截取一个周期长度M=16点，计算其16点FFT，并绘出其幅度谱；

截取序列长度M=10点，计算其16点FFT，绘出其幅度谱，并与〔1〕的结果进行比拟，观察频谱泄漏现象，说明产生频谱泄漏的原因。

实验报告要求

结合实验中所得给定典型序列幅频特性曲线，与理论结果比拟，并分析说明误差产生的原因以及用FFT作谱分析时有关参数的选择方法。

总结实验所得主要结论。

1.计算序列的DTFT和DFT，观察栅栏效应

设，要求用MATLAB实现：

〔1〕计算的傅里叶变换，并绘出其幅度谱；

〔2〕分别计算的4点DFT和8点DFT，绘出其幅度谱。并说明它们和的关系。

代码：

n=0:3;

M=10;

k=0:2*M-1;%观察两个周期

x=u(n)-u(n-4);

X=x*(exp(-j*2*pi/M)).^(n*k);%序列的傅里叶变换

f1=abs(X);

subplot(211);

stem(n,x,fill);

title(R4(n)序列);

xlabel(n);

gridon;

subplot(212);plot(f1);

title(r4(n)的幅度谱);

xlabel(abs(R4(n)));

gridon

〔2〕代码

n=0:3;

M=4;

H=8;

m=100;

k=0:2*M-1;%观察两个周期

k1=0:2*H-1;%观察两个周期

k2=0:2*m-1;%观察两个周期

x=u(n)-u(n-4);

X=x*(exp(-j*2*pi/M)).^(n*k);%序列的傅里叶变换

Y=x*(exp(-j*2*pi/H)).^(n*k1);%序列的傅里叶变换

q=x*(exp(-j*2*pi/m)).^(n*k2);%序列的傅里叶变换

f1=abs(X);

subplot(311);

stem(k,f1,fill);

title(4点DFT幅度谱);

xlabel(abs(R4(n)));

gridon;

f2=abs(Y);

subplot(312);

stem(k1,f2,fill);

title(8点DFT幅度谱);

xlabel(abs(R4(n)));

gridon

subplot(313);

f3=abs(q);

plot(f3);

title(R4(n)DFT幅度谱);

gridon;

根据图像可以看出，4点和8点的DFT是原序列的傅里叶变换在周期上进行等分取值。例如4点〔8点〕就是将其分为4〔8〕等分。

2.计算序列的FFT，观察频谱泄漏

周期为16的信号。

截取一个周期长度M=16点，计算其16点FFT，并绘出其幅度谱；

截取序列长度M=10点，计算其16点FFT，绘出其幅度谱，并与〔1〕的结果进行比拟，观察频谱泄漏现象，说明产生频谱泄漏的原因。

〔3〕

代码：

n=0:15;

x=cos(10*pi/16*n)+cos(12*pi/16*n);

subplot(211);

stem(n,x,fill);

title(x序列图像);

xlabel(n);

gridon;

f2=fft(x);

subplot(212);

stem(n,f2,fill);

title(16点FFT的幅度谱);

xlabel(abs(X(k)));

gridon;

(4)代码：

clc;clear;clf

n=0:15;

m=0:9;

x=cos(10*pi/16*n)+cos(12*pi/16*n);

subplot(211);

stem(n,x,fill);

title(x序列图像);

xlabel(n);

gridon;

x=[cos(10*pi/16*m)+cos(12*pi/16*m),0