[讨论]自己编傅里叶变换和逆变换出现问题

[blackclaw]

问题是这样的 0<t<10^-8时间里 Vi=1，其余时间Vi=0 按照傅里叶变换的定义做它的傅里叶变换 画w vs F(Vi)的图像 再做逆变换，做逆变换后的f(t)图像 f(t)应该跟原函数图像大致相同


首先做傅里叶变换，我在做图时只选择显示f区间为 -5 e+008 到 +5 e+008 得到了比较满意的图像，code如下
syms t w;

f=(-5*10^8:10^6:5*10^8);
w=2*pi*f;

FVi=int(exp(-i.*w*t),t,0,10^-8); %Fourier tranform of Vi(t)

figure(1);
plot(w,subs(FVi))

接下来使用FVi进行逆变换却出现了问题 代码如下：
t=(0:2*10^-8/(length(w)-1):2*10^-8);
a=FVi.*exp(i.*w.*t);
IFVi=1/(2*pi)*int(a,w)

figure(2)
plot(t,subs(IFVi))
title('Inverse Fourier transform of FVi')

提示出错：
??? Input argument "NEWexpr" is undefined.

Error in ==> sym.subs at 104
msg = inputchk(OLDf,OLDexpr,NEWexpr);



请问这样计算思路正确么？ 逆运算时错误在哪里？ 多谢！