MATLAB的离散非周期信号频域分析报告
基于MATLAB的离散信号频域分析、快速傅里叶变换与采样定理
离散信号频域分析
周期离散方波信号频域分析
与周期模拟信号一样,周期离散信号同样可以展开成傅里叶级数形式,并得到离散傅里叶级数(DFS)
X
上式可以看成周期离散信号x(n)的离散傅里叶级数展开。
x
上式是DFS的反变换,记作IDFS并且称X(kΩ)与x(n)构成一对离散傅里叶级数变换对。(以上两式中Ω=2π/N)
在MTALAB中,DFS通过建立周期延拓函数语句实现:
function Xk=DFS(n,x,N)
if N>length(x)
n=0:N-1;
x=[x zeros(1,N-length(x))];
end
k=0:N-1;
WN=exp(-j*2*pi/N);
nk=n'*k;
WNnk=WN.^nk;
Xk=x*WNnk;
end
建立一个离散非周期方波信号
x
R4n通过周期延拓后所得的周期序列利用
clear all;close all;clc;
n=0:3;
x=ones(1,4);
X=fft(x,1024);
Xk1=DFS(n,x,4);
Xk2=DFS(n,x,8);
figure(1);
plot((-1023:2048)/2048*8,[abs(X) abs(X) abs(X)],'--');hol