fourier_seriers.m

matlab从入门到精通第6章书后源码对初学者很有帮助

% fourier_seriers.m
% 创建25个谐波的傅立叶级数
N=25;To=0.2;
% 周期数
n=2*N;
t=linspace(0,To,n+1);
t(end)=[];
f=sawtooth(t,To);
% 复指数形式的傅立叶级数
disp('复指数形式的傅立叶级数：')
Fn=fft(f);
Fn=[conj(Fn(N+1)) Fn(N+2:end) Fn(1:N+1)];
Fn=Fn/n
% 计算三角形式的傅立叶级数
disp('三角形式的傅立叶级数：')
A0=Fn(N+1)
An=2*real(Fn(N+2:end))
Bn=-2*imag(Fn(N+2:end))
% 实际计算结果和转换结果比较
disp('比较：')
idx=-N:N;
Fna=5j./(idx*pi);
Fna(N+1)=5;
Bna=-2*imag(Fna(N+2:end));
Bn_error=(Bn-Bna)./Bna
% 绘制频谱
stem(idx,abs(Fn))
xlabel('Harmonic Index')
title('sawtooth harmonic content')
axis tight;