1主要内容

该程序参考文献《交直流系统潮流计算及相互关联特性分析》，采用5种交直流潮流控制方式：1.定电流定电压 2.定电流定熄弧角 3.定功率定电压 4.定功率定熄弧角 5.定触发角定电流。以9节点系统为例，采用牛拉法进行交直流潮流计算，程序为matlab代码，注释清楚！

程序参考流程图

2 部分代码

%% 计算导纳矩阵
for m=1:line
    a=Line.con(m,1);%首节点a
    b=Line.con(m,2);%末节点b
    z=Line.con(m,8)+1i*Line.con(m,9);%节点ab间的阻抗
    y2=1i*Line.con(m,10)/2;%导纳/2
    if real(z)==0
        k=Line.con(m,7)*Bus.con(b,2)/Bus.con(a,2);%变比K
    else
        k=1;
    end
    Y(a,b)=-1/k/z;%ab间互导纳
    Y(b,a)=Y(a,b);
    Y(a,a)=Y(a,a)+y2+(k-1)/(k*z)+1/k/z;%求自导纳
    Y(b,b)=Y(b,b)+y2+(1-k)/k^2/z+1/k/z;
%     C(a)=data(m,7);%输入节点a,b的补偿电容导纳
%     C(b)=data(m,8);
end
 
%将节点导纳的实虚部分开
B=imag(Y);
G=real(Y);
 
%判断PQ、PV、平衡节点的个数
PQn=size(PQ.con,1);
PVn=size(PV.con,1);
SWn=size(SW.con,1);

DCn=Hvdc.con(:,[1 2])';%直流节点编号
ACn=setdiff(1:n,union(DCn,SW.con(1,1)))';
Ndc=size(DCn,1);%直流节点个数
Nac=size(ACn,1);%交流节点个数
pn=setdiff(1:n,SW.con(1,1))';
qn=setdiff(pn,PV.con(:,1));
Pn=intersect(ACn, pn);%交流节点P方程节点
Qn=intersect(ACn, qn);%交流节点Q方程节点
p=size(Pn,1);   %交流节点P方程个数
q=size(Qn,1);   %交流节点Q方程个数
% YK=zeros(p+q);
DPa=zeros(p,1);
DQa=zeros(q,1);
DPt=zeros(Ndc,1);
DQt=zeros(Ndc,1);
Dd1=zeros(Ndc,1);
Dd2=zeros(Ndc,1);
Dd3=zeros(Ndc,1);
Dd4=zeros(Ndc,1);
Dd5=zeros(1,1);
H=zeros(p+Ndc);
N=zeros(p+Ndc,q+Ndc);
M=zeros(q+Ndc,p+Ndc);
L=zeros(q+Ndc);
JRX=zeros(9,9);
%% 直流迭代初值设置+控制方式
Vd=ones(Ndc,1)*1.1562;%直流电压初值
Id=0.5;
% fai=ones(Ndc,1)*30/180*pi;%功率因数角初值30°
Pd=0.6.*ones(Ndc,1);
Qd=0.3.*ones(Ndc,1);
Sd=0.5.*ones(Ndc,1);
Control_ang=[20/180*pi;17/180*pi];%换流器控制角，第一个为alpha，第二个为gama