本人要用FDTD計(jì)算一個波導(dǎo)縫隙天線的問題。一開始從網(wǎng)上下了一個叫ToyFDTD的開源的c程序做參考。它求解的是一個簡單的矩形波導(dǎo)內(nèi)的場分布，給出了一種最簡單情況下求解FDTD問題的基本框架。然后我就根據(jù)這個例子編自己的問題，加上電壁和吸收邊界條件。
現(xiàn)在講一下我遇到的問題，很多，慢慢講
1 激勵是TE10，采用正弦波和高斯脈沖有什么區(qū)別呢？我看高本慶的書上講用正弦的周期函數(shù)激勵的話一般只要計(jì)算3個周期就可以穩(wěn)定。但是正弦波也是一個隨時間變化的場， 應(yīng)該怎么樣判斷它已經(jīng)穩(wěn)定了呢？
2 激勵的場強(qiáng)如果是1，那么由它產(chǎn)生的場的幅度是不是應(yīng)該都不可能大于1呢？我是這么覺得的，但是ToyFDTD計(jì)算的結(jié)果卻是大于1的？有沒有人能解釋一下呢？
3 關(guān)于吸收邊界的問題。因?yàn)槲疫@個是天線問題，要計(jì)算遠(yuǎn)區(qū)場，所以我的網(wǎng)格空間就不能只包括波導(dǎo)，要用吸收邊界來截?cái)唷Ｎ矣昧?階Mur條件，但是在棱邊上不好處理，不知道有沒有誰有好的方法
4 因?yàn)橐氵h(yuǎn)區(qū)場，覺得近區(qū)遠(yuǎn)區(qū)的轉(zhuǎn)換比較麻煩，有沒有可以參考的程序？

下面是ToyFDTD的源程序，很簡潔，注釋也很詳細(xì)。不過要注意下標(biāo)，因?yàn)樽髡咧皇强紤]理想電壁的情形。
但是從程序上看，他除了一個激勵面之外，其余五個面都是用同樣的方法處理的，感覺有點(diǎn)問題，因?yàn)樵趥鞑シ较蛏喜粦?yīng)該是電壁的。

希望大家和我一起討論討論。
 
1、正弦波的頻譜是一個頻點(diǎn)，而高斯脈沖的頻譜仍然是滿足高斯分布，是一個帶譜。簡諧波的情況一般比較簡單，按照公式或者經(jīng)驗(yàn)選取即可，或者設(shè)一個場觀測點(diǎn)，采兩個周期的場值比較一下就行了。
2、只要輸出能量不大于輸入能量幅度瞬時大于輸入是可以出現(xiàn)的。
3、可以參看西電葛德彪的《電磁場時域有限差分方法》相關(guān)章節(jié)。
4、如果以后要在這個領(lǐng)域安身立命，還是自己編一編，對能力的鍛煉是很重要的。
并沒有看那個程序。不過按一般的算法，由于吸收邊界也只是在無窮遠(yuǎn)才能將入射波衰減為0，而程序只能算有限區(qū)域，所以吸收邊界的邊界也是要截?cái)嗟模詈唵蔚奶幚砭褪窃O(shè)為電壁。不知道是不是和你問的是一個問題，如果說錯了，不要怪我，呵呵。