看论文发现如图公式及其拉普拉斯逆变换 ,
将其带入Matlab R2016b计算验证(win10系统64bit),却计无法复现结果,代码如下
syms a c s y r i
a=c*exp(-s^(1/2)*y)/s^(1/2);
i=ilaplace(a,s,r);
i
输出
c*ilaplace(exp(-s^(1/2)*y)/s^(1/2), s, r)
核心在于,y在此要当做与拉普拉斯逆变换无关的常量。但是一旦上述写法,就无法求出;
反而,将y变为具体数或者赋个值后,就可以计算出结果,比如
y=1时,
i =
(c*exp(-1/(4*r)))/(r^(1/2)*pi^(1/2))
但是根据论文,y实际上属于变换不相干变量,计算过程中相当于常量,但是将其放在e的指数位置似乎行不通。
不知道哪里出了问题,查阅很多资料,没有能解答疑惑的。
如能解惑,感激不尽!