雒向東 趙宇杰 海波 梁曄

【摘? ?要】? ?基于圓柱和球坐標(biāo)系中圓柱和球型矢量波函數(shù)，將自由空間圓柱和球型并矢格林函數(shù)應(yīng)用于介質(zhì)圓柱和球進(jìn)而推出二者的第三類電型并矢格林函數(shù)，并就求解過程特性進(jìn)行比較分析，其結(jié)果可為基于DGF方法處理圓柱和球型介質(zhì)電磁散射等問題提供理論依據(jù)。

【關(guān)鍵詞】? ?柱球坐標(biāo)系;介質(zhì)柱球;矢量波函數(shù);并矢格林函數(shù);比較研究

A Comparative Study of the Characteristics of the Solving Process of

Electric Dyadic Green′s Functions of Medium Column and Sphere

Luo Xiangdong1，2， Zhao Yujie1， Hai Bo1， Liang Ye1

（1. Lanzhou City University， Lanzhou 730070， China;

2. Editorial Department of Journal of Gansu Normal Colleges， Lanzhou 730070， China）

【Abstract】? ? Based on the cylindrical and spherical vector wave functions in the cylindrical and spherical coordinate systems， applying the free space cylindrical and spherical dyadic Green's functions to medium column and sphere， this paper deduces the third kind of electric dyadic Green's functions of the medium cylinder and sphere and makes a comparative analysis of the characteristics of their solving process. The results can provide theoretical basis for dealing with electromagnetic scattering of cylindrical and spherical media with the DGF method.

【Key words】? ? ?cylindrical and spherical coordinate systems; medium cylinder and sphere; vector wave functions; dyadic Green's functions; comparative study

〔中圖分類號〕? O441.4? ? ? ? ? ? ? ? ? ? 〔文獻(xiàn)標(biāo)識碼〕? A ? ? ? ? ? ? ?〔文章編號〕 1674 - 3229（2022）01- 0037 - 04

0? ? ?引言

圓柱體、球體是工程應(yīng)用中常見的電磁輻射與散射模型，這類電磁場邊值問題的求解往往涉及柱型和球型矢量波函數(shù)的構(gòu)建[1-4]。在圓柱和球坐標(biāo)系中分別應(yīng)用圓柱和球型矢量波函數(shù)表示自由空間輻射或散射場極為方便，當(dāng)獲得自由空間圓柱和球型矢量波函數(shù)展開式后，應(yīng)用散射疊加方法可導(dǎo)出其他類型的并矢格林函數(shù)。應(yīng)用并矢格林函數(shù)（簡稱DGF）方法解決電磁理論及工程邊值問題，主要難點就是要構(gòu)建并矢格林函數(shù)。對這一問題國內(nèi)外學(xué)者已做了大量研究[5-14]，如對規(guī)則形狀波導(dǎo)、常用坐標(biāo)系等問題已有一致的結(jié)果，但對介質(zhì)圓柱和球體散射場等問題求解的比較研究未見報道。本文基于圓柱和球坐標(biāo)系中圓柱和球型矢量波函數(shù)，從自由空間構(gòu)建的并矢格林函數(shù)出發(fā)，求解介質(zhì)圓柱和球的第三類電型并矢格林函數(shù)，并就其求解過程特性進(jìn)行比較分析，其結(jié)論可為基于DGF方法解決介質(zhì)圓柱和球形物體散射場等問題提供理論依據(jù)。

1? ? ?圓柱和球型矢量波函數(shù)基本理論

1.1? ?圓柱矢量波函數(shù)基本理論

圓柱矢量波函數(shù)[M、N]定義為[Meonλ、Neonλ]兩組矢量波函數(shù)[8]。

[Meonλ（h）=?×Jn（λr）cossinnφe-ihzez=?nJn（λr）rsincosnφer-?Jn（λr）?rcossinnφeφe-ihz] （1）

[Neonλ（h）=1κλ?×?×Jn（λr）cossinnφe-ihzez=1κλ-ih?Jn（λr）?rcossinnφer±ihnJn（λr）rsincosnφeφ+λ2Jn（λr）cossinnφeze-ihz]（2）

式（1）中雙行符號是同時表示兩個式子的一種符號表示方法，實質(zhì)上是指以下兩式：

[Menλ（h）=?×[Jn（λr）cosnφe-ihzez]=[-nJn（λr）rsinnφer-?Jn（λr）?rcosnφeφ]e-ihzMonλ（h）=?×[Jn（λr）sinnφe-ihzez]=[nJn（λr）rcosnφer-?Jn（λr）?rsinnφeφ]e-ihz]? （3）

文中其他雙行符號式子含義同此，式（3）中下標(biāo)“[e]”表示偶函數(shù)，“[o]”表示奇函數(shù)，式（1）和（2）在整個空間區(qū)域（[∞>r≥0]，2[π][>φ≥0]，[∞>z>-∞]）均有定義，且滿足以下正交歸一性質(zhì)[8]：65E71114-E2A4-42E9-96FB-4A0F5E09D398

[Meonλ（h）?Neon′λ′（-h′）dV=0] （4）

[Meonλ（h）?Meon′λ′（-h）dV=0? ? ? ? ? ? ? ?（n? ?n′）（1+δ0）2π2λδ（λ-λ′）δ（h-h′）? ?（n=n′）]? （5）

[Neonλ（h）?Neon′λ′（-h′）dV=0? ? ? ? ? ? ? ? ? （n? ?n′）（1+δ0）2π2λδ（λ-λ′）δ（h-h′）（n=n′）]? （6）

式（4）-（6）中積分區(qū)域為整個空間，本征值[h]和[λ]為連續(xù)取值。

1.2? ?球型矢量波函數(shù)基本理論

球型矢量波函數(shù)定義為[15]：

[φeomn（κ）=jn（κR）Pmn（cosθ）cossinmφ]? ? ? ? （7）

[Meomn（κ）=?×[φeomnR]=?msinθjn（κR）Pmn（cosθ）sincosmφeθ-jn（κR）?Pmn（cosθ）?θcossinmφeφ] （8）

[Neomn（κ）=1κ?×?×[φeomnR]=n（n+1）κRjn（κR）Pmn（cosθ）cossinmφer+1κR??R[Rjn（κR）]?Pmn（cosθ）?θcossinmφeθ?msinθPmn（cosθ）sincosmφeφ] （9）

球型矢量波函數(shù)正交歸一關(guān)系為：

[Meomn（κ）?Neom′n′（κ′）dV=0] （10）

[Meomn（κ）?Meom′n′（κ′）dV=0? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ?m≠m′，n≠n′（1+δ0）π2n（n+1）（n+m）！κ2（2n+1）（n-m）！δ（κ-κ′） m=m′，n=n′]（11）

[δ0=1m=00? ? ? ? ? ? ? ? ? m≠0 ]? ? ?（12）

[Neomn（κ）?Neom′n′（κ′）dV=0? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? m≠m′，n≠n′（1+δ0）π2n（n+1）（n+m）！κ2（2n+1）（n-m）！δ（κ-κ′）? ? ?m=m′，n=n′] （13）

式（10）（11）和（13）中積分區(qū)域為整個空間區(qū)域。

2? ? ?介質(zhì)圓柱電型并矢格林函數(shù)

設(shè)自由空間圓柱體是各向同性均勻媒質(zhì)構(gòu)成的介質(zhì)圓柱，電型并矢格林函數(shù)在[r=a]界面上滿足電場、磁場在界面切向矢量連續(xù)的邊界條件，這里電型并矢格林函數(shù)為第三類。采用[G（11）e（R，R′）和G（21）e（R，R′）]分別表示圓柱外部和內(nèi)部區(qū)域上的電型并矢格林函數(shù)，源被置于介質(zhì)圓柱外，圓柱內(nèi)外空間對應(yīng)參數(shù)定義為：

[k1=ω（μ0ε0）12，η=（k21-h2）12]? ? [r>a]? ? （14）

[k2=ω（με）12，ξ=（k22-h2）12]? ? ? ?[r

式（14）（15）中[μ、ε]分別表示介質(zhì)柱磁導(dǎo)率和介電常數(shù)，它們可以是實數(shù)也可以是復(fù)數(shù)。基于兩媒質(zhì)分界面存在透射、反射效應(yīng)，第三類電型并矢格林函數(shù)可設(shè)為：

[G（11）e（R，R′）=Ge0（R，R′）+G（11）es（R，R′）? r≥a] （16）

[G（21）e（R，R′）=G（21）es（R，R′）? ? ? ? ?r≤a] （17）

式（16）中[Ge0（R，R′）]是柱坐標(biāo)系中自由空間電型并矢格林函數(shù)，表達(dá)式為[8]：

[Ge0（R，R′）=-1k2ererδ（R-R′）+i8π-∞∞dhn=0∞2-δ0η2M（1）η（h）M′η（-h）+N（1）η（h）N′η（-h）（r>r′）Mη（h）M′（1）η（-h）+Nη（h）N′（1）η（-h） （r

式（18）中上行符號對應(yīng)[r>r′]，下行符號對應(yīng)[r

從物理角度看，[G（11）es（R，R′）]來自介質(zhì)圓柱對[Ge0（R，R′）]的散射，必須滿足無窮遠(yuǎn)處的輻射條件和發(fā)散波的要求;[G（21）e（R，R′）]來自[Ge0（R，R′）]通過圓柱界面的透射，必須滿足在圓柱內(nèi)部有界和駐波特性的要求。因此[G（11）es（R，R′）]和[G（21）es（R，R′）]應(yīng)與[Ge0（R，R′）]有相同的激勵項。考慮到圓柱交界面上滿足的邊界條件，[G（11）es（R，R′）]和[G（21）es（R，R′）]必然有以下的展開式：65E71114-E2A4-42E9-96FB-4A0F5E09D398

[G（11）es（R，R′）=i8π-∞∞dhn=0∞2-δ0η2[AeoηM（1）eoη（h）+BoeηN（1）oeη（h）]M′（1）eoη（-h）+[CeoηN（1）eoη（h）+DoeηM（1）oeη（h）]N′（1）eoη（-h）] （19）

[G（21）es（R，R′）=i8π-∞∞dhn=0∞2-δ0η2[aeoξMeoξ（h）+boeξNoeξ（h）]M′（1）eoη（-h）+[ceoξNeoξ（h）+doeξMoeξ（h）]N′（1）eoη（-h）] （20）

為了求得待定系數(shù)[A、B、C、D和a、b、c、][d，]

應(yīng)用在[r=a]圓柱面上滿足的邊界條件：

[er×[G（11）e（R，R′）-G（21）e（R，R′）]r=a=0] （21）

[er×?×G（11）e（R，R′）μ0-?×G（21）e（R，R′）μr=a=0] （22）

利用這兩式可得到含有待定系數(shù)在內(nèi)的16個線性代數(shù)方程，求解這些代數(shù)方程，便可得到16個待定參數(shù)。將式（19）和（20）代入式（21）得：

[er×[-1k2ererδ（R-R′）+i8π-∞∞dhn=0∞2-δ0η2[Mη（h）M′（1）η（-h）+Nη（h）N′（1）η（-h）]+i8π-∞∞dhn=0∞2-δ0η2[AeoηM（1）eoη（h）+BoeηN（1）oeη（h）]M′（1）eoη（-h）+[CeoηN（1）eoη（h）+DoeηM（1）oeη（h）]N′（1）eoη（-h）-i8π-∞∞dhn=0∞2-δ0η2[aeoξMeoξ（h）+boeξNoeξ（h）]M′（1）eoη（-h）+[ceoξNeoξ（h）+doeξMoeξ（h）]N′（1）eoη（-h）]r=a=0? ? ? ? ? ? ? ? ? ?（23）]

式（23）等效于以下兩個方程：

[er×[（Meoη（h）+AeoηM（1）eoη（h）+BoeηN（1）oeη（h）-aeoξMeoξ（h）-boeξNoeξ（h）]r=a=0]? ? ? ? ? ?（24）

[er×[（Neoη（h）+CeoηN（1）eoη（h）+DoeηM（1）oeη（h）-ceoξNeoξ（h）-doeξMoeξ（h）]r=a=0]? ? ? ? ? ? （25）

對式（24）和（25）兩式求解得：

[-1k1η2H（1）n（ηa）Boeη+1k2ξ2Jn（ξa）boeξ=0-?H（1）n（ηa）?aAeoη?ihnH（1）n（ηa）ak1Boeη+?Jn（ξa）?aaeoξ±ihnJn（ξa）ak2boeξ=?Jn（ηa）?a-η2k1H（1）n（ηa）Ceoη+ξ2k2Jn（ξa）ceoξ=1k1η2Jn（ηa）±ihnH（1）n（ηa）ak1Ceoη-?H（1）n（ηa）?aDoeη?ihnJn（ξa）ak2ceoξ+?Jn（ξa）?adoeξ=?ihnJn（ηa）ak1] （26）

[k1=κηk2=κξ]? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ?（27）

將式（19）和（20）代入式（22）得：

[er×[i8π-∞∞dhn=0∞12-δ0η21μ0[?×Mη（h）′（1）η（-h）+?×Nη（h）′（1）η（-h）]+][i8π-∞∞dhn=0∞2-δ0η21μ0[Aeoη?×M（1）eoη（h）+Beoη?×N（1）oeη（h）]M′（1）eoη（-h）+[?×CeoηN（1）eoη（h）+Doeη?×M（1）oeη（h）]N′（1）eoη（-h）-i8π-∞∞dhn=0∞2-δ0η21μ[aeoξ?×Meoξ（h）+boeξ?×Noeξ（h）]M′（1）eoη（-h）+[ceoξ?×Neoξ（h）+doeξ?×Moeξ（h）]N′（1）eoη（-h）]r=a=0? ? ? ? ?（28）]

式（28）等效于以下兩個方程：

[er×[k1μ0Neoη（h）+Aeoηk1μ0N（1）eoη（h）+Boeηk1μ0M（1）oeη（h）-aeoξNeoξ（h）k2μ-k2μboeξMoeξ（h）]r=a=0]（29）? ? ? ?[er×[Ceoηk1μ0M（1）eoη（h）+Doeηk1μ0N（1）oeη（h）+k2μceoξMeoξ（h）+k2μdoeξNoeξ（h）+k1μ0Meoη（h）]r=a=0] （30）

對式（29）和（30）求解得：

[-η2μ0H（1）n（ηa）Aeoη+ξ2μJn（ξa）aeoξ=η2μ0Jn（ηa）±ihnH（1）n（ηa）aμ0Aeoη-k1μ0?H（1）n（ηa）?aBoeη?ihnJn（ξa）aμaeoξ+k2μ?Jn（ξa）?aboeξ=?ihnJn（ηa）aμ0η2μ0H（1）n（ηa）Doeη+ξ2μJn（ξa）doeξ=0-k1μ0?H（1）n（ηa）?aCeoη?ihnH（1）n（ηa）aμ0Doeη-k2μ?Jn（ξa）?aceoξ?ihnJn（ξa）aμdoeξ=k1μ0?Jn（ηa）?a]? （31）65E71114-E2A4-42E9-96FB-4A0F5E09D398

將式（26）和（31）分別寫成矩陣形式為：

[0-1k1η2H（1）n（ηa）01k2ξ2Jn（ξa）-?H（1）n（ηa）?a?ihnH（1）n（ηa）ak1?Jn（ξa）?a±ihnJn（ξa）ak2-η2μ0H（1）n（ηa）0ξ2μJn（ξa）0±ihnH（1）n（ηa）aμ0-k1μ0?H（1）n（ηa）?a?ihnJn（ξa）aμk2μ?Jn（ξa）?aAeoηBoeηaeoξboeξ=0?Jn（ηa）?aη2μ0Jn（ηa）?ihnJn（ηa）aμ0? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ?（32）]

[-η2k1H（1）n（ηa）? ? ? ? ?0? ? ? ? ? ? ? ? ? ? ? ? ξ2k2Jn（ξa）? ? ? ? 0±ihnH（1）n（ηa）ak1? ? ? ?-?H（1）n（ηa）?a? ? ? ?ihnJn（ξa）ak2? ?Jn（ξa）?a0? ? ? η2μ0H（1）n（ηa）? ? ? ? 0? ? ? ? ? ? ? ? ? ? ? ?ξ2μJn（ξa）-k1μ0?H（1）n（ηa）?a? ??ihnH（1）n（ηa）aμ0? ?-k2μ?Jn（ξa）?a? ??ihnJn（ξa）aμCeoηDoeηceoξdoeξ=1k1η2Jn（ηa）?ihnJn（ηa）ak10k1μ0?Jn（ηa）?a? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? （33）]

求解此矩陣方程就可獲得待定系數(shù)[A、B、C、D和a、b、c、d]的解[16]，將待定系數(shù)代入式（19）和（20），將式（19）和（20）再代入式（16）和（14）即可獲得介質(zhì)圓柱的第三類電型并矢格林函數(shù)。

3? ? ?介質(zhì)球電型并矢格林函數(shù)

介質(zhì)球電型并矢格林函數(shù)為第三類，此時空間被球面分為兩個區(qū)域，即球的內(nèi)部和外部。

假設(shè)：

[k1=ωμ1ε1r≥a]? ? ? ? ? ? ? ? ? ? ? ? ? ?（34）

[k2=ωμ2ε2r≤a]? ? ? ? ? ? ? ? （35）

若介質(zhì)球置于空氣中，則[ε1=ε0，ε2=ε，μ1=μ2][=μ0]，[ε]表示介質(zhì)球介電常數(shù)，一般情況下這些常數(shù)可以是任意的，不受任何限制。在球坐標(biāo)系中自由空間并矢格林函數(shù)[Ge0（R，R′）]本征展開式為[15]：

[Ge0（R，R′）=-1k2RRδ（R-R′）+ik4πm，nCmnM（1）（k1）M′（k1）+N（1）（k1）N′（k1） R>R′M（k1）M′（1）（k1）+N（k1）N′（1）（k1）? ?R

依據(jù)散射基本理論假設(shè)：

[G（11）e（R，R′）=Geo（R，R′）+G（11）es（R，R′）? R≥a]? ? ? （37）

[G（21）e（R，R′）=G（21）es（R，R′）? ? ? ? ? ? ? ? ? R≤a]? ? ? ? ? （38）

因源被置于介質(zhì)球的外部，故散射項可設(shè)為：

[G（11）es（R，R′）=ik14πm，nCmnAnM（1）（k1）M′（1）（k1）+BnN（1）（k1）N′（1）（k1）] （39）

[G（21）es（R，R′）=ik14πm，nCmnCnM（k2）M′（1）（k1）+DnN（k2）N′（1）（k1）] （40）65E71114-E2A4-42E9-96FB-4A0F5E09D398

在球面上第三類電型并矢格林函數(shù)應(yīng)滿足以下邊界條件：

[[er×G（11）e]R=a=[er×G（21）e]R=a]? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ?（41）

[[1μ1er×?×G（11）e]R=a=[1μ2er×?×G（21）e]R=a]? （42）

將式（39）和（40）代入式（41）和（42）得：

[[er×[ik14πm，nCmnM（k1）M′（1）（k1）+N（k1）N′（1）（k1）+ik14πm，nCmnAnM（1）（k1）M′（1）（k1）+BnN（1）（k1）N′（1）（k1）]]R=a=[er×ik14πm，nCmnCnM（k2）M′（1）（k1）+DnN（k2）N′（1）（k1）]R=a]? ? ? ? ? ? ? ? ? ?（43）

[[1μ1er×?×[ik14πm，nCmnM（k1）M′（1）（k1）+N（k1）N′（1）（k1）+ik14πm，nCmnAnM（1）（k1）M′（1）（k1）+BnN（1）（k1）N′（1）（k1）]]R=a=[1μ2er×?×ik14πm，nCmnCnM（k2）M′（1）（k1）+DnN（k2）N′（1）（k1）]R=a? ? ? （44）]? ? ? ? 進(jìn)一步推得：

[[er×[M（k1）+AnM（1）（k1）]]R=a=[er×CnM（k2）]R=a[er×[N（k1）+BnN（1）（k1）]]R=a=[er×DnN（k2）]R=a[k1μ1er×[M（k1）+BnM（1）（k1）]]R=a=[k2μ2er×DnM（k2）]R=a[k1μ1er×[N（k1）+AnN（1）（k1）]]R=a=[k2μ2er×CnN（k2）]R=a]? ? ? ? ? ? ? ? ? ? ?（45）

將[M（k1）、M（1）（k1）、M（k2）、N（k1）、N（1）（k1）、N（k2）]關(guān)系式代入式（45）得：

[jn（ρ1）+Anh（1）n（ρ1）=Cnjn（ρ2）k1μ1[ρ1jn（ρ1）′ρ1+Anρ1h（1）n（ρ1）′ρ1]=k2μ2[Cnρ2jn（ρ2）′ρ2]ρ1jn（ρ1）′ρ1+Bnρ1h（1）n（ρ1）′ρ1=Dnρ2jn（ρ2）′ρ2k1μ1jn（ρ1）+Bnh（1）n（ρ1）=Dnk2μ2jn（ρ2）]? ? ? ? ? ? ? ? ? ? ?（46）

此處，[ρ1=k1a，ρ2=k2a，]上標(biāo)中撇號表示微商運算，上面四個方程可分為兩個代數(shù)方程組，解出系數(shù)便可得到介質(zhì)球的第三類電型并矢格林函數(shù)[G（11）e（R，R′）]和[G（21）e（R，R′）]。

4? ? ?介質(zhì)圓柱與介質(zhì)球電型并矢格林函數(shù)求解比較分析

從式（19）和（20）可以看出，為了滿足介質(zhì)圓柱面上的邊界條件，散射并矢格林函數(shù)包含了待定系數(shù)[Bη、Dη、bξ、dξ]的耦合項。這表明對于介質(zhì)圓柱而言，一個入射的TE波（或TM波）將激勵出TE波和TM波同時存在的散射場。從式（19）和（20）還可看到一個偶性的[M]函數(shù)與一個奇性的[N]函數(shù)在[φ]分量上有相同的角函數(shù)，它們在散射或透射項并矢格林函數(shù)中起類似的作用，反之亦然。對于介質(zhì)球，不必引入TE模（[M]描述）和TM模（[N]描述）之間的耦合項，而這種耦合對介質(zhì)柱是必不可少的。在這兩個不同的邊值問題中，圓柱矢量波函數(shù)的領(lǐng)示矢量為[z]，而球矢量波函數(shù)的領(lǐng)示矢量為[R]，它們的特性是完全不同的。綜上所述，通過比較介質(zhì)球與介質(zhì)圓柱第三類電型并矢格林函數(shù)的求解過程，可以看到球形邊界情形是極其簡單的。

[參考文獻(xiàn)]

[1] 章日榮，邸英杰，李渠塘.均勻柱波導(dǎo)內(nèi)電磁場與并矢格林函數(shù)[J].無線電通信技術(shù)，2000，26（1）：23-25.

[2] 沈忠祥.無輻射介質(zhì)波導(dǎo)中的并矢格林函數(shù)[J].微波學(xué)報，1992，8（2）：37-41.

[3] 劉繼東，張德鋒，朱震，等.圓柱導(dǎo)體電磁散射問題分析[J].電光與控制，2007，14（1）：144-146.

[4] 朱震海，張家宗，洪偉，等.電大尺寸物體電磁輻射與散射特性幾種分析方法的比較研究[J].電波科學(xué)學(xué)報，1995，10（1）：154-162.

[5] Hansen W W.Directional characteristics of any antenna over a plane earth[J].Journalof Applied Physics，1936，7（12）：460-465.

[6] Hansen W W. Transformations useful in certain antenna calculations[J]. Journalof Applied Physics，1937，8（4）：282-286.

[7]JULIUSADAMSSTRATTON.Electromagnetic theory[M].New York：McGraw-Hill，1941：392-423.

[8] 戴振鐸，魯述.電磁理論中的并矢格林函數(shù)[M].武漢：武漢大學(xué)出版社，2005：108-136.

[9] 宋文淼.規(guī)則邊界下電磁場的本征函數(shù)和并矢格林函數(shù)[J].電子科學(xué)學(xué)刊，1986，8（3）：188-195.

[10] LAKHTAKIA A，MACKAY TG.Dyadic green function for an electromagnetic medium inspired by general relativity[J].Chinese physics letters，2006，23（4）：832-833.

[11] 柯亨玉，黃錫文.并矢格林函數(shù)在圓柱坐標(biāo)系中的本征展開及其在幾何繞射理論中的應(yīng)用[J].武漢大學(xué)學(xué)報（自然科學(xué)版），1991（2）：41-52.

[12] 雒向東，張明，海波，等.用磁型格林函數(shù)法構(gòu)建平行板波導(dǎo)波函數(shù)和電場公式[J].貴州師范學(xué)院學(xué)報，2021，37（9）：1-5.

[13] 朱敏.柱形導(dǎo)波結(jié)構(gòu)中場的并矢格林函數(shù)的構(gòu)建[J].真空電子技術(shù)，2004（1）：8-11.

[14] 雒向東，張明，海波，等.幾種格林函數(shù)方法的形成及其應(yīng)用比較研究[J].太原師范學(xué)院學(xué)報（自然科學(xué)版），2021，20（3）：62-69.

[15] CHEN.TO TAI.Dyadic Green′n function in electromagnetic theory[J].New York：IEEE Press，1993：168-170.

[16] 趙利輝，宋紅偉.伴隨矩陣的性質(zhì)及其在解題中的應(yīng)用[J].廊坊師范學(xué)院學(xué)報（自然科學(xué)版），2013，13（4）：24-26.65E71114-E2A4-42E9-96FB-4A0F5E09D398