19、MATLAB PDE 求解与应用实战

MATLAB PDE 求解与应用实战

1. PDE 求解结果与网格数据导出

1.1 导出求解结果

PDE 求解器通常以图形形式展示结果，但实际应用中常需要数值。可将定义的 u 值及获取这些值的网格节点坐标导出到 MATLAB 工作区。操作步骤如下：
1. 点击 Solution 主菜单按钮弹出菜单中的 Export Solution 选项。
2. 在弹出小面板（图 6.19A）中，选择是否更改默认解变量名，然后点击 OK 按钮。
3. 在 MATLAB 工作区窗口，会出现所选名称的变量。一般情况下，u 是列对应当前时间的矩阵；在稳态情况下，u 是包含 69 个 u 值的单列向量。

1.2 导出网格数据

操作步骤如下：
1. 点击 Mesh 主菜单按钮弹出菜单中的 Export Mesh 选项。
2. 在弹出小面板（图 6.19B）中，选择是否更改默认网格变量名，然后点击 OK 按钮。
3. 在 MATLAB 工作区窗口，会出现三个所选名称的变量：
- 两行矩阵 p（点）：第一行是网格点的 x 坐标，第二行是 y 坐标。
- 六行矩阵 e（边）：包含端点索引、起始和结束参数值、边段和子域编号。
- 四行矩阵 t（三角形）：包含角点索引和子域编号。

实际应用中，p 数据通常较重要，但某些命令可能会用到 e 和 t 参数，建议将所有参数都导入工作区。

导出解和网格节点坐标后，可使用以下命令显示其值：

>> disp('x y u'), d