PDE约束优化控制问题-NN结合DAL

已于 2023-02-07 11:42:32 修改
阅读量228 收藏
点赞数
分类专栏： PDE约束优化控制问题 Deep learning 文章标签： 线性代数 概率论 几何学
于 2022-02-16 19:44:23 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/forrestguang/article/details/122970093
版权
Deep learning 同时被 2 个专栏收录
42 篇文章 71 订阅 ¥239.90 ¥399.90
订阅专栏
PDE约束优化控制问题
27 篇文章 4 订阅 ¥259.90 ¥399.90
订阅专栏
本文探讨了偏微分方程（PDE）约束的优化控制问题，通过神经网络（NN）与差分算法（DAL）的结合来求解。文章首先介绍了问题的基本形式，并给出KKT系统的表达式。接着，详细阐述了如何利用PINN解决PDE，并通过三个算例——算例1 Gongwei、算例2 和算例3（Neumann边界条件）展示了方法的应用。最后，利用最小二乘法更新控制变量，确保其在给定的约束范围内。
摘要由CSDN通过智能技术生成

引入算例基本形式
min ⁡ u ∈ U a d , y

了解本专栏 订阅专栏 解锁全文
Galerkin码农选手
关注
专栏目录
订阅专栏
PDE约束优化控制问题-dolfin软件求解
forrestguang的博客
01-28 810
dolfin安装 dolfin是我们发现的一个专门求解PDE约束优化问题的软件，或者说是一个python的库，参考网页添加链接描述，这个软件暂时缺乏用户手册，不过网页上有不少demo，本人就是根据网页上的demo来学习使用dolfin的，具体的原理也请参考网页，本人不做过多解释。这个软件的使用需要结合fenics库，参考fenics官网添加链接描述,fenics有用户手册，自己可以到官网下载。 下面重点讲解dolfin以及相关的库的安装方法，本人都是在linux服务器安装的 首先还是先建一个虚拟环境，本人建
解决PDE的问题：PDE-matlab开发
05-30
各位我正在解这个方程 DuDt+exp(-4x^2)DuDx=0 u(x,0)=atanh(0.1exp(-100x^2)); 对于 matlab，这个 pde 必须很容易学习。 请帮助我可能是我在 mcode 中有错误
Python 相关以及安装 FEniCS
米饭的白色的专栏
06-14 7586
1. Python 相关，安装多版本python和pip共存 这一部分摘抄自：https://blog.csdn.net/hymanjack/article/details/80285400 。 ubuntu18.04自带python3，但是没有python2，pip2，pip3。 sudo apt install python2.7 #安装python2.7 sudo apt inst...
记录 fenics 的使用过程
米饭的白色的专栏
06-25 2498
偶尔要用到 fenics 来计算, 但是各种不熟悉, 每次网上找, 资料很多很乱, 而且 fenics 的版本更新后, 很多函数的用法都报错, 所以这里还是记录一下自己的使用过程 定义网格和空间后, 得到空间的自由度编号 from dolfin import * import numpy as np mesh = UnitSquareMesh(4, 4) # Create a scalar and vector spaces space_S1 = FunctionSpace(mesh, 'CG', 1
fenics入门笔记
热门推荐
qq_41230365的博客
03-18 1万+
fenics名称释义： fe：finite element的简写 cs：computational software的简写 ni：有了fe和cs后，由于最初fenics软件是在芝加哥大学(简称为phoenix）编译的，故而在其间加入ni就很自然而然了。其实就是取个谐音。 fenics有如下几个核心部件： 1. dolfin： dynamic object-oriented lib...
偏微分方程约束下的优化控制问题（PDE-constrained optimal control problems）
forrestguang的博客
02-08 1323
{min⁡(y(x),u(x))∈Y×UJ(y(x),u(x)), s.t. F(y(x),u(x))=0 in Ω, and u(x)∈Uad,\left\{\begin{aligned} &\min _{(y(\mathbf{x}), u(\mathbf{x})) \in Y \times U} J(y(\mathbf{x}), u(\mathbf{x}) ),\\ &\text { s.t. } \ \mathbf{F}(y(\mathbf{x}),u(\mathbf{x
pde.rar_ODE-PDE_PDE MATLAB_mathematics_pde
09-21
标题中的“pde.rar_ODE-PDE_PDE MATLAB_mathematics_pde”暗示了这是一个与偏微分方程（PDE）和常微分方程（ODE）相关的MATLAB工程数学项目。描述中提到的“PDE And ODE in engineering mathematics”进一步确认了这...
Evans-PDE-Solution-Chapter-5-Sobolev.pdf
11-21
根据给定文件的内容，我们可以了解到该文件涉及的是偏微分方程（Partial Differential Equations，简称PDE）领域，特别是Sobolev空间的概念和相关证明。由于文件内容涉及的是Sobolev空间的定理证明，我们将重点分析...
Guy-Gilboa-PDE-based-Image-Filtering.rar_Guy Gilboa_PDE图像_pde 滤波
07-14
Guy Gilboa是该领域的知名专家，他的研究深入探讨了如何利用PDE理论来优化和改善图像质量。本文将详细阐述PDE图像滤波的基本原理、方法及其在实际应用中的价值。 首先，理解PDE滤波的核心在于PDE的特性。PDE是一种...
16年PDE-FIND.rar
07-01
标题中的“16年PDE-FIND.rar”是一个压缩包文件，其中包含了关于如何利用数据来发现和理解偏微分方程（Partial Differential Equations, PDEs）的源代码。PDEs是数学中的一个重要分支，广泛应用于物理、工程、生物等...
pytorch实现NS方程求解-基础PINN
forrestguang的博客
04-14 7481
问题描述Kovasznay problem 参考文献：Dockhorn, T. . “A Discussion on Solving Partial Differential Equations using Neural Networks.” (2019). {−νΔu+uux+vuy+px=f1,−νΔv+uvx+vvy+py=f2,ux+vy=0,u∣∂Ω=u0,v∣∂Ω=v0.\left\{\begin{array}{l} -\nu \Delta u + uu_x + v u_y + p_x = f_
Fouier Net和DeepOnet等求解器算法解读和代码
forrestguang的博客
11-21 3904
首先还是先给定一个泊松方程，其中区域 Ω=[0,1+μ1]×[0,1],μ1∈[1,3]Ω=[0,1+μ_1]×[0,1],μ_1∈[1,3]Ω=[0,1+μ1​]×[0,1],μ1​∈[1,3] ，边界条件 g=1g=1g=1 。 −Δu=f,x∈Ω,u∣∂Ω=g.−Δu=f,x∈Ω,\\ u|∂Ω=g.−Δu=f,x∈Ω,u∣∂Ω=g. 以往的传统数值方法都是针对固定的一个右端项求出一个对应的数值解，这种做法的核心是构造一个输入节点数目为2，输出节点数目为1的模型（有限元句柄函数或者全连接神经网络），基
torch.nn.Parameters用法以及PINN求解PDE和画图
forrestguang的博客
12-14 3382
PINN介绍和代码 以前我的PINN代码一直都是接受(xi,yj)(x_i,y_j)(xi​,yj​)输出u(xi,yj)u(x_i,y_j)u(xi​,yj​) {−Δu=f,inΩ,u=u0,on  ∂Ω.\left\{\begin{array}{l} -\Delta u = f , \quad in \Omega,\\ u=u_0 ,\quad on \;\partial \Omega. \end{array}\right.{−Δu=f,inΩ,u=u0​,on∂Ω.​
基于pytorch的PINN求解含时间项的泊松方程
forrestguang的博客
05-06 3198
方程引出 ut−uxx−uyy−uzz=f,x∈Ω,u_t-u_{xx}-u_{yy}-u_{zz}=f,x\in \Omega,ut​−uxx​−uyy​−uzz​=f,x∈Ω, u(0,x,y,z)=u0,u(0,x,y,z)=u_0,u(0,x,y,z)=u0​, u(t,x,y,z)=h,(x,y,z)∈∂Ωu(t,x,y,z)=h,(x,y,z)\in \partial \Omegau(t,x,y,z)=h,(x,y,z)∈∂Ω 这里同样使用PINN来求解方程。 难点 第一个困难在于我们需要给出几
PFNN两个神经网络组合训练求解泊松方程
forrestguang的博客
09-07 3002
边界条件的自动匹配 仍然考虑这个微分方程 首先给出精确解以及采样过程： def U(X, order,prob): if prob==1: temp = 10*(X[:,0]+X[:,1])**2 + (X[:,0]-X[:,1])**2 + 0.5 if order[0]==0 and order[1]==0: return torch.log(temp) if order[0]==1 and order[1]==0:
PINN求解PDE
forrestguang的博客
04-18 2963
PINN求解Dirichlet泊松方程 泊松方程引入 −Δu=f,x∈Ω=(0,1)2-\Delta u=f,x \in \Omega=(0,1)^2−Δu=f,x∈Ω=(0,1)2 u=g,x∈∂Ωu=g,x \in \partial \Omegau=g,x∈∂Ω import torch import time import numpy as np import matplotlib.pyplot as plt import torch.nn as nn def UU(X, order,prob):
NS方程求解-PointNet和升维思想（效果很差）
forrestguang的博客
04-14 2532
PointNet求解NS方程 参考文献：Kashefi, A. , and T. Mukerji . “Physics-Informed PointNet: A Deep Learning Solver for Steady-State Incompressible Flows and Thermal Fields on Multiple Sets of Irregular Geometries.” (2022). PointNet网络的特点在于中间通过pool提取过一次全局特征，这个特点导致一个问题。
NS方程求解-NSFnet
forrestguang的博客
04-14 2454
NSFnet求解CVD正问题 参考文献： Yang, H. , F. N. Hwang , and X. C. Cai . “NONLINEAR PRECONDITIONING TECHNIQUES FOR FULL-SPACE LAGRANGE-NEWTON SOLUTION OF PDE-CONSTRAINED OPTIMIZATION PROBLEMS *.” Siam Journal on entific Computing (2016). 其中γ=0.01,Re=1,Pr=0.72,Gr=1e
pde约束优化控制问题求解-间断问题
最新发布
09-08
PDE约束优化控制问题是指在优化问题中存在偏微分方程（Partial Differential Equation, PDE）约束的控制问题。常见的PDE约束包括椭圆型、抛物型和双曲型方程等。 而间断问题是指在优化控制问题中，相关参数或变量会突然发生改变，即存在间断点的情况。间断问题的出现常常是由于系统的特性突变、外部扰动或控制策略的调整等原因引起的。 解决PDE约束优化控制问题的间断问题，可以采用以下方法： 1. 近似法：通过使用适当的数值方法，将偏微分方程转化为一系列离散点上的代数方程，然后利用求解代数方程的方法求解。在这种方法下，需要考虑间断点对数值方法的影响，并采取相应的调整措施。 2. 多模型方法：将间断问题转化为多个连续的子问题，每个子问题对应一段不连续的参数或变量。然后在每段连续区间上求解连续的PDE约束优化控制问题，最后将所有段的解整合起来。这种方法的关键是确定每个段的边界条件，以及在段之间实现平滑过渡。 3. 非光滑控制方法：通过使用非光滑控制理论，考虑到间断点的存在，建立非光滑控制模型，并在该模型下进行求解。这种方法的优势是能够处理复杂的间断问题，但需要深入研究非光滑控制理论及其应用。 总之，解决PDE约束优化控制问题求解的间断问题是一个具有挑战性的任务，需要合适的数值方法、数学模型和算法来克服。不同的方法适用于不同的问题，具体的选择应根据具体问题的特点和要求来决定。
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

Galerkin码农选手

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付：¥1
获取中
扫码支付

您的余额不足，请更换扫码支付或充值

打赏作者

实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

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

余额充值