From f9a3b1cd1a2ffda43fe2d46a0d83a88a7cc6e693 Mon Sep 17 00:00:00 2001
From: facat <dugg@21cn.com>
Date: Wed, 10 Apr 2013 17:01:12 +0800
Subject: [PATCH] =?UTF-8?q?=E5=A2=9E=E5=8A=A0=E4=BA=86=E6=94=AF=E8=B7=AF?=
 =?UTF-8?q?=E7=94=B5=E6=B5=81?=
MIME-Version: 1.0
Content-Type: text/plain; charset=UTF-8
Content-Transfer-Encoding: 8bit

Signed-off-by: facat <dugg@21cn.com>
---
 run.asv | 13 +++++++++++++
 run.m   | 31 ++++++++++++++++++++++++++++++-
 2 files changed, 43 insertions(+), 1 deletion(-)
 create mode 100644 run.asv

diff --git a/run.asv b/run.asv
new file mode 100644
index 0000000..8ff07d4
--- /dev/null
+++ b/run.asv
@@ -0,0 +1,13 @@
+addpath('.\Powerflow')
+[~, ~, ~, ~,Volt,Vangle,Y,Yangle,newwordParameter]=pf('ieee4.dat', '0');
+%% 量测量
+% 电压 节点电流 支路 节点功率 支路功率
+%%
+%% 状态量
+% 电压 相角
+%%
+%% 开始生成量测量
+sigma=0.05;% 标准差
+%电压
+mVolt=Volt.*(normrnd(0,sigma,length(Volt),1)+1);
+
diff --git a/run.m b/run.m
index 2125e74..5cfa25b 100644
--- a/run.m
+++ b/run.m
@@ -1,2 +1,31 @@
 addpath('.\Powerflow')
-pf ieee4.dat 0;
+[~, ~, ~, ~,Volt,Vangle,Y,Yangle,r,c,newwordParameter]=pf('ieee4.dat', '0');
+%% 量测量
+% 电压 节点电流 支路电流 节点功率 支路功率
+%%
+%% 状态量
+% 电压 相角
+%%
+%% 开始生成量测量
+sigma=0.05;% 标准差
+%% 电压
+%电压幅值
+rVolt=Volt; %幅值
+mVolt=rVolt.*(normrnd(0,sigma,length(Volt),1)+1);%电压量测量
+%% 电流
+%注入电流
+cmpY=Y.*exp(1j*sparse(r,c,Yangle,length(Y),length(Y)));%复数导纳矩阵
+cmpV=Volt.*exp(1j*Vangle); %复数电压
+cmpI=cmpY*cmpV;% 注入电流
+rI=abs(cmpI); %注入电流量测量要的是电流幅值
+mI=rI.*(normrnd(0,sigma,length(rI),1)+1);%电流量测量
+%% 支路电流
+% 支路电流
+lineI=newwordParameter.line.lineI;
+lineJ=newwordParameter.line.lineJ;
+lineR=newwordParameter.line.lineR;
+lineX=newwordParameter.line.lineX;
+lineB2=newwordParameter.line;
+cmpBranchI=(cmpV(lineI)-cmpV(lineJ))./(lineR+1j*lineX);%复数支路电流
+rBranchI=abs(cmpBranchI);% 支路电流幅值
+mBranchI=rBranchI.*(normrnd(0,sigma,length(rBranchI),1)+1);%支路电流量测量
\ No newline at end of file