主要内容

pdeInterpolant

Interpolant节点数据选择的位置

扩展所有的页面

pdeInterpolant和[p, e t]的代表FEMesh数据不推荐。使用interpolateSolution和evaluateGradient插入一个PDE解决方案及其梯度没有切换到任意点[p, e t]表示。

描述

一个interpolant允许您评估几何中的PDE的解决方案在任何时候。

解决偏微分方程工具箱™返回解决方案值的节点,这意味着网格点。评估一个插值的方案在其他点几何中,创建一个pdeInterpolant对象,然后调用评估函数。

创建

语法

F = pdeInterpolant (p t u)

描述

例子

F = pdeInterpolant (p,t,u)返回一个interpolantF基于数据点p、元素t和数据值的点,u。

使用meshToPet获取p和t数据插值使用pdeInterpolant。

输入参数

全部展开

数据点位置,指定为一个矩阵有两个或三个行。每一列的p是一个2 d或3 d点。有关详细信息,请参见网格数据。

对于二维问题,构造p使用initmesh函数,或出口网菜单的PDE Modeler应用。2 d或3 d几何使用PDEModel对象,获得p使用meshToPet函数model.Mesh。例如,[p, e t] = initmesh (g)或[p, e t] = meshToPet (model.Mesh)。

三角元素,指定为一个矩阵。有关详细信息,请参见网格数据。

对于二维问题,构造t使用initmesh函数,或出口网菜单的PDE Modeler应用。2 d或3 d几何使用PDEModel对象,获得t使用meshToPet函数model.Mesh。例如,[p, e t] = initmesh (g)或[p, e t] = meshToPet (model.Mesh)。

插入数据值,指定为一个向量或矩阵。通常情况下,u返回的是PDE的问题的解决方案assempde,抛物线,双曲,或者另一个解算器。例如,u = assempde (b、p, e、t、c、a、f)。您还可以导出u从解决菜单的PDE Modeler中应用。

矩阵的维度u取决于这个问题。如果np列的数量吗p,N在PDE系统方程的数目,然后呢u有N * np行。第一个np行对应方程1,下一个np行对应方程2,为抛物线或双曲问题,等等。u有一列每个解决方案;否则,u是一个列向量。

对象的功能

评估 插入数据到选定的位置

例子

全部折叠

打开生活的脚本

这个例子显示了如何创建一个pdeInterpolant从一个标量PDE的解决之道。

解方程 - - - - - - Δ u = 1 在单位圆零狄利克雷条件。

g0 = [1, 0, 0, 1];%的圆心与半径1 (0,0)科幻小说=“C1”;g = decsg (g0、科幻小说、科幻小说”);%分解几何矩阵模型= createpde;通用= geometryFromEdges(模型、g);% 0狄利克雷条件applyBoundaryCondition(模型,“边界条件”,…“边缘”(1:gm.NumEdges),…“u”,0);[p, e t] = initmesh(通用);c = 1;= 0;f = 1;u = assempde(模型、p、e t c、a、f);

构造一个interpolant的解决方案。

F = pdeInterpolant (p t u);

评估interpolant在一个正方形的四个角。

撅嘴= [0,1/2,1/2,0;0,0,1/2,1/2);uOut =评估(F,撅嘴)
uOut =4×10.2485 0.1854 0.1230 0.1852

的值uOut (2)和uOut (4)几乎是平等的,他们应该为对称点对称的问题。

版本历史

介绍了R2014b

另请参阅

|

主题