作者doom8199 (～口卡口卡 修～)
看板Grad-ProbAsk
标题Re: [理工] [工数]-高阶非线性O.D.E
时间Tue Dec 8 19:56:41 2009
※ 引述《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
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 140.113.141.151
1F：→ boy210637:看到这题我的直觉反应是用拉式做= = 不过你这样比较快 12/08 20:14
2F：推 iyenn: m(_ _)m 12/08 20:20
3F：推 zendla:请问y'(y"+siny)=0怎麽到[(y')^2/2 - cosy]'=0这一式 12/08 20:26
4F：推 iyenn: 因为D大开天眼看的(逃 12/08 20:49
5F：推 shinyhaung:看来D大已经练到每一题都用正合来解了 超神 12/08 20:57
6F：→ birdhackor:用自变数不出现也可以 12/08 20:59
7F：→ zendla:我前面是令p=y'，pdp/dy=y"，後面就学D大的算法了 12/08 21:03
8F：推 boy210637:D大有轮回眼 XD 12/08 21:06
9F：推 ntust661:拉式不可能 12/08 21:06
10F：→ ntust661:Lsin(y) = ? 12/08 21:07
11F：→ boy210637:我是说直觉上啦 不一定可以做 呵呵 12/08 21:10
12F：推 ntust661:话说今天上德文课我也忽然在脑筋中想到这题XD 12/08 21:11
13F：→ ntust661:结果回来看到这篇文章@_@ 12/08 21:11
14F：→ zendla:我帮你想的具体化了= = 12/08 21:13
15F：→ boy210637:你们有一腿(误) 12/08 21:14
16F：推 birdhackor:我觉得这个难在一阶高次耶...我没想到要半角 12/08 21:15
17F：推 ntust661:看到之後就让我想到单摆XD 12/08 21:16
18F：→ ntust661:1 + cos t = 2cos(t/2)^2 12/08 21:16
19F：推 kagato:d大能提供一下其他作法吗~ 12/08 21:22
20F：推 pigmow:非线性可用taylor级数解 12/08 23:57
21F：→ doom8199:正规解就是上面推文说的，令 y'' = p*(dp/dy) 去解 12/09 01:35
22F：→ doom8199:不过两个好像是一样的东西 OTZ 12/09 01:38