%
% フーリエ級数（方形波）
% https://takuichi.net/hobby/edu/fourier_analysis/fourier.pdf
% のp.24
% MABLAB (by T. Hirano)
%
A=1;    % 振幅
T0=2;   % 周期
tau=1;  % 1の時間
omega0=2*pi/T0; % 基本角周波数
t=-2*T0:(4*T0)/400:2*T0;    % 時間範囲
N=10;   % 級数の和の数
a0=(A*tau)/T0;
y=(t.*0)+a0;   % tと同じ次元のyベクトルを生成し、a0の初期値を代入。
for i=1:N
    ai=(2*A*tau/T0)*sinc(i*omega0*tau/(2*pi));
    dy=ai*cos(i*omega0*t);
    y=y+dy;
end
plot(t,y)