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