16 lines
577 B
Mathematica
16 lines
577 B
Mathematica
|
function Mat_H=FormH(Busnum,Volt,PG,PD,QG,QD,Y,UAngel,r,c,Angle,Loadi)
|
|||
|
|
%%
|
|||
|
|
%QDcos=textread('D:\Project\<EFBFBD><EFBFBD>С<EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD>\<EFBFBD><EFBFBD>С<EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD>\<EFBFBD>ɺ<EFBFBD>919PDQDglys.txt');
|
|||
|
|
%QD(QD~=0)=PD(QD~=0)./tan(QDcos);
|
|||
|
|
%QD(QD_NON_ZERO_IND)=QD_NON_ZERO;
|
|||
|
|
%%
|
|||
|
|
%PD(Loadi)=QD(Loadi)./tan(acos(0.98));
|
|||
|
|
AngleIJ=sparse(r,c,UAngel(r)-UAngel(c)-Angle,Busnum*3,Busnum*3);
|
|||
|
|
% dP=-sparse(Loadi,1,PD,Busnum*3,1)-diag(Volt)*Y.*cos(AngleIJ)*Volt;
|
|||
|
|
% dQ=-sparse(Loadi,1,QD,Busnum*3,1)-diag(Volt)*Y.*sin(AngleIJ)*Volt;
|
|||
|
|
dP=diag(Volt)*Y.*cos(AngleIJ)*Volt;
|
|||
|
|
dQ=diag(Volt)*Y.*sin(AngleIJ)*Volt;
|
|||
|
|
|
|||
|
|
Mat_H=[dP;dQ;];
|
|||
|
|
Mat_H=sparse(Mat_H);
|
|||
|
|
end
|