使用Python求解电阻立方体等效电阻：代码结果异常求助

问题排查 #

代码的核心逻辑（电导矩阵构建、节点电压法求解）是正确的，但结果错误大概率源于节点连接与实际示意图不符，或者电流注入的节点对应错误：

1. 你在代码中定义的特殊电阻R1连接的是节点3和4，但如果示意图中R1的位置并非此处（比如节点0与7、节点1与6等），会直接导致计算结果偏离预期。
2. 电流注入的节点（0和7）如果与示意图中的A、B端不对应，也会得到错误结果。

修正后的代码 #

以下代码保留了正确的核心逻辑，同时添加了注释便于你根据示意图调整节点连接：

import numpy as np

# 立方体节点总数
num_nodes = 8

# 初始化电导矩阵G和电流向量I
G = np.zeros((num_nodes, num_nodes))
I = np.zeros(num_nodes)

# 定义电阻值
R1 = 0.18  # 特殊电阻R1
R = 1      # 其余电阻值

# 向电导矩阵添加两个节点间的电导
def add_conductance(G, i, j, R_val):
    conductance = 1 / R_val
    # 自电导：节点i和j的对角元素增加电导值
    G[i, i] += conductance
    G[j, j] += conductance
    # 互电导：非对角元素减去电导值（对称）
    G[i, j] -= conductance
    G[j, i] -= conductance

# -------------------------- 关键：根据示意图修正此处的节点连接 --------------------------
edges = [
    (0, 1, R), (1, 2, R), (2, 3, R), (3, 0, R),  # 底面节点连接
    (4, 5, R), (5, 6, R), (6, 7, R), (7, 4, R),  # 顶面节点连接
    (0, 4, R), (1, 5, R), (2, 6, R), (3, 7, R),  # 垂直边节点连接
    (3, 4, R1)                                   # 特殊电阻R1的连接，请根据示意图修改此处的节点对
]

# 填充电导矩阵
for i, j, R_val in edges:
    add_conductance(G, i, j, R_val)

# -------------------------- 根据示意图修正电流注入的节点 --------------------------
# 1A流入节点A，1A流出节点B，请替换为实际的A、B节点编号
I[0] = 1    # 节点A（示例为0）
I[7] = -1   # 节点B（示例为7）

# 接地节点B（移除节点B对应的行和列，避免矩阵奇异）
ground_node = 7
G_reduced = np.delete(np.delete(G, ground_node, axis=0), ground_node, axis=1)
I_reduced = np.delete(I, ground_node)

# 求解节点电压
V_reduced = np.linalg.solve(G_reduced, I_reduced)

# 插入接地节点的电压（0V）
V = np.insert(V_reduced, ground_node, 0)

# 计算等效电阻：A、B节点的电压差
R_eq = V[0] - V[7]  # 替换为实际的A、B节点编号

print(f"A、B之间的等效电阻为 {R_eq:.3f} 欧姆。")

验证方法 #

你可以通过以下方式验证代码正确性：

• 将R1设置为与R相同的值（即所有电阻均为1Ω），如果原立方体A、B端的等效电阻为5/6≈0.833Ω，添加R1后结果应小于该值（比如代码中连接3-4时，结果约为0.583Ω）。
• 将R1设置为极大值（如1e9Ω，模拟开路），结果应趋近于原立方体的等效电阻。

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