如何从PDE求解结果中提取等值面？MATLAB实现时内存不足问题咨询

• 你的使用方法确实存在问题：meshgrid的作用是生成规则正交网格矩阵，你将非规则分布的有限元节点坐标直接传入该函数，会生成维度等于节点数三次方的巨型矩阵，自然会出现内存不足的问题。
• 你采用的有限元求解得到的温度场是定义在无结构四面体网格上的，完全不需要转换为规则网格就可以直接提取等值面。

方案1：使用PDE工具箱内置pdeIsosurface函数（最简便，推荐） #

该函数是专门为PDE求解结果设计的等值面提取工具，不需要手动处理网格参数，代码示例如下：

% 直接提取温度为100的等值面
iso = pdeIsosurface(thermalresults, 100);

% 可选：绘制等值面验证结果
patch(iso, 'FaceColor', 'red', 'EdgeColor', 'none')
axis equal
view(3)
camlight
lighting gouraud

返回的iso结构体包含Faces和Vertices字段，兼容MATLAB的patch绘图语法，也可以直接导出为其他格式使用。

方案2：调用无结构网格版本的isosurface函数 #

如果需要更底层的控制，可以直接使用isosurface针对四面体网格的专属语法，不需要生成规则网格：

% 读取网格单元、节点坐标、温度数据
tetra_elements = mesh.Elements;
x = mesh.Nodes(1,:);
y = mesh.Nodes(2,:);
z = mesh.Nodes(3,:);
temperature = thermalresults.Temperature;

% 直接对无结构网格提取100℃等值面
[faces, vertices] = isosurface(x, y, z, tetra_elements, temperature, 100);

% 可选：绘制验证
patch('Faces', faces, 'Vertices', vertices, 'FaceColor', 'blue', 'EdgeColor', 'none')
axis equal
view(3)
camlight
lighting gouraud

该语法直接处理无结构四面体网格，内存占用只和你的原始网格规模相关，完全不会出现内存爆炸的问题。

内容的提问来源于stack exchange，提问作者JackChen