From 3a73dfd84e55e0369d0f5f2b9a0df9dd4a8bb188 Mon Sep 17 00:00:00 2001
From: "dmy@lab" <dmy@lab.com>
Date: Thu, 24 Sep 2015 16:20:22 +0800
Subject: [PATCH] =?UTF-8?q?=E6=B7=BB=E5=8A=A0=E7=94=B5=E5=9C=BA=E8=AE=A1?=
 =?UTF-8?q?=E7=AE=97?=
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 | 31 +++++++++++++++++++------------
 1 file changed, 19 insertions(+), 12 deletions(-)

diff --git a/SelfAdaptSimulation.m b/SelfAdaptSimulation.m
index 7cbf138..b1c0979 100644
--- a/SelfAdaptSimulation.m
+++ b/SelfAdaptSimulation.m
@@ -5,16 +5,16 @@ clear
 % [2].	陈习文, 特高压直流输电线路电磁环境的研究, 2012, 北京交通大学.
 %%
 %设置几个参数
-semi_lineDistance=257;%分裂间距
+semi_lineDistance=457;%分裂间距
 semi_lineCount=4;%分裂数
-ConductorX=[-14500,14500];%导线距地高度
-ConductorY=[16500,16500];%导线间距
+ConductorX=[-13720,0,13720];%导线距地高度
+ConductorY=[20830,20830,20830];%导线间距
 CSM_N=80;%每一个子导线的模拟电荷数
-subconductorR=30;%子导线半径
-phaseN=2;%相数，单回三相
+subconductorR=29.95;%子导线半径
+phaseN=3;%相数，单回三相
 %%
 %设置电压
-Volt_=[500;-500;];
+Volt_=[765/sqrt(3);765/sqrt(3)*exp(1j*4/3*pi);765/sqrt(3)*exp(1j*2/3*pi);];
 Volt=[];
 for vLoop=1:length(Volt_)
     Volt=[Volt;Volt_(vLoop)*ones(CSM_N*semi_lineCount,1);];
@@ -96,12 +96,19 @@ if Loop<maxLoop
     display('Converged.');
 end
 display(Loop);
+%计算场强
+ABCy=imag(repmat(simulationChargeABCPos,1,length(vrfPos)));
+ABCx=real(repmat(simulationChargeABCPos,1,length(vrfPos)));
+y=imag(conj(matVrfPos'));
+x=real(conj(matVrfPos'));
+ERy=sum( ( (ABCy-y)./( (ABCy-y).^2+(ABCx-x).^2 )+(ABCy+y)./( (ABCy+y).^2+(ABCx-x).^2 ) ).*repmat(real(QRI),1,length(vrfPos))./2/pi/eslong,1 );
+EIy=sum( ( (ABCy-y)./( (ABCy-y).^2+(ABCx-x).^2 )+(ABCy+y)./( (ABCy+y).^2+(ABCx-x).^2 ) ).*repmat(imag(QRI),1,length(vrfPos))./2/pi/eslong,1 );  
+ERx=sum( ( (ABCx-x)./( (ABCy-y).^2+(ABCx-x).^2 )+(ABCx-x)./( (ABCy+y).^2+(ABCx-x).^2 ) ).*repmat(real(QRI),1,length(vrfPos))./2/pi/eslong,1 );
+EIx=sum( ( (ABCx-x)./( (ABCy-y).^2+(ABCx-x).^2 )+(ABCx-x)./( (ABCy+y).^2+(ABCx-x).^2 ) ).*repmat(imag(QRI),1,length(vrfPos))./2/pi/eslong,1 );
+E=sqrt(ERy.^2+EIy.^2+ERx.^2+EIx.^2);
+max(E);
+
 scatter(real(simulationChargeABCPos(1:length(simulationChargeABCPos)/1)),imag(simulationChargeABCPos(1:length(simulationChargeABCPos)/1)),[],'r');
 axis equal
 hold on;
-scatter(real(vrfPos),imag(vrfPos),[],'k');
-% scatter(real([simulationChargeAPos(1:10);]),imag([simulationChargeAPos(1:10);]),10,'red');
-% hold;
-% scatter(real([simulationChargeAPos(11:20);]),imag([simulationChargeAPos(11:20);]),10,'blue');
-% scatter(real([simulationChargeAPos(21:30);]),imag([simulationChargeAPos(21:30);]),10,'green');
-% scatter(real([simulationChargeAPos(31:40);]),imag([simulationChargeAPos(31:40);]),10,'yellow');
\ No newline at end of file
+scatter(real(vrfPos),imag(vrfPos),[],'k');
\ No newline at end of file