Re: [工數] Variation of parameters
※ 引述《zirconium (zirconium)》之銘言:
: http://ppt.cc/Da1G
: O'Neil 第六版 82 頁
: 解 Nonhomogeneous Yp 的推導過程
: 裡面如圖中紅框圈起來的
: Simplify this derivative by imposing the condition that
: u1'y1 + u2'y2 = 0
: 不太懂這條件為什麼會成立
: 請板上大大解惑 謝謝
這個方法在high-order的作法也是像這樣一路上都令為0, 最後的再令它等於f(x)
考慮線性微分方程:
(n) (n-1)
y (x) + a(n-1)y (x) + ... + a(0)y(x) = f(x)
y1(x), y2(x), ... yn(x) 是對應的 homogeneous solutions
令yp = u1(x)y1(x) + u2(x)y2(x) + ... + un(x)yn(x)
'
=> yp = (u1y1'+ u2y2'+ ... + unyn') + (u1'y1 + u2'y2 + .... + un'yn)
令 u1'y1 + u2'y2 + .... + un'yn = 0
yp'' = (u1y1''+ u2y2'' + ... + unyn'') + (u1'y1' + u2'y2' + ... + un'yn')
令 u1'y1' + u2'y2' + ... + un'yn' = 0
...
(n) (n) (n) (n)
yp = (u1y1 + u2y2 + ... + unyn ) +
(n-1) (n-1) (n-1)
(u1'y1 + u2'y2 + ... + un'yn )
(n) (n-1)
yp + a(n-1)yp + ... + a(0)yp = f(x) 比較係數
和因為 y1, y2, ... yn 是解
(n-1) (n-1) (n-1)
=> u1y1 + u2y2 + ... + unyn = f(x)
一路下來可以看到, 這樣令它為0的過程可以得到一個規則:
[y1 y2 ... yn ] [u1'] [ 0 ]
[y1'y2'... yn' ] [u2'] [ 0 ]
[ ... ] ... = ...
(n-1) (n-1)
[y1 ... yn ] [un'] [f(x)]
T
wronskian行列式對應的matrix * 參數的一階導數向量 = (0, 0, ..., f)
還蠻漂亮的
而事實上那令它為0是否成立, 就要看u1(x), u2(x), ..., un(x)能不能解的出來
--
※ 發信站: 批踢踢實業坊(ptt.cc)
推
05/15 17:01, , 1F
05/15 17:01, 1F
→
05/15 17:02, , 2F
05/15 17:02, 2F
→
05/15 17:03, , 3F
05/15 17:03, 3F
→
05/15 17:12, , 4F
05/15 17:12, 4F
→
05/15 20:05, , 5F
05/15 20:05, 5F
→
05/15 20:05, , 6F
05/15 20:05, 6F
推
05/16 08:13, , 7F
05/16 08:13, 7F
→
05/16 08:14, , 8F
05/16 08:14, 8F
→
05/16 08:16, , 9F
05/16 08:16, 9F
→
05/16 10:45, , 10F
05/16 10:45, 10F
→
05/16 10:45, , 11F
05/16 10:45, 11F
→
05/16 11:13, , 12F
05/16 11:13, 12F
※ 編輯: yueayase (36.236.229.213), 03/20/2015 21:24:38
→
11/10 11:48, , 13F
11/10 11:48, 13F
→
01/02 15:24,
7年前
, 14F
01/02 15:24, 14F
→
07/07 11:01,
6年前
, 15F
07/07 11:01, 15F
討論串 (同標題文章)
本文引述了以下文章的的內容:
完整討論串 (本文為第 2 之 2 篇):