sanjiao.m

来自「我做了很久的」· M 代码 · 共 12 行

syms t w
ut=sym('Heaviside(t+0.5)-Heaviside(t-0.5)');
Fw=fourier(ut,t,w)
FFw=maple('convert',Fw,'piecewise')
FFFw=FFw*FFw
FFP=abs(FFFw)
subplot(2,1,1)
ezplot(FFFw)
subplot(2,1,2)
ezplot(FFP,[-10*pi 10*pi])
axis([-10*pi 10*pi 0 1])