找回密码
 欢迎注册
查看: 58128|回复: 32

[转载] 一阶非线性微分方程的求解

[复制链接]
发表于 2010-1-6 11:29:01 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?欢迎注册

×
百度贴吧里 有人给出下面的一个微分方程,我给出了解答,比较繁琐,但没人理我,我又不是科班出身,不知道正确不正确,所以贴在这,希望得到澄清~~ $y'(x)=1/(1+x^2)-2y(x)^2,y(0)=0$
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
 楼主| 发表于 2010-1-6 11:33:43 | 显示全部楼层
类似的还有一个: $y'(x)=x^2+y(x)^2,y(0)=0$
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2010-1-6 12:14:01 | 显示全部楼层
你的解法呢? 非线性的很难求解,主要没有固定的求解模式
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2010-1-6 13:06:48 | 显示全部楼层
本帖最后由 数学星空 于 2010-1-6 13:08 编辑 利用软件可以算出: $y'(x)= 1/(1+x^2)-2*y(x)^2$ 得到:$y(x) = (2* c*x+x*arctan(x)+1)/((1+x^2)*arctan(x)+2*c+2*c*x^2+x)$ 又$y(0)=0 $有$c=0$ 即$y(x) = (1+x*arctan(x))/((1+x^2)*arctan(x)+x)$ 对于第二个方程 $y'(x)= x^2+y(x)^2$ 得到:$y(x) = -x*(c*BesselJ(-3/4, (1/2)*x^2)+BesselY(-3/4, (1/2)*x^2))/(c*BesselJ(1/4, (1/2)*x^2)+BesselY(1/4, (1/2)*x^2))$
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
 楼主| 发表于 2010-1-6 13:31:52 | 显示全部楼层
本帖最后由 wayne 于 2010-1-6 13:38 编辑 这是一个Riccati方程:http://en.wikipedia.org/wiki/Riccati_equation 按照上面的链接,作代换y=p'/(2p),则原方程化成 p''=p/(2(1+x^2)). 而这个二阶齐次的我也解不了,于是我就用了Abel的一个技巧:http://mathworld.wolfram.com/AbelsDifferentialEquationIdentity.html 如果知道了二阶齐次微分方程的一个解,那么另一个独立解可以根据阿贝尔恒等式求出, 已知一个解1+x^2,设另一个解为z,那么, z'(1+x^2)-z*(1+x^2)'=C(常数) 所以$\frac{z}{1+x^2}=\int \frac{C}{(1+x^2)^2}dx$ 而后面的也折腾我半天,搞的我是用软件算出来的,最后,算出z来了, $z=c(x+(1+x^2)*arctan(x))$ 所以,$p(x)=C_1(1+x^2)+C_2(x+(1+x^2)*arctan(x))$ 然后,就是确定积分常数了,确定积分常数也很麻烦~~,我几乎是连拼带凑的才得到了最终解: $y(x)=\frac{x}{1+x^2}$
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
 楼主| 发表于 2010-1-6 13:43:23 | 显示全部楼层
本帖最后由 wayne 于 2010-1-6 13:50 编辑 4# 数学星空 谢谢~~ y(0)=0,你把x=0代入你的解,你会发现1/c =0!! 同理,而代入第二个方程,常数没了,式子恒成立, ......... 我猜想这个方程的答案是: $\frac{x BesselJ(\frac{3}{4},\frac{x^2}{2})}{BesselJ(-\frac{1}{4},\frac{x^2}{2})}$
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2010-1-6 14:05:12 | 显示全部楼层
本帖最后由 数学星空 于 2010-1-6 14:09 编辑 算出这一步$p(x)=C_1(1+x^2)+C_2(x+(1+x^2)*arctan(x))$后, 再代入$y={p'}/(2p)$ 可以算出:$y(x)=(c1*x+c2*x*arctan(x)+c2)/(c1+c1*x^2+c2*arctan(x)+c2*arctan(x)*x^2+c2*x)$ 又$y(0)=0$,得到: ${c2}/(c1+c2*arctan(0))=0, 即c1=0 或者c2=0$ 即$y(x)=(1+x*arctan(x))/((1+x^2)*arctan(x)+x)$ (1) $y(x)=x/(1+x^2) $ (2) 将解(1),(2)代入均成立,即有两个解
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
 楼主| 发表于 2010-1-6 14:12:29 | 显示全部楼层
本帖最后由 wayne 于 2010-1-6 14:35 编辑 我想到了一个问题 两个可微的函数有共同的连续区间U,如果存在区间$I\subU$,它们在公共区间I内处处相等,是不是就意味着这两个函数在U内也恒等。 唉,没学数学分析是我的人生一大遗憾,i wonder if i made myself understood,
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
 楼主| 发表于 2010-1-6 14:16:22 | 显示全部楼层
7# 数学星空 好像不满足y(0)=0
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
 楼主| 发表于 2010-1-6 14:29:04 | 显示全部楼层
4# 数学星空 发现你在4#初值代错了,所以产生了这个不正确的解,~~
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
您需要登录后才可以回帖 登录 | 欢迎注册

本版积分规则

小黑屋|手机版|数学研发网 ( 苏ICP备07505100号 )

GMT+8, 2024-11-22 10:21 , Processed in 0.028573 second(s), 17 queries .

Powered by Discuz! X3.5

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表