sanjiao.m

来自「信号与系统时域、频域分析及MATLAB软件的应用及仿真m程序」· M 代码 · 共 8 行

symst  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)
ezplot(FFP,[-10*pi 10*pi])
axis([-10*pi 10*pi 0 1]);