diff --git a/run.m b/run.m
index 8221dd7..726a2d9 100644
--- a/run.m
+++ b/run.m
@@ -426,7 +426,7 @@ SEThreeLoad=SEVoltpABC(:,Loadi).*conj(-SEIpABC);
 % barPlot( phaseLoadQError, 10,'相对误差%','分布密度','无功负荷误差');
 %% 计算统计量
 %目标函数均值
-rVoltABCV=conj([VoltpA';VoltpC';VoltpC']);
+rVoltABCV=conj([VoltpA';VoltpB';VoltpC']);
 mVoltABCV=conj([mVoltpA';mVoltpB';mVoltpC']);
 rPD3P=real([ phaseASpotLoadP'+1j*phaseASpotLoadQ';
     phaseBSpotLoadP'+1j*phaseBSpotLoadQ';