[交流] 师兄师姐们，帮小弟看看这个 Laplace变换后的二阶常微分方程怎么解啊已有3人参与

» 猜你喜欢

🌟 比利时新鲁汶大学（UCLouvain）诚邀CSC博士加入Pascal Gehring教授团队已经有0人回复
第一性原理计算方向2026级博士申请 PRB*1，四级484 已经有1人回复
物理学I论文润色/翻译怎么收费? 已经有182人回复
求助VISSIM破解版软件已经有0人回复
求2026年在台湾举行的物理和材料领域国际学术会议信息已经有0人回复
求国际会议网站已经有1人回复
求取一些关于纳米材料和纳米技术相关的英文PPT。已经有0人回复
【复旦大学】二维材料方向招收2026年博士研究生1名已经有0人回复
北京纳米能源与系统研究所王中林院士/曹南颖研究员课题组2026级硕/博/博后招生已经有10人回复
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生已经有5人回复

» 本主题相关价值贴推荐，对您同样有帮助:

请求高手给小弟指点一下这个反应！！急！！急！！已经有6人回复

1楼2014-07-17 11:45:43

吕春莹

铜虫 (初入文坛)

应助: 5 (幼儿园)
金币: 82.8
红花: 1
帖子: 24
在线: 19.8小时
虫号: 2815782
注册: 2013-11-20
性别: MM
专业: 核技术及其应用

★
小木虫: 金币+0.5, 给个红包，谢谢回帖

引用回帖:

9楼: Originally posted by 月只蓝 at 2014-07-20 19:32:16
不是哈，是Maple。...

好吧，一开始我们老师也是教我们maple，后来就用matlab了，加上要用神经网络什么的，后来我看到另外一个老师用mathmatic的积分和求导等式子的形式漂亮多了，但是找不到安装包

10楼2014-07-20 19:54:58

查看全部 24 个回答

吕春莹

铜虫 (初入文坛)

应助: 5 (幼儿园)
金币: 82.8
红花: 1
帖子: 24
在线: 19.8小时
虫号: 2815782
注册: 2013-11-20
性别: MM
专业: 核技术及其应用

★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
fegg7502: 金币+1, 鼓励交流 2014-07-21 08:21:23

首先，这不是二阶常微分方程，“常”表示y''，y'，y的系数都是常数，形式为ay''+by'+cy=0,这类问题可以用特征方程求解。你的问题是二阶非线性微分方程。
  我想了下也不知道怎么解，用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)/(2*a), (s^(1/2)*x)/a^(1/2)))/(-1)^((a + 1)/(2*a))

>> pretty(Y)

                                             a + 1
                                             -----       /          1/2 \
                                             2 a       | a + 1  s x |
   a + 1                               C2 x    besseli| - -----, ------ |
   -----       /          1/2 \                   | 2 a    1/2  |
   2 a       | a + 1 s x |                   \          a    /
  C3 x    besselk| - -----, - ------ | + ------------------------------------
               | 2 a    1/2  |                   a + 1
               \          a    /                   -----
                                                         2 a
                                                      (-1)
>> 解中由用到bessel函数

3楼2014-07-20 11:01:11

吕春莹

铜虫 (初入文坛)

应助: 5 (幼儿园)
金币: 82.8
红花: 1
帖子: 24
在线: 19.8小时
虫号: 2815782
注册: 2013-11-20
性别: MM
专业: 核技术及其应用

★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
fegg7502: 金币+1, 鼓励交流 2014-07-22 08:04:26

引用回帖:

3楼: Originally posted by 吕春莹 at 2014-07-20 11:01:11
首先，这不是二阶常微分方程，“常”表示y''，y'，y的系数都是常数，形式为ay''+by'+cy=0,这类问题可以用特征方程求解。你的问题是二阶非线性微分方程。
我想了下也不知道怎么解，用matlab算:
>> syms a ...

但是，反过来 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)/(2*a), -(s^(1/2)*x)/a^(1/2)))/a^(1/2)))/a^(1/2) + (besselk(-(a + 1)/(2*a), -(s^(1/2)*x)/a^(1/2))*(a + 1))/(2*a*x^2) - ((a + 1)*((s^(1/2)*besselk(1 - (a + 1)/(2*a), -(s^(1/2)*x)/a^(1/2)))/a^(1/2) - (besselk(-(a + 1)/(2*a), -(s^(1/2)*x)/a^(1/2))*(a + 1))/(2*a*x)))/(2*a*x)) + (C2*x^((a + 1)/(2*a))*((s^(1/2)*((besseli(1 - (a + 1)/(2*a), (s^(1/2)*x)/a^(1/2))*((a + 1)/(2*a) - 1))/x + (s^(1/2)*besseli(-(a + 1)/(2*a), (s^(1/2)*x)/a^(1/2)))/a^(1/2)))/a^(1/2) + (besseli(-(a + 1)/(2*a), (s^(1/2)*x)/a^(1/2))*(a + 1))/(2*a*x^2) - ((a + 1)*((s^(1/2)*besseli(1 - (a + 1)/(2*a), (s^(1/2)*x)/a^(1/2)))/a^(1/2) - (besseli(-(a + 1)/(2*a), (s^(1/2)*x)/a^(1/2))*(a + 1))/(2*a*x)))/(2*a*x)))/(-1)^((a + 1)/(2*a)) + (C3*x^((a + 1)/(2*a) - 1)*(a + 1)*((s^(1/2)*besselk(1 - (a + 1)/(2*a), -(s^(1/2)*x)/a^(1/2)))/a^(1/2) - (besselk(-(a + 1)/(2*a), -(s^(1/2)*x)/a^(1/2))*(a + 1))/(2*a*x)))/a + (C3*x^((a + 1)/(2*a) - 2)*((a + 1)/(2*a) - 1)*besselk(-(a + 1)/(2*a), -(s^(1/2)*x)/a^(1/2))*(a + 1))/(2*a) + (C2*x^((a + 1)/(2*a) - 1)*(a + 1)*((s^(1/2)*besseli(1 - (a + 1)/(2*a), (s^(1/2)*x)/a^(1/2)))/a^(1/2) - (besseli(-(a + 1)/(2*a), (s^(1/2)*x)/a^(1/2))*(a + 1))/(2*a*x)))/((-1)^((a + 1)/(2*a))*a) + (C2*x^((a + 1)/(2*a) - 2)*((a + 1)/(2*a) - 1)*besseli(-(a + 1)/(2*a), (s^(1/2)*x)/a^(1/2))*(a + 1))/(2*(-1)^((a + 1)/(2*a))*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) - (besselk(-(a + 1)/(2*a), -(s^(1/2)*x)/a^(1/2))*(a + 1))/(2*a*x)) + (C2*x^((a + 1)/(2*a))*((s^(1/2)*besseli(1 - (a + 1)/(2*a), (s^(1/2)*x)/a^(1/2)))/a^(1/2) - (besseli(-(a + 1)/(2*a), (s^(1/2)*x)/a^(1/2))*(a + 1))/(2*a*x)))/(-1)^((a + 1)/(2*a)) + (C3*x^((a + 1)/(2*a) - 1)*besselk(-(a + 1)/(2*a), -(s^(1/2)*x)/a^(1/2))*(a + 1))/(2*a) + (C2*x^((a + 1)/(2*a) - 1)*besseli(-(a + 1)/(2*a), (s^(1/2)*x)/a^(1/2))*(a + 1))/(2*(-1)^((a + 1)/(2*a))*a))/x - s*(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)/(2*a), (s^(1/2)*x)/a^(1/2)))/(-1)^((a + 1)/(2*a)))

这回我不明白了

4楼2014-07-20 11:09:10

孙小辉

铜虫 (小有名气)

应助: 0 (幼儿园)
金币: -10.3
散金: 59
红花: 1
帖子: 108
在线: 63.8小时
虫号: 3310688
注册: 2014-07-07
性别: GG
专业: 石油、天然气地质学

引用回帖:

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)/( ...

您好，能不能发个截图上去，后面这个是它的解析解吗？

5楼2014-07-20 15:14:34