24小时热门版块排行榜    

CyRhmU.jpeg
查看: 316  |  回复: 2
当前主题已经存档。

fengligang

金虫 (小有名气)

[交流] 【求助】一个微分方程的求解

在用分离变量法求解一个二阶线性偏微分方程时,得到下面这样一个二阶常微分方程:
x^2*y''+2*x*y'+(a*x^2+b)*y=0
其中a,b为常数。
怎样解这个微分方程呢?好像求解球体的径向振动时也涉及这个方程。
各位虫友帮忙看一下,谢谢了。
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fspdlh

金虫 (正式写手)

★ ★
点点风(金币+2,VIP+0):感谢参与讨论,加分鼓励! 3-6 16:42
>> y=dsolve('x^2*D2y+2*x*Dy+(a*x^2+b)*y=0','x')

y =

C1/x^(1/2)*besselj(1/2*(1-4*b)^(1/2),a^(1/2)*x)+C2/x^(1/2)*bessely(1/2*(1-4*b)^(1/2),a^(1/2)*x)


>> simplify(y)

ans =

(C1*besselj(1/2*(1-4*b)^(1/2),a^(1/2)*x)+C2*bessely(1/2*(1-4*b)^(1/2),a^(1/2)*x))/x^(1/2)
2楼2009-03-06 13:24:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fengligang

金虫 (小有名气)

引用回帖:
Originally posted by fspdlh at 2009-3-6 13:24:
>> y=dsolve('x^2*D2y+2*x*Dy+(a*x^2+b)*y=0','x')

y =

C1/x^(1/2)*besselj(1/2*(1-4*b)^(1/2),a^(1/2)*x)+C2/x^(1/2)*bessely(1/2*(1-4*b)^(1/2),a^(1/2)*x)


>> simplify(y)

ans ...

非常感谢你,用matlab或mathmatic可以很方便的算出来。
但具体的求解过程是怎样的呢?

[ Last edited by fengligang on 2009-3-6 at 16:47 ]
3楼2009-03-06 16:26:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 fengligang 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见