MATLAB通过输入文件进行水准网平差的代码

以下是MATLAB实现水准网平差的代码示例：

% 读取数据文件，数据格式为 station_name, height, observation_name, measured_height, standard_deviation
data = readtable('data.txt');

% 构建A、L矩阵
num_stations = numel(unique(data.station_name));
num_observations = size(data, 1);
A = zeros(num_observations, num_stations);
L = zeros(num_observations, 1);
station_names = unique(data.station_name);
for i = 1:num_observations
    station_index = find(strcmp(station_names, data.station_name{i}));
    A(i, station_index) = 1;
    L(i) = data.measured_height(i) - data.height(i);
end

% 计算系数矩阵和常数向量
P = diag(1 ./ (data.standard_deviation .^ 2));
ATPA = A' * P * A;
ATPL = A' * P * L;

% 解方程 Ax = b
x = ATPA \ ATPL;

% 计算高程值
station_heights = x + mean(data.height);

其中，data.txt为输入文件，第一列为测站名称，第二列为测站高程，第三列为观测点名称，第四列为观测点高程，第五列为观测点高程的标准差。执行完上述代码后，station_heights为求解得到的各个测站高程值。