24小时热门版块排行榜    

查看: 2559  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 资源与环境 调剂申请(333分) +7 holy J 2026-03-21 7/350 2026-03-24 17:24 by xiaohai104
[考研] 材料考研调剂生 +3 黄粱一梦千年 2026-03-24 3/150 2026-03-24 17:00 by barlinike
[考研] 材料292调剂 +8 橘颂思美人 2026-03-23 8/400 2026-03-24 16:33 by laoshidan
[考研] 材料专硕331求调剂 +4 鲜当牛 2026-03-24 4/200 2026-03-24 15:58 by JourneyLucky
[考研] 一志愿华东理工大学081700,初试分数271 +5 kotoko_ik 2026-03-23 6/300 2026-03-24 10:29 by 学术搬砖er
[考研] 一志愿哈工大,085400,320,求调剂 +3 gdlf9999 2026-03-24 3/150 2026-03-24 10:08 by 搏击518
[考研] 291求调剂 +8 hhhhxn.. 2026-03-23 8/400 2026-03-23 23:15 by peike
[考研] 384求调剂 +3 子系博 2026-03-22 6/300 2026-03-23 21:45 by 子系博
[考研] 考研化学308分求调剂 +7 你好明天你好 2026-03-23 8/400 2026-03-23 18:39 by macy2011
[考研] 317求调剂 +12 申子申申 2026-03-19 18/900 2026-03-22 22:23 by luoyongfeng
[考研] 269专硕求调剂 +6 金恩贝 2026-03-21 6/300 2026-03-22 14:31 by ColorlessPI
[考研] 材料学硕301分求调剂 +7 Liyouyumairs 2026-03-21 7/350 2026-03-21 22:31 by peike
[考研] 0703化学调剂 +4 妮妮ninicgb 2026-03-21 4/200 2026-03-21 18:39 by 学员8dgXkO
[考研] 297求调剂 +11 戏精丹丹丹 2026-03-17 12/600 2026-03-21 17:47 by ColorlessPI
[考研] 一志愿重庆大学085700资源与环境总分308求调剂 +7 墨墨漠 2026-03-20 7/350 2026-03-21 16:36 by barlinike
[考研] 材料学学硕080502 337求调剂-一志愿华中科技大学 +4 顺顺顺mr 2026-03-18 5/250 2026-03-21 10:22 by luoyongfeng
[考研] 求调剂 +6 Mqqqqqq 2026-03-19 6/300 2026-03-21 08:04 by JourneyLucky
[考研] 330求调剂 +4 小材化本科 2026-03-18 4/200 2026-03-20 23:13 by JourneyLucky
[考研] 288求调剂 +16 于海海海海 2026-03-19 16/800 2026-03-20 22:28 by JourneyLucky
[考研] 一志愿 南京航空航天大学大学 ,080500材料科学与工程学硕 +5 @taotao 2026-03-20 5/250 2026-03-20 20:16 by JourneyLucky
信息提示
请填处理意见