[已解决] Matlab 拉普拉斯逆变换问题

[复制链接]
我打这里走过 发表于 2022-5-8 22:28:12
看论文发现如图公式及其拉普拉斯逆变换 屏幕截图 2022-05-08 214921.png ,将其带入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的指数位置似乎行不通。
不知道哪里出了问题,查阅很多资料,没有能解答疑惑的。

如能解惑,感激不尽!

最佳答案


maple1314168 发表于 2022-5-8 23:11:50
syms a c s y r i
assume(y>0)
a=c*exp(-s^(1/2)*y)/s^(1/2);
i=ilaplace(a,s,r);
回复此楼

4 条回复


maple1314168 发表于 2022-5-8 23:11:50
syms a c s y r i
assume(y>0)
a=c*exp(-s^(1/2)*y)/s^(1/2);
i=ilaplace(a,s,r);
回复此楼

我打这里走过 发表于 2022-5-9 10:11:47
maple1314168 发表于 2022-5-8 23:11
syms a c s y r i
assume(y>0)
a=c*exp(-s^(1/2)*y)/s^(1/2);

您好,为什么y>0呢,可否解释下呀。因为后续把字母作为常量参加运算的情况很多,比如里面还会蹦出e f g h,搞不懂如何约束。

maple1314168 发表于 2022-5-10 21:33:32
我打这里走过 发表于 2022-5-9 10:11
您好,为什么y>0呢,可否解释下呀。因为后续把字母作为常量参加运算的情况很多,比如里面还会蹦出e f g h ...

这个是对符号工具箱的要求。知道的信息多,最好让软件知道。
否则,很难简化。最简单的开平方,不知道正负号如何展开?不展开可能面对复杂的计算。如果没有已知信息的话,可以先试试mathematica或maple。

我打这里走过 发表于 4 天前
maple1314168 发表于 2022-5-10 21:33
这个是对符号工具箱的要求。知道的信息多,最好让软件知道。
否则,很难简化。最简单的开平方,不知道正 ...

恩恩,感谢~
您需要登录后才可以回帖 登录 | 注册

本版积分规则

相关帖子
相关文章
热门教程
站长推荐
快速回复 返回顶部 返回列表