簡(jiǎn)介這份《電磁場(chǎng)與電磁波矢量分析》PPT課件源自電子科技大學(xué)編寫(xiě)、高等教育出版社與高等教育電子音像出版社2005年出版的教材面向電子信息、通信工程等專業(yè)本科生及考研復(fù)習(xí)者用于夯實(shí)電磁場(chǎng)理論的數(shù)學(xué)基礎(chǔ)。課件系統(tǒng)梳理矢量代數(shù)、三種常用正交曲線坐標(biāo)系、標(biāo)量場(chǎng)的梯度、矢量場(chǎng)的通量與散度、環(huán)流與旋度、無(wú)旋場(chǎng)與無(wú)散場(chǎng)、拉普拉斯運(yùn)算以及格林定理和亥姆霍茲定理并配有矢量加減、點(diǎn)積、叉積、標(biāo)量三重積與矢量三重積的坐標(biāo)分量推導(dǎo)便于理解電場(chǎng)與磁場(chǎng)性質(zhì)及波動(dòng)方程求解。資源包共1個(gè)PPT文件約3.56MB以幻燈片形式呈現(xiàn)公式推導(dǎo)與圖示適合課堂同步學(xué)習(xí)或考前快速回顧。目前已有114人學(xué)習(xí)可作為電磁場(chǎng)與電磁波課程入門(mén)的配套參考。1. 電磁場(chǎng)與電磁波矢量分析.ppt一份課件背后藏著的場(chǎng)論入門(mén)硬骨頭如果你正在搜「電磁場(chǎng)與電磁波矢量分析.ppt」大概率不是想找一份能直接交差的幻燈片而是被梯度、散度、旋度這三個(gè)詞卡住了——課本翻了三遍公式抄了一整頁(yè)合上書(shū)還是不知道哈密頓算子到底在算什么。這份課件標(biāo)題指向的是電磁場(chǎng)理論里最底層也最容易被跳過(guò)的一關(guān)用矢量分析把「場(chǎng)」這件事從幾何直覺(jué)翻譯成可計(jì)算的微分方程。它解決的不是某一道題而是讓你后面看麥克斯韋方程組時(shí)不再靠死記硬背。適合兩類人正在學(xué)電磁場(chǎng)與電磁波、信號(hào)與系統(tǒng)相關(guān)課程的學(xué)生以及需要補(bǔ)數(shù)學(xué)基礎(chǔ)再上手天線、微波、電磁兼容仿真的工程師。矢量分析不過(guò)關(guān)后面全是玄學(xué)。2. 矢量分析到底在算什么從標(biāo)量場(chǎng)到矢量場(chǎng)的三個(gè)算子2.1 為什么電磁場(chǎng)必須先過(guò)矢量分析這一關(guān)電磁場(chǎng)和電路最大的區(qū)別在于電路里電壓電流是標(biāo)量節(jié)點(diǎn)和支路是離散的而電磁場(chǎng)里電場(chǎng)強(qiáng)度 E、磁感應(yīng)強(qiáng)度 B 都是空間每一點(diǎn)都有大小和方向的矢量函數(shù)。你沒(méi)法用「一個(gè)數(shù)」描述它必須用「一個(gè)矢量場(chǎng)」描述它。這就帶來(lái)三個(gè)基本問(wèn)題場(chǎng)在某個(gè)方向變化多快場(chǎng)從某點(diǎn)往外是發(fā)散還是匯聚場(chǎng)在局部有沒(méi)有旋轉(zhuǎn)趨勢(shì)這三個(gè)問(wèn)題分別對(duì)應(yīng)梯度、散度和旋度。很多課件把這三個(gè)算子直接甩公式學(xué)生背完就忘。真正要建立的是幾何圖像梯度是標(biāo)量場(chǎng)最陡上升方向散度是矢量場(chǎng)通量密度旋度是環(huán)量密度。電磁場(chǎng)里電場(chǎng)由電位負(fù)梯度得到磁場(chǎng)是無(wú)散場(chǎng)變化的電場(chǎng)和磁場(chǎng)互相激發(fā)靠的是旋度方程。你后面看到的 ?×E -?B/?t 和 ?×H J ?D/?t本質(zhì)就是旋度算子在描述「場(chǎng)怎么卷起來(lái)」。所以矢量分析不是數(shù)學(xué)課的前置知識(shí)它就是電磁場(chǎng)語(yǔ)言本身。跳過(guò)它后面所有推導(dǎo)都是黑匣子。2.2 哈密頓算子的三種用法與坐標(biāo)系選擇哈密頓算子 ? 在不同坐標(biāo)系下形式不同這是初學(xué)者翻車的重災(zāi)區(qū)。直角坐標(biāo)系下它寫(xiě)成三個(gè)偏導(dǎo)分量看起來(lái)最干凈但一旦問(wèn)題有球?qū)ΨQ或柱對(duì)稱還硬用直角坐標(biāo)就是自找麻煩。坐標(biāo)系梯度 ?φ散度 ?·A旋度 ?×A直角(?φ/?x, ?φ/?y, ?φ/?z)?Ax/?x?Ay/?y?Az/?z行列式展開(kāi)柱坐標(biāo)(?φ/?ρ, (1/ρ)?φ/?φ, ?φ/?z)(1/ρ)?(ρAρ)/?ρ…含 1/ρ 因子球坐標(biāo)(?φ/?r, (1/r)?φ/?θ, …)(1/r2)?(r2Ar)/?r…含 1/(r2sinθ) 因子選坐標(biāo)系的判斷標(biāo)準(zhǔn)很簡(jiǎn)單邊界形狀決定坐標(biāo)系。平行板用直角同軸線用柱坐標(biāo)點(diǎn)電荷和天線輻射用球坐標(biāo)。我一般會(huì)先畫(huà)邊界邊界跟哪個(gè)坐標(biāo)面重合就用哪個(gè)。提示柱坐標(biāo)和球坐標(biāo)的散度、旋度公式不要死背記住「面積因子和體積因子」的來(lái)路現(xiàn)場(chǎng)推比背錯(cuò)強(qiáng)。2.3 用 Python 驗(yàn)證散度與旋度的幾何意義光看公式?jīng)]感覺(jué)可以用數(shù)值方式把散度和旋度「畫(huà)」出來(lái)。下面這段代碼在直角坐標(biāo)系下構(gòu)造一個(gè)簡(jiǎn)單矢量場(chǎng)計(jì)算散度和旋度并驗(yàn)證「散度為零」和「旋度為零」分別對(duì)應(yīng)什么幾何形態(tài)。import numpy as np # 構(gòu)造網(wǎng)格 x, y, z np.meshgrid(np.linspace(-2, 2, 21), np.linspace(-2, 2, 21), np.linspace(-2, 2, 21), indexingij) h x[1,0,0] - x[0,0,0] # 網(wǎng)格步長(zhǎng) # 場(chǎng)1徑向發(fā)散場(chǎng) A (x, y, z)散度應(yīng)為3旋度應(yīng)為0 Ax, Ay, Az x, y, z # 數(shù)值散度中心差分 div_A (np.gradient(Ax, h, axis0) np.gradient(Ay, h, axis1) np.gradient(Az, h, axis2)) # 數(shù)值旋度分量 curl_x np.gradient(Az, h, axis1) - np.gradient(Ay, h, axis2) curl_y np.gradient(Ax, h, axis2) - np.gradient(Az, h, axis0) curl_z np.gradient(Ay, h, axis0) - np.gradient(Ax, h, axis1) print(場(chǎng)1 散度均值:, div_A.mean()) # 接近3 print(場(chǎng)1 旋度模均值:, np.sqrt(curl_x**2curl_y**2curl_z**2).mean()) # 接近0 # 場(chǎng)2旋轉(zhuǎn)場(chǎng) A (-y, x, 0)散度應(yīng)為0旋度z分量應(yīng)為2 Ax2, Ay2, Az2 -y, x, np.zeros_like(z) div_A2 (np.gradient(Ax2, h, axis0) np.gradient(Ay2, h, axis1) np.gradient(Az2, h, axis2)) curl_z2 np.gradient(Ay2, h, axis0) - np.gradient(Ax2, h, axis1) print(場(chǎng)2 散度均值:, div_A2.mean()) # 接近0 print(場(chǎng)2 旋度z均值:, curl_z2.mean()) # 接近2這段代碼的關(guān)鍵在于用中心差分替代解析偏導(dǎo)網(wǎng)格步長(zhǎng) h 越小越接近理論值。運(yùn)行后你會(huì)看到徑向場(chǎng)的散度接近 3旋度接近 0旋轉(zhuǎn)場(chǎng)的散度接近 0旋度 z 分量接近 2。這就是散度和旋度最直白的幾何含義——一個(gè)管「發(fā)不發(fā)散」一個(gè)管「轉(zhuǎn)不轉(zhuǎn)」。參數(shù)說(shuō)明網(wǎng)格范圍 -2 到 2 是經(jīng)驗(yàn)值太小看不出差分誤差太大邊界效應(yīng)明顯21 個(gè)點(diǎn)保證步長(zhǎng) 0.2精度夠用又不至于太慢。如果你把 h 改成 0.05散度會(huì)更接近 3但計(jì)算量上升。3. 把課件里的公式變成可復(fù)現(xiàn)的推導(dǎo)從梯度到麥克斯韋3.1 梯度、散度、旋度的組合恒等式怎么記矢量分析里有一組恒等式課件上通常列一頁(yè)但真正在電磁場(chǎng)里反復(fù)用的是兩條?×(?φ)0 和 ?·(?×A)0。前者說(shuō)「梯度的旋度為零」對(duì)應(yīng)靜電場(chǎng)可以寫(xiě)成電位梯度后者說(shuō)「旋度的散度為零」對(duì)應(yīng)磁場(chǎng)可以寫(xiě)成矢量磁位旋度。記憶方法不是背而是理解梯度場(chǎng)沒(méi)有渦旋旋度場(chǎng)沒(méi)有源。這兩條恒等式直接決定了電磁場(chǎng)里哪些量可以引入勢(shì)函數(shù)。靜電場(chǎng)引入電位 φ因?yàn)?E-?φ磁場(chǎng)引入矢量磁位 A因?yàn)?B?×A。沒(méi)有這兩條恒等式勢(shì)函數(shù)根本沒(méi)法定義。我一般會(huì)讓學(xué)生自己用行列式展開(kāi)驗(yàn)證一遍比看十遍結(jié)論管用。3.2 用 SymPy 做符號(hào)推導(dǎo)驗(yàn)證 ?×(?φ)0數(shù)值驗(yàn)證只能看趨勢(shì)符號(hào)推導(dǎo)才能確認(rèn)恒等式。下面用 SymPy 定義一般標(biāo)量函數(shù)直接算梯度的旋度看是否恒為零。import sympy as sp x, y, z sp.symbols(x y z) phi sp.Function(phi)(x, y, z) # 任意標(biāo)量函數(shù) # 梯度 grad_phi [sp.diff(phi, x), sp.diff(phi, y), sp.diff(phi, z)] # 旋度?×F (?Fz/?y - ?Fy/?z, ?Fx/?z - ?Fz/?x, ?Fy/?x - ?Fx/?y) curl_grad [ sp.diff(grad_phi[2], y) - sp.diff(grad_phi[1], z), sp.diff(grad_phi[0], z) - sp.diff(grad_phi[2], x), sp.diff(grad_phi[1], x) - sp.diff(grad_phi[0], y) ] # 化簡(jiǎn) for i, comp in enumerate(curl_grad): print(f分量{i}:, sp.simplify(comp))運(yùn)行結(jié)果三個(gè)分量全為 0。這不是巧合而是偏導(dǎo)交換次序的結(jié)果。SymPy 的 simplify 會(huì)幫你把混合偏導(dǎo)合并。參數(shù)上phi 定義成一般函數(shù)而不是具體表達(dá)式才能證明恒等式對(duì)任意標(biāo)量場(chǎng)成立。同樣的方法可以驗(yàn)證 ?·(?×A)0把 A 定義成三個(gè)一般函數(shù)即可。這兩個(gè)驗(yàn)證做完你對(duì)勢(shì)函數(shù)的來(lái)路就清楚了。3.3 從矢量恒等式到麥克斯韋方程組的推導(dǎo)路徑課件里通常直接給麥克斯韋方程組但如果你知道矢量分析可以自己推出來(lái)。路徑是這樣的從法拉第電磁感應(yīng)定律的實(shí)驗(yàn)事實(shí)出發(fā)寫(xiě)成積分形式再用斯托克斯定理把線積分轉(zhuǎn)成面積分得到 ?×E -?B/?t。安培環(huán)路定理加上位移電流修正用同樣方法得到 ?×H J ?D/?t。高斯定理用散度定理轉(zhuǎn)成微分形式得到 ?·D ρ 和 ?·B 0。這里用到的兩個(gè)積分定理——斯托克斯定理和高斯散度定理——本質(zhì)上是旋度和散度在宏觀尺度上的體現(xiàn)。斯托克斯定理說(shuō)「旋度的面積分等于環(huán)量」高斯定理說(shuō)「散度的體積分等于通量」。沒(méi)有矢量分析這兩個(gè)定理沒(méi)法表達(dá)麥克斯韋方程組只能停留在積分形式邊界條件也沒(méi)法推。所以一份好的矢量分析課件應(yīng)該把這條推導(dǎo)鏈講清楚而不是只列公式。你拿到課件后可以自己按這個(gè)路徑推一遍推完再看麥克斯韋方程組感覺(jué)完全不一樣。4. 避坑與排查矢量分析學(xué)習(xí)中最容易翻車的五個(gè)點(diǎn)4.1 現(xiàn)象柱坐標(biāo)和球坐標(biāo)散度公式總是記混原因柱坐標(biāo)散度里有一項(xiàng) (1/ρ)?(ρAρ)/?ρ球坐標(biāo)里有一項(xiàng) (1/r2)?(r2Ar)/?r形式相似但因子不同。死記硬背必然混。解決記住「面積因子」來(lái)路。柱坐標(biāo)中 ρ 方向的面積微元是 ρ dρ dφ所以散度里出現(xiàn) 1/ρ 和 ρ 的乘積球坐標(biāo)中 r 方向的面積微元是 r2 sinθ dθ dφ所以出現(xiàn) 1/r2 和 r2 的乘積。每次用之前花十秒推一下面積因子比背錯(cuò)強(qiáng)。4.2 現(xiàn)象數(shù)值計(jì)算旋度時(shí)結(jié)果和理論值差很多原因差分步長(zhǎng)選得太大或者邊界處理沒(méi)做。中心差分在邊界只能退化成單側(cè)差分誤差會(huì)累積。解決把網(wǎng)格步長(zhǎng) h 減小到 0.05 以下邊界點(diǎn)單獨(dú)用單側(cè)差分公式。如果還差很多檢查場(chǎng)函數(shù)定義有沒(méi)有寫(xiě)錯(cuò)分量順序。我見(jiàn)過(guò)最常見(jiàn)的錯(cuò)誤是把 curl_x 寫(xiě)成 ?Az/?y - ?Ay/?z符號(hào)搞反。4.3 現(xiàn)象?×(?φ)0 在數(shù)值上不為零原因數(shù)值差分本身有誤差離散旋度和離散梯度不滿足交換律。解決這是正?，F(xiàn)象不是代碼錯(cuò)。把步長(zhǎng)減小誤差會(huì)按 h2 下降。如果要做嚴(yán)格驗(yàn)證用 SymPy 符號(hào)推導(dǎo)不要用數(shù)值。4.4 現(xiàn)象分不清梯度、散度、旋度的輸入輸出類型原因三個(gè)算子名字太像容易混。解決記住輸入輸出類型。梯度標(biāo)量場(chǎng)進(jìn)矢量場(chǎng)出。散度矢量場(chǎng)進(jìn)標(biāo)量場(chǎng)出。旋度矢量場(chǎng)進(jìn)矢量場(chǎng)出。用一句話記梯度升維散度降維旋度保維。電磁場(chǎng)里電位是標(biāo)量電場(chǎng)是矢量E-?φ 就是梯度升維磁感應(yīng)強(qiáng)度是矢量磁通量是標(biāo)量Φ∫B·dS 就是散度降維的積分版。4.5 現(xiàn)象看課件覺(jué)得懂了做題就卡原因矢量分析是操作型技能不是知識(shí)型技能?？炊茖?dǎo)不等于會(huì)算。解決每學(xué)一個(gè)算子至少手算三個(gè)不同坐標(biāo)系的例子。直角坐標(biāo)算一遍柱坐標(biāo)算一遍球坐標(biāo)算一遍。算完再用 Python 數(shù)值驗(yàn)證。手算和數(shù)值對(duì)上了才算真會(huì)。5. 進(jìn)階技巧用矢量分析快速判斷場(chǎng)的形態(tài)學(xué)完基礎(chǔ)之后有一個(gè)很實(shí)用的技巧拿到一個(gè)矢量場(chǎng)表達(dá)式先算散度和旋度就能判斷它像什么場(chǎng)。散度為零、旋度不為零像磁場(chǎng)散度不為零、旋度為零像靜電場(chǎng)兩者都為零是調(diào)和場(chǎng)無(wú)源無(wú)旋兩者都不為零是一般時(shí)變場(chǎng)。這個(gè)判斷在仿真前處理里特別有用。比如你拿到一個(gè)天線近場(chǎng)數(shù)據(jù)先算散度和旋度分布就能知道哪些區(qū)域是輻射場(chǎng)、哪些是感應(yīng)場(chǎng)。輻射場(chǎng)區(qū)域旋度占主導(dǎo)感應(yīng)場(chǎng)區(qū)域散度不可忽略。再進(jìn)一步可以用矢量分析檢查邊界條件。理想導(dǎo)體表面切向電場(chǎng)為零法向磁場(chǎng)為零這些條件寫(xiě)成散度和旋度形式后可以直接在數(shù)值代碼里做投影。我一般會(huì)在 FDTD 或 FEM 代碼里加一段診斷每步算一下散度殘差如果殘差突然變大說(shuō)明網(wǎng)格或邊界出了問(wèn)題。最后說(shuō)個(gè)血淚經(jīng)驗(yàn)矢量分析不要等到用的時(shí)候再補(bǔ)它是電磁場(chǎng)所有后續(xù)內(nèi)容的語(yǔ)法。我見(jiàn)過(guò)太多人卡在麥克斯韋方程組回頭補(bǔ)矢量分析結(jié)果發(fā)現(xiàn)不是方程難是算子沒(méi)學(xué)透。如果你現(xiàn)在正在看這份課件建議把梯度、散度、旋度各手推三遍再用代碼驗(yàn)證一遍后面會(huì)省很多后悔藥。希望幫到你。本文還有配套的精品資源點(diǎn)擊獲取