untitled.m

来自「利用矩阵分析的方法对线性方程组进行求解分析」· M 代码 · 共 7 行

%画原方程为一个三角波
%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)