一区新风向:PINN+贝叶斯!!

今天给大家推荐一个好发一区TOP的思路:PINN+贝叶斯!

贝叶斯的引入,为PINN的改进提供了新视角,在提高预测准确性、可靠性、计算效率和规避过拟合方面,展现了强大的优势。比如模型BPIELM便在准确性提高30%的同时,训练速度提高70倍!

主要在于:PINN利用深度学习的强大计算能力,能够高效地处理复杂的数据集和模型结构,而贝叶斯方法则提供了有效的模型选择和参数估计手段,两者相辅相成。尤为特殊的是,该结合能在不确定性条件下提供更丰富的信息输出,如后验标准差等,使模型更好地处理复杂数据集和场景,这是其他方法做不到的。

此外,该思路目前还在蓝海期,创新空间很大!为让大家能紧跟领域前沿,找到更多idea启发,我给大家准备了8种创新思路及开源代码!

论文原文+开源代码需要的同学看文末

论文:Bayesian Physics-Informed Extreme Learning Machine for Forward and Inverse PDE Problems with Noisy Data
内容

该论文介绍了一种名为贝叶斯物理信息极学习机(BPIELM)的新方法,用于解决带有噪声数据的正向和逆向偏微分方程(PDE)问题。BPIELM结合了极学习机(ELM)和贝叶斯方法,通过引入先验概率分布并使用贝叶斯方法估计参数的后验分布来量化噪声数据带来的不确定性。与物理信息神经网络(PINN)和极学习机(PIELM)相比,BPIELM在预测准确性和计算成本方面表现出显著优势,尤其是在噪声场景下。

论文:Randomized Physics-Informed Neural Networks for Bayesian Data Assimilation
内容

该论文提出了一种名为随机化物理信息神经网络(rPINN)的方法,用于在存在噪声数据的逆偏微分方程(PDE)问题中进行不确定性量化。该方法通过在物理信息神经网络(PINN)的损失函数中引入随机噪声,并通过对不同噪声实现的最小化问题的解来近似后验分布,从而克服了传统贝叶斯PINN(BPINN)方法在非线性逆PDE问题中汉密尔顿蒙特卡洛(HMC)收敛失败的问题。

论文:Physics-Informed Bayesian Optimization of Variational Quantum Circuits
内容

该论文提出了一种新的贝叶斯优化方法,用于变分量子电路(VQE)的物理信息贝叶斯优化。该方法通过在VQE的目标函数中引入一个与量子电路结构相匹配的VQE核,显著降低了后验不确定性。同时,文章提出了一种名为EMICoRe的新型采集函数,能够主动利用VQE核的归纳偏差,将预测不确定性低的区域视为间接“观测”到的区域。这种方法仅需要在搜索域中的三个点进行观测,就能确定整个一维子空间中的完整目标函数。

论文:Correcting model misspecification in physics-informed neural networks
内容

该论文提出了一种新方法来纠正物理信息神经网络(PINNs)在学习控制方程时可能出现的模型错误规范问题。研究者们通过引入另一个深度神经网络(DNN)来模拟不完美模型与观测数据之间的差异,从而纠正PINNs中的错误规范模型。

 关注下方《人工智能学起来》

回复“BPINN”获取全部论文+开源代码

码字不易,欢迎大家点赞评论收藏

### 使用 MATLAB 实现物理信息神经网络(PINN) #### 构建 PINN 的基本框架 为了在 MATLAB 中实现物理信息神经网络(PINN),需要定义个多层感知器(MLP)模型,该模型能够接受时间和空间坐标作为输入,并输出相应的解。对于特定的物理问题,比如无压含水层瞬变流或是 Burgers 方程,可以通过调整损失函数中的物理约束项来适应不同的应用场景[^1]。 ```matlab % 定义 MLP 结构 layers = [ featureInputLayer(2,'Name','input') % 输入层, 接受 t 和 x fullyConnectedLayer(20,'Name','fc1') reluLayer('Name','relu1') fullyConnectedLayer(20,'Name','fc2') reluLayer('Name','relu2') fullyConnectedLayer(1,'Name','fc3') regressionLayer('Name','output')]; % 输出层, 对应 u(t,x) % 创建 DNN 模型 net = trainNetwork(...); ``` #### 设定物理约束条件 为了让神经网络理解并遵循给定的物理规律,在训练过程中应当引入额外的惩罚项至总损失函数之中。这些惩罚项代表了待解决问题所特有的微分方程及其边界条件或初始状态等先验知识[^2]。 ```matlab function loss = computeLoss(net,params,t,x,u_true) % 计算预测值 u_pred = forward(net,[t(:),x(:)]); % 物理残差计算 res = ...; % 根据具体 PDE 表达式编写 % 总损失=数据拟合误差+物理致性罚项 lambda_data = 1; lambda_phys = 1e-4; loss = mean((u_pred-u_true).^2)*lambda_data + ... sum(res.^2)/numel(res)*lambda_phys; end ``` #### 处理数值稳定性问题 针对某些复杂场景下可能出现的梯度爆炸或者消失现象,可以考虑采用改进策略如自适应权重缩放机制、特征标准化处理等方式提高算法收敛性能以及泛化能力[^3]。 ```matlab % 应用批量归化改善优化过程 layers = [ ... batchNormalizationLayer('Name','bn1') ... ]; ```
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值