作者WRO (效應)
看板comm_and_RF
標題[問題] 請問ADI-FDTD如何建立其矩陣方程式
時間Sun Apr 28 10:21:17 2013
大家好,上來打擾大家一個問題
我最近在讀ADI-FDTD,
主要是參考
"A New FDTD Algorithm Based on Alternating-Direction
Implicit Method"這篇文章
但對於它要如何建立矩陣方程式有個疑問百思不得解
依照文中的方法,以TE為例子,Ey和Ex無法直接像Yee的FDTD一樣直接解出來
而是要先創造出足夠多的方程式,並建立成矩陣形式之後一次解出來
我把我的理解畫成下圖(圖中是以解 TE mode 的 Ey 為例子)
http://homepage.ntu.edu.tw/~d97942009/1.bmp
圖中是FDTD的網格,黑框所框起來的代表每一"條"方程式會牽扯到的變數
按照我的理解,第二條方程式則是像下圖的棕色框一樣建立
http://homepage.ntu.edu.tw/~d97942009/2.bmp
根據該文中(2')式的描述,每一條方程式的未知數一定是由緊鄰的三個Ey所組成
簡言之就是 Ey(i-1) + Ey(i) + Ey(i+1) = 已知數
所以要創建夠多的方程式,就只能像圖中一樣,緊鄰的三格逐一向下推動
如果只框未知數,就像下圖
(當然只有Ey是未知,Hz是因為畫圖沒辦法而被框進去的)
http://homepage.ntu.edu.tw/~d97942009/0.bmp
我的問題是,如果按照這種規則,假設共有 N 個 Ey 待解,
那把這個框框推到底,應該也只能創造出 N-2 條方程式
怎麼能解出 N 個未知數??
我想很久還是無法補齊殘缺的兩條方程式,
所以想上來問問有經驗的先進,我的誤解在哪,是不是能有什麼建議?
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 114.37.189.166
1F:→ mp19990920:boundary condition? 111.67.53.154 04/28 15:08