commit 0e7f2af5d1d0ed4bd24086566eae06a86ff170b6
Author: dmy@lab <dmy@lab.com>
Date:   Thu Aug 20 08:53:06 2015 +0800

    已完成。双回路220kV。
    
    Signed-off-by: dmy@lab <dmy@lab.com>

diff --git a/.gitignore b/.gitignore
new file mode 100644
index 0000000..21229f7
--- /dev/null
+++ b/.gitignore
@@ -0,0 +1 @@
+*.asv
\ No newline at end of file
diff --git a/Field.m b/Field.m
new file mode 100644
index 0000000..b2b253c
--- /dev/null
+++ b/Field.m
@@ -0,0 +1,51 @@
+%%
+%���þ��񷨼���ܿ��߳�ǿ
+%���� Τ�� ��� ��ͬ�˲��ܶ�����·��ĵ糡ǿ�Ⱥ͸�Ӧ��ѹ�� 
+%%
+clc
+clear
+%��ѹ�����������ϵ��£�������ң�
+%˫��
+Volt=[220; 220*exp(1j*240*pi/180); 220*exp(1j*120*pi/180);220; 220*exp(1j*240*pi/180); 220*exp(1j*120*pi/180) ];
+%�γɵ�λϵ������
+H=diag([28.5, 21, 27.5, 28.5, 21, 27.5 ]);
+r=diag([0.2,0.2,0.2,0.2,0.2,0.2]);
+Dcomma=[0,58.5,55,62.42,58.98,55.5;
+        58.5,0,51.5,58.98,55.5,52.1;
+        55, 51.5 0,55.5,52.1 48.7;
+        62.42,58.98,55.55,0,58.5,55;
+        58.98,55.55,52.13,58.5,0,51.5;
+        55.55,52.13,48.73,55,51.5,0];
+D=[   0  3.51 7.02 7.2 8.27 10.4;
+    3.51 0 3.51 8.27 7.8 8.8
+    7.02 3.51 0 10.4 8.8 8.4
+    7.2 8.27 10.48 0 3.5 7.02
+    8.27 7.8 8.8 3.51 0 3.51
+    10.4 8.8 8.4 7.02 3.51 0];
+eslong=1;
+P=1/pi/eslong*log(2*H./r);
+P(isnan(P))=0;
+P2=1/pi/eslong*log(Dcomma./D);
+P2(isnan(P2))=0;
+P=P+P2;
+%����
+QRI=inv(P2)*Volt;
+%��������
+ABCx=[-3.6; -3.9; -4.2; 3.6; 3.9; 4.2];
+ABCy=[31; 27.5; 24; 31; 27.5; 24];
+%����泡ǿ
+plotY=zeros(60,1);
+plotX=zeros(60,1);
+for I=-30:30
+    x=I;
+    y=0;
+    ERy=sum( (ABCy-y)./( (ABCy-y).^2+(ABCx-x).^2 ).*real(QRI)./2/pi/eslong );
+    EIy=sum( (ABCy-y)./( (ABCy-y).^2+(ABCx-x).^2 ).*imag(QRI)./2/pi/eslong );
+    
+    ERx=sum( (ABCx-x)./( (ABCy-y).^2+(ABCx-x).^2 ).*real(QRI)./2/pi/eslong );
+    EIx=sum( (ABCx-x)./( (ABCy-y).^2+(ABCx-x).^2 ).*imag(QRI)./2/pi/eslong );
+    E=sqrt(ERy.^2+EIy.^2+ERx.^2+EIx.^2);
+    plotY(I+31)=E;
+    plotX(I+31)=x;
+end
+plot(plotX,plotY);
\ No newline at end of file