24小时热门版块排行榜    

查看: 1863  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 BnSlaJi8z5Vo 2026-08-17 3/150 2026-08-17 17:38 by CBFiACZS7HAK
[考博] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 i7NFEVbjQMM5 2026-08-16 5/250 2026-08-17 16:38 by CBFiACZS7HAK
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +34 医学老男孩 2026-08-13 81/4050 2026-08-17 16:20 by yudaoqian88
[基金申请] 时间戳变了,能看出什么问题? +4 基诺咪客 2026-08-17 4/200 2026-08-17 16:10 by Vivilian
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 EWi2p09MOOv8 2026-08-16 6/300 2026-08-17 14:02 by EYK67XJMj64E
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 6GojgJvkDudM 2026-08-16 4/200 2026-08-17 13:49 by EYK67XJMj64E
[基金申请] 今天系统多次维护,明天很可能放榜! +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 i7NFEVbjQMM5 2026-08-16 4/200 2026-08-17 12:10 by 4GBAYCdVQoK3
[考研] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 6GojgJvkDudM 2026-08-16 4/200 2026-08-17 11:46 by 4GBAYCdVQoK3
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 6GojgJvkDudM 2026-08-16 3/150 2026-08-17 09:49 by xLVPIRuSvCUe
[论文投稿] 岩土工程学报什么时候才能终审完呐 25+3 yeager111 2026-08-10 4/200 2026-08-16 20:42 by tfang
[基金申请] 2027广东省杰青 +3 奶牛小黑 2026-08-15 6/300 2026-08-16 20:07 by 奶牛小黑
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +10 majunge000 2026-08-11 12/600 2026-08-16 08:18 by xli1984
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 是这周出结果还是下周出结果? +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
信息提示
请填处理意见