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



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）;%余弦項(xiàng)矩陣，大小為n*size（t）
AN=（2*E）*sin（N*pi*t0/T1）./（N*pi）;%余弦項(xiàng)的系數(shù)，是一個(gè)大小為1*n的矩陣

AN_1=meshgrid（ANones（size（t）））‘;%將余弦系數(shù)擴(kuò)展成一個(gè)大小為n*si