应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 软件使用与操作 » 师兄师姐们，帮小弟看看这个 Laplace变换后的二阶常微分方程怎么解啊
242/3返回列表上一页123下一页
查看: 2248  |  回复: 23
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

月只蓝

主管区长 (职业作家)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
10楼: Originally posted by 吕春莹 at 2014-07-20 19:54:58
好吧，一开始我们老师也是教我们maple，后来就用matlab了，加上要用神经网络什么的，后来我看到另外一个老师用mathmatic的积分和求导等式子的形式漂亮多了，但是找不到安装包...

安装包小木虫里就有的吧。
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
11楼2014-07-20 20:00:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

孙小辉

铜虫 (小有名气)
引用回帖:
6楼: Originally posted by 月只蓝 at 2014-07-20 16:21:44
结果见附图。

附图1.png

第一个解的时候，这种明显不是贝塞尔方程的形式（一次项分母有x），这位师兄，能不能和小弟我说一下，第一种解形式是怎么得来的……在线等……

[ 发自小木虫客户端 ]
赞一下 回复此楼
12楼2014-07-21 07:46:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

孙小辉

铜虫 (小有名气)
引用回帖:
7楼: Originally posted by 吕春莹 at 2014-07-20 19:08:10
我在3楼写的是用matlab算
syms a s y x
Y=dsolve('a*D2y-Dy/x-s*y=0','x')

解析解是
Y =
C3*x^((a + 1)/(2*a))*besselk(-(a + 1)/(2*a), -(s^(1/2)*x)/a^(1/2)) + (C2*x^((a + 1)/(2*a))*besseli(-(a + 1) ...

嗷嗷，太谢谢了，不过是什么原因呢，为什么重新代入就不行了，，，，
赞一下 回复此楼
13楼2014-07-21 08:32:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

孙小辉

铜虫 (小有名气)
引用回帖:
11楼: Originally posted by 月只蓝 at 2014-07-20 20:00:59
安装包小木虫里就有的吧。...

我刚才表达的不准确，师兄，为什么正常解法和laplace解法得到的解析解形式不一样，还有，第一种解析解的求解过程，师兄您有相关的资料吗
赞一下 回复此楼
14楼2014-07-21 08:35:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
14楼: Originally posted by 孙小辉 at 2014-07-21 08:35:30
我刚才表达的不准确，师兄，为什么正常解法和laplace解法得到的解析解形式不一样，还有，第一种解析解的求解过程，师兄您有相关的资料吗...

这是软件解出来的。具体解的过程涉及的专业知识，我并不是很了解。
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
15楼2014-07-21 08:41:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

孙小辉

铜虫 (小有名气)
引用回帖:
6楼: Originally posted by 月只蓝 at 2014-07-20 16:21:44
结果见附图。

附图1.png

这位师兄，送佛送到西啊。得到Y（x)的解析解之后，用Maple软件进行：
（1）比如由Y(0)=1;Y(1)=2确定出Y（x)的C1,C2
（2）画出0<x<10之间的Y（x）曲线。
谢谢啦。。。
赞一下 回复此楼
16楼2014-07-21 16:28:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

孙小辉

铜虫 (小有名气)
引用回帖:
4楼: Originally posted by 吕春莹 at 2014-07-20 11:09:10
但是，反过来 ans=a*diff(Y,x,2)-diff(Y,x)/x-s*Y

ans =

a*(C3*x^((a + 1)/(2*a))*((s^(1/2)*((besselk(1 - (a + 1)/(2*a), -(s^(1/2)*x)/a^(1/2))*((a + 1)/(2*a) - 1))/x + (s^(1/2)*besselk(-(a + 1)/( ...

这位师姐，得到Y（x)的解析解之后，用Matlab进行：
（1）比如由Y(0)=1;Y(1)=2确定出Y（x)的C1,C2
（2）画出0<x<10之间的Y（x）曲线。
谢谢啦。。。
                                   求截图
赞一下 回复此楼
17楼2014-07-21 16:29:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
16楼: Originally posted by 孙小辉 at 2014-07-21 16:28:09
这位师兄，送佛送到西啊。得到Y（x)的解析解之后，用Maple软件进行：
（1）比如由Y(0)=1;Y(1)=2确定出Y（x)的C1,C2
（2）画出0<x<10之间的Y（x）曲线。
谢谢啦。。。...

Y(0)=1,Y(1)=2？
有没有Y'(0)=的数据？
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
18楼2014-07-21 16:49:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
16楼: Originally posted by 孙小辉 at 2014-07-21 16:28:09
这位师兄，送佛送到西啊。得到Y（x)的解析解之后，用Maple软件进行：
（1）比如由Y(0)=1;Y(1)=2确定出Y（x)的C1,C2
（2）画出0<x<10之间的Y（x）曲线。
谢谢啦。。。...

Y(0)=1;Y(1)=2的条件的话，就是边值问题，哪有（或者很难）有解析解。
Y(0)=,Y'(0)= 的条件是初值问题，可以试着找解析解。
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

孙小辉
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
19楼2014-07-21 16:54:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

吕春莹

铜虫 (初入文坛)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
17楼: Originally posted by 孙小辉 at 2014-07-21 16:29:55
这位师姐，得到Y（x)的解析解之后，用Matlab进行：
（1）比如由Y(0)=1;Y(1)=2确定出Y（x)的C1,C2
（2）画出0<x<10之间的Y（x）曲线。
谢谢啦。。。
                                   求截图...

你的方程x在分母位置，x不可以为0，所以你的条件Y（0）=1需要修改
用matlab解方程组用
x=1;y1=subs(Y);
x=2;y2=subs(Y);
[C1,C2]= solve（‘y1=1’，‘y2=2’）
解出C1，C2，你再代回去
y=subs(Y)
x=[1:0.2:10]
plot(Y)
给我红花吧^^
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

孙小辉
20楼2014-07-21 18:52:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 孙小辉 的主题更新
242/3返回列表上一页123下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 娱乐 +3 Tide man 2026-08-03 3/150 2026-08-04 09:58 by gltch
[基金申请] 面上提前没消息，有中的吗 +14 archvillain 2026-08-02 17/850 2026-08-04 09:46 by archvillain
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +3 f1dmPx051StA 2026-08-03 7/350 2026-08-04 09:31 by v1TJC4qJauAB
[博后之家] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 jRl3mE6ddZGq 2026-08-03 5/250 2026-08-04 09:27 by v1TJC4qJauAB
[教师之家] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +5 cu9Nq1xK233Z 2026-08-03 10/500 2026-08-04 09:15 by v1TJC4qJauAB
[基金申请] 系统今天又提示维护了，估计离放榜不远了 +12 winnerche 2026-07-29 16/800 2026-08-04 08:18 by Equinoxhua
[有机交流] 一个有机合成实验室都需要哪些设备？ 50+3 kf2781974 2026-07-31 10/500 2026-08-04 08:05 by 88817753
[找工作] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 9/450 2026-08-04 07:37 by k7OM8YghWbkC
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +4 cu9Nq1xK233Z 2026-08-03 13/650 2026-08-04 07:30 by k7OM8YghWbkC
[考博] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 cu9Nq1xK233Z 2026-08-03 6/300 2026-08-04 06:10 by k7OM8YghWbkC
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 jRl3mE6ddZGq 2026-08-03 6/300 2026-08-04 04:47 by k7OM8YghWbkC
[公派出国] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 7/350 2026-08-04 04:46 by k7OM8YghWbkC
[考研] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 8/400 2026-08-04 03:21 by k7OM8YghWbkC
[基金申请] 什么时候能放榜呀？ +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
[基金申请] 面上再次挂了，太难了，躺也躺不了，倦也卷不过，小学校之殇！ +19 低垂的野花 2026-07-31 25/1250 2026-08-03 09:25 by gy116024
[基金申请] 2026年国自然面上资助率 +16 布布和一二 2026-07-30 22/1100 2026-08-03 09:20 by gy116024
[基金申请] 微信指数没变化，科研之友没阅读 +15 wangze12014 2026-07-28 19/950 2026-08-02 20:05 by 蔡棒棒菂
[考博] 2027年申博 50+3 射雕英雄胜 2026-07-30 3/150 2026-08-02 09:36 by lfy8008
[高分子] HXDI做水性聚氨酯乳液，是不是特别容易出渣 15+3 yuyusuv 2026-07-29 3/150 2026-07-31 09:21 by huizingga
[基金申请] 你们的时间戳变了吗 +3 archvillain 2026-07-30 4/200 2026-07-30 18:53 by levinzhwen
信息提示
关闭
请填处理意见
关闭