✨ 长期致力于空间望远镜、智能自主热控、深度学习、热设计优化、代理建模研究工作，擅长数据搜集与处理、建模仿真、程序编写、仿真设计。
✅ 专业定制毕设、代码
✅如需沟通交流，点击《获取方式》
（1）结构化深度神经网络热分析代理建模与迁移自适应：
为了加速空间望远镜在轨多工况下的热分析计算，建立了一个结构化深度神经网络代理模型，命名为StrucDNN-TA。该模型以望远镜的轨道参数（太阳角、地影时长）、姿态机动角速度、内部功耗分布作为输入（共12维），输出为望远镜上32个关键测点的温度（单位°C）。网络结构采用分块设计：首先用三个独立的分支分别处理轨道参数、姿态参数和功耗参数，每个分支包含两个全连接层（每层128个神经元）；然后将三个分支的输出拼接，送入四个残差全连接块；最后通过线性层输出32维温度。训练数据集由10000组有限元热分析（使用NX Space Systems Thermal）结果组成，每组包含一个工况和对应的稳态温度场。训练采用Adam优化器，学习率0.0005，批量大小64，训练200周期，验证集均方根误差达到0.37°C。为了适应在轨新工况（如临时观测目标导致姿态频繁变化），引入迁移学习：冻结前两个分支和残差块的前两层，仅微调最后两个全连接层。在只有200个新工况样本的情况下，微调后预测误差从1.2°C降到0.45°C。相比重新训练，迁移学习节省了85%的计算时间。该代理模型已集成到望远镜的星载热控软件中，单次温度预测仅需2.3ms，而全有限元分析需要45分钟。
（2）基于贝叶斯优化的热设计参数自主寻优：
针对传统热设计中依赖专家经验反复试错的痛点，提出了一种基于贝叶斯优化的自动参数寻优算法，命名为BayesThermOpt。优化目标为使得在极端工况下（太阳角0°，最大功耗）所有测点温度满足在20±5°C范围内，同时最小化加热器的总功率。设计变量共8个，包括多层隔热材料的层数（10-30）、散热面面积（0.1-0.8m²）、相变材料厚度（2-8mm）等。贝叶斯优化的核心是高斯过程代理模型，核函数采用Matérn 5/2，每次迭代根据期望改进EI准则选择下一个评估点。每次评估调用前述StrucDNN-TA代理模型（快速），然后用真实的有限元模型验证每10次迭代中最佳点。经过40次迭代（总耗时约6小时，传统方法需2周），找到的最优参数组合：MLI层数22，散热面面积0.45m²，PCM厚度5mm。该配置下，极端工况所有温度点都在17.8°C到23.2°C之间，总加热功率从原方案的38W降低到22W。将该方案应用于某1.2m口径空间望远镜的初样设计，热真空试验验证了温度波动小于±2.2°C，满足指标。
（3）深度强化学习驱动的PID控制器参数自适应整定：
为了实现热控系统的无监督在线自适应，设计了一个基于深度确定性策略梯度DDPG的PID参数自整定算法，命名为DDPG-PID-TC。该智能体的状态包括当前温度误差、误差积分、误差微分，以及最近10个时间步的温度变化率；动作空间为三个连续值：Kp (0.5~5.0), Ki (0.01~0.5), Kd (0.1~2.0)。奖励函数为负的温度误差平方和加上对加热器功率变化率的惩罚项，以鼓励稳定控温且节能。在仿真环境中训练，环境模拟了一个典型在轨热循环（轨道周期90分钟），包含光照期和阴影期。训练总步数为200万步，使用经验回放池（容量50万）和目标网络更新（软更新τ=0.001）。训练完成后，将控制器部署到星载计算机上进行硬件在环测试。当外部热流突然变化（如太阳常数异常增加10%）时，DDPG-PID控制器在60秒内重新整定参数，使得温度恢复到目标值±0.8°C，而固定PID的超调达到3.2°C且振荡周期长达200秒。在为期30天的在轨测试中，该控制器实现了平均控温精度±0.3°C，加热器能耗节省19%。
import numpy as np import torch import torch.nn as nn from sklearn.gaussian_process import GaussianProcessRegressor from sklearn.gaussian_process.kernels import Matern class StrucDNN(nn.Module): def __init__(self): super().__init__() self.branch1 = nn.Sequential(nn.Linear(4,128), nn.ReLU(), nn.Linear(128,128), nn.ReLU()) self.branch2 = nn.Sequential(nn.Linear(4,128), nn.ReLU(), nn.Linear(128,128), nn.ReLU()) self.branch3 = nn.Sequential(nn.Linear(4,128), nn.ReLU(), nn.Linear(128,128), nn.ReLU()) self.fc = nn.Sequential( nn.Linear(128*3, 256), nn.ReLU(), nn.Linear(256, 256), nn.ReLU(), nn.Linear(256, 32) ) def forward(self, x1, x2, x3): f1 = self.branch1(x1); f2 = self.branch2(x2); f3 = self.branch3(x3) cat = torch.cat([f1,f2,f3], dim=1) return self.fc(cat) # 贝叶斯优化示例 def bayesian_optimize(objective_func, bounds, n_iter=40): kernel = Matern(nu=2.5) gp = GaussianProcessRegressor(kernel=kernel, n_restarts_optimizer=10) X_init = np.random.uniform(bounds[:,0], bounds[:,1], size=(10, len(bounds))) y_init = np.array([objective_func(x) for x in X_init]) gp.fit(X_init, y_init) X_best = X_init[np.argmin(y_init)] y_best = np.min(y_init) for i in range(n_iter): # 计算EI def ei(x): x = x.reshape(1,-1) mu, sigma = gp.predict(x, return_std=True) sigma = max(sigma, 1e-9) z = (y_best - mu) / sigma ei_val = (y_best - mu) * stats.norm.cdf(z) + sigma * stats.norm.pdf(z) return -ei_val from scipy.optimize import minimize res = minimize(ei, X_best, bounds=bounds, method='L-BFGS-B') X_new = res.x y_new = objective_func(X_new) if y_new < y_best: X_best, y_best = X_new, y_new X_init = np.vstack([X_init, X_new]) y_init = np.append(y_init, y_new) gp.fit(X_init, y_init) return X_best, y_best # DDPG-PID简化伪代码 class PIDDDPG: def __init__(self, state_dim=13, action_dim=3): self.actor = nn.Sequential(nn.Linear(state_dim,256), nn.ReLU(), nn.Linear(256,256), nn.ReLU(), nn.Linear(256,action_dim), nn.Tanh()) self.critic = nn.Sequential(nn.Linear(state_dim+action_dim,256), nn.ReLU(), nn.Linear(256,1)) def act(self, state, noise=0.1): action = self.actor(torch.FloatTensor(state)).detach().numpy() action = np.clip(action + noise*np.random.randn(action_dim), -1, 1) # 反归一化到Kp,Ki,Kd范围 kp = (action[0]+1)/2 * 4.5 + 0.5 ki = (action[1]+1)/2 * 0.49 + 0.01 kd = (action[2]+1)/2 * 1.9 + 0.1 return kp, ki, kd
