应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 专业学科区 » 数学 » 工科数学 » 【求助】求解形如y''+a*y'^2+by'+c=0 二阶微分方程
151/2返回列表12下一页
查看: 2472  |  回复: 14
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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大小判断的问题，不知能否求出适当的解析解，还请各位帮助
回复此楼

» 猜你喜欢
1楼 2010-08-31 10:25:22
已阅   回复此楼   关注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人) 回复此楼
明德厚学、求是创新
2楼2010-08-31 20:04:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fyq98

木虫 (正式写手)
★
xu17211(金币+1):谢谢参与
这类问题的解析解可以通过幂级数法求得，即假设解为幂级数，代入方程，根据对应次幂系数相等，得到一些代数方程，解之即得
赞一下(2人) 回复此楼
克难奋进
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:
这类问题的解析解可以通过幂级数法求得，即假设解为幂级数，代入方程，根据对应次幂系数相等，得到一些代数方程，解之即得

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

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

实际应用数值计算居多
赞一下(2人) 回复此楼
长期招收博士生，参见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时: ......
赞一下(1人) 回复此楼
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):谢谢参与
看起简单，真是挺难。一定要求解么？不能用常微分方程的定性理论的一些东西来分析解的性质等等？
赞一下(1人) 回复此楼
拳打脚踢
10楼2010-09-02 11:50:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xu17211 的主题更新
151/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭