作者fonlintw0621 (fonlintw0621)
看板Grad-ProbAsk
标题Re: [理工] [工数]-高阶非线性O.D.E
时间Tue Dec 8 22:32:47 2009
※ 引述《doom8199 (～口卡口卡 修～)》之铭言：
: ※ 引述《zendla (夏夜薄荷)》之铭言：
: : 出现三角的题目，感觉不是很好做啊，请大大们教一下
: : Find the solution of the nonlinear ordinary differential equation
: : y" + siny = 0
: : with initial conditions: y(0) = 0, y'(0) = 2. [清大工科]
: : y y
: : ans:ln|sec── + tan──| = x
: : 2 2
: ---
: 感觉这题作法很多样化
: 我是这样想:
: y" + siny = 0
: → y'(y'' + siny) = 0 for y'≠0
: (y')^2
: → [ ______ - cosy ]' = 0
: 2
: → (y')^2 - 2cosy = c1
: 把 x=0 带入:
: 2^2 - 2cos0 = c1 → c1 = 2
: 所以 (y')^2 = 2cosy + 2
: → y' = ± 2cos(y/2)
: → ∫ (1/2)sec(y/2) dy = ∫ ±1 dx
: → ln｜sec(y/2) + tan(y/2)｜ = ±x + c2
: 再由 y(0)=0 可解出 c2=0
: y'(0)=2 可知 ln｜sec(y/2) + tan(y/2)｜ = -x 不合
: 因此 ln｜sec(y/2) + tan(y/2)｜ = x
y(0) = 0
y'(0) = 2
y'' = -sin y = 0
y''' = -cos y * y' = -2
2
y'''' = sin y * y' - cos y * y'' = 0
3
y''''' = cos y * y' + sin y * 2 y' y'' + siny y' * y'' - cosy * y'''
= 10
看到题目不会作
直接泰勒展开 送他
y'(0)
y (x) = y(0) + -------- x .........................
1 !
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 114.36.212.94
1F：推 hihaka2001:这题好像在级数解有看到ㄟXD说真的变线性真技巧有点难 12/08 23:40
2F：推 pigmow:刚刚没看到这篇 我直觉也是用这招 非线性本来就没固定方法 12/09 00:01
3F：→ birdhackor:如果可以求出一般通解 用级数展开应该会扣点分吧 12/09 00:02
4F：→ birdhackor:我猜的...不保证... 12/09 00:03