fourier_seriers.asv

来自「matlab从入门到精通第6章书后源码对初学者很有帮助」· ASV 代码 · 共 21 行

% 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('')