From 286caf83a4c5b70b5bdf6cc7939cbf079fe14a6c Mon Sep 17 00:00:00 2001
From: "dmy@lab" <dmy@lab.lab>
Date: Tue, 7 Apr 2015 10:39:04 +0800
Subject: [PATCH] =?UTF-8?q?=E5=8A=A0=E4=B8=80=E4=B8=AA=E5=B0=8F=E7=9A=84?=
 =?UTF-8?q?=E7=90=86=E8=AE=BA=E9=AA=8C=E8=AF=81=E4=BB=A3=E7=A0=81=E7=89=87?=
 =?UTF-8?q?=E6=AE=B5?=
MIME-Version: 1.0
Content-Type: text/plain; charset=UTF-8
Content-Transfer-Encoding: 8bit

Signed-off-by: dmy@lab <dmy@lab.lab>
---
 theory/t.m | 23 +++++++++++++++++++++++
 1 file changed, 23 insertions(+)
 create mode 100644 theory/t.m

diff --git a/theory/t.m b/theory/t.m
new file mode 100644
index 0000000..92644e5
--- /dev/null
+++ b/theory/t.m
@@ -0,0 +1,23 @@
+clc
+clear
+
+P=[12;16;15];
+mP=P.*(1+normrnd(0,0.03,length(P),1));
+Q=[7;15;10];
+mQ=Q.*(1+normrnd(0,0.03,length(Q),1));
+rVolt=[1;
+    1*exp(1j*-123/180*pi);
+    1*exp(1j*119/180*pi);
+    ];
+mVolt=[1*(1+normrnd(0,0.03));
+    1*(1+normrnd(0,0.03))*exp(1j*-120/180*pi);
+    1*(1+normrnd(0,0.03))*exp(1j*120/180*pi);
+    ];
+a=exp(1j*2*pi/3);
+Tp2f=1/3*[1 1 1;
+    1 a a^2;
+    1 a^2 a];
+rCurr=(P+1j*Q)./rVolt;
+mCurr=(mP+1j*mQ)./mVolt;
+rI=Tp2f*rCurr
+mI=Tp2f*mCurr
\ No newline at end of file