iT邦幫忙

第 11 屆 iThome 鐵人賽

DAY 21
0
AI & Data

連接數學與現實世界的橋樑 -- 數學建模系列 第 21

Day 21 : 動力系統的圖解 -- 向量場圖

  • 分享至 

  • xImage
  •  

今天就來介紹如何透過圖像的方式來描述動力系統在平衡點附近的性質,並用該訊息與向量場圖一起獲得在整個狀態空間上動態行為的圖形描述。

範例

考慮下圖的電路圖,
https://ithelp.ithome.com.tw/upload/images/20190922/20119600PcyYw1upAj.jpg
其中https://chart.googleapis.com/chart?cht=tx&chl=v_c表示電容上的電壓,https://chart.googleapis.com/chart?cht=tx&chl=i_R表示經過電阻的電流,而https://chart.googleapis.com/chart?cht=tx&chl=f(x)為電阻的https://chart.googleapis.com/chart?cht=tx&chl=v-i特徵。給定https://chart.googleapis.com/chart?cht=tx&chl=L%20%3D%201, https://chart.googleapis.com/chart?cht=tx&chl=C%20%3D%20%5Cfrac%7B1%7D%7B3%7D, https://chart.googleapis.com/chart?cht=tx&chl=f(x)%20%3D%20x%5E3%20%2B%204x,確定這個電路隨時間變化的行為。

  1. 提出問題
    根據克希荷夫電流定律,進入一個節點的電流和等於流出的電流和;克希荷夫電壓定律,閉路上的所有電壓之合為0,可以寫出我的們的假設
    變量:
    https://chart.googleapis.com/chart?cht=tx&chl=v_c%20%3D%20電容上的電壓
    https://chart.googleapis.com/chart?cht=tx&chl=i_c%20%3D%20過電容的電流
    https://chart.googleapis.com/chart?cht=tx&chl=v_R%20%3D%20電阻上的電壓
    https://chart.googleapis.com/chart?cht=tx&chl=i_R%20%3D%20過電阻的電流
    https://chart.googleapis.com/chart?cht=tx&chl=v_L%20%3D%20電感上的電壓
    https://chart.googleapis.com/chart?cht=tx&chl=i_L%20%3D%20過電感的電流
    假設:
    https://chart.googleapis.com/chart?cht=tx&chl=C%20%5Cfrac%7Bdv_c%7D%7Bdt%7D%20%3D%20i_c
    https://chart.googleapis.com/chart?cht=tx&chl=v_R%20%3D%20f(i_R)
    https://chart.googleapis.com/chart?cht=tx&chl=L%20%5Cfrac%7Bdi_L%7D%7Bdt%7D%20%3D%20v_L
    https://chart.googleapis.com/chart?cht=tx&chl=i_c%20%3D%20i_R%20%3D%20i_L
    https://chart.googleapis.com/chart?cht=tx&chl=v_c%20%2B%20v_R%20%2B%20v_L%20%3D%200
    https://chart.googleapis.com/chart?cht=tx&chl=L%20%3D%201
    https://chart.googleapis.com/chart?cht=tx&chl=C%20%3D%20%5Cfrac%7B1%7D%7B3%7D
    https://chart.googleapis.com/chart?cht=tx&chl=f(x)%20%3D%20x%5E3%20%2B%204x
    目標:
    確定六個變量隨時間變化的行為
  2. 選擇建模方法
    利用連續時間動力系統為此問題建模,並畫出向量場圖。誠如先前討論的,假設有一個動力系統在平衡點的鄰域內有一階連續偏導數,而https://chart.googleapis.com/chart?cht=tx&chl=A是在平衡點的一階偏導數矩陣。當https://chart.googleapis.com/chart?cht=tx&chl=x在平衡點附近時,該動力系統的行為會與線性系統相似,即https://chart.googleapis.com/chart?cht=tx&chl=x'%20%3D%20F(x)%20%5Capprox%20A(x-%20x_0),對此我們可以透過兩者的向量場圖來進行驗證。
    這在我們要介紹一個新名詞,同胚(homeomorphism),定義為具有連續逆映射的連續函數,https://chart.googleapis.com/chart?cht=tx&chl=G%20%20%3A%20%20%5Cmathbb%20R%5E2%20%5Crightarrow%20%5Cmathbb%20R%5E2
    如果https://chart.googleapis.com/chart?cht=tx&chl=A的特徵值的實部都不為0,則存在一個同胚https://chart.googleapis.com/chart?cht=tx&chl=G,它把系統https://chart.googleapis.com/chart?cht=tx&chl=x'%20%3D%20Ax的向量場圖映射為https://chart.googleapis.com/chart?cht=tx&chl=x'%20%3D%20F(x)的向量場圖,並滿足https://chart.googleapis.com/chart?cht=tx&chl=G(0)%20%3D%20x_0。這也就是說,除了某些扭曲外,https://chart.googleapis.com/chart?cht=tx&chl=x'%20%3D%20F(x)在平衡點https://chart.googleapis.com/chart?cht=tx&chl=x_0附近的向量場看起來像線性系統的向量場。它意味著,我們可以利用線性逼近的方式,來得到非線性動力系統在平衡點附近的行為。
  3. 推導模型的數學表達式
    https://chart.googleapis.com/chart?cht=tx&chl=x_1%20%3D%20i_R%20%3D%20i_L%20%3D%20i_c, https://chart.googleapis.com/chart?cht=tx&chl=x_2%20%3D%20v_c,則可得
    https://chart.googleapis.com/chart?cht=tx&chl=%5Cfrac%7Bx_%7B2%7D'%7D%7B3%7D%20%3D%20x_1
    https://chart.googleapis.com/chart?cht=tx&chl=v_R%20%3D%20x_%7B1%7D%5E3%20%2B%204%20x_1
    https://chart.googleapis.com/chart?cht=tx&chl=x_%7B1%7D'%20%3D%20v_L
    https://chart.googleapis.com/chart?cht=tx&chl=x_2%20%2B%20v_R%20%2B%20v_L%20%3D%200
    整理後得到
    https://chart.googleapis.com/chart?cht=tx&chl=x_%7B1%7D'%20%3D%20f_%7B1%7D(x_1%2C%20x_2)%20%3D%20-x_%7B1%7D%5E3%20-%204%20x_1%20-%20x_2
    https://chart.googleapis.com/chart?cht=tx&chl=x_%7B2%7D'%20%3D%20f_%7B2%7D(x_1%2C%20x_2)%20%3D%203%20x_1
  4. 求解模型
    從上述的方程組可得到向量場圖如下,
    https://ithelp.ithome.com.tw/upload/images/20190923/20119600IdRkHMJbb7.jpg
    從圖中可知平衡點為https://chart.googleapis.com/chart?cht=tx&chl=(0%20%2C%20%200),但難以判定是否穩定,因此對方程組求偏導數,得到
    https://chart.googleapis.com/chart?cht=tx&chl=%5Cfrac%7B%5Cpartial%20f_1%7D%7B%5Cpartial%20x_1%7D%20%3D%20-3%20x_%7B1%7D%5E2%20-%204
    https://chart.googleapis.com/chart?cht=tx&chl=%5Cfrac%7B%5Cpartial%20f_1%7D%7B%5Cpartial%20x_2%7D%20%3D%20-1
    https://chart.googleapis.com/chart?cht=tx&chl=%5Cfrac%7B%5Cpartial%20f_2%7D%7B%5Cpartial%20x_1%7D%20%3D%203
    https://chart.googleapis.com/chart?cht=tx&chl=%5Cfrac%7B%5Cpartial%20f_2%7D%7B%5Cpartial%20x_2%7D%20%3D%200
    計算偏導數矩陣得,https://chart.googleapis.com/chart?cht=tx&chl=A%20%3D%20%5Cleft%20(%20%5Cbegin%7Barray%7D%20%7Bccc%7D%20%20-4%20%26%20-1%20%5C%5C%203%20%26%200%20%5Cend%7Barray%7D%5Cright)
    求奇特徵值,https://chart.googleapis.com/chart?cht=tx&chl=%5Clambda%5E2%20%2B%204%20%5Clambda%20%2B%203%20%20%3D0,得https://chart.googleapis.com/chart?cht=tx&chl=%5Clambda%20%3D%20-3%20%2C%20-1
    因為兩者的實部都是負的,所以平衡態是穩定的。
    進一步求其特徵向量,當https://chart.googleapis.com/chart?cht=tx&chl=%5Clambda%20%3D%20-3,對應的特徵向量為https://chart.googleapis.com/chart?cht=tx&chl=%5Cleft(%20%5Cbegin%7Barray%7D%7Bccc%7D%20%20%20-1%20%5C%5C%20%201%20%20%5Cend%7Barray%7D%5Cright)%20e%5E%7B-3t%7D
    https://chart.googleapis.com/chart?cht=tx&chl=%5Clambda%20%3D%20-1,對應的特徵向量為https://chart.googleapis.com/chart?cht=tx&chl=%5Cleft(%20%5Cbegin%7Barray%7D%7Bccc%7D%20-1%20%5C%5C%203%20%5Cend%7Barray%7D%5Cright)%20e%5E%7B-t%7D
    則線性系統的通解可以寫成,
    https://chart.googleapis.com/chart?cht=tx&chl=%5Cleft(%20%5Cbegin%7Barray%7D%7Bccc%7D%20x_1%20%20%5C%5C%20%20x_2%20%5Cend%7Barray%7D%5Cright)%20%3D%20c_1%20%5Cleft(%20%5Cbegin%7Barray%7D%7Bccc%7D%20-1%20%20%5C%5C%20%201%20%5Cend%7Barray%7D%5Cright)%20e%5E%7B-3t%7D%20%2B%20c_2%20%5Cleft(%20%5Cbegin%7Barray%7D%7Bccc%7D%20-1%20%20%5C%5C%20%203%20%5Cend%7Barray%7D%5Cright)%20e%5E%7B-t%7D,其中https://chart.googleapis.com/chart?cht=tx&chl=c_1%2C%20c_2為任意實常數。根據此通解畫出來的向量場如下,
    https://ithelp.ithome.com.tw/upload/images/20190923/20119600IDwCBfkfKt.jpg
    綜合非線性與線性的向量場圖得到原非線性系統的完整向量場圖如下,
    https://ithelp.ithome.com.tw/upload/images/20190923/201196004ZMkNr7mar.jpg
  5. 表達分析結果
    不論電路的初始狀態如何,電流與電容上的壓降最終都會趨近於0。

上一篇
Day 20 : 離散系統的特徵值方法
下一篇
Day 22 : 動態系統的模擬求解
系列文
連接數學與現實世界的橋樑 -- 數學建模30
圖片
  直播研討會
圖片
{{ item.channelVendor }} {{ item.webinarstarted }} |
{{ formatDate(item.duration) }}
直播中

尚未有邦友留言

立即登入留言