如何用matlab实现：假设三角网的节点坐标存储在一个 Nx2 的矩阵 node 中，每一行代表一个节点的横纵坐标。三角形的拓扑关系存储在一个 Mx3 的矩阵 tri 中，每一行代表一个三角形的三个节点的编号。 % 假设三角网的真实高程数据存储在一个长度为 N 的一维数组 ztrue 中，三角网的计算高程数据存储在一个长度为 N 的一维数组 zcalc 中。

以下是一个简单的 Matlab 函数，用于计算三角网的均方根误差（RMSE）：

function rmse = triangle_rmse(node, tri, ztrue, zcalc)
% 计算三角网的均方根误差（RMSE）
% 输入参数：
% node: Nx2 的矩阵，每一行代表一个节点的横纵坐标
% tri: Mx3 的矩阵，每一行代表一个三角形的三个节点的编号
% ztrue: 长度为 N 的一维数组，代表三角网的真实高程数据
% zcalc: 长度为 N 的一维数组，代表三角网的计算高程数据
% 输出参数：
% rmse: 三角网的均方根误差

% 计算每个三角形的面积
tri_area = zeros(size(tri, 1), 1);
for i = 1:size(tri, 1)
    node1 = node(tri(i, 1), :);
    node2 = node(tri(i, 2), :);
    node3 = node(tri(i, 3), :);
    tri_area(i) = 0.5 * abs(det([node2-node1; node3-node1]));
end

% 计算三角网的均方根误差
rmse = sqrt(mean((zcalc - ztrue).^2));
end

该函数首先计算每个三角形的面积，然后计算三角网的均方根误差。使用该函数可以如下调用：

% 假设三角网的节点坐标存储在一个 Nx2 的矩阵 node 中，每一行代表一个节点的横纵坐标。
% 三角形的拓扑关系存储在一个 Mx3 的矩阵 tri 中，每一行代表一个三角形的三个节点的编号。
% 假设三角网的真实高程数据存储在一个长度为 N 的一维数组 ztrue 中，三角网的计算高程数据存储在一个长度为 N 的一维数组 zcalc 中。

rmse = triangle_rmse(node, tri, ztrue, zcalc);

其中，参数 node、tri、ztrue、zcalc 分别对应三角网的节点坐标、拓扑关系、真实高程数据和计算高程数据。函数返回的 rmse 即为三角网的均方根误差。