作者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