一种极端天气微波链路畸变时空校准方法与流程

本发明涉及微波通信监测与气象观测，尤其涉及一种极端天气微波链路畸变时空校准方法。

背景技术：

1、随着城市通信网络和移动回传系统的快速建设，大量微波链路持续输出接收信号电平、载波相位、交叉极化鉴别度和多普勒频移等电磁特征。降雨、冰雹、冻雨及其复合天气会通过吸收、散射、去极化、多径传播和天线罩覆水或覆冰改变上述特征。因此，利用现有微波链路感知气象信息，能够在不大规模增设专用传感器的情况下获得路径积分观测，可用于短时强降水监测、灾害天气预警和通信链路运行状态评估。

2、现有微波链路气象反演方法，通常以晴空接收电平为参考，根据降雨比衰减与降雨率之间的幂函数关系计算降雨量，或采用经验回归模型对信号衰减进行补偿。该类方法在小雨、中雨及水凝物沿程分布较均匀时具有一定适用性，但在暴雨核、冰雹带、冻雨覆冰和雨雹混合等极端天气中，链路沿程介质呈现强非均匀性，不同电磁特征之间又存在非线性耦合，使得单一衰减特征和统一补偿模型难以保持校准精度。

3、第一，现有方法多采用固定晴空均值或简单滑动平均作为基线，难以将设备温漂、发射功率控制变化、天线老化等慢变分量与极端天气造成的快变畸变分离，基线偏移会使附加衰减、相位偏移和极化偏移被系统性高估或低估。

4、第二，微波链路、气象雷达、卫星遥感和地面气象站的采样周期、传输时延及空间支撑域并不一致。现有的最近时刻、最近站点或链路中点匹配方式，未考虑天气单体移动方向和速度所产生的平流时滞，在暴雨核快速移入或移出链路时容易造成时间与空间的同时错配。

5、第三，暴雨、冰雹和冻雨对微波信号的作用机理存在显著差异。暴雨主要引起吸收与散射衰减，冰雹还会产生非瑞利散射、去极化和短时闪烁，冻雨则会形成具有积累与融化迟滞特性的覆冰附加衰减。采用单一回归模型处理所有工况，难以在天气类型切换及复合事件中正确分配各物理机制的贡献。

6、第四，现有物理模型多以整条链路的平均衰减为输出，难以描述局部暴雨核或冰雹带沿链路的分段分布及传播关系；纯数据驱动模型虽可以拟合非线性残差，但在训练样本未覆盖的极端工况中，可能产生衰减为负、路径分段积分与链路总衰减不一致等违反传播规律的结果。

7、第五，现有校准方法通常只输出单一点估计，缺乏与观测完整率、多源匹配质量及模型外推程度相关联的不确定性评价。当雷达遮挡、地面站缺测或天气突变时，校准值仍会以确定结果输出，难以为气象预警和链路运维提供可靠的风险边界。

8、第六，现有在线校准通常采用统一学习率更新基线、传播系数和非线性校正参数。该方式既可能使快速变化的天气参数跟踪滞后，也可能使本应稳定的晴空基线和物理传播系数被单次异常事件改写，导致长期参数漂移。

9、因此，如何将电磁特征基线分离、异步多源观测时空匹配、不同极端天气传播机理、路径分段非线性校准、不确定性评价以及在线参数更新统一到同一闭环中，已成为提高极端天气条件下微波链路电磁特征校准精度、可解释性和工程可用性的关键技术问题。

技术实现思路

1、为了解决上述问题，本发明的目的在于提供一种极端天气微波链路畸变时空校准方法，通过采集微波链路收发电平、衰减量、相位及误码率等电磁特征数据，并融合雷达、雨量、温湿度、风速和覆冰等多源气象观测数据；采用快慢时间尺度分解构建能够随设备状态变化自适应更新的动态晴空基线，实现设备漂移与天气畸变的分离；根据气象单元移动速度、传播方向、链路位置和有效覆盖长度，建立传播时延感知的时空匹配模型，实现异步气象观测与链路传播路径的准确对应；利用极端天气物理衰减规律约束混合专家模型，对暴雨、冰雹、冻雨、大风及其复合天气进行识别，并确定不同天气因素对链路畸变的贡献程度；进一步将传播路径划分为多个具有空间关联关系的链路区段，通过路径图残差校准模型描述局部畸变在相邻区段间的传播、叠加和衰减过程，获得链路级畸变重构值；在此基础上，对校准结果进行不确定性量化，输出校准值及其可信区间，并在观测缺失或模型可信度不足时执行安全回退；最后采用变点触发的快慢参数协同更新机制，分别更新天气响应参数和设备基线参数，形成“数据采集-时空匹配-天气判型-路径重构-可信校准-在线更新”的闭环校准过程。本发明能够提高极端及复合天气条件下微波链路电磁特征校准的准确性、稳定性和可解释性，为链路状态监测、异常诊断及通信质量保障提供可靠依据。

2、为了实现上述目的，本发明提供的一种极端天气微波链路畸变时空校准方法是这样实现的：

3、一种极端天气微波链路畸变时空校准方法，包括多源数据同步采集模块、双时间尺度鲁棒基线分解模型、传播路径分段模型、气象观测单元标准化模型、时空非平衡最优传输对齐模型、物理约束混合专家模型、传播路径图残差校准模型、物理一致性约束模型、校准可信度评估模型、安全回退模型以及快慢参数在线更新模型；所述多源数据同步采集模块用于获取微波链路发射端和微波链路接收端输出的多频点和/或双极化电磁特征序列，并获取气象雷达、卫星遥感、地面气象站及链路位置数据，依据统一授时信号为各类数据写入采集时间戳，所述电磁特征至少包括接收信号电平、载波相位、极化角、交叉极化鉴别度和多普勒频移中的两种；所述双时间尺度鲁棒基线分解模型用于对经过质量控制的电磁特征序列进行分解，获得晴空基线分量、极端天气畸变分量、设备慢漂移分量和随机噪声分量，并以不随单个极端天气事件变化的晴空基线分量构建反事实晴空参考；所述传播路径分段模型用于根据微波链路收发端的三维坐标、天线挂高、数字高程模型、有效地球半径和大气折射参数重建链路传播路径，并按照传播环境变化和空间分辨率将所述传播路径划分为多个路径单元；所述气象观测单元标准化模型用于将不同来源、不同采样周期和不同空间分辨率的气象观测转换为同时具有时间支撑域、空间支撑域和质量置信度的观测单元；所述时空非平衡最优传输对齐模型用于根据路径单元与观测单元之间的采样时间差、天气单体平流时滞、三维空间距离、观测支撑域重叠度和物理特征差异建立联合代价矩阵，通过带边缘松弛项的熵正则非平衡最优传输求解路径单元与观测单元之间的耦合矩阵，并根据所述耦合矩阵将异步多源气象观测对齐至相应路径单元；所述物理约束混合专家模型用于接收对齐后的气象观测、去基线电磁畸变、多频差分特征、双极化差分特征和观测质量特征，通过门控单元确定不同极端天气类型的初始权重，并根据各天气专家输出结果相对于物理传播规律的违背量修正所述初始权重，再由暴雨专家、冰雹专家、冻雨专家和复合事件专家分别输出候选校准量并进行加权融合，以获得各路径单元的分段物理校准结果；所述传播路径图残差校准模型用于以路径单元为节点，并根据路径单元之间的相邻关系、空间距离、天气单体移动方向和传播时滞构建具有有向边和边权的传播路径图，对所述分段物理校准结果进行非线性残差修正，以表征局部电磁畸变在传播路径上的迁移、累积和衰减关系；所述物理一致性约束模型用于对各路径单元的校准结果施加路径积分守恒约束，并将修正结果投影至由电磁传播机理和气象参数取值范围共同确定的物理可行域，得到校准后的电磁特征及对应气象参数；所述校准可信度评估模型用于根据校准样本的非一致性分数计算有限样本分位数，为校准后的电磁特征或气象参数生成可信区间；所述安全回退模型用于在校准结果可信度低于设定阈值时抑制传播路径图残差校准模型的非线性残差输出，切换至基于物理传播关系的保守校准模式，并输出低可信度标记；所述快慢参数在线更新模型用于对校准残差进行在线变点检测，在变点统计量超过设定阈值时快速更新门控参数、图残差参数和时延参数，在变点统计量于稳定窗口内持续低于设定阈值时，以低于快参数更新学习率的学习率更新晴空基线参数和物理传播参数，从而形成极端天气快速响应参数与设备长期演化参数相分离的闭环时空校准机制。

4、本发明的双时间尺度鲁棒基线分解模型为：

5、对第条链路在时刻采集的电磁特征定义为

6、（1）

7、式（1）中，表示第条微波链路在时刻采集的观测电磁特征向量；表示与日周期、季节变化及正常大气折射条件相关的晴空基线分量；表示暴雨、冰雹、冻雨或复合极端天气导致的稀疏或分段稀疏畸变分量；表示设备温漂、发射功率控制变化或天线老化导致的慢漂移分量；表示随机噪声分量；表示微波链路索引，，表示参与校准的微波链路总数；表示离散采样时刻；上标表示观测量；黑体符号表示由不同频点、极化方式或电磁特征类型组成的向量。

8、将滑动窗口内的多链路观测组织为矩阵，质量掩码矩阵记为，通过下式求解晴空基线、天气畸变和设备漂移：

9、（2）

10、式（2）中，表示滑动窗口内由多条链路、多频点和多极化电磁特征组成的观测矩阵；、和分别表示待求解的晴空基线矩阵、极端天气畸变矩阵和设备慢漂移矩阵；表示与维数相同的质量掩码矩阵，有效观测位置取相应质量权重，缺测位置取零；表示对应元素相乘；表示frobenius范数；表示核范数，用于保持相邻链路及双极化通道的低秩共同基线；和分别表示沿时间方向的一阶差分算子和二阶差分算子；表示基线时间差分的范数，用于保持晴空基线缓慢变化；表示组稀疏范数，用于使同一极端天气事件在相关电磁通道上成组出现；表示设备慢漂移二阶差分的平方frobenius范数，用于限制慢漂移曲率；、、和分别表示低秩约束、基线时间平滑约束、天气组稀疏约束和设备漂移平滑约束的非负权重系数；表示在、和上使目标函数取得最小值。分解后，以作为若无当前极端天气时的反事实晴空参考，其中，为晴空基线分量的估计值，用于表示微波链路在不受极端天气影响时的正常电磁特征基准。为设备慢漂移分量的估计值，用于表示器件老化、温度变化或系统状态变化造成的缓慢偏移。以作为后续校准的天气畸变输入，其中表示第条微波链路在时刻的去基线电磁畸变向量；表示第条微波链路在时刻的观测电磁特征向量；表示由式（2）估计得到的晴空基线向量；表示由式（2）估计得到的设备慢漂移向量；表示微波链路索引；表示离散采样时刻；变量上方的尖帽符号表示相应参量的估计值。作为后续时空匹配、天气判型和路径分段校准的天气畸变输入。由此，当前天气事件不会被写入晴空参考，设备慢漂移也不会被误识别为气象畸变。

11、本发明的传播时延感知非平衡最优传输匹配模型为：

12、链路发射端和接收端的三维坐标分别为和。考虑有效地球半径及折射修正后，将链路划分为个路径单元，第个路径单元中心坐标为

13、（3）

14、式（3）中，表示第个路径单元中心的三维坐标向量；和分别表示微波链路发射端和接收端的三维坐标向量；表示路径单元索引，；表示传播路径划分得到的路径单元总数；表示第个路径单元中心在收发端直线路径上的归一化位置；表示第个路径单元的地球曲率修正向量；表示第个路径单元的大气折射修正向量。对第个多源气象观测单元，设其中心时刻为、中心坐标为、天气单体速度为。由观测单元到路径单元的平流时滞估计为：

15、（4）

16、式（4）中，表示第个气象观测单元所表征的天气单体平流至第个路径单元所需的估计时滞；表示待搜索的候选时滞；表示允许搜索的最大时滞绝对值；表示使目标函数取得最小值时对应的候选时滞；表示第个路径单元中心坐标；表示第个气象观测单元中心坐标；表示第个观测单元对应天气单体的三维平流速度向量；表示天气单体经过时长后的预测位置；表示欧氏范数；和分别表示路径单元索引和气象观测单元索引。

17、联合匹配代价定义为

18、（5）

19、式（5）中，表示第个路径单元与第个气象观测单元之间的联合匹配代价；、、、和分别表示时间差异项、三维空间距离项、支撑域不重叠项、物理特征差异项和观测质量项的非负权重；表示第个路径单元对应的链路电磁特征采样时刻；表示第个观测单元的中心时刻；表示由式（5）得到的平流时滞；和分别表示时间差及空间距离的尺度归一化参数；表示两坐标点之间的三维空间距离函数；、和分别表示路径单元中心坐标、观测单元中心坐标和天气单体平流速度；表示路径单元影响域与观测单元空间支撑域的交并比；表示第个路径单元的电磁畸变描述向量；表示第个观测单元的气象物理描述向量；表示物理特征差异的尺度或协方差矩阵；表示以加权的马氏距离平方；表示第个观测单元的质量置信度；表示防止对数自变量为零的正数；表示自然对数。由此，时间接近但空间不相交的观测，或者空间接近但已被天气单体移动错开的观测，均不会获得不合理的高权重。

20、设路径单元质量向量为，观测单元质量向量为，耦合矩阵通过下式获得：

21、（6）

22、式（7）中，表示最优耦合矩阵；表示待求解的非负耦合矩阵，其元素表示第个气象观测单元向第个路径单元分配的质量；表示由式（5）构成的联合代价矩阵；表示两矩阵对应元素乘积之和；表示熵正则系数；表示kullback-leibler散度；表示路径单元的归一化质量向量；表示气象观测单元的归一化质量向量；表示与的外积；上标表示矩阵转置；表示元素均为1且维数与相应矩阵乘法匹配的列向量；和分别表示耦合矩阵的行边缘质量和列边缘质量；和分别表示路径单元边缘松弛项和观测单元边缘松弛项的权重。后两个边缘松弛项允许存在雷达遮挡、地面站缺测及链路质量降级，而不强制错误的一一对应。对齐到路径单元的气象向量为:

23、（7）

24、式（7）中，表示对齐到第个路径单元的融合气象观测向量；表示最优耦合矩阵中第行第列的元素；表示第个气象观测单元的气象要素向量，所述气象要素至少包括降雨率、雷达反射率、温度、湿度、风速、冰雹概率或覆冰量中的一种；式中的求和运算表示对全部候选气象观测单元求和；表示防止分母为零的正数；和分别表示路径单元索引和气象观测单元索引；变量上方的波浪符号表示经时空耦合加权对齐后的参量。

25、本发明的物理约束混合专家校准模型为：

26、极端天气物理基线由路径单元上的比衰减积分得到。暴雨衰减为

27、（8）

28、式（8）中，表示暴雨在整条微波链路上产生的路径积分衰减；表示路径单元总数；表示微波链路工作频率；表示水平极化、垂直极化或其他设定极化方式；和分别表示由与共同确定的降雨比衰减系数和幂指数；表示第个路径单元的等效降雨率；表示第个路径单元的有效降雨覆盖长度，单位为km；式中的求和运算表示对第1个至第个路径单元的降雨衰减贡献求和。

29、冰雹等效散射衰减为

30、（9）

31、式（9）中，表示冰雹粒子散射和吸收造成的整链路等效衰减；表示冰雹粒子等效直径；和分别表示参与计算的冰雹粒子最小直径和最大直径；表示第个路径单元内粒径为的冰雹粒子数浓度谱分布；表示直径为的冰雹粒子在工作频率下的消光效率因子；表示冰雹粒子的等效几何截面积；表示对冰雹粒径进行积分的微元；表示第个路径单元的有效冰雹覆盖长度；式中的积分运算表示对至粒径范围内的消光贡献积分，式中的求和运算表示对第1个至第个路径单元求和。

32、冻雨覆冰附加衰减为

33、（10）

34、式（10）中，表示时刻冻雨及覆冰引起的附加衰减，单位为db；表示由工作频率、极化方式和环境温度决定的路径覆冰单位长度衰减系数；表示时刻的链路有效覆冰长度；表示工作频率下天线罩单位覆冰厚度的附加衰减系数；表示时刻天线罩的等效覆冰厚度；表示路径或天线罩附近的环境温度。

35、将经时空非平衡最优传输对齐模型对齐至第i个路径单元的多源气象观测，与该路径单元对应的去基线电磁畸变、多频差分特征、双极化差分特征和观测质量特征进行组合，构成门控输入zi。第k个天气专家的物理违背量为:

36、（11）

37、式（12）中，表示第个天气专家对第个路径单元输出结果的物理违背量；表示天气专家索引；表示物理约束项索引；表示物理约束项总数；表示第类天气对应的第项物理约束函数，所述物理约束至少包括气象参数取值范围、衰减非负、极化响应和路径积分一致性中的两项；表示第个天气专家对第个路径单元输出的候选校准向量；表示仅保留大于零的约束违背部分；表示路径衰减一致性项的非负权重系数；表示由链路观测分配或换算到第个路径单元的观测衰减；表示第个天气专家依据相应物理传播模型计算的第个路径单元物理衰减；表示绝对值。门控权重为

38、（12）

39、式（12）中，表示第个路径单元属于或采用第个天气专家的归一化门控权重；表示门控单元根据输入向量计算的第个天气专家的原始得分；表示由对齐气象观测、去基线电磁畸变、多频差分特征、双极化差分特征和观测质量特征组成的第个路径单元门控输入向量；表示物理违背惩罚系数；表示由式（11）计算的物理违背量；表示控制权重分布平滑程度的门控温度参数；表示天气专家总数；表示分母求和所用的天气专家索引；表示自然指数函数；分母用于使同一路径单元各天气专家权重之和等于1。候选校准量由各专家输出，融合结果为

40、（13）

41、式（13）中，表示物理约束混合专家模型对第个路径单元输出的融合候选校准向量；上标表示混合专家融合结果；表示第个天气专家在第个路径单元上的门控权重；表示第个天气专家的候选校准映射函数；表示第个路径单元的门控及专家输入向量；表示第个天气专家的待训练参数集合；表示天气专家索引；表示天气专家总数；式中的求和运算表示对第1个至第个天气专家的候选输出进行加权求和；变量上方的尖帽符号表示估计值。

42、当暴雨与冰雹同时存在时，暴雨专家与冰雹专家可同时获得非零权重；物理违背量会抑制与观测不相容的专家，从而避免强制单标签导致的复合事件混淆。

43、本发明的路径分段图残差校准模型为：

44、以路径单元为节点构建图。节点i的初始状态为。从节点到节点i的有向边权为

45、（14）

46、式（14）中，表示传播路径图中由节点指向节点i的归一化有向边权；表示节点i与节点对应路径单元中心之间的空间距离；表示空间距离衰减尺度；表示天气单体由节点对应位置传播至节点对应位置的时滞；表示传播时滞衰减尺度；表示天气移动方向一致性项的权重；表示由节点指向节点的连线方向角；表示天气单体移动方向角；表示节点连线方向与天气移动方向的一致程度；表示节点的入邻居节点集合；表示对邻居节点求和时使用的索引；和分别表示节点与邻居节点之间的空间距离和传播时滞；表示自然指数函数；分母用于使节点的全部入边权之和等于1。第层消息更新为

47、（15）

48、式（15）中，表示第层消息传递后节点的状态向量；和分别表示第层节点和其邻居节点的状态向量；表示图消息传递层索引；表示非线性激活函数；表示第层作用于节点自身状态的权重矩阵；表示第层作用于邻居消息的权重矩阵；表示由式（15）得到的节点指向节点的有向边权；表示节点的入邻居节点集合；表示第层偏置向量；式中的求和运算表示对节点的全部入邻居消息求和。

49、路径图模型只学习物理基线未能解释的残差，最终校准量为

50、（16）

51、式（16）中，表示第个路径单元最终校准后的电磁特征或气象参数向量；表示由式（13）得到的混合专家物理候选校准向量；表示路径图残差模型对第个路径单元输出的非线性残差修正向量；表示由衰减非负、气象参数范围、极化响应及路径积分关系共同确定的物理可行域；表示将括号内向量投影至物理可行域的投影算子；表示路径单元索引。

52、模型训练目标为

53、（17）

54、式（17）中，表示路径图残差校准模型的总训练损失；表示由校准输出与共址地面站、雷达定量产品或通信设备自检参考值之间差异构成的数据拟合损失；、和分别表示路径积分一致性项、图总变分平滑项和物理约束项的非负权重；表示整条微波链路的实际观测总衰减；表示第个路径单元校准后的比衰减估计值；表示第个路径单元的有效长度；表示欧氏范数；表示传播路径图的有向边集合；表示连接节点与节点的一条边；表示相应有向边权；和分别表示节点和节点的校准输出向量；表示范数；表示作用于节点校准输出的第项物理约束；表示物理约束项索引；表示仅惩罚大于零的约束违背部分。式中第二项保证路径分段校准结果积分后与链路实际总衰减一致，第三项在天气单体连续区域内抑制非物理振荡，第四项限制输出落入物理可行域。

55、本发明的可信区间生成与安全回退模型为：

56、在独立校准样本上定义第个样本的非一致性分数

57、（18）

58、式（18）中，表示第个独立校准样本的非一致性分数；表示独立校准样本索引；表示第个样本的参考真值或经质量审核的参考量；表示第个样本的模型校准估计值；表示校准估计值与参考量之间的绝对误差；表示模型对第个样本给出的标准差或尺度不确定性估计值；表示防止分母为零并表征不可约噪声下限的正数；变量上方的尖帽符号表示估计值。

59、取有限样本修正后的分位数

60、（19）

61、式（19）中，表示置信水平为时经有限样本修正得到的非一致性分数分位数；表示目标失覆盖率，且；表示目标可信区间覆盖率；表示独立校准样本总数；表示由式（18）计算得到的全部独立校准样本非一致性分数；表示样本集合的第分位数；表示向上取整运算；表示考虑有限样本数量后的修正分位水平。

62、则实时样本的可信区间为

63、（20）

64、式（20）中，表示实时样本在时刻的可信区间；表示时刻的校准估计值；表示由式（19）获得的有限样本修正分位数；表示模型对时刻输出的标准差或尺度不确定性估计值；表示不可约噪声下限；区间左端点和右端点分别表示校准结果的可信下界和可信上界；表示实时采样时刻。定义可信度，其中，表示时刻的综合可信度；表示可信区间宽度；表示被校准电磁特征或气象参数的有效量程；表示防止分母为零的正数；表示由观测完整率、最优传输耦合质量及传感器质量共同确定的时刻数据质量系数；表示自然指数函数；表示安全回退触发阈值。当时，将图残差系数逐步衰减至零，并采用由式（8）至式（10）得到的物理基线输出，以防止低质量数据驱动不受约束的外推。

65、本发明的快慢参数分离在线更新模型为：

66、设实时残差为，在长度为的参考窗口上计算变点统计量

67、（21）

68、式（21）中，表示时刻的在线变点统计量；表示截至时刻的最近个采样周期内校准残差向量的样本均值；表示时刻的实时残差向量，满足；表示时刻的参考电磁特征或参考气象参数向量；表示时刻的校准估计向量；表示变点检测参考窗口长度；表示稳定状态下残差向量的基准均值；表示稳定状态下残差向量的基准协方差矩阵；表示所述基准协方差矩阵的逆矩阵；上标表示向量转置；表示实时采样时刻。越大，表示当前窗口残差分布偏离稳定状态越明显。

69、当时，判定天气类型、移动速度或数据时延发生突变，冻结慢参数，按较大学习率更新快参数：

70、（22）

71、式（22）中，和分别表示第次更新前和第次更新后的快参数向量；所述快参数至少包括混合专家门控参数、路径图残差参数和传播时延参数；表示快参数学习率；表示时刻根据实时校准样本计算的损失函数；表示损失函数关于快参数向量的梯度；表示在线参数更新次数。当式（21）得到的大于变点阈值时执行式（22），其中，表示判定天气类型、移动速度或数据时延发生突变的统计阈值。

72、当持续个采样周期后，再更新慢参数：

73、（23）

74、式（23）中，和分别表示第次更新前和第次更新后的慢参数向量；所述慢参数至少包括晴空基线参数和物理传播参数；表示慢参数学习率；表示式（22）中的快参数学习率，且；表示对慢参数向量求梯度；表示时刻的实时校准损失；表示历史参数保持项的非负权重；表示慢参数分量索引；表示稳定历史事件对第个慢参数的重要度；表示当前慢参数向量中的第个分量；表示所述第个慢参数在当前稳定更新阶段开始前保存的历史值；式中的求和运算表示对全部慢参数分量求和；表示允许慢参数更新前，变点统计量持续不高于阈值所需的稳定采样周期数。由此，在突发极端天气中快速跟踪门控、图残差和传播时滞，同时避免传播系数及晴空基线被单个事件破坏。

75、由于本发明采用双时间尺度鲁棒基线分解、时延感知时空匹配、物理约束混合专家、传播路径图残差校准、可信区间与安全回退以及快慢参数分离更新等技术手段，构建微波链路电磁特征畸变的闭环时空校准机制，从而具有以下有益效果：

76、1.本发明通过双时间尺度鲁棒基线分解，将极端天气畸变与设备慢漂移相分离；结合平流时滞、三维空间距离、支撑域重叠度和观测质量，实现异步多源气象观测与链路路径单元的准确匹配；进一步利用物理约束混合专家和传播路径图残差模型，分别表征不同极端天气及复合事件的传播特性和路径分段畸变规律，从而提高极端天气条件下电磁特征校准的准确性、物理一致性和稳定性。

77、2.本发明通过有限样本可信区间对校准结果的不确定性进行量化，并在观测缺失、模型外推或可信度不足时自动切换至物理基线保守校准模式；同时采用变点触发的快慢参数分离更新机制，使天气响应参数快速更新、晴空基线参数和物理传播参数缓慢更新，从而兼顾突发天气跟踪能力与长期运行稳定性，并可将校准结果用于链路信号恢复、路径气象参数反演和低可信度告警。