facat
/
potential
1
0
Fork 0
Code Issues Pull Requests Releases Wiki Activity

1.加了计算完成提示 

2.改了误差统计方法

Signed-off-by: dmy@lab <dmy@lab.com>
This commit is contained in:
dmy@lab 2015-09-23 18:21:22 +08:00
parent 798c0e1142
commit d17ed3632d
1 changed files with 7 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

12
SelfAdaptSimulation.m
View File

@ -18,9 +18,9 @@ arc=2*pi/semi_lineCount;
%子导线中心到导线中心的距离 %子导线中心到导线中心的距离
R=semi_lineDistance/2/sin(arc/2); R=semi_lineDistance/2/sin(arc/2);
%计算模拟电荷的位置 %计算模拟电荷的位置
r1=0; r1=-10;
error=10000; error=10000;
for I=1:100 for Loop=1:10000
simulationChargePos=ones(semi_lineCount,1); simulationChargePos=ones(semi_lineCount,1);
for I=1:semi_lineCount for I=1:semi_lineCount
simulationChargePos(I)=exp(1j*((I-1)*arc+arc/2))*(R+r1);%逆时针转一个角度 simulationChargePos(I)=exp(1j*((I-1)*arc+arc/2))*(R+r1);%逆时针转一个角度
@ -51,7 +51,7 @@ for I=1:100
break; break;
end 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); vrfRelPos=exp(1j*vrfRelA)*r(1);
%移动坐标，使验证的子导线中心和实际子导线中心重合。 %移动坐标，使验证的子导线中心和实际子导线中心重合。
@ -63,7 +63,9 @@ for I=1:100
Pij=1/2/pi/eslong*log(vrf2MirrorDistance./vrf2ConductorDistance); Pij=1/2/pi/eslong*log(vrf2MirrorDistance./vrf2ConductorDistance);
%计算电压 %计算电压
V=Pij*QRI; V=Pij*QRI;
error=sum(abs(V-500)./500); error=sum(abs(V-500)./500)/length(V);
r1=r1+r(1)/20; r1=r1+r(1)/100;
end end
display('Finished.');
display(Loop);
% scatter(real([simulationChargeAPos;vrfPos]),imag([simulationChargeAPos;vrfPos])); % scatter(real([simulationChargeAPos;vrfPos]),imag([simulationChargeAPos;vrfPos]));