离散化傅里叶级数编程求各次谐波含量

 

%x=sin(w.*t)+0.3.*sin(3.*w.*t)+0.1.*sin(5.*w.*t)+0.08.*sin(7.*w.*t);

tic%开始计时

%取样
fs=256;%采样频率
N=256;%采样点数
T=N./fs;
dett=1./fs;
t=0:dett:T-dett;

f=10;%频率 hz
w=2.*pi.*f;
X=sin(w.*t)+0.3.*sin(3.*w.*t)+0.1.*sin(5.*w.*t)+0.08.*sin(7.*w.*t);
subplot(2,1,1)
stem(t,X)
title('取样点');


Xn=ones(1,N);
Z=0;
for k=0:N-1
    for n=0:N-1
        Y=X(n+1).*exp(-j.*2.*pi./N).^(n.*k);
        Z=Z+Y;
        Xn(k+1)=Z;
    end
    Z=0;
end

%Xn=fft(X);
Xn=Xn./N.*2;
Xn(1)=Xn(1)./2;
Xn(N)=Xn(N)./2;

subplot(2,1,2)
stem(t,abs(Xn))
title('离散傅里叶变换')

toc%结束计时