看起来你在构建无向管网的加权邻接矩阵时遇到了维度对齐的问题，咱们先拆解一下你现有代码里的问题，再给出正确的实现方案。

现有代码的核心问题分析 #

第一段代码的问题 #

• weighted_length = abs(weighted_length - weighted_length.'); 逻辑有误：你本来想让矩阵对称，但这个操作会让对称位置的值变成原差值的绝对值，而非直接复制权重到对称位置，属于绕弯路且易出错的写法。
• 最后一行 find(adj) == find(weighted_length); 仅做了无意义的比较——既没有保存结果，且因find返回的索引顺序差异，这个比较根本无法验证边的位置是否对应。
• 潜在风险：如果Length的长度和tril(adj)中非零元素数量不匹配（比如邻接矩阵不是严格对称的无向图结构），会直接抛出维度错误。

第二段代码的问题 #

• 循环内的adj(i)==1是致命错误：i是循环索引，adj(i)取的是邻接矩阵第i个线性元素，而非positions(i)对应的边的位置，这会导致赋值逻辑完全混乱。
• weighted_length = weighted_length - weighted_length'; 会让对称位置的值变成相反数，完全不符合无向图加权矩阵的要求（无向图的(i,j)和(j,i)权重应当一致）。

修正后的加权邻接矩阵实现（针对管道长度） #

首先明确前提：你的管网是无向图，936个节点对应936×936的对称邻接矩阵，1274条管道对应邻接矩阵中1274个唯一的边（下三角或上三角各1274个非零元素），Length数组长度为1274，对应每条唯一管道的长度。

function [weighted_matrix] = CreateWeightedAdjacencyMatrix(edge_weights, adj)
    % edge_weights: 长度为1274的数组，对应每条唯一管道的权重（比如长度、直径组合值）
    % adj: 936×936的邻接矩阵，无向对称
    
    % 初始化加权矩阵
    weighted_matrix = zeros(size(adj));
    
    % 找到下三角的非零位置（对应唯一的边）
    [r, c] = find(tril(adj));
    
    % 提前校验权重数组长度是否匹配，避免维度错误
    if length(edge_weights) ~= length(r)
        error('edge_weights的长度与邻接矩阵中唯一边的数量不匹配！');
    end
    
    % 给无向边的两个对称位置赋值相同权重
    for k = 1:length(r)
        weighted_matrix(r(k), c(k)) = edge_weights(k);
        weighted_matrix(c(k), r(k)) = edge_weights(k);
    end
end

代码说明 #

1. 用tril(adj)提取下三角非零元素，得到每条唯一边的行、列索引，确保数量为1274，和edge_weights（如Length）长度严格匹配。
2. 循环给每条边的两个对称位置赋相同权重，保证矩阵对称，贴合无向管网的特性。
3. 添加了长度校验逻辑，提前抛出错误，方便快速定位维度问题。

扩展到多因子加权（结合管道+节点属性） #

如果要融合管道长度、直径，以及节点高程、压力差等多因子生成权重，可以先计算每条边的综合权重，再传入上述函数。示例如下：

% 假设你已准备好以下变量：
% Length: 1274×1的管道长度数组
% Diameter: 1274×1的管道直径数组
% NodeElevation: 936×1的节点高程数组
% NodePressure: 936×1的节点压力数组

% 获取每条唯一边的节点对
[r, c] = find(tril(adj));

% 计算综合权重（公式可根据业务需求调整）
elevation_diff = abs(NodeElevation(r) - NodeElevation(c));
pressure_diff = abs(NodePressure(r) - NodePressure(c));
composite_weights = (Length ./ Diameter) + elevation_diff + pressure_diff;

% 生成多因子加权邻接矩阵
weighted_matrix = CreateWeightedAdjacencyMatrix(composite_weights, adj);

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