TimeDependentResults

随时间变化的偏微分方程解和导出量

展开全部页面

描述

一个TimeDependentResults对象以便于绘图和后处理的形式包含偏微分方程及其梯度的解。

一个TimeDependentResults对象包含解及其在三角形或四面体网格节点上计算的梯度，由generateMesh．
节点上的解决方案值显示在NodalSolution财产。
的中显示了求解时间SolutionTimes财产。
节点上解值梯度的三个分量出现在XGradients，YGradients,ZGradients属性。
的数组尺寸NodalSolution，XGradients，YGradients,ZGradients使您能够为指定的时间指标和偏微分方程系统中的方程指标提取解和梯度值。

将解或其梯度插值到自定义网格(例如，由meshgrid),使用interpolateSolution或evaluateGradient．

创建

有几种方法可以创建TimeDependentResults对象:

方法解决与时间相关的问题solvepde函数。这个函数返回一个PDE解TimeDependentResults对象。这是推荐的方法。
方法解决与时间相关的问题抛物线或双曲函数。然后使用createPDEResults函数来获取TimeDependentResults返回的PDE解决方案抛物线或双曲．请注意,抛物线而且双曲都是遗留函数。它们不推荐用于解决PDE问题。

属性

全部展开

`网`- - - - - -有限元网格
`FEMesh`对象

有限元网格，返回为FEMesh属性对象。

`NodalSolution`- - - - - -节点上的解决方案值
向量|数组

节点上的解决方案值，作为向量或数组返回。的尺寸NodalSolution,请参阅解、梯度和通量的维度金宝搏官方网站．

数据类型:双
复数支持:金宝app是的

`SolutionTimes`- - - - - -解决方案*
真正的向量

解时间，作为实向量返回。SolutionTimes是一样的tlist输入solvepde，或tlist对遗留的输入抛物线或双曲解决者。

数据类型:双

`XGradients`- - - - - -x-节点上的梯度分量
向量|数组

x-节点上梯度的分量，作为向量或数组返回。的尺寸XGradients,请参阅解、梯度和通量的维度金宝搏官方网站．

数据类型:双
复数支持:金宝app是的

`YGradients`- - - - - -y-节点上的梯度分量
向量|数组

y-节点上梯度的分量，作为向量或数组返回。的尺寸YGradients,请参阅解、梯度和通量的维度金宝搏官方网站．

数据类型:双
复数支持:金宝app是的

`ZGradients`- - - - - -z-节点上的梯度分量
向量|数组

z-节点上梯度的分量，作为向量或数组返回。的尺寸ZGradients,请参阅解、梯度和通量的维度金宝搏官方网站．

数据类型:双

对象的功能

`evaluateCGradient`	评估PDE溶液的通量
`evaluateGradient`	评估任意点上PDE解的梯度金宝搏官方网站
`interpolateSolution`	插值PDE解到任意点

例子

全部折叠

一个抛物问题的解法

打开实时脚本

用二维几何解一个抛物线问题。

创建并查看几何图形:一个正方形和一个圆形子域。

以(1,1)为圆心的正方形%Rect1 = [3;4;0;2;2;0;0;0;圆心为(1.5,0.5)的圆Circ1 = [1;1.5;.75;0.25];在圆上附加额外的0Circ1 = [Circ1;zero (length(rect1)-length(Circ1)，1)];Gd = [rect1,circ1];Ns = char(“通过rect1”，“circ1”）;Ns = Ns ';科幻小说=“通过rect1 + circ1”；[dl,bt] = decsg(gd,sf,ns);pdegplot (dl,“EdgeLabels”，“上”，“FaceLabels”，“上”)轴平等的ylim ([-0.1, 2.1])

图中包含一个轴对象。axis对象包含11个类型为line, text的对象。

在PDE模型中包含几何图形。

Model = createpde();geometryFromEdges(模型、dl);

设置边界条件，上面和左边边缘温度为10。

applyBoundaryCondition(模型,“边界条件”，“边缘”(2、3),“u”10);

设定初始条件，正方形区域温度为0，圆形区域温度为100。

setInitialConditions(模型中,0);100年setInitialConditions(模型,“脸”2);

定义模型系数。

specifyCoefficients(模型,“米”0,' d ', 1“c”, 1“一个”0,“f”, 0);

求解从0到1/2的问题，步骤为0.01。

generateMesh(模型,“Hmax”, 0.05);Tlist = 0:0.01:0.5;Results = solvepde(model,tlist);

画出乘以0.02、0.04、0.1和0.5的解。

sol = results. nodesolution;次要情节(2 2 1)pdeplot(模型,“XYData”索尔(:3)标题(“时间0.02”) subplot(2,2,2)“XYData”索尔(:5)标题(“时间0.04”) subplot(2,2,3)“XYData”索尔(:11)标题(“时间0.1”) subplot(2,2,4)“XYData”索尔(:51)标题(“时间0.5”）

图中包含4个轴对象。标题为Time 0.02的坐标轴对象1包含一个patch类型的对象。标题为Time 0.04的Axes对象2包含一个patch类型的对象。标题为Time 0.1的坐标轴对象3包含一个patch类型的对象。标题为Time 0.5的Axes对象4包含一个patch类型的对象。