简介：这是一份结合数据驱动与模型驱动方法的机械工程专业资料，面向具备机械工程或热力学背景的研究人员、工程师，尤其适合从事主轴轴承系统热特性分析与设计优化的专业人士。文档基于论文方法，完整实现了混合驱动框架，涵盖热网络建模、SIAN与Sobol全局灵敏度分析、粒子滤波与稀疏识别结合的多层粒子滤波（MLPF）算法，能够帮助读者在仅有少量温度传感器条件下预测不同工况的温度场。包内为1个docx文件，共64KB，包含详细的理论推导、Python可运行代码及逐步解释，并讨论了热源计算、热阻建模、热-机耦合效应等关键技术点。读者可通过实际代码复现混合驱动框架，掌握参数敏感性分析与多层粒子滤波的落地实现思路，为开发更精准的热特性预测工具、优化主轴轴承设计提供直接参考。该资源已有66人学习下载，内容紧凑、实践性强，适合用于课题研究、算法验证及工程应用对照。
1. 混合驱动框架：主轴轴承热分析里“模型算不准、数据喂不够”的折中方案
主轴轴承的温度场预测，在机床热误差补偿和轴承故障诊断里一直是绕不过去的坎。纯模型方法把热网络搭得很细，但接触热阻、油膜摩擦热这些参数受装配预紧力和热-变形耦合影响，出厂实测值常常跟试验对不上；纯数据驱动方法又依赖传感器密度，而轴承内部能做测点的位置只有那么几个，数据量根本不足以支撑高维回归。这篇论文提出的混合驱动框架，把热网络模型当骨架，用Sobol全局灵敏度分析筛出真正影响输出的不确定参数，再用粒子滤波加稀疏识别实时修正这些参数，最后只靠三五个温度测点重构整个温度场。适合做主轴系统热特性研究、机床热误差补偿或轴承状态监测的工程师和研究生——资源里附带的Python代码把建模、滤波、验证的主流程都串起来了，可以直接改着跑。
2. 先搭物理骨架：热网络模型与温度场仿真的可运行实现
2.1 为什么用热网络而不是有限元
主轴轴承的热分析核心要回答两个问题：热量从哪里来，热量往哪里去。摩擦热源集中在滚动体与滚道的接触区，热量从内圈经过滚动体传到外圈，再经轴承座散到壳体和环境中；轴端同时有润滑油带走一部分热，转速和载荷一变，热源和热阻都会跟着变。
有限元能把这个过程算得很精细，但代价是建模周期长、在线预测一次耗时几十秒到几分钟，而且接触热阻在有限元里依然只能靠经验值去填空。热网络方法则把主轴、内圈、滚动体、外圈、轴承座各离散成一个节点，节点之间用热阻连接，每个节点配一个热容方程，最后形成一组常微分方程。用scipy.integrate.odeint求解，单步计算只要几毫秒，完全跟得上工况变化，这是论文选择热网络作为模型驱动部分的关键原因。
另外一点实际工程的考量是，热网络的节点数和结构可以直接映射到可测的位置。论文里传感器放在哪几个节点，粒子滤波的观测函数就提取哪几个节点的温度，物理结构、模型节点、传感器位置三者一一对应，排查问题时会非常直观。
2.2 热平衡方程的代码落地与参数解释
资源里的ThermalNetworkModel实现了最基础的热平衡方程。下面这段是在原代码基础上补充了环境热阻初始化的可运行版本：
import numpy as np from scipy.integrate import odeint class ThermalNetworkModel: def __init__(self, nodes, resistances, heat_sources, R_env=None): self.nodes = nodes self.R = resistances self.Q = heat_sources self.T_ambient = 25.0 if R_env is None: self.R_env = np.full(nodes, 50.0) # 默认环境热阻，单位 K/W else: self.R_env = R_env def heat_balance(self, T, t): dTdt = np.zeros(self.nodes) for i in range(self.nodes): for j in range(self.nodes): if i != j and self.R[i, j] > 0: dTdt[i] += (T[j] - T[i]) / self.R[i, j] dTdt[i] += self.Q[i] dTdt[i] += (self.T_ambient - T[i]) / self.R_env[i] return dTdt def simulate(self, T0, t): return odeint(self.heat_balance, T0, t)逻辑说明：heat_balance是核心方法，它对每个节点写出能量守恒式：相邻节点温差除以热阻得到热流，热源项直接叠加，最后用(T_ambient - T[i]) / R_env[i]表示该节点与环境之间的散热。simulate把初始温度数组T0和时间点数组t交给odeint求解常微分方程组。注意我在__init__里补了R_env的默认值——原始样板代码里heat_balance用到了这个属性但构造函数没赋值，直接跑会报AttributeError，这是复现这个资源时第一个会遇到的问题。
参数说明：
nodes：节点数量。论文示例用 5 个节点，实际工程我建议按 内圈、滚动体、外圈、轴承座、轴端、壳体 拆 6~10 个节点。节点分得越细，热源位置描述越准，但待估计参数也越多，粒子滤波的负担会明显上升。刚开始跑通流程时别贪多，5 到 6 个节点足够验证算法。resistances：热阻矩阵，单位 K/W。物理上热阻必须满足互易性，即R[i][j]与R[j][i]相等。示例代码里先随机生成再(R + R.T) / 2对称化、对角线清零，这个顺序不能乱。heat_sources：热源向量，单位 W。示例里Q[1]=1.0、Q[2]=0.5表示两个主要发热位置。实际数值从摩擦热公式算出来填进去。R_env：节点到环境的散热热阻，默认给 50 K/W，实际根据壳体散热面积和风冷条件调整，数量级一般在 10~200 K/W。
我自己跑这个模型的经验是，先用所有热阻取同一个数量级做“等热阻测试”。如果稳态温度场出现某个节点明显不合理的高温，先检查热阻矩阵对称化和环境散热阻，而不是急着调热源大小。热源偏差只会让整体温度水平平移，热阻不对称却会在节点之间制造假温差，两种现象完全不同。
2.3 摩擦热源计算：从 Palmgren 公式到热-变形耦合修正
热网络模型里的Q向量是直接给定的，但实际工程中，轴承摩擦热源要由转速、载荷和润滑状态算出来。论文给了改进的BearingHeatSourceCalculator，用 Palmgren 公式计算基础摩擦扭矩，再叠加热-变形耦合修正：
class BearingHeatSourceCalculator: def __init__(self, bearing_params): self.params = bearing_params self.deformation_coeff = 0.02 # 热变形耦合系数，需根据实验标定 def calculate_heat_source(self, rpm, load, temp_diff): n = rpm / 60 # 转/秒 dm = self.params['pitch_diameter'] fr = load M0 = 10e-7 * self.params['f0'] * (n * dm)**3 if n > 0 else 0 M1 = self.params['f1'] * fr * dm torque = M0 + M1 heat = 2 * np.pi * n * torque # 摩擦热功率，单位 W heat *= (1 + self.deformation_coeff * temp_diff) # 热变形耦合修正 return heat逻辑说明：M0对应与转速和润滑剂粘度相关的粘性摩擦扭矩，M1对应与载荷相关的载荷摩擦扭矩，两者相加得到总摩擦扭矩。摩擦功率等于扭矩乘以角速度，即2 * np.pi * n * torque，作为热源功率注入热网络对应节点。最后一行是关键改进：温升导致轴承零件热膨胀、接触压力变化、摩擦热增加，形成正反馈，deformation_coeff * temp_diff就是这个反馈环的一阶近似。
参数说明：
f0、f1：轴承类型相关的摩擦系数，查轴承手册获取，示例中10e-7是量级参考。pitch_diameter：轴承节圆直径，单位米，这是热源计算里最基本的几何参数。deformation_coeff：热变形耦合系数，论文没有给通用值，我建议先在 0.01~0.05 范围做敏感性扫描，再根据稳态温度实验反推。
这个热源计算器与热网络模型的衔接方式是：每个时间步先用当前转速、载荷和上一时间步的温度差算出Q，再代入heat_balance求解下一步温度。这就是后面混合驱动框架里“控制输入是转速、模型预测受工况驱动”的直接体现。
3. 参数敏感性分析：用 Sobol 指标筛出“真正需要被估计”的未知参数
3.1 为什么要做全局灵敏度分析而不是逐个试参
热网络模型一搭起来就会面临“参数太多、数据太少”的矛盾：热阻矩阵里的每个非零元素、每个节点背后的热源系数、环境散热阻，加起来十几二十个参数，而温度测点往往只有两三个。把所有参数都丢给粒子滤波去估计，状态维度能冲到十几维，粒子数量呈指数级膨胀，收敛速度和精度都会崩掉。
论文的工程思路值得借鉴：先用 SIAN 结构可识别性分析加 Sobol 全局灵敏度指标做筛选，只把“对输出影响大、且具备可辨识性”的参数留作不确定参数，其余参数固定为标称值。Sobol 方法的核心是方差分解：把模型输出的总方差分解为每个参数单独变化的贡献以及参数交互作用的贡献。一阶敏感性指数S_i表示参数X_i单独变化时对输出方差的贡献占比；总敏感性指数ST_i额外统计该参数与其他参数交互作用的贡献。示例代码里给了一个简化版，用中位数分组近似条件方差。严格做法是用 Saltelli 采样方案生成两组样本矩阵来计算一阶与总阶指标。参数不多时两者排序结果通常一致，用简化版完全够用。
温度模型还有一个特点值得单独说明：参数间往往存在强相关性。比如内圈热阻和外圈热阻对测点温度的贡献方向相似时，即使 Sobol 指标都高，两者也无法同时被独立估计。SIAN 结构可识别性分析就是用来抓这种相关性的——它从模型结构上判断哪些参数组合能由观测数据唯一确定。两个指标配合使用的逻辑是：Sobol 回答“参数影响大不大”，SIAN 回答“参数能不能被唯一确定”，只有两个条件同时满足，这个参数才值得放进粒子滤波的状态向量。
3.2 基于蒙特卡洛采样的简化 Sobol 实现
def sobol_sensitivity_analysis(model, param_ranges, n_samples=1000): n_params = len(param_ranges) samples = np.zeros((n_samples, n_params)) for i in range(n_params): samples[:, i] = np.random.uniform( param_ranges[i][0], param_ranges[i][1], n_samples) outputs = np.array([model(sample) for sample in samples]) var_total = np.var(outputs) sensitivities = [] for i in range(n_params): # 按第 i 个参数的中位数把样本切成两半，比高半区方差占比 cond_var = np.var(outputs[samples[:, i] > np.median(samples[:, i])]) / var_total sensitivities.append(cond_var) return sensitivities逻辑说明： 第一步，在每个参数的取值范围内做均匀随机采样，得到n_samples组参数组合；第二步，逐组调用模型函数得到对应输出；第三步，对每个参数，把样本按该参数值的中位数切成上下两半，计算上半区输出方差占总方差的比例，比值越高说明输出对该参数越敏感。简化版的代价是灵敏度绝对值不可靠，但用于参数排序筛选完全够用。
参数说明：
param_ranges：每个参数的取值范围列表，例如[(0.01, 0.1), (1.0, 5.0)]。取值范围直接影响排序可靠性——范围过窄会把真正敏感的参数低估，过宽又会放大无关参数的波动。n_samples：采样数。5 个参数以内建议 1000~2000 组，模型本身是 ODE 求解，计算量能接受；如果换成有限元模型做同样的敏感性分析，就不现实了。model：一个函数，输入参数数组、返回标量输出。示例里用“第一个节点在 t=10s 时的温度”作为输出，你也可以换成稳态温度或温度变化速率。
实际使用中，我一般把敏感性计算结果排序后画个条形图，把排名前 3~5 的参数标记为“待估计”，其余的固定为标称值。这里特别注意：敏感性高不等于可辨识性好。如果两个参数的 Sobol 指标都很高、但物理上作用方向相同，粒子滤波会在这两个参数之间反复横跳，最终收敛到一组“看起来合理、实际是错误组合”的值。此时要结合 SIAN 或者干脆做一次参数相关性分析，把线性相关的参数合并掉一个。
3.3 敏感性结果如何指导粒子滤波的状态维度设计
粒子滤波的状态维度就是待估计参数的数量。把热阻矩阵所有非零元素都放进状态，维度可能冲上 15 维，1000 个粒子都未必收敛。做一轮 Sobol 筛选后，往往只把 3~5 个热阻或热源系数留作状态维，其余按名义值固定，粒子的搜索空间立刻被压缩到可管理的范围。
论文的逻辑到这里就很清晰了：热网络模型先按标称参数跑一个初始温度场，再把 Sobol 筛选出的不确定参数注入粒子群，让粒子在这些参数的小邻域内扩散搜索。传感器只需要覆盖敏感节点对应的物理位置——通常是内圈或外圈附近——就能把参数辨识出来。这种做法直接决定了后面“只用少量传感器就能预测整个温度场”这个结论能否成立，因为状态维度低了，稀疏观测才有意义。
一个小技巧是，在做敏感性分析时不要只用一个输出指标。比如你可以同时计算“10 秒时的温度”和“稳态温度”两个输出的 Sobol 指数，如果一个参数对稳态温度不敏感、但对瞬态响应敏感，它同样值得保留在状态里。这样筛选出的参数集合对工况变化和时间尺度都更加稳健。
4. 把模型和数据拧在一起：粒子滤波、稀疏识别与混合预测的实现细节
4.1 粒子滤波三步骨架：预测、更新、重采样
粒子滤波解决的核心问题是：模型参数不确定，观测数据只有零星几个温度测点，如何在参数空间中维持一组带权重的假设样本并随时间更新。标准的三步是预测、更新、重采样，示例如下：
class ParticleFilter: def __init__(self, n_particles, state_dim, process_noise, obs_noise): self.n_particles = n_particles self.state_dim = state_dim self.process_noise = process_noise self.obs_noise = obs_noise self.particles = np.random.randn(n_particles, state_dim) self.weights = np.ones(n_particles) / n_particles def predict(self, transition_func, control_input=None): for i in range(self.n_particles): if control_input is not None: self.particles[i] = transition_func(self.particles[i], control_input) else: self.particles[i] = transition_func(self.particles[i]) self.particles[i] += np.random.multivariate_normal( np.zeros(self.state_dim), self.process_noise) def update(self, observation, observation_func): likelihoods = np.zeros(self.n_particles) for i in range(self.n_particles): predicted_obs = observation_func(self.particles[i]) residual = observation - predicted_obs exponent = -0.5 * residual.T @ np.linalg.inv(self.obs_noise) @ residual likelihoods[i] = np.exp(exponent) self.weights *= likelihoods self.weights /= np.sum(self.weights) effective_n = 1.0 / np.sum(self.weights**2) if effective_n < self.n_particles / 2: self.resample() def resample(self): indices = np.random.choice( range(self.n_particles), size=self.n_particles, p=self.weights) self.particles = self.particles[indices] self.weights = np.ones(self.n_particles) / self.n_particles def estimate(self): mean = np.average(self.particles, weights=self.weights, axis=0) cov = np.cov(self.particles, rowvar=False, aweights=self.weights) return mean, cov逻辑说明：predict把每个粒子按照状态转移函数推进一步，再叠加过程噪声，用来表达参数本身的慢漂移和未建模扰动。update用高斯似然函数，根据传感器位置的预测残差给每个粒子重新加权。effective_n是有效粒子数，当它低于粒子总数一半时触发重采样——淘汰权重过小的粒子，复制权重大的粒子，随后权重被重新均分。
参数说明：
n_particles：粒子数量，示例用 100。状态维度 3~5 时 200 个粒子基本够用；维度超过 10 就要加到 1000 以上并提前评估计算量。process_noise：过程噪声协方差，物理含义是“参数随时间变化的自由程度”。设得太大，参数会到处跳，跟踪曲线噪声很大；设得太小，跟不上转速突变。论文场景下我建议在 0.005~0.02 区间调。obs_noise：观测噪声协方差，取温度传感器实际噪声方差的 1.5~3 倍较合适。设得太小，似然函数过于尖锐，会出现权重退化。
重采样是粒子滤波最容易出问题的地方。如果观测噪声过小，只有极少数粒子拿到有效权重，重采样后粒子多样性基本耗尽，粒子群坍缩到同一点，后续观测对估计的修正效果趋近于零。为了验证滤波器是否健康，我通常在每一轮观测更新后打印effective_n，如果它长期低于粒子总数的三分之一，就要回头检查观测噪声的设置，而不是急着增加粒子数。
4.2 多层粒子滤波 MLPF 到底“多”在哪里
单层粒子滤波在参数估计上有个短板：当待估参数之间存在相关性，或系统同时存在多个时间尺度的变化——比如快速转速变化叠加慢速热积累——单层滤波很难兼顾快速跟踪和稳定收敛。论文提出的 MLPF 思路是：把参数按响应时间尺度分层，慢变参数放在慢层滤波，快变参数放在快层滤波，层与层之间用估计结果互相衔接；每层内部结合稀疏识别约束参数演化方向，避免盲目搜索。
工程实现上，MLPF 可以理解为一个滤波器组：底层就是上面这个ParticleFilter类，各层共享观测函数，但过程噪声和时间步长不同。慢层用较小的过程噪声做慢速精细搜索，快层用较大的过程噪声做快速追踪。每一轮循环先把快层估计结果传给慢层作为参考均值，慢层再把平滑后的参数反馈给快层做约束。真正复现时不需要发明新算法，关键是设置好层间接口和不同尺度的噪声参数。
4.3 SINDy 稀疏识别：让数据反哺一个解释得动的模型
数据驱动部分论文用的是稀疏识别方法。SINDy 的核心是：对状态时间序列构造候选特征库，通过稀疏回归选出系数非零的少数项，得到一个稀疏、可解释的微分方程。用pysindy库实现很简洁：
import pysindy as ps def sparse_identification(X, t): model = ps.SINDy( optimizer=ps.STLSQ(threshold=0.1), # 稀疏阈值，低于此值的系数归零 feature_library=ps.PolynomialLibrary(degree=2), differentiation_method=ps.FiniteDifference(), ) model.fit(X, t=t) model.print() return model逻辑说明：X是被识别变量的时间序列，每行是一个时刻的状态；t是采样时间点。SINDy 内部先用有限差分计算时间导数，再在多项式特征库上做稀疏回归。STLSQ(threshold=0.1)表示绝对值小于 0.1 的系数会被置零，保留下来的项构成近似物理模型。
参数说明：
degree=2：特征库最高二次多项式。热网络这类物理过程二次项已经足够，更高次项通常在拟合噪声而不是物理规律。threshold：稀疏阈值。越小保留项越多，拟合精度高但解释性差；越大越稀疏，但可能丢掉真实存在的弱物理项。建议打印结果看系数量级后再调。
SINDy 在这个框架里不是用来替代热网络模型的，它有两个实际用途：一是对测点温度做局部动态识别，提取数据中主导的演化特征；二是为粒子滤波提供参数的先验轨迹——如果 SINDy 识别出某个热阻系数的时变规律，可以把对应状态的过程噪声设得更小，让粒子沿着轨迹搜索，避免在物理上不可能的区间浪费粒子。
4.4 混合预测策略：加权合并与测量更新
def hybrid_prediction(self, current_state, control_input): # 基于热网络模型的预测：取 simulate 最后一步温度 model_pred = self.thermal_model.simulate(current_state, [0, 1])[-1] if self.data_driven_model is not None: features = np.concatenate([current_state, [control_input]]) data_correction = self.data_driven_model.predict([features])[0] hybrid_pred = model_pred * 0.7 + data_correction * 0.3 else: hybrid_pred = model_pred return hybrid_pred逻辑说明：hybrid_prediction先用热网络模型基于当前状态外推一个单位时间步，再用随机森林等数据驱动模型对同一时刻的输出做偏差修正，最后按 0.7/0.3 权重合并。之所以保留 30% 的数据驱动分量，是为了在模型标定误差明显时保留一定的纠错能力；保留 70% 的模型分量，是为了保证在未见过的工况下预测结果仍然服从物理规律。
权重是固定的 0.7/0.3 吗？我的建议是把它当初始值，跑一轮之后根据残差分布调整。平稳工况下固定权重问题不大，但转速突变瞬间随机森林可能给出极端修正值，这时固定权重的混合预测会出现尖峰抖动。更好的做法是滑动窗口计算两路预测残差方差，让权重随残差倒变化——这一段在第 5 章配合具体故障场景展开。
测量更新的衔接方式是，update_with_measurements把传感器位置的温度作为观测值喂给粒子滤波器，由滤波器输出对状态向量的修正。整体循环是：模型预测 → 数据驱动修正 → 粒子滤波校正 → 进入下一步。
5. 避坑指南：主轴轴承热分析实战中的五条踩坑记录
这一章是实打实复现这个项目时记录下来的问题，每条都是现象、原因、解决的完整链路。
5.1 一跑 simulate 就 AttributeError，温度场完全出不来
现象：把热网络模型实例化后调用simulate，立刻抛出AttributeError: 'ThermalNetworkModel' object has no attribute 'R_env'。
原因：原始示例代码的heat_balance方法里使用了self.R_env，但构造函数只初始化了R和Q，漏掉了R_env这个属性。这是资源本身的遗漏，不是使用问题。
解决：在__init__中显式初始化self.R_env = np.full(nodes, 50.0)，再按实际散热条件调整到 10~200 K/W 之间。从那以后我每次拿到别人的热网络代码，第一件事就是在构造函数里排查所有属性是否都有初值，这是最廉价的自检习惯。
5.2 温度场出现非物理解：热阻矩阵不对称导致的假温差
现象：模拟结果中，外圈节点温度长期高于内圈节点，但实际物理过程是内圈发热更集中，温度理应更高。调整热源数值后仍然无法消除这个反常温差。
原因：热阻矩阵只随机填充了上三角，没有做对称化，或者手动修改了R[i][j]而忘记同步R[j][i]。热流计算依赖成对热阻，不对称会让同一个温差在两个方向算出不同热流，相当于凭空制造了热汇或热源。
解决：封装一个辅助函数统一处理热阻赋值，保证成对写入；同时做等温场检验——所有节点初始温度相同，如果模拟过程中温度场还出现节点间温差，说明基础配置有误。这是成本最低的模型自检手段。
5.3 有效粒子数骤降，重采样后粒子群完全重复
现象：跑十几步后打印effective_n，发现经常小于 10，重采样后粒子矩阵的多行完全相同，继续更新但估计值不变化。
原因：观测噪声协方差obs_noise设得太小，似然函数过于尖锐，少数粒子吃掉几乎全部权重。重采样时反复复制这几个粒子，多样性耗尽。
解决：把obs_noise放大到传感器实际噪声方差的两到三倍，并将重采样触发阈值从n_particles / 2放宽到n_particles / 3。如果重采样后仍然大量重复，改用正则化重采样，在复制后的粒子上叠加微小的多维高斯噪声。
5.4 SINDy 识别出一堆高次交叉项，物理上无法解释
现象：model.print()输出的方程里出现x0^2*x1这类交叉项，系数还不小，但温度动态模型不应当包含这种组合。
原因：特征库选了二次多项式，但实际上系统主导的动态可能只是一次项加常数项；稀疏阈值 0.1 又偏低，数值差分产生的噪声被当作有效信号纳入了回归。
解决：先把threshold提升到 0.2 或 0.3，观察残差变化；如果结果仍然复杂，把degree降到 1，从线性模型开始识别。识别结果必须能被物理过程解释，在这样的框架里，白盒的可解释性比训练集拟合精度更重要。
5.5 转速突变时混合预测出现尖峰抖动
现象：转速从 3000 r/min 跳到 6000 r/min 时，预测温度先冲高再回落，出现一个明显尖峰。
原因：固定权重 0.7/0.3 中，随机森林修正项在训练数据里没有覆盖足够的突变样本，转速跳变瞬间给出极端修正值；权重又无法动态降低这个噪声的影响。
解决：改成自适应权重，用滑动窗口统计最近 N 步两路预测的残差方差，权重随残差方差倒数变化。窗口长度取 10~20 步即可，转速突变瞬间数据驱动残差方差飙升，权重自动回落，混合预测回到模型主导，尖峰消失。
6. 用少量传感器验证温度场预测的可靠度：三张图加一个交叉验证
把这个框架用在实际轴承实验台上时，传感器一般只布两三个位置，比如外圈外表面和轴承座边缘。如何证明其他位置的预测温度可信？我的验证流程固定是下面四步。
第一张图，真值与估计值对比温度曲线。把所有有测点位置的模拟温度和粒子滤波估计值画在一起，确认趋势一致，再计算测点位置的 RMSE。经验值上是，稳态 RMSE 在 1℃ 以内算合格；如果超过 3℃，先检查热源计算参数f0、f1的量级，再观察测点布置位置是否避开了强非线性区。
第二张图，残差频谱图。把每个时间的预测残差做 FFT，看到明显的周期性峰说明存在未建模的周期扰动，可能是油膜振荡或冷却液周期性启停。这时需要让 MLPF 的多层结构把快周期和慢周期分开处理，或者在热网络里增加对应的周期热源项。
第三张图，参数收敛曲线。这条曲线画的是被粒子滤波当作状态量的热阻或热源系数随时间的变化。判断健康与否的标准是：曲线应当在若干次观测更新后趋于平滑，而不是每次新观测都剧烈摆动。如果参数持续漂移、始终找不到稳定平台，说明问题已经超出滤波器调参能解决的范围，需要回头检查热网络模型结构本身。
第四个环节是留一测点交叉验证，这是我最推荐也最容易被忽略的方法：把 N 个测点中的一个在滤波过程的观测更新中剔除，只用剩余测点做更新，然后检查“被剔除位置”的估计误差。如果保留测点位置的 RMSE 明显低于未保留位置，说明温度场的可观测性偏弱，传感器布点需要调整，或者需要增设测点。很多实际工况里加传感器不现实，但这个方法至少能告诉你当前布点方案下预测值在空间上哪些区域可信。从那以后，我每次搭混合驱动热分析流程都会强制先跑一遍留一测点交叉验证，再确定最终布点方案——这个习惯帮我躲过了多次实验台上的无效数据，希望帮到你。
本文还有配套的精品资源，点击获取
