请问各位，我在进行Laplace逆变换正弦荷载时，在一定时间（10^10 s）以后结果不好，是为什么？谢谢！
代码如下：
<< NumericalInversion.m
P[s_] := Qs = LaplaceTransform[100*(1 + Sin[2 \[Pi]*10^-6*t]), t, s]
time = Table[10^i, {i, 0, 10, 0.01}]
ft = Crump[P[s], s, time]
fig = ListPlot[Transpose[{Log10[time], ft}], PlotRange -> All,
Joined -> True, InterpolationOrder -> 3]

图是这样的

确切的说是10的7次方以后就不行了，按理说正弦荷载是不会收敛的

从图上看大概是数值误差的问题，把定义改成
P[s_] = LaplaceTransform[100*(1 + Sin[2 \[Pi]*10^-6*t]), t, s];
试下
如果还不行的话再提高time的精度，比如time = SetPrecision[Table[10^i, {i, 0, 10, 0.01}], 30]

问题解决了，这个是程序包的问题，这个问题用InverseLaplaceTransform就可以很好的解决了，谢谢各位！

……就我的使用经验来说，FT函数和GWR函数应该算是目前已有的数值拉普拉斯反变换程序包里性能比较好的，当然了，在你这个问题上它们也跪了。（10^10这个参数过于极端。）考虑到数值拉普拉斯反变换本身就是病态问题，我并不期待有哪个纯数值方法可以扛得住任意极端情况。靠符号计算来规避这一困难对这个问题固然是可行的，不过万一出现了找不到符号解的情况嘛……这个等有人问了再说吧。
顺便楼主所说的程序包，以及FT和GWR，都可以比较容易地在网上找到的，但是这不是你不给链接的理由：
http://library.wolfram.com/infocenter/MathSource/2691/
http://library.wolfram.com/infocenter/MathSource/4738/
http://library.wolfram.com/infocenter/MathSource/5026/

顺便发一下我用InverseLaplaceTransform搞的程序：
<< NumericalInversion.m
P[s_] := Qs = LaplaceTransform[100*(1 + Sin[2 \[Pi]*10^-6*t]), t, s]
time = Table[10^i, {i, 0, 10, 0.002}]
ft = InverseLaplaceTransform[P[s], s, time]
fig = ListPlot[Transpose[{Log10[time], ft}], PlotRange -> All,
Joined -> True, InterpolationOrder -> 3]