Re: [工數] 一題ODE

看板Math作者 (趙哥)時間14年前 (2011/08/20 18:06), 編輯推噓5(5013)
留言18則, 3人參與, 最新討論串3/3 (看更多)
※ 引述《handsboy (夠夠拋兒瑞久)》之銘言: : dy 3 2 -1 : ---- = x (y-x) + x y , x>0 : dx : 想不出來-.- : 有人會解嗎 此為Riccati ODE 按照解題步驟慢慢算即可.. 1. y=x為一特解 令y=x+u 帶入ODE (會變成Bernoulli ODE) 帶入得 1+u'=x^3*u^2 + (x+u)/x u'=x^3*u^2 + u/x (有合併對象 用觀察法較易..) u^-2 du = (1/x)(1/u)dx + x^3dx d(1/u) + (1/x)(1/u)dx = - x^3dx , I= exp(lnx)=x d(x/u) = -x^4dx , xu = (-x^5+c)/5 u = 5x/(c-x^5) , y = x+u = x + 5x/(c-x^5) ===============另解==================== 令y-x=z 整理得 z'=x^3*z^2 + z/x xdz-zdx = x^4*z^2 dx d(zx^-1)/(x^-2) = x^4*z^2 dx 1/(x^-2*z^2) = x^4 dx ,積分得 -1/z*x^-1 = (x^5+c)/5 z = -5x/(x^5+c) , y = x - 5x/(x^5+c) -- ※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 140.113.122.118

08/20 18:11, , 1F
哦哦哦哦哦 原來是RICCATI -.-
08/20 18:11, 1F

08/20 18:12, , 2F
不過我不太懂另一個大大怎樣看出變異係數法可以解的
08/20 18:12, 2F

08/20 18:12, , 3F
@@
08/20 18:12, 3F

08/20 18:14, , 4F
打錯 參數變異法 (好像一樣? XD
08/20 18:14, 4F

08/20 18:17, , 5F
參數變易法是哪個阿 哈哈
08/20 18:17, 5F

08/20 18:19, , 6F
沒事 我想錯了@@
08/20 18:19, 6F

08/20 18:20, , 7F
我剛看了一下 用z=y-x後 觀察法應該也可以解
08/20 18:20, 7F
※ 編輯: doublewhi 來自: 140.113.122.118 (08/20 18:27)

08/20 18:23, , 8F
太久沒看一階ODE 就..變的廢廢的 -.-
08/20 18:23, 8F

08/20 18:23, , 9F
令y=x+ 1/u 帶回好像比較快 直接解一階線性
08/20 18:23, 9F

08/20 18:28, , 10F
可以問如何看出y=x+ 1/u嗎 @@
08/20 18:28, 10F
※ 編輯: doublewhi 來自: 140.113.122.118 (08/20 18:28)

08/20 19:21, , 11F
就把你的u 令成我的1/u 哈哈 算出來應該一樣
08/20 19:21, 11F

08/20 22:21, , 12F
u直接改另1/u算是直接跳兩步 跟教科書的算法相同
08/20 22:21, 12F

08/20 22:23, , 13F
可化成參數變異法則是將y另為-(1/R(x)u)(du/dx)
08/20 22:23, 13F

08/20 22:24, , 14F
就可以升階成二階ODE
08/20 22:24, 14F

08/20 23:55, , 15F
樓上可以詳細介紹下升階的做法和想法嗎?
08/20 23:55, 15F

08/21 14:09, , 16F
升階我是在書上看到的 不過我推倒不出來所以沒用過這
08/21 14:09, 16F

08/21 14:10, , 17F
種解法 這種解法優點是可以不用假設特解救可求通解
08/21 14:10, 17F

08/21 14:11, , 18F
化減後會變u"-(Q+R'/R)u'+PRu=0
08/21 14:11, 18F
文章代碼(AID): #1EJuSTB- (Math)
討論串 (同標題文章)
本文引述了以下文章的的內容:
工數
2
6
完整討論串 (本文為第 3 之 3 篇):
工數
2
6
文章代碼(AID): #1EJuSTB- (Math)