24小时热门版块排行榜    

查看: 2551  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085700资源与环境308求调剂 +3 墨墨漠 2026-03-18 3/150 2026-03-18 22:35 by bingxueer79
[考研] 一志愿吉林大学材料学硕321求调剂 +4 Ymlll 2026-03-18 6/300 2026-03-18 22:15 by li123456789.
[考研] 321求调剂 +3 何润采123 2026-03-18 3/150 2026-03-18 21:27 by li123456789.
[考研] 295求调剂 +3 一志愿京区211 2026-03-18 5/250 2026-03-18 17:03 by zhaoqian0518
[考研] 08工科 320总分 求调剂 +5 梨花珞晚风 2026-03-17 5/250 2026-03-18 14:49 by haxia
[考研] 0703化学求调剂 总分331 +3 ZY-05 2026-03-13 3/150 2026-03-18 10:58 by macy2011
[考研] 0703化学调剂 +3 妮妮ninicgb 2026-03-17 3/150 2026-03-18 10:29 by macy2011
[考研] 085600材料与化工求调剂 +5 绪幸与子 2026-03-17 5/250 2026-03-17 16:40 by laoshidan
[考研] 26考研求调剂 +6 丶宏Sir 2026-03-13 6/300 2026-03-17 16:13 by 醉在风里
[考研] 278求调剂 +3 Yy7400 2026-03-13 3/150 2026-03-17 08:24 by laoshidan
[考研] 11408 一志愿西电,277分求调剂 +3 zhouzhen654 2026-03-16 3/150 2026-03-17 07:03 by laoshidan
[考研] 070300化学学硕求调剂 +6 太想进步了0608 2026-03-16 6/300 2026-03-16 16:13 by kykm678
[考博] 东华理工大学化材专业26届硕士博士申请 +6 zlingli 2026-03-13 6/300 2026-03-15 20:00 by ryzcf
[考研] 070305求调剂 +3 mlpqaz03 2026-03-14 4/200 2026-03-15 11:04 by peike
[考研] 289求调剂 +4 这么名字咋样 2026-03-14 6/300 2026-03-14 18:58 by userper
[考研] 招收0805(材料)调剂 +3 18595523086 2026-03-13 3/150 2026-03-14 00:33 by 123%、
[考研] 26调剂/材料科学与工程/总分295/求收留 +9 2026调剂侠 2026-03-12 9/450 2026-03-13 20:46 by 18595523086
[考研] 295求调剂 +3 小匕仔汁 2026-03-12 3/150 2026-03-13 15:17 by vgtyfty
[考研] 290求调剂 +7 ADT 2026-03-12 7/350 2026-03-13 15:17 by JourneyLucky
[考研] 274求调剂 +3 S.H1 2026-03-12 3/150 2026-03-13 15:15 by JourneyLucky
信息提示
请填处理意见