1 / 72

# 4-4 Undetermined Coefficients – Superposition Approach - PowerPoint PPT Presentation

4-4 Undetermined Coefficients – Superposition Approach. This section introduces some method of “ guessing ” the particular solution. 4-4-1 方法適用條件. (1). (2). Suitable for linear and constant coefficient DE.

I am the owner, or an agent authorized to act on behalf of the owner, of the copyrighted work described.

## PowerPoint Slideshow about ' 4-4 Undetermined Coefficients – Superposition Approach' - ezra-perry

Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author.While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server.

- - - - - - - - - - - - - - - - - - - - - - - - - - E N D - - - - - - - - - - - - - - - - - - - - - - - - - -
Presentation Transcript

This section introduces some method of “guessing” the particular solution.

4-4-1 方法適用條件

(1)

(2)

Suitable for linear and constant coefficient DE.

(3) g(x), g'(x), g'' (x), g'''(x), g(4)(x), g(5)(x), ………contain finite number of terms.

4-4-2 方法

g(x) 長什麼樣子，particular solution 就應該是什麼樣子.

(計算時要把 A, B, C, … 這些 unknowns 解出來)

(from text page 143)

It comes from the “form rule”. See page 198.

yp = ?

yp = ?

yp = ?

Example 2 (text page 141)

Step 1: find the solution of the associated homogeneous equation

Guess

Step 2: particular solution

A = 6/73, B = 16/73

Step 3: General solution:

Example 3 (text page 142)

Step 1: Find the solution of

Step 2: Particular solution

guess

guess

Step 3: General solution

4-4-4 方法的解釋

Form Rule:yp should be a linear combination of g(x), g'(x),

g'' (x), g'''(x), g(4)(x), g(5)(x), …………….

Why? 如此一來，在比較係數時才不會出現多餘的項

When g(x) = xn

When g(x) = cos kx

When g(x) = exp(kx)

When g(x) = xnexp(kx)

:

:

Example 4

(text page 142)

Particular solution guessed by Form Rule:

(no solution)

Why?

Glitch condition 1: The particular solution we guess belongs to the complementary function.

For Example 4

Complementary function

(text page 144)

From Form Rule, the particular solution is Aex

Example 8 (text page 145)

Step 1

Step 2

Step 3

Step 4 Solving c1 and c2 by initial conditions

(最後才解 IVP)

Example 11 (text page 146)

From Form Rule

yp只要有一部分和 yc相同就作修正

If we choose

If we choose

A = 1/6, B = 1/3, C = 3, E = 12

Glitch condition 2:g(x), g'(x), g'' (x), g'''(x), g(4)(x), g(5)(x), ……………

contain infinite number of terms.

If g(x) = ln x

If g(x) = exp(x2)

:

:

4-4-6 本節需要注意的地方

(1) 記住 Table 4.1 的 particular solution 的假設方法(其實和 “form rule” 有相密切的關聯)

(2) 注意 “glitch condition”

(3) 所以要先算 complementary function，再算 particular solution

(4) 同樣的方法，也可以用在 1st order 的情形

(5) 本方法只適用於 linear, constant coefficient DE

For a linear DE:

Annihilator Operator:

4-5-1 方法適用條件

(1) Linear， (2) Constant coefficients

(3) g(x), g'(x), g'' (x), g'''(x), g(4)(x), g(5)(x), ………contain finite number of terms.

Example 1: (text page 150)

annihilator: D4

annihilator: D + 3

annihilator: (D − 2)2

(D − 2)2 = D2 − 4D + 4

(x − 2)2 = x2 − 4x + 4

 (D − 2)2 = D2 − 4D + 4

If

then the annihilator is

If

b1 0 or b2 0

then the annihilator is

Example 2: (text page 151)

annihilator

Example 5: (text page 154)

annihilator

Example 6: (text page 155)

annihilator

If g(x) = g1(x) + g2(x) + …… + gk(x)

Lh[gh(x)] = 0 but Lh[gm(x)]  0 if m h,

then the annihilator of g(x) is the product of Lh (h = 1 ~ k)

Proof:

(因為 L1, L2為 linear DE with constant coefficient,

L1L2 = L2L1 )

:

:

Therefore,

Example 7 (text page 154)

annihilator: D3

annihilator: D − 5

annihilator: (D− 2)3

annihilator of g(x): D3 (D− 2)3 (D − 5)

Step 2-1

Find the annihilator L1 of g(x)

Step 2-2

(homogeneouslinear & constant coefficient DE)

then

Use the method in Section 4-3 to find the solution of

Step 2-4

Find the particular solution.

The particular solution yp is a solution of

but not a solution of

(Proof):

Since , if g(x)  0, should be nonzero.

Moreover, .

Step 2-5

Solve the unknowns

particular solution yp

solutions of

particular solution yp solutions of

 solutions of

Example 3 (text page 152)

Step 1: Complementary function

(solution of the associated homogeneous function)

Step 2-1: Annihilation: D3

Step 2-2:

Step 2-3:

auxiliary function

roots: m1 = m2 = m3 = 0, m4 = −1, m5 = −2

Solution for :

Step 2-4 : particular solution

Step 2-5:

Step 3:

Example 4 (text page 153)

Step 1: Complementary function

From auxiliary function, m2− 3m = 0, roots: 0, 3

Step 2-1: Find the annihilator

D− 3 annihilate but cannot annihilate

(D2 + 1) annihilate but cannot annihilate

(D − 3)(D2 + 1) is the annihilator of

Step 2-2:

auxiliary function:

solution of :

Step 2-4: particular solution

Step 2-5:

Step 3: general solution

4-5-5 本節要注意的地方

(1) 所以要先算 complementary function，再算 particular solution

(2) 若有兩個以上的 annihilator，選其中較簡單的即可

(3) 計算 auxiliary function 時有時容易犯錯

(4) 的解和 的解不一樣。

(5) 這方法，只適用於 constant coefficient linear DE

(因為，還需借助 auxiliary function)

The thing that can be done by the annihilator approach can always be done by the “guessing” method in Section 4-4, too.

4-6-1 方法的限制

The method can solve the particular solution for any linear DE

(1) May not have constant coefficients

(2) g(x) may not be of the special forms

4-6-2 Case of the 2 nd order linear DE

associated homogeneous equation:

Suppose that the solution of the associated homogeneous equation is

Then the particular solution is assumed as:

(方法的基本精神)

zero

zero

| |: determinant

4-6-3 Process for the 2 nd Order Case

Step 2-1 變成standard form

Step 2-2

Step 2-3

Step 2-4

Step 2-5

Example 1 (text page 159)

Step 1: solution of

Step 2-2:

Step 2-3:

Step 2-5:

Step 3:

Example 2 (text page 159)

Step 1: solution of

Step 2-1: standard form:

Step 2-2:

Step 2-3:

Step 2-4:

(未完待續)

Step 3:

Note: 課本 Interval (0, /6) 改為(0, /3)

Example 3 (text page 160)

Note: 沒有 analytic 的解

Solution of the associated homogeneous equation:

The particular solution is assumed as:

W k: replace the kth column of W by

For example, when n = 3,

Step 2-1 變成 standard form

Step 2-2 Calculate W, W1, W2, …., Wn (see page 237)

Step 2-3 ………

Step 2-4 …….

Step 2-5

4-6-7 本節需注意的地方

• 養成先解 associated homogeneous equation 的習慣

• 記熟幾個重要公式

• 這裡 | | 指的是 determinant

• (4) 算出 u1(x) 和 u2(x) 後別忘了作積分

• (5) f(x) = g(x)/an(x) (和 1st order 的情形一樣，使用 standard form)

• (6) 計算 u1'(x) 和 u2'(x) 的積分時，+ c可忽略

• 因為 我們的目的是算particular solution yp

• yp是任何一個能滿足原式的解

• (7) 這方法解的範圍，不包含 an(x) = 0 的地方

4-7-1 解法限制條件

k

not constant coefficients

but the coefficients of y(k)(x) have the form of

ak is some constant

associated homogeneous equation

particular solution

4-7-2 解法

Associated homogeneous equation of the Cauchy-Euler equation

Guess the solution as y(x) = xm , then

4-7-3 For the 2nd Order Case

auxiliary function:

roots

[Case 1]: m1 m2 and m1, m2 are real

two independent solution of the homogeneous part:

[Case 2]: m1= m2

Use the method of reduction of order

Note 1: 原式

Note 2: 此時

If y2(x) is a solution of a homogeneous DE

then c y2(x) is also a solution of the homogeneous DE

If we constrain that x > 0, then

[Case 3]: m1 m2 and m1, m2 are the form of

two independent solution of the homogeneous part:

Example 1 (text page 163)

Example 2 (text page 164)

Example 3 (text page 165)

4-7-4 For the Higher Order Case

Process:

auxiliary function

Step 1-1

roots

n independent solutions

Step 1-2

solution of the nth order associated homogeneous equation

Step 1-3

(1) 若 auxiliary function在 m0的地方只有一個根

(2) 若 auxiliary function在 m0的地方有 k個重根

(3) 若 auxiliary function在 + j和 −j的地方各有一個根

(未出現重根)

(4) 若 auxiliary function在 + j和 −j的地方皆有 k個重根

Example 4 (text page 166)

auxiliary function

4-7-5 Nonhomogeneous Case

To solve the nonhomogeneous Cauchy-Euler equation:

Method 1: (See Example 5)

(1) Find the complementary function (general solutions of the associated homogeneous equation) from the rules on pages 244-247, 251-252.

(2) Use the method in Sec.4-6 (Variation of Parameters) to find the particular solution.

(3) Solution = complementary function + particular solution

Method 2: See Example 6，重要

Set x = et, t = ln x

Example 5 (text page 166, illustration for method 1)

Step 1 solution of the associated homogeneous equation

auxiliary function

Step 2-2 Particular solution

Step 2-3

Step 2-5

Step 3

Example 6 (text page 167, illustration for method 2)

Set x = et, t = ln x

(chain rule)

Therefore, the original equation is changed into

( 別忘了 t = ln x要代回來)

Note 1: 以此類推

Note 2: 簡化計算的小技巧：結合兩種解 nonhomogeneous Cauchy-Euler equation 的長處

4-7-6 本節要注意的地方

(1) 本節公式記憶的方法：

(2) 如何解 particular solution?

Variation of Parameters 的方法

(3) 解的範圍將不包括 x = 0 的地方 (Why?)

(1) numerical approach (Section 4-9-3)

(2) using special function (Chap. 6)

(3) Laplace transform and Fourier transform (Chaps. 7, 11, 14)

(4) 查表 (table lookup)

(1) 即使用了 Section 4-7 的方法，大部分的 DE還是沒有辦法解

(2) 所幸，自然界真的有不少的例子是 linear DE

Section 4-4 5, 6, 14, 17, 18, 24, 26, 32, 33, 39, 42

Section 4-5 2, 7, 13, 18, 31, 45, 62, 69, 70

Section 4-6 4, 5, 8, 13, 14, 17, 18, 21, 24, 25, 30

Section 4-7 11, 17, 18, 20, 21, 24, 32, 34, 36, 37, 42

Review 4 2, 19, 20, 23, 25, 26, 27, 28, 30, 31, 33, 38