Re: [理工] [工數]-高階非線性O.D.E

看板Grad-ProbAsk作者 (~口卡口卡 修~)時間14年前 (2009/12/08 19:56), 編輯推噓11(11011)
留言22則, 9人參與, 最新討論串2/3 (看更多)
※ 引述《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

12/08 20:14, , 1F
看到這題我的直覺反應是用拉式做= = 不過你這樣比較快
12/08 20:14, 1F

12/08 20:20, , 2F
m(_ _)m
12/08 20:20, 2F

12/08 20:26, , 3F
請問y'(y"+siny)=0怎麼到[(y')^2/2 - cosy]'=0這一式
12/08 20:26, 3F

12/08 20:49, , 4F
因為D大開天眼看的(逃
12/08 20:49, 4F

12/08 20:57, , 5F
看來D大已經練到每一題都用正合來解了 超神
12/08 20:57, 5F

12/08 20:59, , 6F
用自變數不出現也可以
12/08 20:59, 6F

12/08 21:03, , 7F
我前面是令p=y',pdp/dy=y",後面就學D大的算法了
12/08 21:03, 7F

12/08 21:06, , 8F
D大有輪迴眼 XD
12/08 21:06, 8F

12/08 21:06, , 9F
拉式不可能
12/08 21:06, 9F

12/08 21:07, , 10F
Lsin(y) = ?
12/08 21:07, 10F

12/08 21:10, , 11F
我是說直覺上啦 不一定可以做 呵呵
12/08 21:10, 11F

12/08 21:11, , 12F
話說今天上德文課我也忽然在腦筋中想到這題XD
12/08 21:11, 12F

12/08 21:11, , 13F
結果回來看到這篇文章@_@
12/08 21:11, 13F

12/08 21:13, , 14F
我幫你想的具體化了= =
12/08 21:13, 14F

12/08 21:14, , 15F
你們有一腿(誤)
12/08 21:14, 15F

12/08 21:15, , 16F
我覺得這個難在一階高次耶...我沒想到要半角
12/08 21:15, 16F

12/08 21:16, , 17F
看到之後就讓我想到單擺XD
12/08 21:16, 17F

12/08 21:16, , 18F
1 + cos t = 2cos(t/2)^2
12/08 21:16, 18F

12/08 21:22, , 19F
d大能提供一下其他作法嗎~
12/08 21:22, 19F

12/08 23:57, , 20F
非線性可用taylor級數解
12/08 23:57, 20F

12/09 01:35, , 21F
正規解就是上面推文說的,令 y'' = p*(dp/dy) 去解
12/09 01:35, 21F

12/09 01:38, , 22F
不過兩個好像是一樣的東西 OTZ
12/09 01:38, 22F
文章代碼(AID): #1B7Zxy3i (Grad-ProbAsk)
文章代碼(AID): #1B7Zxy3i (Grad-ProbAsk)