diff --git a/IPMLoop.m b/IPMLoop.m
index 364f926..84fae39 100644
--- a/IPMLoop.m
+++ b/IPMLoop.m
@@ -2,8 +2,8 @@ function [ V1r,V1i,I1r,I1i,isConverged ] = IPMLoop(V1measurement,wV1r,wV1i,I1mea
 %把每个序的循环写在这个函数中。其实也就是内点法循环。
 V1r=1*ones(busNum,1);
 V1i=0*ones(busNum,1);
-I1r=0.0*ones(length(Loadi),1);%注入电流，相当于放电机电流，与负荷电流负荷相反
-I1i=0.0*ones(length(Loadi),1);%注入电流，相当于放电机电流，与负荷电流负荷相反
+I1r=0.1*ones(length(Loadi),1);%注入电流，相当于放电机电流，与负荷电流负荷相反
+I1i=0.1*ones(length(Loadi),1);%注入电流，相当于放电机电流，与负荷电流负荷相反
 KK=0;
 plotGap=zeros(1,60);
 %初始化
diff --git a/feeder13/data1.txt b/feeder13/data1涓夌浉涓嶅钩琛�.txt
similarity index 100%
rename from feeder13/data1.txt
rename to feeder13/data1涓夌浉涓嶅钩琛�.txt
diff --git a/run.m b/run.m
index d69405e..9a5ef51 100644
--- a/run.m
+++ b/run.m
@@ -269,15 +269,33 @@ for II=1:3
     
     if II==1
         fprintf('正序\n');
+        wV1r(1:2)=1e10;
+        wV1i(1:2)=1e10;
+        wI1r(1:2)=1e10;
+        wI1i(1:2)=1e10;
+        tic
         [ V1r,V1i,I1r,I1i,isConverged1 ]=IPMLoop(V1measurement,wV1r,wV1i,I1measurement,wI1r,wI1i,BalI1r,BalI1i,busNum,Loadi,fsY11,Balance,1 );
+        toc
     end
     if II==2
         fprintf('负序\n');
+        wV2r(1:2)=1e10;
+        wV2i(1:2)=1e10;
+        wI2r(1:2)=1e10;
+        wI2i(1:2)=1e10;
+        tic
         [ V2r,V2i,I2r,I2i,isConverged2 ]=IPMLoop(V2measurement,wV2r,wV2i,I2measurement,wI2r,wI2i,BalI2r,BalI2i,busNum,Loadi,fsY22,Balance,0 );
+        toc
     end
     if II==3
         fprintf('零序\n');
+        wV0r(1:2)=1e10;
+        wV0i(1:2)=1e10;
+        wI0r(1:2)=1e10;
+        wI0i(1:2)=1e10;
+        tic
         [ V0r,V0i,I0r,I0i,isConverged0 ]=IPMLoop(V0measurement,wV0r,wV0i,I0measurement,wI0r,wI0i,BalI0r,BalI0i,busNum,Loadi,fsY00,Balance,0 );
+        toc
     end
     
 end