CST如何利用探針求出S參數(shù)
求解器:時(shí)域求解器
入射:平面波 正入射
邊界:open 和 周期
結(jié)構(gòu):3000nm 空氣(1)和 3000nm 玻璃(1.5)
理論上可以算出,電場(chǎng)振幅透過(guò)率平方為0.96
但是如果通過(guò)放探針探測(cè)電場(chǎng)振幅除以入射振幅求出來(lái)的數(shù)值完全不對(duì),請(qǐng)問(wèn)是什么原因,正確的做法是什么?
而且,對(duì)于同一均勻介質(zhì),入射平面波設(shè)置1V/m, 在多個(gè)波長(zhǎng)外放致探針,如果改變介質(zhì)的折射率,探針檢測(cè)到的電場(chǎng)振幅居然會(huì)變,比如介質(zhì)折射率為1,探測(cè)到的是1,折射率為1.5,探測(cè)的就是0.84了。不知道是什么原理。
真實(shí)奇怪的結(jié)果,而且在頻域求解器中也存在同樣的問(wèn)題,通過(guò)探針的方式得出的投射譜和計(jì)算結(jié)果是不一樣的。
嘿嘿。用探針測(cè)s參數(shù),我之前做天線仿真的時(shí)候也試過(guò),我在之前的帖子里面也談過(guò),這是一個(gè)很蛋疼的問(wèn)題。
首先,你用測(cè)到的電場(chǎng)幅度去直接除激勵(lì),如果你除以的是理論的激勵(lì)值的化,我基本可以確定你得不到cst直接計(jì)算給你的值,因?yàn)閏st在計(jì)算中使用了一個(gè)隱形的歸一話,你加的電場(chǎng)是1v/m,實(shí)際計(jì)算中可能是0。8v/m,至于為什么這么做我也不知道。
還有一點(diǎn),如果是需要計(jì)算s參數(shù)的話,因?yàn)閰?shù)S_{11},、S_{12},、S_{21},和S_{22},給出的是每一端口的反射功率波與入射功率波的關(guān)系。所以應(yīng)該是反射場(chǎng)除入射場(chǎng),比如說(shuō)我們假設(shè)全場(chǎng)為Etotal,入射場(chǎng)為Einc.則有反射場(chǎng)Ediff = Etotal-Einc。你用探針測(cè)得的是Etotal,所以正確的算法應(yīng)該是FT{(Etotal-Einc)/Einc}.其中Einc應(yīng)該是歸一化以后,在cst計(jì)算中真實(shí)的入射電場(chǎng)值,及你定義的理論值乘上一個(gè)系數(shù)。我用這個(gè)辦法算過(guò)天線的s參數(shù),得出來(lái)的結(jié)果和cst直接計(jì)算出來(lái)的結(jié)果完全一致,我不太清楚在你這個(gè)問(wèn)題里面是否適用。
關(guān)于第二個(gè)問(wèn)題,也需要注意你探車到的是全場(chǎng),而不只是入射場(chǎng)或者是放射場(chǎng)。如果你已經(jīng)做了減法,則需要注意你減的入射場(chǎng)要是cst真實(shí)計(jì)算中的數(shù)值。用商業(yè)軟件就是這點(diǎn)不好,因?yàn)樗泻芏嘤?jì)算過(guò)程是隱性的,計(jì)算中的物理參數(shù)并非完全受你控制,所以你要確定你設(shè)定的物理?xiàng)l件是否就是cst計(jì)算中的真實(shí)條件,商業(yè)軟件在計(jì)算的時(shí)候因?yàn)閿?shù)學(xué)需要常常會(huì)做一些隱性的歸一化或者等效變換,這樣就會(huì)導(dǎo)致方針的最終結(jié)果是對(duì)的,但是如果我門提取一些中間結(jié)果,然后用自己理解的物理模型來(lái)計(jì)算,就可能導(dǎo)致最終結(jié)果和軟件輸出的不一樣。這也就是我現(xiàn)在為什么要苦13的自己編寫(xiě)仿真代碼的原因。
我搞電磁學(xué)算是半路出家,本科和碩士階段都沒(méi)怎么系統(tǒng)的學(xué)過(guò),用cst也不到半年,如果有說(shuō)的不對(duì)的還望各位大大出來(lái)指正,以免誤導(dǎo)他人。
好的,非常感謝!
利用DFT變換算反射沒(méi)問(wèn)題,透射問(wèn)題還是一樣的??磥?lái)對(duì)于多介質(zhì)仿真探針不能在FDTD求解器下做定量分析了。