function?err=fs（n）
%?n=100;
%?n為傅里葉級數
E=1;%E為方波信號的高度，
T1=0.5;%T1為方波的周期;
Fn=2000;%Fn為方波的采樣頻率
t0=T1/2;%t0為脈沖的寬度
T=2;%T為方波的截取的時間長度
t?=?-T/2:1/Fn:T/2;%t為橫軸的時間坐標
squ?=?E*（1+square（2*pi/T1*t））/2;%產生方波脈沖信號，周期是T1，高度是E
squ=cat（2squ（T1/4*Fn:end）squ（1:T1/4*Fn-1））;%把方波信號移位1/4周期，形成偶對稱的脈沖



a0=E*t0/T1;
N=1:n;
%?AN=（2*E*t0/T1）.*sin（N*pi*t0/T1）./（N*pi*t0/T1）;
COSIN=cos（2*pi/T1*N‘*t）;%余弦項矩陣，大小為n*size（t）
AN=（2*E）*sin（N*pi*t0/T1）./（N*pi）;%余弦項的系數，是一個大小為1*n的矩陣

AN_1=meshgrid（ANones（size（t）））‘;%將余弦系數擴展成一個大小為n*si