Section 7-4 Operational Properties II

Slides:



Advertisements
Similar presentations
1 Lecture 5 Properties of LTI Systems The Solution of LCCDE.
Advertisements

Differentiation 微分 之二 以公式法求函數的微分. Type 函數形式 Function f (x) Derivative d f (x) /d x c=constant 常數 c0 Power of x xaxa a x a-1 Trigonometric 三角函數 sin x cos.
663 Chapter 14 Integral Transform Method Integral transform 可以表示成如下的積分式的 transform  kernel Laplace transform is one of the integral transform 本章討論的 integral.
楊學成 老師 Chapter 1 First-order Differential Equation.
工職數學 第四冊 第一章 導 數 1 - 1 函數的極限與連續 1 - 2 導數及其基本性質 1 - 3 微分公式 1 - 4 高階導函數.
數位訊號處理 第4章 離散時間訊號與LTI系統之傅利葉分析

5.1 自然對數函數:微分 5.2 自然對數函數:積分 5.3 反函數 5.4 指數函數:微分與積分 5.5 一般底數的指數函數和應用 5.6 反三角函數:微分 5.7 反三角函數:積分 5.8 雙曲函數.
Differential Equations (DE)
XI. Hilbert Huang Transform (HHT)
工程數學 Chapter 12 PDE 楊學成 老師.
Chapter three the Z Transform Z 变换
3-3 Modeling with Systems of DEs
Euler’s method of construction of the Exponential function
Signals and Systems Lecture 28
IV. Implementation IV-A Method 1: Direct Implementation 以 STFT 為例
期末考的範圍遠遠多於期中考,要了解的定理和觀念也非常多
Linear Programming: Introduction and Duality
V. Homomorphic Signal Processing
數 學 基 礎 2 ※ 本章主要目的 1. 介紹拉氏轉換的基本理論。 2. 舉例說明應用拉氏轉換求解線性常微分方程式的方法。
課程大綱 第一章 Laplace 變換 1.1 基本概念與定理 1.2 常係數之線性微分方程式的 Laplace 變換解
Chapter 7 The Laplace Transform
附錄一 Methods of Solving the First Order Differential Equation
Section 7-2 Inverse Transforms and Transforms of Derivatives
Differential Equations (DE)
Differentiation 微分 之二 以公式法求函數的微分.
Chapter 12 Boundary-Value Problem in Rectangular Coordinates
Chapter 11 Orthogonal Functions and
On Some Fuzzy Optimization Problems
期末考的範圍遠遠多於期中考,要了解的定理和觀念也非常多
期末考的範圍遠遠多於期中考,要了解的定理和觀念也非常多
Ch.2 Modeling in the Frequency Domain
1.1 線性方程式系統簡介 1.2 高斯消去法與高斯-喬登消去法 1.3 線性方程式系統的應用
Chapter 14 Integral Transform Method
第二十九單元 方向導數與梯度.
II. Short-time Fourier Transform
Randomized Algorithms
4-4 Undetermined Coefficients – Superposition Approach
VI. Brief Introduction for Acoustics
4-5 Undetermined Coefficients – Annihilator Approach
第 一 單 元 不定積分.
2-4 Exact Equations 方法的條件 任何 first order DE 皆可改寫成 的型態 (1) 當 成立時,
Chapter 5 Modeling with Higher Order Differential Equations
Chapter 2 聯立線性方程式與矩陣 授課教師:李金鳳(Amy Lee)
义务教育课程标准实验教科书 九年级 上册 28.1 锐角三角函数(第4课时) 人民教育出版社.
Advanced Digital Signal Processing 高等數位訊號處理
第三章 付里叶分析 离散付氏级数的数学解释(The Mathematical Explanation of DFS)
XIV. Orthogonal Transform and Multiplexing
VII. Data Compression (A)
Summary for Chapters 24 摘要: 24章
Definition of Trace Function
4-7 Cauchy-Euler Equation
96學年度第一學期電機系教學助理課後輔導進度表(一)
工程數學 Chapter 10 Fourier Series , Integrals , and Transforms 楊學成 老師.
weihuang[AT]mail.ustc.EDU.cn Summer 2018, Hefei math/phys
4-1 Linear Differential Equations: Basic Theory
第4章 连续时间傅立叶变换 The Continuous-Time Fourier Transform
Q & A.
F F F F F F F 第二章 连续时间信号与系统的时域分析 本章要点 常用典型信号 连续时间信号的分解 连续时间系统的数学模型
96學年度第二學期電機系教學助理課後輔導進度表(二)(查堂重點)
96學年度第二學期電機系教學助理課後輔導進度表(三)(查堂重點)
IV. Implementation IV-A Method 1: Direct Implementation 以 STFT 為例
96學年度第二學期電機系教學助理課後輔導進度表(一)(查堂重點)
本講義為使用「訊號與系統,王小川編寫,全華圖書公司出版」之輔助教材
補充 數值方法 數值方法.
II. Short-time Fourier Transform
Principle and application of optical information technology
Computer Architecture
Presentation transcript:

Section 7-4 Operational Properties II 7-4-1 Derivatives of Transforms 比較: Laplace 微分 乘 sn Laplace 乘 tn 微分

Proof of the Theorem of Derivatives of Transforms:

Example 1 (text page 283) 練習:為何

7-4-2 Convolution (旋積) Definition of convolution: (標準定義) When f(t) = 0 for t < 0 and g(t) = 0 for t < 0 , 上方的式子可以簡化為下方的式子

Convolution 的物理意義 (重要) 當 Input f() 對 output y(t) 的影響為 g(t−) g(t−) 只和 t 與  之間的差有關 Input f() 對 output y(t) 的影響,決定於 t 與  之間的差 例如: f() 是在  這個時間點上太陽照射到某個地方的熱量 g(t-) 可想像成是經過了t- 的時間之後,還未幅射回外太 熱量比例 y(t) 可想像成是溫度

7-4-3 Convolution Theorem Multiplication Convolution Proof: note (A) 見後頁說明 令 t =  +  note (B ) 見後頁說明

note (A) 定理: note (B) 積分範圍的改變:

Example 3 (text page 284) Example 4 (text page 285)

7-4-4 Integration (想成 “負一次微分”)

Example: Example:

7-4-5 Transform of a Periodic Function Theorem 7.4.3 When f(t + T) = f(t) then Proof: 令 f1(t) = f (t) when 0  t < T f1(t) = 0 otherwise f (t) = f1(t) + f1(t − T) + f1(t − 2T) + f1(t − 3T) + ……………. = f1(t) + f1(t −T)u(t − T) + f1(t −2T) u(t −2T ) + f1(t −3T) u(t −3T ) + …………….

L{f (t)} = L{f1(t)} + L{f1(t −T)u(t − T)} +L{ f1(t − 2T) u(t − 2T)} :

Example 7 (text page 288) Square Wave (方波) 的例子 for 0  t < 1 for 1  t < 2

Example 8 (text page 288) E(t) 為 page 482 之方波

長除法 k = 0, 1, 2, 3, …..

先使用 的公式 再算出 註:雖然也可以用 來算 但是較麻煩且容易出錯

7-4-6 Section 7.4 要注意的地方 (1) 注意代公式的順序 (例:Page 487 例子) (2) 熟悉 convolution (3) 變成積分時,別忘了加上 initial value 如 (4) 一定要背熟幾個重要的 properties (7 大性質)

Section 7.5 The Dirac Delta Function 7-5-1 Unit Impulse for t < t0−a or t > t0+a for t0−a  t  t0+a 稱作 unit impulse 高= 1/2a 面積 = 1 t0−a t0+a t-axis

7-5-2 Dirac Delta Function for t = t0 for t  t0 Fig. 7-5-2

7-5-3 Properties of the Dirac Delta Function (1) Integration (2) Sifting when t0  [p, q] Proof: 當 a 很小的時候,f(t) = f(t0) for t0−a  t  t0+a

(3) Laplace transform of (t − t0) when t0 > 0 Proof: (4) Relation with the unit step function

7-5-4 Example Example 1(a) (text page 293) 0  t < 2 t  2

7-5-5 幾個名詞 where (1) w(t) = L−1{W(s)} 稱作 weight function 或 impulse response Note: When Q(s) = 0 (no initial condition) and G(s) = 1 (g(t) = (t)), Y(s) = W(s), y(t) = w(t). (2) 許多文獻把 Dirac delta function (t − t0) 亦稱作 delta function , impulse function,或 unit impulse function

7-5-6 本節要注意的地方 (1) Dirac delta function 不滿足 Theorem 7.1.3 7-5-6 本節要注意的地方 (1) Dirac delta function 不滿足 Theorem 7.1.3 (2) 幾個定理記熟,本節即可應付自如

Section 7-6 Systems of Linear Differential Equations Chapter 7 的應用題 比較:類似的問題,也曾經在 Section 4-8 出現過 7-6-1 雙彈簧的例子

7-6-2 電路學的例子 (由第2, 3 個式子) Fig. 3-3-4 Fig. 7-6-2

Example 2 (text page 297) E(t) = 60 V, L = 1 h, R = 50 , C = 10−4 f, i1(t) = i2(t) = 0 ……….. (式1) …... (式2) (式1) × 1 + (式2) × s

複習:分子如何算出? 將 代入式 (1)

7-6-3 Double Pendulum (雙單擺) 的例子

Example 3 (text page 298) m1 = 3, m2 = 1, l1 = 12 = 16, 1(0) = 1, 2(0) = 1, '1(0) = 0, '2(0) = 0 Laplace

………….. (式1) ………….. (式2) (式1) × (16s2 + g)  (式2) × 4s2

將 代入 (式2) 直接用之前的式子

分式分解快速驗算技巧 將 s = 0, s = 1, 或其他的值代入,看等號是否成立

7-6-4 本節需要注意的地方 (1) 正負號勿寫錯 (2) 要熟悉聯立方程式的變數消去法 (3) 多學習,甚至多「研發」簡化計算的技巧

Review of Chapter 7 (1) Laplace transform 定義 Inverse Laplace transform If and f(t) is piecewise continuous of exponential order then (2) 7 大transform pairs (看講義 page 420)

transform pairs 補充 f(t) F(s) tsin(kt) tcos(kt) tsinh(kt) tcosh(kt) u(ta) f (t) = f (t+2a) 1 a b

(3) 7 大 properties input Laplace transform (1) Differentiation (Sec 7-2) (2) Multiplication by t (Sec7-4) (3) Integration (Sec 7-4) (續)

input Laplace transform (4) Multiplication by exp (Sec7-3) (5.1) Translation (Sec 7-3) (5.2) Translation (Sec 7-3) (6) Convolution (Sec 7-4) (7) Periodic Input (Sec 7-4) f(t) = f(t + T)

Properties 補充 input Laplace transform Scaling f(t / a) aF(as) Multiple Integrations Integration for s f(t) / t

(4) 簡化運算的方法 分式分解 (see pages 439-444) Initial conditions (see pages 452, 453) (5) Delta function 的四大性質 Pages 491, 492

(6) General solutions Laplace transform 的 general solution,可以用 initial conditions 來表示。 例子: 用 Laplace transform :

和 Section 4-3 的解相比較 將 代入

Exercise for practice Sec. 7-1: 5, 8, 9, 18, 32, 33, 38, 41, 46, 48 Sec. 7-2:   11, 20, 23, 27, 30, 36, 37, 41, 42, 43 Sec. 7-3:   10, 16, 19, 20, 24, 34, 44, 56, 62, 70, 74, 83 Sec. 7-4:   8, 13, 25, 28, 30, 43, 52, 53, 54, 59, 61 Sec. 7-5:   5, 8, 11, 12, 15 Sec. 7-6:   8, 11, 12, 14, 15 Review 7:  12, 24, 25, 29, 36, 37, 40, 41, 42 Maclaurin seris = Taylor series Homework 3 (due: 12/2) (1) Sec. 6-1 31 (2) Sec. 6-1 34 (3) Sec. 6-2 19 (4) Sec. 6-2 32 (5) Sec. 6-3 6 (6) Review 6 6 (7) Sec. 7-1 7 (8) Sec. 7-1 42 (9) Sec. 7-2 24 (10) Sec. 7-3 46

Chapter 8 Systems of Linear First-Order Differential Equations 另一種解「聯立微分方程式」的方法 (1) Section 4.8: (2) Chapter 7: (3) Chapter 8: Using matrix operations 註:本章這學期只教不考

比較 (1) 這 3 種方法都只適用於 linear & constant coefficients 的情形 註:其實 Laplace transform 可用來解 nonlinear & non-constant coefficient DEs, 但過程頗為複雜 (2) Laplace transform 的方法優於 Section 4-8 的方法的地方, 在於可以輕易的解決 initial condition 的問題 注意:但是,若 boundary conditions 不是在 t = 0 的地方,用 Laplace transform 需要花一番功夫。

(3) 無論是 Section 4-8 的方法,還是 Laplace transform, 運算量皆不少 Chapter 8 的方法可以減少 1st order 聯立微分方程式的運算量 但 2nd order 以上反而比 Laplace transform 麻煩

Section 8.1 Preliminary Theory 方法的限制: (a) linear, (b) 1st order DEs (c) full rank (n 個 dependent variable 需要 n 個式子) 名詞: linear system (pp. 521) homogeneous, nonhomogeneous (pp. 522) solution vector (pp. 522) fundamental set of solutions (pp. 526) general solution (pp. 526) complementary function (pp. 530) particular solution (pp. 530) 本節學習秘訣:和 Section 4-1 相比較

8-1-1 表示法和名詞 假設有 n 個 dependent variables x1(t), x2(t), ……., xn(t), 8-1-1 表示法和名詞 假設有 n 個 dependent variables x1(t), x2(t), ……., xn(t), n 個只有針對其中一個dependent variable 作微分的 linear DEs : : 稱作 linear system

Matrix form of a linear system fn(t) = 0 for all n homogeneous linear system otherwise nonhomogeneous linear system solution vector

其中 可改寫成 Example 2 (text page 306) 皆為 的解

If x1(t0) = r1, x2(t0) = r2, ………., xn(t0) = rn, linear system 可寫成 subject to

8-1-2 基本定理 將 Section 4-1 的幾個定理改成 vector 和 matrix 的型態 8-1-2 基本定理 將 Section 4-1 的幾個定理改成 vector 和 matrix 的型態 [Theorem 8.1.1] If the entries of A and F are continuous on a common interval that contains the point t0, then the initial value problem on the previous page has a unique solution on this interval. (比較 Theorem 4.1.1, page 137)

[Theorem 8.1.2] For the homogeneous linear system if X1, X2, …., Xk are the solution of then is also a solution of [Definition 8.1.3 and Theorem 8.1.5] If the size of A is n × n and X1, X2, …., Xn are the linearly independent solutions of , then X1, X2, …., Xn are said to be a fundamental set of solutions. Then, the general solution of is c1, c2, ….., cn are arbitrary constants (比較 Theorem 4.1.5, page 144)

[Theorem 8.1.3] Linearly dependent / independent 判斷方式 (課本用 | | 來表示 det)

Either W(X1, X2, …., Xn)  0 for every t or W(X1, X2, …., Xn) = 0 linearly independent dependent (比較 Wronskian, page 147)

Example 4 (text page 308) determinant

[Theorem 8.1.6] General solution for nonhomogeneous system subject to 稱作為 complementary function particular solution (比較講義 page 149)

8-1-3 本節要注意的地方 (1) 大部分的定理和 Section 4-1 相似 8-1-3 本節要注意的地方 (1) 大部分的定理和 Section 4-1 相似 (2) 當一個式子出現 2 個 dependent variable 的微分時 先化成講義 page 521 linear system 的型態

Section 8.2 Homogeneous Linear Systems 8-2-1 本節摘要 (A) 解法的限制: 同講義 page 520 ,但多了二個限制 (d) homogeneous (e) 最好是 constant coefficients

(B) 解法 (constant coefficients) size of A: n × n 假設解為 a = 1, 2, …., n 其中 a: A 的 eigenvalue Ka: A 的 eigenvector (AKa = Ka) 證明見講義 529 頁 General solution:

(C) 三種情形 Case 1: A has distinct eigenvalues: 解法如前一頁 Case 2: A has repeated eigenvalues 當 a 的 multiplicities 為 m Case 2.1 可以找到 a 的 m 個 linearly independent eigenvectors 解法同前一頁 Case 2.2 無法找到 a 的 m 個 linearly independent eigenvectors 若只有 1 個 linearly independent eigenvector,將解表示成 :

注意: : Case 2.3 無法找到 a 的 m 個 linearly independent eigenvectors 有超過 1 個 linearly independent eigenvector 其實,也可以用類似方法求解,但較為複雜

Case 3 若 a =  + j 為 A 的eigenvalues, A 為 real matrix b =  − j 必為 A 的eigenvalues 若 Ka = B1 + jB2 為 a 所對應的 eigenvector Kb = B1 − jB2 必為 b 所對應的 eigenvector 此時,可將解改寫成 (D) 名詞與其他  phase portrait (見 535 頁,注意其畫法和觀察法) trajectory (534 頁) phase plane (536 頁) multiplicity (540 頁)

8-2-2 方法 假設解為 (和 Section 4-3 相似)

的問題變成 (和 linear algebra 當中解 eigenvector, eigenvalue 的問題相同)

 是 A 的 eigenvalue K 是 A 的 eigenvector 當  算出後,K 為使得 由 det(A −I) = 0 算出 (稱作 characteristic equation) 成立 的任一個滿足 K  0 的解

Example 1 (text page 313)  = −1, 4 (i) When  = −1 k2 = − k1 設 k1 = 1, k2 = −1

(ii) When  = 4 k2 = 2 k1/3 設 k1 = 3, k2 = 2

trajectory Fig. 8.2.1

c1 < 0, c2 > 0 phase portrait c1 < 0, c2 = 0 c1 = 0, c2 > 0 c1 = 0, c2 < 0 c1 > 0, c2 > 0 c1 > 0, c2 = 0 Fig. 8.2.2

phase plane: 即前頁的 x-y plane repeller attractor

8-2-3 Case 1: Distinct Eigenvalues 根據 eigenvalues ,分成 3 cases Case 1: Distinct eigenvalues Case 2: Repeated eigenvalues Case 3: Complex eigenvalues

Example 2 (text page 314)  = −3, −4, 5 (distinct)

When  = −3 3rd row: k2 = 0 1st row: −k1 + k2 + k3 = −k1 + k3 = 0, k1 = k3 When  = −4,  = 5 (自己練習解解看)

8-2-4 Case 2: Repeated Eigenvalues 有時, det(A − I) 會出現 ( − a)m a 被稱作 eigenvalue of multiplicity m (這種情形較複雜,但也是本節的重點)

Case 2.1 當 a 的 multiplicity 為 m (m > 1) 時,有的時候可以將 m 個 linearly independent eigenvectors 全部找出來。 此時,solutions 解法和 Case 1 相同 注意: 當 A = AT 時, 若 a 的 multiplicity 為 m ,一定可以找到 a 所對應的 m 個 linearly independent eigenvectors Example 3 (text page 316)

(i) 當  = −1 row operation new 2nd row = old 2nd row + 1st row new 3rd row = old 3rd row − 1st row 3 個 variables, 1 個式子 3 − 1 = 2 2 個 linearly independent solutions

2 個 linearly independent solutions (第一個 solution) 設 k1 = 0, k2 = 1 k3 = 1 (第二個 solution) 設 k1 = 1, k2 = 0 k3 = −1 Check: 的確互為 linearly independent 為 A 在  = −1 時的 eigenvectors 小技巧: 任意給定其他 n-1 個 unknowns 的值 再將最後一個 unknown 的值算出來 通常可以得到一個新的 independent solution (但是也有的時候得到的解不為 independent, 所以要 check)

(ii) 當  = 5 算出來的 eigenvector 為 General solution for Example 3:

Case 2.2 當 a 的 multiplicity 為 m (m > 1) 時,有的時候只能找出 1 個 linearly independent eigenvector Ka,1。 將a 所對應的 m 個解表示成 : Ka,1: 唯一滿足 AKa,1 = aKa,1 的 eigenvector Ka,q (q  1) 的求法如後頁

554 當 (p = 1, 2, …, m) 由 (q = 0, 1, …., m−1) 的係數得出 比較 當中 由 Ka,1 求出 Ka,2 由 Ka,2 求出 Ka,3 : 由 Ka,p-1 求出 Ka,p :

註:(1) 課本 page 316 頁中 K11 = K21 = …. = Km1 K22 = K32 = …. = Km2 K33 = K43 = …. = Km3 : (2) 經常有多個 linearly independent 解 在這種情形下,我們只需找出其中一個解即個 (但是必需以可以繼續解下去為條件,如 page 557)

Example 5 (text page 319) eigenvalues: 2, 2, 2 only one independent solution:

選擇其中一個 solution: 注意: 若選擇 Ka,2 為其他的值,最後的解還是一樣的

其中一個 solution: General solution of Example 5

Case 2.3 當 a 的 multiplicity 為 m (m > 1) 時,有的時候只能找出 2 ~ m −1 個 linearly independent eigenvectors 。 例子:Section 8-2 Exercises 31 and 50 three independent solutions:

Set k2 = 1, k4 = 0 Choose 但 無解

Set 無解 的解為 Set (當一個解無法繼續算時,嚐試由其他的解來算)

General solution for Exercises 31 and 50

8-2-5 Case 3: Complex Conjugated Eigenvalues 其實和 Case 1 (distinct eigenvalues) 相同 只是用不同的方式來表示 solutions 當 a =  + j 和b =  − j (,  為 real) 皆為 A 的 eigenvector 且 A 為 real matrix 若 Ka = B1 + jB2 (B1, B2 為 real)是 a 所對應的 eigenvector 則 Kb = B1 − jB2 必為是 b 所對應的 eigenvector Proof:

此時,可將解改寫成 (證明如後)

因此,兩個 linearly independent solutions 可改寫為

Example 6 (text page 322) 已知  = 2i 為其中一個 eigenvalue,所對應的 eigenvector 為 可以迅速判斷 2 個 independent solutions 為

8-2-6 高階線性聯立微分方程的解法 解法:將問題變成 1st order DE

8-2-7 Section 8-2 要注意的地方 (1) 方法適用的情形 (a) linear, (b) 1st order DEs, (c) full rank (n 個 dependent variable 需要 n 個式子), (d) homogeneous, (e) constant coefficients (2) 複習並熟悉算 eigenvector 的方法 (可以研究快速法) (我們只要得出任何一個 eigenvector 或任何一組 linearly independent eigenvectors 即可,因此可以選擇當中較簡單的 ) (3) Case 2 比較複雜,要多加練習 (4) 注意 page 551 找 independent solution 的小技巧 (5) Case 2.2 選擇其中一組解即可 (但是要可以繼續解下去)

(6) 計算前,確定 的係數皆為 1 (standard form) (7) 熟悉原理,才不會背錯公式

Section 8.3 Nonhomogeneous Linear Systems 的 particular solution 本節討論如何找 (方法 1) undetermined coefficients 猜 particular solutions,類似 Section 4-4 (方法 2) variation of parameters ,類似 Section 4-6 (t): fundamental matrix,定義見 page 579

(方法 2) variation of parameters ,with initial conditions 名詞: fundamental matrix

8.3.2 方法一: Undetermined Coefficients 和 Section 4.4 的方法相似 根據 F(t) 來「猜」 particular solution 複習講義 page 191 (1) 出現 tn (2) 出現 cos(at) (3) 出現 exp(bt)

(4) 出現 「綜合」 (5) 只要 F(t) 其中有一個 entry 有某一項 則 particular solution 其他每一個 entry 都要根據這一項來猜 particular solution 的型態 (見 567 頁的注意) (這一點和 Section 4.4 稍有所不同) (6) 和 homogeneous solution 有重覆時,不只乘 t,原來的 term 也保留 (見 569 頁) (這一點也和 Section 4.4 有所不同)

Example 3 (text page 328) solving the complementary function eigenvalues of A: 2, 4 corresponding eigenvectors for  = 2: corresponding eigenvectors for  = 4: complementary function

解 particular solution 因為 ,所以假設 particular solution 為  注意,每一個 entry 皆有 1, t, e−t From

補充的範例 (Example 1 in text page 327 的變型) 注意,這一項要保留 設 particular solution 為 乘 t

choose

8.3.3.1 方法二: Variation of Parameters 先找 complementary function (solution of the associated homogeneous DE) fundamental matrix

令 由於 每個 column 都是 associated homogeneous DE 的解

some constants

Example 4 (text page 331) eigenvalues of A:  = −2, −5 eigenvectors of A : fundamental matrix

8.3.3.2 和 initial value problems 相結合 在此時,可改寫成定積分的型態 Since thus

8.3.4 Section 8.3 需要注意的地方 (1) 2 × 2 matrix 的 eigenvector 快速算法 (2) 注意 undetermined coefficient 的方法和 Section 4.4 異同處 (3) Variation of parameters 的部分,關鍵在是否能將公式背起來 (4) 通常 undetermined coefficient 的方法會比較容易解 而 variation of parameters 較複雜,但適用於任何情形 (5) 同樣記得先算 complementary function (homogeneous 部分 的 solution),再算 particular solution

Section 8.4 Matrix Exponential 把 linear system 當成一般 1st order DE 來解 (比較 Section 2-3) (1) (2) (3) (With initial condition X(t0) = X0) 其中 可以由 Laplace transform (see page 589) 或 eigenvector-eigenvalue decomposition (see page 591) 算出

8.4.2 For Homogeneous Systems solution: 的定義 h = k − 1

8.4.3 For Nonhomogeneous Systems solution: 或 比較: with initial conditions

constant column vector 8.4.4 Computation of 一定有這樣的 initial condition constant column vector 令 X(s) 為 X(t) 的 Laplace transform

Example 1 (text page 335) Determine

殺雞焉用牛刀…… 複習 linear algebra 當中, 的求法 (1) eigenvector-eigenvalue decomposition for A e1, e2, e3, ……, en 皆為 A 的eigenvectors, 皆為 n × 1 的 column 1, 2, 3 , …., n 為 e1, e2, e3, ……, en 所對應的eigenvalues

例如, Example 1 當中

8.4.5 注意 (1) 本節可以解的問題,用 Sections 8-2, 8-3 的方法也可以解 (2) 熟悉公式和 的計算 8.4.5 注意 (1) 本節可以解的問題,用 Sections 8-2, 8-3 的方法也可以解 (2) 熟悉公式和 的計算 (3) 使用 eigenvalue-eigenvector decomposition 的方法時,別忘了將  變成

Chapter 11 Orthogonal Functions and Fourier Series 複習: linear algebra 關於 orthogonal (正交) basis 的介紹 在 linear algebra 當中 (1) inner product (2) orthogonal (3) 若 f1, f2, …., fN 為 orthogonal set, 則

例如 在只有三個 entry 的情形下 是一組 orthogonal set f3 f2 f1 問題:在 continuous 當中該如何定義 orthogonal?

Section 11.1 Orthogonal Functions 11.1.1 綱要:熟悉幾個重要定義 (7) normalize (pp. 604) (1) inner product (pp. 597) (8) complete (pp. 605) (2) orthogonal (pp. 599) (9) orthogonal series expansion (pp. 606) (3) orthogonal set (pp. 600) (10) generalized Fourier series (pp. 606) (4) square norm (pp. 602) (11) weight function (pp. 608) (5) norm (pp. 602) (6) orthonormal set (pp. 602) 學習方式:(1) 可以多和 linear algebra 當中的定義多比較 (2) 複習三角函式的公式 (see pp. 611-612)

11.1.2 定義 (1) inner product on an interval [a, b] (f1, f2 為 real 時) 11.1.2 定義 (1) inner product on an interval [a, b] (f1, f2 為 real 時) 比較: discrete case 補充:more standard definition for inner product with conjugation

Inner product 性質 (a) (f1, f2) = (f2, f1)* (b) (k f1, f2) = k (f1, f2), k 為 scalar (或稱為constant) (c) (f, f) = 0 if and only if f = 0, (f, f) > 0 if and only if f  0, (d) (f1 + f2 , g) = (f1, g) + (f2 , g) discrete case 亦有這些性質

(2) orthogonal on an interval [a, b] (f1, f2 為 real 時) (more standard definition) 或 比較: discrete case 例子: 當 [a, b] = [-1, 1], 1 和 xk (k 為奇數) 互為 orthogonal 注意:任何 even function 和任何 odd function 在 [-a, a] 之間必為 orthogonal

(3) orthogonal set 有一組 functions 0(x), 1(x), 2(x), 3(x), ……….. 若 for m  n 則 0(x), 1(x), 2(x), 3(x), ……….. 被稱作 orthogonal set on an interval [a, b]

Example 1 (text page 399) Show that the set {1, cosx, cos2x, cos3x, …..} is an orthogonal set on the interval [−, ] when one of the functions is 1 when both the two functions are not 1

(4) square norm 比較: discrete case (5) norm (6) orthonormal set 對一個 orthogonal set, 若更進一步的滿足 for all n 則被稱為 orthonormal set

Example 2 (text page 400) Calculate the norms of {1, cosx, cos2x, cos3x, …..} 運用三角函式公式 {1, cosx, cos2x, cos3x, …..} normalization as a orthonormal set

(7) normalize 將 norm 變為 1 (x) 注意,此時 可藉由 normalization, 將 orthogonal set 變成 orthonormal set

(8) complete 若在 interval [a, b] 之間,任何一個 function f(t) 都可以表示成 0(x), 1(x), 2(x), 3(x), ………..的 linear combination 則 0(x), 1(x), 2(x), 3(x), ……….. 被稱作 complete 比較:在 linear algebra 當中 對 3-tuple vector 而言 e1 = [1, 0, 0], e2 = [0, 1, 0], e3 = [0, 0, 1] 為 complete

(9)(10) 若 0(x), 1(x), 2(x), 3(x), ………..為complete 可將 f(x) 表示成 被稱作 (9) orthogonal series expansion 當 0(x), 1(x), 2(x), 3(x), ……….. 不為 orthogonal, cn 不易算 當 0(x), 1(x), 2(x), 3(x), ……….. 為 orthogonal cn 被稱作 (10) generalized Fourier series

當 0(x), 1(x), 2(x), 3(x), ……….. 為 orthonormal

11.1.3 Orthogonal with Weight Function (11) inner product with weight function 其中 w(x) 被稱作 weight function 加上了 weight function 後 (11-1) orthogonal 的定義改成 for m  n (11-2) square norm 的定義改成

(11-3) norm 的定義改成 (11-4) orthonormal 的定義改成 for m  n (11-5) normalize 的算法改成

(11-6) orthogonal series expansion of f(x) 以及 generalize Fourier series 的算法改成

11.1.4 三角函數表 (要複習) cos(a + b) cosa cosb − sina sinb sin(a + b) 11.1.4 三角函數表 (要複習) cos(a + b) cosa cosb − sina sinb sin(a + b) sina cosb + cosa sinb cos(a − b) cosa cosb + sina sinb sin(a − b) sina cosb − cosa sinb cosa cosb [cos(a + b) + cos(a − b)]/2 sina sinb [cos(a − b) − cos(a + b)]/2 sina cosb [sin(a + b) + sin(a − b)]/2

cos(2a) cos2a − sin2a or 1 − 2sin2a or 2cos2a − 1 sin(2a) 2sin a cos a cos2a [cos(2a) + 1]/2 sin2a [1 − cos(2a)]/2

11.1.5 Section 11.1 需要注意的地方 (1) Norm 和 square of norm 要分清楚 作 normalization 時,要除以 norm (2) 熟悉三角函數的公式 (i) 記住幾個,其他的就不難推算出來 (ii) 許多公式可以由 導出來

複習: Legendre polynomials 是一種 orthogonal set if m  n 其他常用的 orthogonal set Hermite polynomials (with weight function) (補充) Chebyshev polynomials (with weight function) (補充) Cosine series Sine series Fourier series