24小时热门版块排行榜    

查看: 2563  |  回复: 14

xu17211

铜虫 (初入文坛)

[交流] 【求助】求解形如y''+a*y'^2+by'+c=0 二阶微分方程 已有2人参与

求解形如y''+a*y'^2+by'+c=0 二阶微分方程,其中a,b,c随输入而变化,会出现b^2-4*a*c>0或<0的情况,a后项为一阶微分平方项。请高手指教,谢谢
若使用matlab直接求解会出现不合适实际的b^2-4*a*c大小判断的问题,不知能否求出适当的解析解,还请各位帮助
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

xiegangmai

版主 (职业作家)

我没头衔

优秀版主优秀版主优秀版主


xu17211(金币+1):谢谢参与
首先将y''+a*y'^2+by'+c=0 二阶微分方程变换为一阶微分方程组的形式
y(1) = y(2);
y(2) = -a * y(1) ^ 2 - b * y(1) - c;

给定具体a、b、c值,可用ode45可求解。
用MATLAB求不了解析解。
明德厚学、求是创新
2楼2010-08-31 20:04:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fyq98

木虫 (正式写手)


xu17211(金币+1):谢谢参与
这类问题的解析解可以通过幂级数法求得,即假设解为幂级数,代入方程,根据对应次幂系数相等,得到一些代数方程,解之即得
克难奋进
4楼2010-08-31 23:26:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

onesupeng

金虫 (职业作家)


xu17211(金币+1):谢谢参与
引用回帖:
Originally posted by fyq98 at 2010-08-31 23:26:56:
这类问题的解析解可以通过幂级数法求得,即假设解为幂级数,代入方程,根据对应次幂系数相等,得到一些代数方程,解之即得

怎么可能?你说这个特征值接法,线性方程才可以用。

至于非线性,一般意义上不会有解析解,特殊情况有。但是对于一些问题,采用渐进分析这一套可能会有用,例如摄动法,同论摄动这些。

实际应用数值计算居多
长期招收博士生,参见http://fsl-unsw.com
5楼2010-09-01 00:17:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

xu17211

铜虫 (初入文坛)

如果不应用matlab进行求解,请问大家能给出合理的解析解的思路吗?
3楼2010-08-31 22:33:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fyq98

木虫 (正式写手)

引用回帖:
Originally posted by onesupeng at 2010-09-01 00:17:33:


怎么可能?你说这个特征值接法,线性方程才可以用。

至于非线性,一般意义上不会有解析解,特殊情况有。但是对于一些问题,采用渐进分析这一套可能会有用,例如摄动法,同论摄动这些。

实际应用数值计算 ...

why not try?
克难奋进
6楼2010-09-01 08:31:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

onesupeng

金虫 (职业作家)

引用回帖:
Originally posted by fyq98 at 2010-09-01 08:31:12:

why not try?

你try try看看
长期招收博士生,参见http://fsl-unsw.com
7楼2010-09-01 12:48:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xxppyy

木虫 (小有名气)


xu17211(金币+1):谢谢参与
y'=p
dp/dt+ap^2+bp+c=0
dp/(ap^2+bp+c)=-dt
两边积分
b^2-4ac<0时:


b^2-4ac>0时: ......
8楼2010-09-01 23:21:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiegangmai

版主 (职业作家)

我没头衔

优秀版主优秀版主优秀版主

用MATLAB的dsolve命令求解

syms a b c
dsolve('D2y + a * Dy ^ 2 + b * Dy + c = 0')

ans =

-(log(b^2 - 4*a*c)/2 - log(C6 + C5*cosh(t*(b^2 - 4*a*c)^(1/2)) + C5*sinh(t*(b^2 - 4*a*c)^(1/2))) - log(a) + (b*t)/2 + (t*(b^2 - 4*a*c)^(1/2))/2)/a

不知道结果对不
明德厚学、求是创新
9楼2010-09-02 10:03:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

scuxchen

金虫 (著名写手)


xu17211(金币+1):谢谢参与
看起简单,真是挺难。一定要求解么?不能用常微分方程的定性理论的一些东西来分析解的性质等等?
拳打脚踢
10楼2010-09-02 11:50:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xu17211 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 086000生物与医药292求调剂 +4 小小陈小小 2026-03-22 7/350 2026-03-25 19:07 by 星空星月
[考研] 材料调剂 +3 iwinso 2026-03-23 3/150 2026-03-25 11:29 by greychen00
[考研] 271求调剂 +4 生如夏花… 2026-03-22 4/200 2026-03-25 11:25 by userper
[考研] 求调剂323材料与化工 +4 1124361 2026-03-24 4/200 2026-03-25 11:19 by shulmg
[考研] 311求调剂 +3 冬十三 2026-03-24 3/150 2026-03-24 21:31 by peike
[考研] 【考研调剂】化学专业 281分,一志愿四川大学,诚心求调剂 +15 吃吃吃才有意义 2026-03-19 16/800 2026-03-24 19:51 by 了了了了。。
[考研] 321求调剂 +4 Ymlll 2026-03-24 4/200 2026-03-24 14:44 by sprinining
[考研] 一志愿吉大化学322求调剂 +4 17501029541 2026-03-23 6/300 2026-03-24 10:21 by 戴围脖的小蚊子
[考研] 341求调剂(一志愿湖南大学070300) +5 番茄头--- 2026-03-22 6/300 2026-03-23 23:45 by Txy@872106
[考研] 327求调剂 +5 prayer13 2026-03-23 5/250 2026-03-23 22:11 by 星空星月
[考研] 336化工调剂 +4 王大坦1 2026-03-23 5/250 2026-03-23 18:32 by allen-yin
[考研] 接收2026硕士调剂(学硕+专硕) +4 allen-yin 2026-03-23 6/300 2026-03-23 15:04 by 汪!?!
[考研] 328求调剂 +4 LHHL66 2026-03-23 4/200 2026-03-23 14:55 by lbsjt
[考研] 263求调剂 +6 yqdszhdap- 2026-03-22 9/450 2026-03-23 12:57 by yqdszhdap-
[考研] 298求调剂 +8 上岸6666@ 2026-03-20 8/400 2026-03-23 11:02 by laoshidan
[考研] 298求调剂一志愿211 +3 上岸6666@ 2026-03-20 3/150 2026-03-22 15:50 by ColorlessPI
[考研] 一志愿南大,0703化学,分数336,求调剂 +3 收到VS 2026-03-21 3/150 2026-03-21 18:42 by 学员8dgXkO
[考研] 求调剂 +3 .m.. 2026-03-21 4/200 2026-03-21 16:25 by barlinike
[考研] 一志愿武汉理工材料工程专硕调剂 +9 Doleres 2026-03-19 9/450 2026-03-20 22:36 by JourneyLucky
[考研] 一志愿 南京航空航天大学大学 ,080500材料科学与工程学硕 +5 @taotao 2026-03-20 5/250 2026-03-20 20:16 by JourneyLucky
信息提示
请填处理意见