简介：基于ITU-R P.676-11建议书开发的MATLAB复现代码，面向研究大气气体对无线电波衰减的科研人员、通信工程师及电磁场专业学生。运行main函数即可直接输出完整衰减曲线，覆盖大气吸收损耗、毫米波传播衰减等典型场景；所有源码行均带注释，便于理解算法流程和调节频率、温湿度等参数。压缩包共6个文件，包括5个.m脚本与1个txt说明，整体仅3KB，轻量精简；源码按主程序加gfun、gamo、gamw等子函数模块组织，结构清晰，既适合新手逐行研读，也方便有经验者直接调用函数嵌入自有仿真框架。txt说明对运行环境和文件用途做了简要交代，能帮助减少上手成本。目前已有417人学习下载，适合无线通信与电磁场方向入门到进阶的实践者。借助这套代码，可对照ITU-R P.676-11建议书中的气体衰减模型快速验证计算结果，并为卫星通信、雷达链路等工程中的大气衰减评估提供基础算法支持与可扩展参考。

1. 毫米波大气衰减：一张曲线卡住多少链路预算

做毫米波通信或雷达链路的人，迟早会撞上同一道坎：电磁波在大气里跑一公里，被氧气和水汽吃掉多少dB。频率上了10 GHz以后，这不是可以忽略的背景噪声，而是直接决定发射功率、天线口径和可用带宽的硬约束。P.676-11 建议书就是用来算这笔账的标准算法——它把氧气和水汽几十条吸收谱线、连续吸收项、温度和气压修正全部整理成可编程的数学形式。工程上常说的“大气吸收损耗”“电磁波的大气衰减”“毫米波大气衰减曲线”，基本都是用这套方法算出来的。这套MATLAB复现代码就是围绕该建议书实现的可运行版本，输入频率、气压、温度和水汽密度，输出1 GHz到1000 GHz全程的比衰减曲线。适合正在做星地链路、无线回传、遥感辐射计选频，或者对“22 GHz为何是水汽谷、60 GHz为何是氧气峰”追根问底的人。

2. P.676-11 建议书的核心：谱线参数表与两套衰减算法，先把物理模型对齐

2.1 为什么衰减曲线是一堆吸收峰的叠加：氧和水汽的谱线结构

大气气体对电磁波的衰减不是单调递增的。氧分子在50 GHz到70 GHz之间有一整簇磁偶极子吸收谱线，中心大致在60 GHz附近，再加上118.75 GHz那条独立谱线；水汽分子则有22.235 GHz、183.31 GHz、325.15 GHz、380.20 GHz等多个转动吸收峰。每一个峰都对应某一个具体的量子跃迁，衰减曲线本质上是这些洛伦兹线型叠加的结果。

所以复现代码的第一步，不是直接背公式，而是搞清楚谱线表长什么样。氧气谱线有三十几条，水汽谱线在P.676-11里也列了几十条，每条谱线都有中心频率、谱线强度、线宽参数，而线宽参数又要随温度和气压做二次修正。常见误区是把衰减当成频率的简单多项式，或者只拿“每公里衰减随频率升高单调变大”来估算。真实曲线在22 GHz、60 GHz、183 GHz附近会出现明显的驼峰，在峰与峰之间反而有相对低损耗的窗口。做毫米波频段规划的人，恰恰要利用这些低谷窗口，而不是避开所有峰。

2.2 分离谱线法 vs 近似法：精度、适用范围和参数边界

P.676-11 给出了两种计算路径。第一种是分离谱线法，把每一条吸收谱线单独计算贡献再求和，适用频率范围覆盖1 GHz到1000 GHz，在峰值附近表现好，但计算量比较大。第二种是近似法，用一组多项式拟合替代逐条谱线求和，适合快速估算和低频段使用，但在60 GHz氧气带和几个强水汽峰附近误差会明显放大。

工程选型时，我的习惯是这样：如果只是做30 GHz以下频段的宏观覆盖估算，近似法足够了；如果要评估E频段（71到86 GHz）、120 GHz或者更高频段的点对点链路，必须上分离谱线法。两套方法用到的气象输入完全一样，只是内部对谱线的处理方式不同。复现的资源包里通常会把两套方法都写进去，方便做交叉比照。下表是两套算法的典型差异：

2.3 复现前必须确认的气象学输入量：温度、气压、水汽密度如何进入公式

衰减公式的输入不是单一的“天气状况”，而是三个物理量：温度T（单位K，也有版本用摄氏度）、气压P（单位hPa或百帕）、水汽密度ρ（单位g/m³）。三者共同决定分子碰撞频率和谱线展宽程度。水汽密度尤其关键，它直接控制22 GHz、183 GHz这些水汽峰的强度；空气湿度越大，衰减峰越高，而氧气峰几乎不受湿度影响。

如果你手里只有相对湿度RH，需要先把RH换算成绝对湿度。常见做法是先查当前温度下的饱和水汽压，再用RH乘以饱和水汽密度得到ρ。这一步经常成为复现时数值对不上的主要来源——有人直接把RH当ρ丢进公式，结果22 GHz峰值小了一个数量级。建议在代码入口统一转化：

% 由相对湿度 RH(%) 和温度 T(K) 估算水汽密度 rho(g/m^3) % 采用常见近似公式，适合 0~40°C 范围 T_C = T_K - 273.15; es = 6.112 * exp(17.62 * T_C / (243.12 + T_C)); % 饱和水汽压 hPa rho = RH * es / (461.5 * T_K) * 100; % 由气体方程换算密度

逻辑说明：先求出当前温度下饱和水汽压es，再利用水汽气体常数 Rs=461.5 J/(kg·K)，把分压力换算成绝对湿度。乘以100是处理 hPa 与 Pa 的单位差。很多误差其实是这一步埋下的，所以先把输入量烫平再谈算法。

3. 用MATLAB把大气损耗曲线跑出来：模块拆解、参数设置与绘制输出

3.1 主程序框架：扫频范围、气象场景与输出存根

复现的第一件事，是把建议书里的数学公式拆成几个独立的MATLAB函数，而不是把几百行公式堆在一个脚本里。我习惯分成三层：主脚本负责扫频和组织输出，中间层负责调用氧气和水汽两个谱线模块，最底层是单条谱线的线型计算。下面是一个可用的主程序骨架：

%% main_p676_demo.m % P.676-11 大气气体衰减复现，扫频得到比衰减曲线 f_list = [1:0.25:40, 40.5:0.05:80, 80.5:0.25:200, 200.1:0.5:400]; % GHz % 频率分布设计：10GHz 以下粗扫，60GHz 氧带细扫，高频段中粒度 P_hPa = 1013.25; % 气压 hPa，海平面标准大气 T_K = 288.15; % 温度 K，对应 15°C rho = 7.5; % 水汽密度 g/m^3，中纬度夏季典型值 gamma = zeros(size(f_list)); for idx = 1:numel(f_list) gamma(idx) = atten_gas(f_list(idx), P_hPa, T_K, rho); end

逻辑说明：f_list 不是等间距的，因为60 GHz附近氧气谱线极密，步长太大会把峰切掉；低频段谱线稀少，可以放宽步长。循环内调用的 atten_gas 是统一入口，它内部会把氧气和水汽的贡献相加。参数说明里最关键的是单位——频率一律用GHz，气压用hPa，否则公式里的常数对不上。初次跑的时候建议先跑频率上限到400 GHz，确认无报错再扩展到1000 GHz。

3.2 氧吸收谱线模块：谱线表的组织方式与求和循环

氧气谱线表是复现里最容易出错的地方。P.676-11 中氧气有三十多条谱线，每条包含中心频率、谱线强度等参数。MATLAB代码里一般把谱线表存成矩阵，行是谱线，列是参数。下面是简化的代码结构，实际谱线数据要按建议书完整录入：

function att_o = atten_oxygen(f_GHz, P_hPa, T_K) % 氧气谱线表占位：行 = [中心频率GHz, 谱线强度参数, 线宽参数...] % 实际使用时替换为 P.676-11 Table 中的完整数据 o2_lines = [ 50.474, 0.94, 1.0; 55.400, 0.63, 1.0; 59.591, 2.30, 1.0; 60.306, 2.30, 1.0; 61.150, 2.30, 1.0; 62.384, 2.30, 1.0; 64.678, 1.45, 1.0; 66.600, 0.94, 1.0; 68.500, 0.63, 1.0; 118.75, 1.80, 1.0; ]; theta = 300 / T_K; % 温度比倒数的标准变形 att_o = 0; for row = 1:size(o2_lines, 1) f0 = o2_lines(row, 1); % 单线贡献：压力与温度修正的洛伦兹线型 df = 0.9 * P_hPa / 1013.25 * (300 / T_K)^0.7; % 简化线宽 S_i = o2_lines(row, 2) * P_hPa / 1013.25 * theta^2.5; % 简化强度修正 delta = f_GHz - f0; line_shape = (df / (delta^2 + df^2)) / pi; att_o = att_o + S_i * line_shape * f_GHz / f0; end att_o = att_o * f_GHz * 1e-3; % 粗略单位校准，正式版按建议书公式走 end

逻辑说明：这段代码是结构示范，不是对建议书的逐字翻译。重点在于每条谱线的贡献由强度、线宽、频率偏离度共同决定，线宽随气压升高而变宽、随温度升高而变窄，强度则与气压成正比。参数说明里，theta 是 300/T_K，这是毫米波谱线公式里反复出现的标准量；如果你把 T_K 的单位误写成摄氏度，theta 会变成负数，衰减直接翻车。单位校准系数要根据P.676-11的具体公式重算，这里只是先把结构跑通。

3.3 水汽吸收谱线模块与连续吸收项处理

水汽衰减和氧气最大的不同，在于除了22 GHz、183 GHz这些离散谱线外，还有一个不可忽略的连续吸收项。这个连续项来自水汽分子碰撞导致的远翼吸收，它会填平谱线之间的低谷，尤其对20 GHz到40 GHz的平坦段影响明显。忽略连续项是复现结果“看起来峰对了、谷不对”的头号原因。

function att_w = atten_wv(f_GHz, P_hPa, T_K, rho) % 水汽谱线头部参数表（简化示意） wv_lines = [ 22.235, 0.11; 67.000, 0.02; 119.995, 0.02; 183.310, 0.23; 325.153, 1.50; 380.197, 0.35; 439.000, 0.10; ]; att_w = 0; for row = 1:size(wv_lines, 1) f0 = wv_lines(row, 1); S_i = wv_lines(row, 2) * rho * (300 / T_K)^3.5; % 水汽强度与湿度成正比 df = 0.8 * P_hPa / 1013.25 * (300 / T_K)^0.8; delta = f_GHz - f0; line_shape = (df / (delta^2 + df^2)) / pi; att_w = att_w + S_i * line_shape * f_GHz / f0; end % 连续吸收项：与水汽密度平方和频率平方近似成正比 att_c = 2.5e-4 * rho^2 * (f_GHz / 22)^2 / (P_hPa / 1013.25); att_w = att_w + att_c; end

逻辑说明：水汽谱线强度随湿度线性增加，上一段代码里的 S_i 乘了 rho；连续项则更复杂，它近似与 rho 的平方和频率平方相关，气压升高时还会被压缩。参数说明里那个 2.5e-4 是先验估算值，严谨复现时要以建议书原文里的连续吸收系数公式为准。这里保留它是因为很多简化版本把连续项漏掉，导致183 GHz以上高频段曲线整体偏低。

3.4 结果物化：绘制标准曲线、导出CSV，为链路预算提供输入

跑出数值只是第一步，工程上还要把曲线画出来、把数据导出去给链路仿真用。P.676-11 的标准画法通常用对数纵轴，因为从1 GHz到400 GHz，衰减值跨越了三个数量级：低频段不到0.1 dB/km，60 GHz峰能到15 dB/km上下。如果还用线性纵轴，低频段会被压成一条直线。

figure('Color', 'w'); semilogy(f_list, gamma, 'b-', 'LineWidth', 1.5); xlabel('Frequency (GHz)'); ylabel('Specific attenuation (dB/km)'); title('P.676-11 Atmospheric Gas Attenuation (15C, 7.5g/m3)'); grid on; xlim([1, 400]); % 导出 CSV，供外部链路仿真读取 out_tbl = table(f_list', gamma', ... 'VariableNames', {'freq_GHz', 'gamma_dB_per_km'}); writetable(out_tbl, 'p676_attenuation.csv');

逻辑说明：semilogy 是对数纵轴绘图，能把22 GHz水汽峰、60 GHz氧带、183 GHz次峰同时显示清楚。table 和 writetable 的组合会把结果写成两列CSV，第一列频率、第二列比衰减。参数说明里的CSV文件名可以自己改，但要留意频率单位—其他仿真软件读入时往往默认Hz，如果直接用GHz数据会导致链路预算偏差几十dB，这也是常見的交接事故点。

4. 复现时最容易翻车的五个坑：数值偏差、图表错位与性能卡顿排查

4.1 现象：22 GHz水汽峰高度比参考曲线低约30%

原因：水汽密度单位被误用。有人把相对湿度RH（比如50%）直接当作水汽密度，而P.676-11公式需要绝对湿度g/m³。在海平面15°C时，50% RH对应的绝对湿度大约6.5 g/m³，如果直接把50塞进去，强度被放大了近8倍；但更多人是把 ρ 写成kg/m³，数值又缩小了1000倍。解决：在入口统一做单位换算，并用“干空气ρ=0时水汽衰减必须归零”做自检。

4.2 现象：60 GHz氧气带出现双峰错位或峰值凹陷

原因：扫频步长太大或者谱线表合并错误。氧气谱线在59到62 GHz之间极为密集，相邻谱线间距只有几十MHz量级，如果用1 GHz步长扫描，峰值正好被采样点错过。解决：在50到70 GHz区间加密步长到0.05 GHz甚至更细；同时检查谱线表中有没有把相邻两条谱线误删或合并成一条。60 GHz带对谱线完整性极其敏感，差一条线峰值形状都会变。

4.3 现象：曲线在低频段出现负衰减

原因：线型函数中的分母出现零值或者单位混乱。当频率f远小于第一条谱线中心频率时，f/f0项很小，如果线宽被设成负值或者温度单位被误填成摄氏度，计算结果会反号。另一个常见原因是MATLAB里指数或平方运算溢出。解决：先固定一条谱线单测，检查f=1 GHz时氧衰减是否为正值，再逐步打开其他谱线。所有线宽和强度的温度修正项都要以建议书原文为准，不要用简化公式里的占位常数硬算。

4.4 现象：与旧版本或网络上的参考值对不上

原因：P.676-11 相比旧版本对水汽连续吸收系数和部分氧气谱线线宽公式做过调整。如果拿老版本的建议书编号下的参数表混用，22 GHz谷值和80到120 GHz段的数值都会偏。解决：先确认参数表版本，再按建议书附录给出的算例对点校准。最稳妥的做法是对照官方算例里的3个校验频点（典型的是22.5 GHz、60 GHz、183.3 GHz），偏差控制到1%以内再往下走。

4.5 现象：扫频点一多，脚本卡到怀疑人生

原因：双层循环逐点逐线求和，还没做向量化。MATLAB对纯标量循环不友好，如果频率点2000个、谱线60条、三层嵌套，单次运行可以卡到几分钟。解决：把谱线表做成矩阵，用bsxfun或隐式扩展一次性算出所有频率点的线型贡献，再用sum按频率维求和。下面的优化片段把主循环消掉，时间能降两个数量级：

% 向量化谱线求和示例 [f_mesh, f0_mesh] = meshgrid(f_list, o2_lines(:,1)); df_mesh = 0.9 * P_hPa / 1013.25 * (300/T_K)^0.7; delta_mesh = f_mesh - f0_mesh; line_mesh = (df_mesh ./ (delta_mesh.^2 + df_mesh.^2)) / pi; att_o = sum(o2_lines(:,2) .* line_mesh .* (f_mesh ./ f0_mesh), 1);

逻辑说明：meshgrid 把频率列表和谱线中心频率展开成二维矩阵，每一行对应一条谱线、每一列对应一个频率点；line_mesh 一次性算出整个频率网格的线型值。sum 按第一维求和，得到整个频率列表上的氧气贡献。参数说明里，o2_lines(:,2) 是谱线强度列向量，MATLAB直接用它和矩阵做点乘，隐式扩展把列向量广播到每个频率点。这一步跑完后，主循环就可以彻底删掉了。

5. 快速验证你的曲线是否可信：官方算例匹配与物理自洽性检查

5.1 第一步：对官方算例，而不是凭感觉看曲线形状

曲线形状整体差不多，不代表数值准确。P.676-11 附录中给出了一些频点和气象条件下的比衰减算例，复现后第一件事是把这些点逐一对齐。做这种事别用“肉眼看柱状图差不多”来蒙混，直接把MATLAB的计算结果与文档值做差，打印成表：

check_freq = [22.5, 60.0, 183.3]; % GHz 校验频点 check_ref = [0.18, 14.8, 10.7]; % 参考值量级，具体以文档为准 calc_val = arrayfun(@(f) atten_gas(f, 1013.25, 288.15, 7.5), check_freq); disp(table(check_freq', calc_val', check_ref', ... calc_val' ./ check_ref' * 100, ... 'VariableNames', {'fGHz','Calc','Reference','RatioPercent'}));

逻辑说明：arrayfun 对每个校验频点调用一次核心函数，得到比衰减；ratioPercent 列显示计算值与参考值的百分比。这项检查的价值在于把“大概吻合”变成“可接受偏差”。参数说明里 check_ref 的数值是量级示意，实际复现时应该从建议书附录中读取，不同气象条件的参考值不同。我一般要求偏差小于1%，达不到就先排查谱线表录入和单位，而不是去调拟合系数。

5.2 第二步：干空气归零试验与极端气象检查

物理自洽性是最快的排错手段。把水汽密度设成0，水汽模块输出必须为0，总衰减应等于纯氧气衰减；把水汽密度设成50 g/m³，22 GHz峰不能出现负值或震荡。极端检查比日常频点更早暴露公式里的隐藏bug。

我习惯跑三个场景：干空气、高湿度、低温低压。低温低压模拟高空环境，气压降到500 hPa时，氧气峰应该明显变矮变窄，因为谱线展宽变小、强度降低。如果高空曲线反而比地面高，说明温度修正符号写反了。这类问题靠对照算例表格是查不出来的，因为算例只覆盖地面标准气象。

5.3 第三步：与近似公式和在线参考背靠背比较

分离谱线法和近似法在同一次复现里可以互相校验。在低频段，两者偏差应该在10%以内；在60 GHz处方波应该拉大。如果低频段偏差也很大，大概率是分离谱线法里连续项算错了。做这一步不需要额外下载数据，直接在代码里调用两个函数对比即可。

6. 把这张曲线用起来：频率选带、链路预算接入与批量扫描技巧

6.1 从比衰减到链路损耗：一个简单的EIRP算例

比衰减单位是dB/km，乘上链路距离就是气体带来的附加损耗。做地面点对点链路时，我习惯把这条曲线和雨衰曲线放一起看。气体衰减是“晴天也有”的底噪，雨衰是“下雨才慌”的增量。比如12 km的点对点链路，在22 GHz附近气体衰减按0.2 dB/km量级算，总附加损耗约2.4 dB；在60 GHz氧带中心按15 dB/km量级算，总损耗高达180 dB，基本宣告链路不可用。这个对比能直接解释为什么E频段通信通常避开60 GHz。

6.2 批量扫描具体技巧：锁定低损耗频窗而不是只看峰值

高频应用选频不能只看衰减曲线的谷值点，还要考虑可用带宽和邻近峰的影响。我通常写一个批量脚本，把1到100 GHz按200 MHz步长扫一遍，找出连续2 GHz窗口内平均衰减最小的频段作为候选。注意不要只看单点最低值——单点谷值往往夹在两个峰之间，带宽不够复用。

% 滑动平均找低损耗频窗 win_size = 200; % 频点数对应 2 GHz avg_gamma = movmean(gamma, win_size); [min_avg, min_idx] = min(avg_gamma); best_freq = f_list(min_idx); fprintf('最佳中心频点: %.2f GHz\n', best_freq);

逻辑说明：movmean 对衰减曲线做滑动平均，窗口宽度对应实际需要的带宽；min 找出平均衰减最小的中心频点。参数说明里 win_size 与频率步长相关，如果步长是0.01 GHz，200个点就是2 GHz。这个技巧能避免被人眼峰值误导。

6.3 海拔与气候变化场景的快速摸排

同一条链路，夏天和冬天气体衰减能差出几个dB，原因主要是水汽密度大幅变化。我习惯在脚本里加一组场景循环：海平面标准、低湿度高原、高湿热带，每个场景输出一条曲线，并打印选频窗口的变化。这一组结果直接放进链路预算文档，比单独一条标准曲线有说服力得多。

这份复现代码里，最值钱的不是那条标准曲线本身，而是你能把曲线拆开、按自己的场景重新拼装的能力。从那以后，我每次拿到任何一款大气衰减模型，都强制走一遍“对算例、干空气归零、低频与近似法交叉验证”这三板斧，再也不敢直接拿网上的参数表硬算。希望帮到你。

本文还有配套的精品资源，点击获取