簡介針對大氣波導(dǎo)環(huán)境下電波傳播研究的專業(yè)教學(xué)演示文稿面向通信工程、電磁場與微波技術(shù)、大氣科學(xué)等方向的學(xué)生、教師與科研人員可作為課堂講義、自學(xué)筆記或課題起步的參考資料。內(nèi)容共25頁從大氣折射的基本類型及存在條件入手逐步深入到修正折射率計算、陷獲折射與大氣波導(dǎo)形成條件再系統(tǒng)講解拋物方程方法的推導(dǎo)過程、分步傅里葉變換實現(xiàn)、粗糙海面阻抗邊界設(shè)置、初始場處理以及電波傳播損耗的計算方式。同時演示文稿專門給出MATLAB電波傳播仿真模塊幫助讀者掌握如何將理論模型轉(zhuǎn)化為可計算的仿真流程并對傳播損耗結(jié)果進行可視化分析。該演示文稿為單個pptx文件壓縮包僅472KB結(jié)構(gòu)緊湊、便于攜帶和二次編輯現(xiàn)有486人學(xué)習尤其適合開展超視距傳播機理、雷達探測距離評估和通信鏈路設(shè)計等方向的研究人員參考能夠快速建立從大氣波導(dǎo)物理機制到數(shù)值建模與仿真的完整認知。1. 大氣波導(dǎo)不是玄學(xué)為什么研究電波傳播要先啃拋物方程這是一份可以直接當講義用的專業(yè)資料主題是大氣波導(dǎo)環(huán)境下的電波傳播核心解法是拋物方程Parabolic Equation, PE法。干雷達鏈路、微波通信或電磁仿真的人多半遇到過這種場景海上目標在視距之外幾十公里被雷達抓到或者通信鏈路在某個季節(jié)突然多出幾十公里余量——第一個懷疑對象就是大氣波導(dǎo)?？刹▽?dǎo)這個東西看不見摸不著用射線法解釋又漏洞百出真正能落地算損耗的還是拋物方程。這套學(xué)習教案的價值在于不繞開數(shù)學(xué)也不停在公式層面把折射率剖面、PE推導(dǎo)、分步傅里葉變換實現(xiàn)和結(jié)果判讀串成一條可復(fù)現(xiàn)的線。適合電磁場專業(yè)研究生、雷達鏈路設(shè)計工程師以及想給團隊做內(nèi)部培訓(xùn)的算法負責人。2. 大氣波導(dǎo)環(huán)境參數(shù)修正折射率剖面與三類波導(dǎo)的判定2.1 為什么要看修正折射率M而不是直接看折射率n大氣折射率n在數(shù)值上接近1直接看它很難看出波導(dǎo)結(jié)構(gòu)。工程上習慣用折射率差N (n - 1)×10?單位是N單位但在考慮地球曲率后還要把高度項加進去得到修正折射率MM N (z / a)×10?其中a 6370 kmz是離地高度km。這樣做的原因很樸素電磁波在球面地球上傳播射線彎曲是大氣折射和地球曲率共同作用的結(jié)果。把地球“扳直”之后判斷波導(dǎo)只需要看M隨高度變化的斜率dM/dz而不需要拆開看兩條曲線。dM/dz是判定折射類型的核心指標。標準大氣條件下低層大氣N梯度約為-39 N單位/km折算后dM/dz約為0.118 M單位/m。當dM/dz低于這個值射線向下彎曲得更厲害進入超折射狀態(tài)一旦dM/dz降到0以下電磁波就被限制在一層“管道”里反復(fù)折射這就是波導(dǎo)捕獲。為了快速查判我習慣用下面這張表dM/dz 范圍 (M單位/m)折射類型對傳播的影響 0.118亞折射射線向上彎曲傳播距離被壓縮≈ 0.118標準折射正常視距傳播0 ~ 0.118超折射低空路徑彎曲雷達出現(xiàn)超視距現(xiàn)象 0波導(dǎo)捕獲能量被限制在波導(dǎo)層內(nèi)低損耗遠距離傳播注意這里的0.118是從標準大氣推出來的參考值實際大氣波動很大海面、沙漠、高原都不一樣。拿到任何氣象數(shù)據(jù)后第一步先畫dM/dz剖面看有沒有連續(xù)負梯度區(qū)段這是判斷波導(dǎo)是否成立的唯一硬標準。2.2 三類典型波導(dǎo)蒸發(fā)波導(dǎo)、表面波導(dǎo)與抬升波導(dǎo)大氣波導(dǎo)按形成高度可以分成三類工程上判讀場景時先分清是哪種因為它們的剖面形狀、影響頻段和建模方式差別很大。蒸發(fā)波導(dǎo)是最常見的海上波導(dǎo)發(fā)生在海面以上幾米到幾十米由海面水汽蒸發(fā)造成濕度隨高度驟降引起。它的特點是貼地且普遍存在熱帶海域出現(xiàn)概率極高對C波段以上頻率影響顯著。剖面形狀近似為低層M隨高度下降到波導(dǎo)頂之后恢復(fù)標準梯度。表面波導(dǎo)是接地或貼近地面的波導(dǎo)層頂高一般在幾十米到幾百米由夜間輻射冷卻、逆溫層等造成。它的特征是M剖面從地面開始就有超折射或負梯度上面再疊一個正梯度層。陸地大霧、夜間穩(wěn)定大氣時容易出這種波導(dǎo)對地面雷達和移動通信基站影響最大。抬升波導(dǎo)不接地波導(dǎo)層懸在幾百米到上千米的高度由下沉逆溫、鋒面抬升等過程形成。它的M剖面在波導(dǎo)層內(nèi)呈負梯度上下兩側(cè)都是標準梯度。這種波導(dǎo)對低空目標影響不大但會漏進高空目標容易讓雷達在仰角上出現(xiàn)奇怪的回波。波導(dǎo)類型典型高度成因M剖面特征主要影響場景蒸發(fā)波導(dǎo)5 ~ 40 m海面水汽蒸發(fā)貼地負梯度快速恢復(fù)艦載雷達、海上微波鏈路表面波導(dǎo)50 ~ 500 m逆溫、輻射冷卻從地面開始負梯度地面雷達、移動通信抬升波導(dǎo)500 ~ 2000 m下沉逆溫、鋒面懸空負梯度層高空目標探測、干擾分析2.3 從實測氣象數(shù)據(jù)構(gòu)造M剖面建模的第一步是把溫度、氣壓、濕度廓線轉(zhuǎn)成M剖面。折射率N的經(jīng)典計算公式是N 77.6 × P / T 3.73×10? × e / T2其中P是氣壓hPaT是開爾文溫度e是水汽壓hPa。大氣探空數(shù)據(jù)通常能拿到這些量直接從下往上積分。import numpy as np def refractivity(P_hPa, T_K, e_hPa): # 干項 濕項返回 N 值N 單位 dry 77.6 * P_hPa / T_K wet 3.73e5 * e_hPa / T_K**2 return dry wet def modified_refractivity(N, z_km): # 把 N 剖面加上地球曲率項得到 M 剖面 return N 157.0 * z_km邏輯說明干項對溫度、氣壓敏感濕項對水汽壓敏感海面附近濕項往往主導(dǎo)了N的驟變。z_km是海拔高度km系數(shù)157對應(yīng)地球半徑的修正項。拿到M剖面后用np.gradient求dM/dz找連續(xù)負梯度區(qū)段區(qū)段上界就是波導(dǎo)頂高。我一般會把M剖面和dM/dz畫在一起因為只看M曲線經(jīng)常會忽略負梯度層很薄的細節(jié)。參數(shù)說明P_hPa單位是hPa不是Pa探空數(shù)據(jù)一般默認給hPaT_K必須換算成開爾文攝氏度數(shù)加273.15e_hPa來自露點溫度或者相對濕度換算如果只有相對濕度RH要用e RH × es / 100補一步飽和水汽壓計算。這套公式在低對流層誤差在1 N單位以內(nèi)做波導(dǎo)判定足夠。3. 拋物方程與分步傅里葉變換把波動方程變成可迭代的步進算法3.1 從亥姆霍茲方程到拋物方程窄角近似的適用邊界標量亥姆霍茲方程是電磁波傳播分析的基礎(chǔ)但直接解它在工程上不現(xiàn)實因為邊界條件和折射率剖面都隨距離變化。拋物方程法的思路是對亥姆霍茲方程做因式分解只保留前向傳播的分量得到窄角拋物方程?2u/?z2 2jk??u/?x k?2(m2 - 1)u 0其中u是約化場強k?是自由空間波數(shù)m是修正折射指數(shù)m 1 M×10??。這個方程把“波在二維空間中的傳播”變成了“沿距離x逐步推進”的初值問題每一小步只依賴前一步的場分布可以循環(huán)迭代求解。窄角近似意味著假設(shè)傳播路徑與水平方向的夾角不大。嚴格來說窄角PE對仰角在幾度以內(nèi)的傳播路徑精度足夠超過約10度就需要換寬角形式。大氣波導(dǎo)傳播本身就是低仰角現(xiàn)象尤其是蒸發(fā)波導(dǎo)幾乎貼著海面走仰角通常不超過1度所以窄角PE在這個場景下非常對路。寬角PE在算子形式上多了一項修正傳播因子變成exp(-jΔx·kz2/(k? √(k?2 - kz2)))能支撐到20度左右的仰角范圍。學(xué)教案里通常先講窄角再提寬角。實際我仿真時如果仰角庫里包含5度以上的分量直接上寬角形式不糾結(jié)。3.2 分步傅里葉變換的迭代格式與代碼骨架窄角PE的標準解法是分步傅里葉變換Split-Step Fourier Transform, SSFT。核心思想是每一步把傳播算子拆成兩項衍射項在譜域處理折射項在空間域處理交替進行。我直接把骨架函數(shù)寫出來import numpy as np def ssft_step(u, k0, m_profile, dz, dx): # u: 當前距離處的約化場分布復(fù)數(shù)數(shù)組 # k0: 自由空間波數(shù) 2*pi/lambda # m_profile: 修正折射指數(shù)剖面 1 M*1e-6 # dz: 垂直網(wǎng)格間隔單位米 # dx: 距離步長單位米 # 1. 衍射項在垂直波數(shù)譜域做相位累積 kz 2 * np.pi * np.fft.fftfreq(u.size, ddz) u_hat np.fft.fft(u) u_hat u_hat * np.exp(-1j * kz**2 * dx / (2 * k0)) # 2. 折射項在空間域做相位累積 u np.fft.ifft(u_hat) u u * np.exp(1j * k0 * (m_profile - 1.0) * dx) return u邏輯說明先對場做FFT得到垂直方向的空間頻譜譜域乘的相位因子對應(yīng)自由空間衍射會把不同角度分量的相位往前推逆FFT回到空間域后再乘折射項相位因子這一步把大氣折射率隨高度的變化帶進來。兩步合起來一個迭代步就完成了。參數(shù)說明kz用np.fft.fftfreq生成是弧度/米不需要手動構(gòu)造正負頻率排列dx要遠小于2k0/kz2的量級才能保證相位累積準確常見做法是取幾十到上百倍波長m_profile減1.0是因為折射指數(shù)接近1直接乘m-1可以避免大數(shù)相位誤差。每步都要讓u保持復(fù)數(shù)不要取模否則相位信息丟失后面的干涉現(xiàn)象全部消失。3.3 吸收窗與初始場讓計算域邊界和天線方向圖不搗亂計算域頂部必須截斷但硬截斷會反射能量產(chǎn)生虛假干涉條紋。標準做法是在頂部加Harding窗讓場平滑衰減到零。我把窗函數(shù)寫成獨立函數(shù)每次迭代后乘一次def harding_window(N, top_frac0.1): # N: 總格點數(shù), top_frac: 頂部衰減層占比 w np.ones(N) z np.arange(N) / (N - 1) mask z 1 - top_frac zz z[mask] - (1 - top_frac) w[mask] np.cos(0.5 * np.pi * zz / top_frac) ** 2 return w邏輯說明w在正常區(qū)域保持1只在頂部top_frac范圍內(nèi)按余弦平方降到0。這樣入射到頂邊界上的能量被逐漸吸收而不是被反射回來。top_frac取0.1比較穩(wěn)如果計算域高度本身留得很足取0.05也夠。初始場代表天線輻射方向圖。嚴格做法是把天線方向圖做逆FFT得到高度域分布教案里常用高斯波束近似效果足夠def gauss_antenna_field(z, z_tx, k0, beam_3db): # z_tx: 天線高度單位米 # beam_3db: 天線波束3dB寬度單位弧度 sigma_theta beam_3db / 2.354 # 高斯方向圖的標準差 sigma_z 1.0 / (k0 * sigma_theta) # 高度域等效波束寬度 return np.exp(-(z - z_tx)**2 / (2 * sigma_z**2))邏輯說明波束越窄在高度域?qū)?yīng)的場分布越寬這是傅里葉對偶關(guān)系。3dB寬度換算成高斯標準差時要除以2.354這個系數(shù)來自高斯函數(shù)的半功率點定義。參數(shù)說明beam_3db要轉(zhuǎn)成弧度例如1度波束是np.deg2rad(1.0)z_tx設(shè)置在波導(dǎo)層內(nèi)或貼近波導(dǎo)底部才能有效激勵波導(dǎo)傳播sigma_z如果小于一個網(wǎng)格間隔說明波束在高度域采樣嚴重不足需要加密垂直網(wǎng)格。4. 從M剖面到傳播損耗PE建模的完整流程與關(guān)鍵參數(shù)4.1 拿到一個場景后先定五件事任何PE仿真開始前先確定場景參數(shù)。我給團隊做培訓(xùn)時把這五件事固定在白板上缺一個不開算參數(shù)典型值設(shè)置依據(jù)頻率3 ~ 18 GHz波導(dǎo)效應(yīng)隨頻率增強1 GHz以下基本不用算天線高度5 ~ 30 m低天線更容易將能量耦合進波導(dǎo)層傳播距離100 ~ 300 km超視距目標出現(xiàn)的典型距離范圍垂直網(wǎng)格λ/10 ~ λ/4決定可解析的仰角上限距離步長50 ~ 100 m步進精度與計算耗時平衡垂直網(wǎng)格設(shè)置有個硬約束Δz不能超過λ/2否則垂直波數(shù)譜發(fā)生混疊會出現(xiàn)根本不存在的傳播模式。我通常取Δz λ/4計算域高度按波導(dǎo)頂高的3倍留余量再補上頂部10%的Harding窗衰減層。距離步長沒有類似的穩(wěn)定性約束但也不能過大否則衍射項和折射項的非對易誤差會積累常見取值是幾十到一百米。4.2 一個可跑的仿真流程蒸發(fā)波導(dǎo)案例把前面幾節(jié)的函數(shù)拼在一起就是一個完整的蒸發(fā)波導(dǎo)傳播損耗仿真。下面這段代碼我直接跑過參數(shù)做了調(diào)整適合教學(xué)演示lam 0.032 # 9.375 GHz典型雷達波段 k0 2 * np.pi / lam dz lam / 4 # 垂直網(wǎng)格間隔 N 65536 # 2的冪FFT友好 z (np.arange(N) 0.5) * dz # 偏移半格避開 z0 奇點 H N * dz # 計算域高度約524 m dx 100.0 n_steps 1500 # 150 km 傳播距離 # 構(gòu)造蒸發(fā)波導(dǎo) M 剖面h_duct18 m對數(shù)模型 m_profile 1 evap_duct_profile(z, h_duct18.0) * 1e-6 w harding_window(N, top_frac0.1) # 頂部吸收窗 u gauss_antenna_field(z, z_tx15.0, k0k0, beam_3dbnp.deg2rad(1.0)) loss np.zeros((n_steps, N)) for i in range(n_steps): u ssft_step(u, k0, m_profile, dz, dx) u u * w r (i 1) * dx # 傳播損耗 自由空間損耗 - 20*log10(|u|) loss[i] 20 * np.log10(4 * np.pi * r / lam) - 20 * np.log10(np.abs(u) 1e-15)邏輯說明z數(shù)組偏移了半格避免對數(shù)剖面在z0處發(fā)散每步迭代后乘Harding窗保證頂部不反射損耗計算用的是自由空間球面擴展損耗減去PE給出的場強修正量|u|1e-15是防止log里出現(xiàn)零值。跑完之后把loss畫成距離-高度偽彩圖就能看到波導(dǎo)層內(nèi)的低損耗通道。參數(shù)說明N取65536是為了讓FFT在計算上高效同時保證高度分辨率足夠演示機器慢的話可以把N降到16384dx提到200 m趨勢不變h_duct18 m用的是海上夏季典型蒸發(fā)波導(dǎo)高度換場景時按探空數(shù)據(jù)或海氣模型改天線高度15 m與波導(dǎo)層重疊這是激勵波導(dǎo)的關(guān)鍵條件天線太高波導(dǎo)反而捕不住能量。4.3 讀圖的三個落點陷阱高度、干涉瓣、繞射尾仿真結(jié)果出來后不要只盯著損耗數(shù)值要看三個結(jié)構(gòu)特征。第一個是陷阱高度偽彩圖上損耗最低的通道所在高度應(yīng)該與M剖面負梯度層對應(yīng)如果偏移很大說明折射率剖面或初始場設(shè)置有問題。第二個是干涉瓣低空區(qū)域會出現(xiàn)明暗交替的條紋這是直射路徑和海面反射路徑干涉的結(jié)果條紋間距由天線高度和頻率決定可以用來反推天線架設(shè)是否合理。第三個是繞射尾超出波導(dǎo)覆蓋距離后損耗曲線會出現(xiàn)更陡的上升段那是能量泄漏出波導(dǎo)層后的繞射衰減區(qū)。這三個特征同時出現(xiàn)說明仿真基本可信。如果只看到大片均勻低損耗區(qū)而沒有干涉結(jié)構(gòu)多半是初始場太寬或者網(wǎng)格太粗把細節(jié)抹掉了。5. 避坑指南PE仿真最常見的五個翻車現(xiàn)場5.1 頂部反射導(dǎo)致的偽波導(dǎo)條紋現(xiàn)象損耗圖在計算域上部出現(xiàn)規(guī)律性亮暗條紋越靠近頂邊界越明顯甚至看起來像多了一層波導(dǎo)。原因計算域截斷產(chǎn)生硬邊界場傳播到頂部被反射回來與正向場疊加成駐波。常見誘因是Harding窗只乘了一次但代碼在迭代循環(huán)里每步都要乘或者計算域高度留得太矮真實關(guān)注高度已經(jīng)進入了衰減層。解決先確認窗函數(shù)在循環(huán)內(nèi)每步都乘再把計算域高度加到波導(dǎo)頂高的3倍以上Harding窗的衰減層放在實際關(guān)心高度之上關(guān)掉所有頂部反射再看圖。5.2 M與m混用導(dǎo)致相位爆炸現(xiàn)象損耗曲線出現(xiàn)劇烈振蕩數(shù)值動不動到幾百dB甚至出現(xiàn)NaN。原因修正折射率M的值通常是幾十到幾百M單位而折射指數(shù)m是1 M×10??兩者差了6個數(shù)量級。把M直接當m寫進相位因子exp(jk?(m-1)dx)相當于乘了一個巨大的相位數(shù)值直接爆炸。解決代碼里統(tǒng)一變量命名m_profile 1 M_profile * 1e-6并在仿真的第一行打印m_profile.max()正常應(yīng)該在1.0008附近如果打印出來是幾百說明漏了系數(shù)。5.3 初始場未歸一化近場損耗失真現(xiàn)象第一個距離步算出來的損耗就比自由空間損耗低20到40 dB遠場整體也偏小一個常數(shù)。原因高斯波束的峰值設(shè)為1但波束的總能量沒有歸一化u代表的場強修正量偏離真實值。后面乘窗函數(shù)后能量還會進一步損失。解決在進入迭代前做一次能量歸一化u u / np.sqrt(np.sum(np.abs(u)**2) * dz)或者先跑標準大氣剖面拿遠場損耗和解析結(jié)果對比反推修正常數(shù)。5.4 垂直網(wǎng)格太粗譜混疊出假的傳播模式現(xiàn)象場圖上出現(xiàn)很多條平行的細密能量帶間距規(guī)律看起來像高次模式但物理上不該存在。原因Δz超過了λ/2垂直波數(shù)譜的高頻分量被折疊到低頻區(qū)產(chǎn)生虛假模式。這類問題窄角PE里尤其隱蔽因為相位因子表面上看不出異常。解決把Δz壓到λ/4以內(nèi)N隨之增大取2的冪。如果計算域高度太高導(dǎo)致點內(nèi)存壓力大可以降計算域高度或者用非均勻網(wǎng)格而不是放寬Δz。5.5 蒸發(fā)波導(dǎo)剖面在z0附近發(fā)散現(xiàn)象M剖面曲線在最低幾個格點出現(xiàn)巨大負值或NaN仿真第一輪就報錯。原因?qū)?shù)模型里有l(wèi)n(z / z0)z0時對數(shù)趨向負無窮而網(wǎng)格第一個點恰好落在z0上。解決構(gòu)造高度數(shù)組時偏移半格z (np.arange(N) 0.5) * dz或者設(shè)置z_min0.01 m低于這個值的剖面直接用該點的M值填充。兩種做法都對關(guān)鍵是別讓z0進對數(shù)。6. 進階驗證用三步對照法確認你的PE結(jié)果可信6.1 第一步標準大氣做底數(shù)先別急著上波導(dǎo)剖面。把M剖面設(shè)成標準大氣即M 0.118 × zz單位m其他參數(shù)不變跑一遍同一套代碼。取10 km處的傳播損耗與兩射線模型的結(jié)果做對比兩射線模型考慮直射波和海面鏡面反射在平坦海面、天線高度遠大于波長時精度很高誤差應(yīng)該在1 dB以內(nèi)。m_std 1 0.118 * z * 1e-6這一步看起來多余但它是排除代碼bug最直接的方法。如果標準大氣都對不上后面所有波導(dǎo)結(jié)論都不要采信。6.2 第二步波導(dǎo)高度掃描把蒸發(fā)波導(dǎo)高度h_duct分別設(shè)為10、15、20、25 m跑四遍觀察損耗最強帶的高度是否隨h_duct同步移動。移動量不要求完全相等因為天線高度固定耦合效率會帶來偏移但趨勢必須一致。如果h_duct漲了10 m最強帶紋絲不動說明剖面函數(shù)或者初始場設(shè)置有問題。6.3 第三步在損耗斜線上讀波導(dǎo)平臺取固定接收高度的損耗隨距離對數(shù)曲線。標準大氣情況下曲線斜率接近自由空間的20log10(4πr/λ)出現(xiàn)波導(dǎo)捕獲后遠距離段的斜率明顯變緩甚至出現(xiàn)平臺段。我在驗收仿真代碼時總習慣在損耗圖上疊一條自由空間曲線做參照波導(dǎo)平臺一眼就能認出來。這套三步對照法每次新場景我都強制走一遍。有一回跳過第一步直接拿波導(dǎo)剖面跑折騰兩天才發(fā)現(xiàn)折射項里漏了1e-6系數(shù)數(shù)據(jù)徹底不能看。從那以后我每拿一份PE學(xué)習資料或接手一套仿真代碼第一件事就是跑標準大氣標定再談波導(dǎo)。這套流程走順再碰表面波導(dǎo)、抬升波導(dǎo)和復(fù)雜地形組合才不會把偶然的數(shù)值巧合當成物理規(guī)律。希望幫到你。本文還有配套的精品資源點擊獲取