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