%畫原方程為一個三角波
%syms t w y p;
%y=(sin((1*t)/2)./t).^2
%p=fourier(y,t,w)
w=linspace(-pi/2,pi/2,10000) ;
p=-1/4*pi*(w-1).*(Heaviside(-w,-1)-Heaviside(w,+1))+1/2*pi*w.*(Heaviside(-w,0)-Heaviside(w,0))-1/4*pi*(w+1).*(Heaviside(-w,+1)-Heaviside(w,-1))
plot(w,p)