diff --git a/run.m b/run.m index 05910a9..648c014 100644 --- a/run.m +++ b/run.m @@ -181,8 +181,8 @@ mIf2=-mIf2; %% 设定量测误差 %% 先算正序的 %平衡节点电流 -BalI1r=real(-sum(rIf1)); -BalI1i=imag(-sum(rIf1)); +BalI1r=real(-sum(mIf1)); +BalI1i=imag(-sum(mIf1)); %电压 %制作量测量 % mfsY11=fsY11; @@ -210,8 +210,8 @@ wI1i=abs(imag(I1measurement)).*sigma; % f=sum(([real(I1measurement);imag(I1measurement)]-[I1r;I1i]).^2)+sum((real(rV1)-V1r).^2)+sum((imag(rV1)-V1i).^2); % fprintf('目标值 %f\n',full(f)); %% 算负序的 -BalI2r=real(-sum(rIf2)); -BalI2i=imag(-sum(rIf2)); +BalI2r=real(-sum(mIf2)); +BalI2i=imag(-sum(mIf2)); %电压 %制作量测量 % mfsY22=fsY22; @@ -238,8 +238,8 @@ wI2i=abs(imag(I2measurement)).*sigma; % f=sum(([real(I2measurement);imag(I2measurement)]-[I2r;I2i]).^2)+sum((real(rV2)-V2r).^2)+sum((imag(rV2)-V2i).^2); % fprintf('目标值 %f\n',full(f)); %% 算零序 -BalI0r=real(-sum(rIf0)); -BalI0i=imag(-sum(rIf0)); +BalI0r=real(-sum(mIf0)); +BalI0i=imag(-sum(mIf0)); %电压 %制作量测量 % mfsY00=fsY00;