From 4daca01ca06da3c2a300b0fa781d157d7cfbaac2 Mon Sep 17 00:00:00 2001
From: "dugg@lab-desk" <dugg@lab-desk>
Date: Thu, 9 Oct 2014 09:05:33 +0800
Subject: [PATCH] =?UTF-8?q?=E5=87=86=E5=A4=87=E6=8D=A2=E7=94=A8=E5=86=85?=
 =?UTF-8?q?=E7=82=B9=E6=B3=95=E6=9D=A5=E5=81=9A?=
MIME-Version: 1.0
Content-Type: text/plain; charset=UTF-8
Content-Transfer-Encoding: 8bit

Signed-off-by: dugg@lab-desk <dugg@lab-desk>
---
 run.m | 18 +++++++++++-------
 1 file changed, 11 insertions(+), 7 deletions(-)

diff --git a/run.m b/run.m
index 6feda3e..c2bf845 100644
--- a/run.m
+++ b/run.m
@@ -92,22 +92,26 @@ while(k<=kmax && maxD> EPS)
     %形成负荷序电流的测量值
     mIf0=If0;
     mIf1=If1;
-    mIf1(1)=0.02;
-    mIf1(3)=0.03;
+    mIf1(3)=-mIf1(2);
+%     mIf1(1)=0.02;
+%     mIf1(3)=0.03;
     Loadi=[ 1 2 3];
     mIf2=If2;
     %计算
     fsY11=fsY11+sparse(Balance,Balance,ones(length(Balance),1),busNum,busNum);%这里要置0，置1，否则是奇异的
     %%做最小二乘法
     %先做正序的
-    Z=[%这里要加-号，因为用的Z是注入电流
-        -real(mIf1(Loadi));
-        -imag(mIf1(Loadi));
+    Z=[
+        -real(mIf1(Loadi));%这里要加-号，因为用的Z是注入电流
+        -imag(mIf1(Loadi));%这里要加-号，因为用的Z是注入电流
+        %加电压量测
+        [1 1 1]';
         ];
     H=fsY11(Loadi,:);
     H=[
-        -real(H),imag(H);
-        -imag(H),-real(H);
+        -real(H),imag(H),;
+        -imag(H),-real(H),;
+        eye(3),eye(3);
         ];
     %     J=fsY11(Loadi,:);
     %     J=[