From d17ed3632dde98e4b1e273515ebf96201ed282e4 Mon Sep 17 00:00:00 2001
From: "dmy@lab" <dmy@lab.com>
Date: Wed, 23 Sep 2015 18:21:22 +0800
Subject: [PATCH] =?UTF-8?q?1.=E5=8A=A0=E4=BA=86=E8=AE=A1=E7=AE=97=E5=AE=8C?=
 =?UTF-8?q?=E6=88=90=E6=8F=90=E7=A4=BA=202.=E6=94=B9=E4=BA=86=E8=AF=AF?=
 =?UTF-8?q?=E5=B7=AE=E7=BB=9F=E8=AE=A1=E6=96=B9=E6=B3=95?=
MIME-Version: 1.0
Content-Type: text/plain; charset=UTF-8
Content-Transfer-Encoding: 8bit

Signed-off-by: dmy@lab <dmy@lab.com>
---
 SelfAdaptSimulation.m | 12 +++++++-----
 1 file changed, 7 insertions(+), 5 deletions(-)

diff --git a/SelfAdaptSimulation.m b/SelfAdaptSimulation.m
index 90a4393..3fcb2f8 100644
--- a/SelfAdaptSimulation.m
+++ b/SelfAdaptSimulation.m
@@ -18,9 +18,9 @@ arc=2*pi/semi_lineCount;
 %子导线中心到导线中心的距离
 R=semi_lineDistance/2/sin(arc/2);
 %计算模拟电荷的位置
-r1=0;
+r1=-10;
 error=10000;
-for I=1:100
+for Loop=1:10000
     simulationChargePos=ones(semi_lineCount,1);
     for I=1:semi_lineCount
         simulationChargePos(I)=exp(1j*((I-1)*arc+arc/2))*(R+r1);%逆时针转一个角度
@@ -51,7 +51,7 @@ for I=1:100
         break;
     end
     %选检验导线上一个角度
-    vrfRelA=[0.0;2*pi/6;2*pi/3;2*pi/1];%vrf=verify
+    vrfRelA=linspace(0,2*pi)';%vrf=verify
     %计算检验点相对于子导线的位置
     vrfRelPos=exp(1j*vrfRelA)*r(1);
     %移动坐标，使验证的子导线中心和实际子导线中心重合。
@@ -63,7 +63,9 @@ for I=1:100
     Pij=1/2/pi/eslong*log(vrf2MirrorDistance./vrf2ConductorDistance);
     %计算电压
     V=Pij*QRI;
-    error=sum(abs(V-500)./500);
-    r1=r1+r(1)/20;
+    error=sum(abs(V-500)./500)/length(V);
+    r1=r1+r(1)/100;
 end
+display('Finished.');
+display(Loop);
 % scatter(real([simulationChargeAPos;vrfPos]),imag([simulationChargeAPos;vrfPos]));
\ No newline at end of file