24小时热门版块排行榜    

查看: 1950  |  回复: 8

swustxxl

金虫 (正式写手)


[交流] 【求助】拉普拉斯求解~~~

鄙人非数学专业,想用拉普拉斯变换解一个微分方程组,由于太过复杂常用的拉普拉斯变换表查不到对应的象函数,上网查一些非常用的象函数也没找到,所以被卡住了~~~请问这题该怎么解,不胜感激~~~



原题已用附件上传~~~

[ Last edited by swustxxl on 2011-2-10 at 16:21 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

swustxxl

金虫 (正式写手)


没人回答呀
2楼2011-02-10 17:06:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

leedobb

金虫 (正式写手)


其实很简单啊,就是解线性方程组吗,
上式可进一步简化为

sI* V = K* V + B* V_0

其中V=(L(x1),L(x2),...)向量,I为单位矩阵,K为系数矩阵,V_0为初值向量

这样可以解得
V= (sI-K)^(-1)*B*V_0

然后用matlab的符号计算算算结果,就行了
3楼2011-02-15 17:29:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

swustxxl

金虫 (正式写手)


引用回帖:
Originally posted by leedobb at 2011-02-15 17:29:13:
其实很简单啊,就是解线性方程组吗,
上式可进一步简化为

sI* V = K* V + B* V_0

其中V=(L(x1),L(x2),...)向量,I为单位矩阵,K为系数矩阵,V_0为初值向量

这样可以解得
V= (sI-K)^(-1)*B*V_0

然后 ...

等了许久了。。。还是先谢谢
你可能没搞清楚我的问题,你说的这个线性方程我当然会解,不用matlab,我手工都能。问题是之后怎么办?拉普拉斯变换中最重要的不是求L(X1),L(X2),...而是怎么把象函数还原为原函数,即怎么求得X1,X2,X3...

现在我可以求出L(X),但因为在拉普拉斯变换表中找不到对应的原函数,所以才发帖求助
4楼2011-02-16 18:31:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

leedobb

金虫 (正式写手)


那你把解得的结果拿出来看看,一般而言变回去最重要的一个技巧是把
两次方以上的多项式化成若干个
a/(bS+c)的分式相加最可以了。
5楼2011-02-16 19:49:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

信彼南山

木虫 (著名写手)


太长了
matlab里面不是有ilaplace么
6楼2011-02-16 20:46:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

swustxxl

金虫 (正式写手)


引用回帖:
Originally posted by 信彼南山 at 2011-02-16 20:46:14:
太长了
matlab里面不是有ilaplace么

结果很复杂,没法化成加和形式,我不会用matlab,通常拉普拉斯变换我都是手工做,如果你会的话帮我解下好吗?我愿意多出金币

[ Last edited by swustxxl on 2011-2-16 at 22:50 ]
7楼2011-02-16 22:46:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

信彼南山

木虫 (著名写手)


这个我感觉不是特别难,就是麻烦点
而且不知道你为什么要求这个解析解,那些k应该都是具体数值吧?
如果有具体的数值,直接用simulink做一下就可以了。

如果非要求这个表达式的话,相对麻烦一点
8楼2011-02-16 23:06:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ghw_nit

铁杆木虫 (正式写手)


swustxxl(金币+10): 用mathematica也可以解,但解析形式太复杂,没什么实际意义,不过还是多谢 2011-02-19 09:34:18
你到控制的版面吧,这是一个传递函数求解的问题,其实很简单的
9楼2011-02-17 08:53:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 swustxxl 的主题更新
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 上海工程技术大学 激光智能制造课题组 2027级博士研究生招生 +5 oxidpl 2026-09-20 7/350 2026-09-27 10:43 by Equinoxhua
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 ViqFlLSxHDGt 2026-09-26 3/150 2026-09-27 10:42 by 4IW0sJvtEqX8
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 ViqFlLSxHDGt 2026-09-26 3/150 2026-09-27 10:21 by 4IW0sJvtEqX8
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 ViqFlLSxHDGt 2026-09-26 3/150 2026-09-27 10:15 by 4IW0sJvtEqX8
[考博] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 ViqFlLSxHDGt 2026-09-26 3/150 2026-09-27 10:03 by 4IW0sJvtEqX8
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 oeyPlfmMUpOC 2026-09-26 3/150 2026-09-27 09:13 by 4IW0sJvtEqX8
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 GxESXJmptOtQ 2026-09-25 4/200 2026-09-27 06:42 by XgCC8uTcwILl
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 GxESXJmptOtQ 2026-09-25 4/200 2026-09-27 06:39 by XgCC8uTcwILl
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 GxESXJmptOtQ 2026-09-25 3/150 2026-09-27 06:32 by XgCC8uTcwILl
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 vQZoDrrm7mUF 2026-09-25 3/150 2026-09-27 06:08 by XgCC8uTcwILl
[公派出国] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 vQZoDrrm7mUF 2026-09-25 3/150 2026-09-27 06:05 by XgCC8uTcwILl
[论文投稿] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 vQZoDrrm7mUF 2026-09-25 4/200 2026-09-27 05:56 by XgCC8uTcwILl
[找工作] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +5 vQZoDrrm7mUF 2026-09-25 5/250 2026-09-27 05:53 by XgCC8uTcwILl
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 vQZoDrrm7mUF 2026-09-25 4/200 2026-09-27 05:50 by XgCC8uTcwILl
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 vQZoDrrm7mUF 2026-09-25 4/200 2026-09-27 05:44 by XgCC8uTcwILl
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 vQZoDrrm7mUF 2026-09-25 3/150 2026-09-27 05:44 by XgCC8uTcwILl
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 vQZoDrrm7mUF 2026-09-25 6/300 2026-09-27 05:42 by XgCC8uTcwILl
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 PGTSwSC3F6nQ 2026-09-24 5/250 2026-09-25 11:19 by yhq9807
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 PGTSwSC3F6nQ 2026-09-24 7/350 2026-09-25 11:18 by yhq9807
[基金申请] 师弟论文见刊大半年才想起来申请专利,还能抢救一下吗? +4 13108017953 2026-09-21 4/200 2026-09-24 10:06 by kudofaye
信息提示
请填处理意见