facat
/
distributionnetwork-power2current
1
0
Fork 0
Code Issues Pull Requests Releases Wiki Activity

添加把结果转换会三相值 

Signed-off-by: dugg@lab-desk <dugg@lab-desk>
This commit is contained in:
dugg@lab-desk 2014-10-29 21:26:39 +08:00
parent 098fe0e7c1
commit 7409723aae
1 changed files with 12 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

13
run.m
View File

@ -298,8 +298,8 @@ sigmaI0=normrnd(0,sigma,length(Loadi),1);
I0measurement=mIf0(Loadi).*(1+sigmaI0);%测量值是等效发电机电流 I0measurement=mIf0(Loadi).*(1+sigmaI0);%测量值是等效发电机电流
wI0r=abs( real(mIf0(Loadi).*sigmaI0) ); wI0r=abs( real(mIf0(Loadi).*sigmaI0) );
wI0i=abs( imag(mIf0(Loadi).*sigmaI0) ); wI0i=abs( imag(mIf0(Loadi).*sigmaI0) );
[ V0r,V0i,I0r,I0i ]=IPMLoop(V0measurement,wV0r,wV0i,I0measurement,wI0r,wI0i,BalI0r,BalI0i,busNum,Loadi,fsY00,Balance,0 ); [ V0r,V0i,I0r,I0i ]=IPMLoop(V0measurement,wV0r,wV0i,I0measurement,wI0r,wI0i,BalI0r,BalI0i,busNum,Loadi,fsY00,Balance,0 );
%%
% f=sum(([real(Imeasurement);imag(Imeasurement)]-[-I1r;-I1i]).^2)+sum((real(Vmeasurement)-V0r(Loadi)).^2)+sum((imag(Vmeasurement)-V0i(Loadi)).^2); % f=sum(([real(Imeasurement);imag(Imeasurement)]-[-I1r;-I1i]).^2)+sum((real(Vmeasurement)-V0r(Loadi)).^2)+sum((imag(Vmeasurement)-V0i(Loadi)).^2);
% fprintf('目标值 %f\n',full(f)); % fprintf('目标值 %f\n',full(f));
%状态量 %状态量
@ -319,3 +319,14 @@ wI0i=abs( imag(mIf0(Loadi).*sigmaI0) );
% SEQD=sparse(zeros(busNum,1)); % SEQD=sparse(zeros(busNum,1));
%检查目标函数 %检查目标函数
%% 把三序合成三相
% 三相电压
SEVf0=V0r+1j*V0i;
SEVf1=V1r+1j*V1i;
SEVf2=V2r+1j*V2i;
SEVoltpABC=Tf2p*conj([ SEVf0'; SEVf1'; SEVf2']);
% 三序电流
SEIf0=I0r+1j*I0r;
SEIf1=I1r+1j*I1r;
SEIf2=I2r+1j*I2r;
SEIpABC=Tf2p*conj([SEIf0';SEIf1';SEIf2']);