0
收藏
微博
微信
复制链接
【PCB考证】 零基础学PCB 嵌入式干货 100G电子设计资料 阻抗下载神器 学员评价 广告合作 校企合作申请

基于交替迭代法的交直流混合系统潮流计算

电磁MATLAB
2024-07-25 16:34
364

由于受端负荷量持续增大，直流输电受端系统的电压稳定性能自然引起了人们的关注。目前，交直流混合电网的电压稳定分析大致分为两类，即静态电压稳定分析和动态电压稳定分析，本聚焦于前者。传统分析中，有效短路比评价交直流电网的强度，ＶＳＩ指标评价受端电网电压稳定性，但两者均未能揭示受端电网电压失稳的机理。

1 交直流系统的等值模型

如图所示为交直流系统等值的示意

72f04c16214a17ae6692ee4a7fad06.jpg

其基本特性方程可表示为式

9243c47769b4e220f5340ec3466a73.jpg

a7675d61ad03f46b7f715170a45758.jpg

cf07bcd687448be6f32ef8ea766e6e.jpg


1.1 整流器定电流逆变器定电压（ＣＣ－ＣＵ）

b165775fa28bb4c1e06dccde3cb0ee.jpg


1.2 整流器定电流逆变器定关断角（ＣＣ－ＣＩＡ）

d2e2cc8904a9be38c0ad9f23c7b36d.jpg


1.3 整流器定功率逆变器定电压（ＣＰ－ＣＵ）

895e49696757f574c15df6c7c32108.jpg


1.4 整流器定功率逆变器定关断角（ＣＰ－ＣＩＡ）

8bc9f453071ec7f3f756a6f41b9446.jpg


2 iEEE9节点系统算例

7与8节点改为直流支路。

5addd31656f557dca2c41eda196181.jpg


3 matlab程序运行结果

a52188aeaa177a15eb8a9a970b543b.jpg


4 matlab程序

1）主函数


%%  5中控制模式，直流未知数为UD ID SD PD QD Contrl_ang 参考王云鹏文章附件控制方式
%function flow_calculationclear;clc;close allticeval('case9')T=0.00001;%% 参数kk=1;%迭代次数line=size(Line.con,1);%支路数目n=max(Bus.con(:,1));%找出第1,2列的最大值，即节点数TT=ones(n,1)*T;%收敛误差设置%% 创建导纳矩阵YY=zeros(n);%创建导纳矩阵YU=ones(n,1);dU=zeros(n,1);U(PV.con(:,1))=PV.con(:,5);U(SW.con(:,1))=SW.con(:,4);th=zeros(n,1);dth=zeros(n,1);Pg=zeros(n,1);Pl=zeros(n,1);Ql=zeros(n,1);Qg=zeros(n,1);Pg(PV.con(:,1))=PV.con(:,4);Pl(PQ.con(:,1))=PQ.con(:,4);Ql(PQ.con(:,1))=PQ.con(:,5);
%% 计算导纳矩阵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;间互导纳    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% Control_ang(2)=30/180*pi;W=cos(Control_ang);Kt=ones(Ndc,1);%直流支路两端变压器变比值设置
%% 直流参数设置% ku=0.995;%换相效应常数，王锡凡181页Rdc=0.0625;%直流电阻Xc=[0.1345;0.1257];%整流及逆变侧变压器电抗，写入真实值disp('请输入控制方式:')disp('方式1：整流侧定电流，逆变侧定电压')disp('方式2：整流侧定电流，逆变侧定熄弧角')disp('方式3：整流侧定功率，逆变侧定电压')disp('方式4：整流侧定功率，逆变侧定熄弧角')disp('方式5：整流侧定触发角，逆变侧定电流')Controlmode=input('方式：');
%% 交直流计算主循环for mm=1:1000        H=zeros(p Ndc);    N=zeros(p Ndc,q Ndc);    M=zeros(q Ndc,p Ndc);    L=zeros(q Ndc);
%求ΔPa，ΔQa,ΔPt,ΔQtPu=0;Qu=0;for m=1:p    for k=1:n        Pu=Pu U(Pn(m))*U(k)*(G(Pn(m),k)*cos(th(Pn(m))-th(k)) B(Pn(m),k)*sin(th(Pn(m))-th(k)));    end    DPa(m)=Pg(Pn(m))-Pl(Pn(m))-Pu;    Pu=0;endfor m=1:q    for k=1:n        Qu=Qu U(Qn(m))*U(k)*(G(Qn(m),k)*sin(th(Qn(m))-th(k))-B(Qn(m),k)*cos(th(Qn(m))-th(k)));    end    DQa(m)=Qg(Qn(m))-Ql(Qn(m))-Qu;    Qu=0;end for m=1:Ndc    for k=1:n        Pu=Pu U(DCn(m))*U(k)*(G(DCn(m),k)*cos(th(DCn(m))-th(k)) B(DCn(m),k)*sin(th(DCn(m))-th(k)));        Qu=Qu U(DCn(m))*U(k)*(G(DCn(m),k)*sin(th(DCn(m))-th(k))-B(DCn(m),k)*cos(th(DCn(m))-th(k)));    end    if any(DCn(m)==Hvdc.con(:,1))        DPt(m)=Pg(DCn(m))-Pl(DCn(m))-Pu-Pd(m);        DQt(m)=Qg(DCn(m))-Ql(DCn(m))-Qu-Qd(m);    else        DPt(m)=Pg(DCn(m))-Pl(DCn(m))-Pu Pd(m);        DQt(m)=Qg(DCn(m))-Ql(DCn(m))-Qu-Qd(m);%%再仔细想想直流无功的正负号    end    Pu=0;    Qu=0; end
%矩阵H的形成for x=1:(p Ndc)    for y=1:(p Ndc)        if pn(x)==pn(y)            for m=1:n                H(x,x)=H(x,x) U(pn(x))*U(m)*(G(pn(x),m)*sin(th(pn(x))-th(m))-B(pn(x),m)*cos(th(pn(x))-th(m)));            end                H(x,x)=H(x,x)-U(pn(x))*U(pn(x))*(G(pn(x),pn(x))*sin(th(pn(x))-th(pn(x)))-B(pn(x),pn(x))*cos(th(pn(x))-th(pn(x))));        else            H(x,y)=-U(pn(x))*U(pn(y))*(G(pn(x),pn(y))*sin(th(pn(x))-th(pn(y)))-B(pn(x),pn(y))*cos(th(pn(x))-th(pn(y))));        end    endendHaa=H(Pn-1,Pn-1);Hat=H(Pn-1,DCn-1);Hta=H(DCn-1,Pn-1);Htt=H(DCn-1,DCn-1);H=[Haa,Hat;Hta,Htt];
%矩阵N的形成for x=1:(p Ndc)    for y=1:(q Ndc)        if pn(x)==qn(y)            for m=1:1:n                N(x,y)=N(x,y)-U(pn(x))*U(m)*(G(pn(x),m)*cos(th(pn(x))-th(m)) B(pn(x),m)*sin(th(pn(x))-th(m)));            end                N(x,y)=N(x,y) U(pn(x))*U(pn(x))*(G(pn(x),pn(x))*cos(th(pn(x))-th(pn(x))) B(pn(x),pn(x))*sin(th(pn(x))-th(pn(x))))-2*U(pn(x))^2*G(pn(x),pn(x));        else            N(x,y)=-U(pn(x))*U(qn(y))*(G(pn(x),qn(y))*cos(th(pn(x))-th(qn(y))) B(pn(x),qn(y))*sin(th(pn(x))-th(qn(y))));        end    endend% N(pn==DCn(1),qn==DCn(1))=N(pn==DCn(1),qn==DCn(1))-3*sqrt(2)/pi*cos(Control_ang(1))*Id;% N(pn==DCn(2),qn==DCn(2))=N(pn==DCn(2),qn==DCn(2))-3*sqrt(2)/pi*cos(Control_ang(2))*Id;Naa=N(Pn-1,Qn-Qn(1) 1);Nat=N(Pn-1,DCn-Qn(1) 1);Nta=N(DCn-1,Qn-Qn(1) 1);Ntt=N(DCn-1,DCn-Qn(1) 1);N=[Naa,Nat;Nta,Ntt];%矩阵M的形成for x=1:(q Ndc)    for y=1:(p Ndc)        if qn(x)==pn(y)            for m=1:1:n                M(x,y)=M(x,y)-U(qn(x))*U(m)*(G(qn(x),m)*cos(th(qn(x))-th(m)) B(qn(x),m)*sin(th(qn(x))-th(m)));            end                M(x,y)=M(x,y) U(qn(x))*U(qn(x))*(G(qn(x),qn(x))*cos(th(qn(x))-th(qn(x))) B(qn(x),qn(x))*sin(th(qn(x))-th(qn(x))));        else            M(x,y)=U(qn(x))*U(pn(y))*(G(qn(x),pn(y))*cos(th(qn(x))-th(pn(y))) B(qn(x),pn(y))*sin(th(qn(x))-th(pn(y))));        end    endendMaa=M(Qn-Qn(1) 1,Pn-1);Mat=M(Qn-Qn(1) 1,DCn-1);Mta=M(DCn-Qn(1) 1,Pn-1);Mtt=M(DCn-Qn(1) 1,DCn-1);M=[Maa,Mat;Mta,Mtt];%矩阵L的形成for x=1:(q Ndc)    for y=1:(q Ndc)        if qn(x)==qn(y)            for m=1:1:n                L(x,y)=L(x,y)-U(qn(x))*U(m)*(G(qn(x),m)*sin(th(qn(x))-th(m))-B(qn(x),m)*cos(th(qn(x))-th(m)));            end                L(x,y)=L(x,y) U(qn(x))*U(qn(x))*(G(qn(x),qn(x))*sin(th(qn(x))-th(qn(x)))-B(qn(x),qn(x))*cos(th(qn(x))-th(qn(x)))) 2*U(qn(x))^2*B(qn(x),qn(x));        else            L(x,y)=-U(qn(x))*U(qn(y))*(G(qn(x),qn(y))*sin(th(qn(x))-th(qn(y)))-B(qn(x),qn(y))*cos(th(qn(x))-th(qn(y))));        end    endend% L(qn==DCn(1),qn==DCn(1))=N(qn==DCn(1),qn==DCn(1))-3*sqrt(2)/pi*cos(Control_ang(1))*Id*tan(fai(1));% L(qn==DCn(2),qn==DCn(2))=N(qn==DCn(2),qn==DCn(2))-3*sqrt(2)/pi*cos(Control_ang(2))*Id*tan(fai(2));Laa=L(Qn-Qn(1) 1,Qn-Qn(1) 1);Lat=L(Qn-Qn(1) 1,DCn-Qn(1) 1);Lta=L(DCn-Qn(1) 1,Qn-Qn(1) 1);Ltt=L(DCn-Qn(1) 1,DCn-Qn(1) 1);L=[Laa,Lat;Lta,Ltt];
%求Δd1----Δd3 for m=1:Ndc     Dd1(m)=Vd(m)-3*sqrt(2)/pi*U(DCn(m))*W(m) 3/pi*Xc(m)*Id;     Dd2(m)=Pd(m)-Vd(m)*Id;     Dd3(m)=Sd(m)-3*sqrt(2)/pi*U(DCn(m))*Id;     Dd4(m)=Sd(m).^2-Pd(m).^2-Qd(m).^2;          endDd5=(1/Rdc*Vd(1)-1/Rdc*Vd(2))-Id; %% JPX JQX JRV JRX
switch Controlmode    case 1JRX_1=eye(2);

JRX_2=[-3*sqrt(2)/pi*Kt(1)*U(DCn(1)),0;0,-3*sqrt(2)/pi*Kt(2)*U(DCn(2))];JRX_3=zeros(2);

JRX=[[1;0],JRX_2,zeros(2,6);[-Id;0],JRX_3,JRX_1,JRX_3,JRX_3;[0;0],JRX_3,JRX_3,JRX_1,JRX_3;...    [0;0],JRX_3,[-2*Pd(1),0;0,-2*Pd(2)],[2*Sd(1),0;0,2*Sd(2)],[-2*Qd(1),0;0,-2*Qd(2)];1/Rdc,zeros(1,7),0];%%为使雅可比矩阵不奇异，将JRX对角线元素为0设为1

JPX=zeros(2,9);JPX(1,4)=-1;JPX(2,5)=1;
JQX=zeros(2,9);JQX(1,8)=-1;JQX(2,9)=-1;
JRV_1=[-3*sqrt(2)/pi*W(1)*U(DCn(1)),0;0,-3*sqrt(2)/pi*W(2)*U(DCn(2))];JRV_2=[-3*sqrt(2)/pi*Id*U(DCn(1)),0;0,-3*sqrt(2)/pi*Id*U(DCn(2))];JRV=[JRV_1;zeros(2);JRV_2;zeros(3,2)];
J=[H,N,[zeros(6,9);JPX];M,L,[zeros(4,9);JQX];zeros(9,8),zeros(9,4),JRV,JRX];DD=-inv(J)*[DPa;DPt;DQa;DQt;Dd1;Dd2;Dd3;Dd4;Dd5]; %将增量加到U th上
   dth(Pn)=DD(1:6);    dth(DCn)=DD(7:8);    dU(Qn)=DD(9:12).*U(Qn);    dU(DCn)=DD(13:14).*U(DCn);    dVd=DD(15);    dW=DD(16:17);    dPd=DD(18:19);    dSd=DD(20:21);    dQd=DD(22:23);            th=th dth;    U=U dU;    Vd(1)=Vd(1) dVd;    W=W dW;    Pd=Pd dPd;    Sd=Sd dSd;    Qd=Qd dQd;   case 2JRX_1=eye(2);

JRX_2=[-3*sqrt(2)/pi*Kt(1)*U(DCn(1));0];JRX_3=zeros(2);

JRX=[JRX_1,JRX_2,zeros(2,6);[-Id,0;0,-Id],[0;0],JRX_1,JRX_3,JRX_3;[0;0],JRX_3,JRX_3,JRX_1,JRX_3;...    [0;0],JRX_3,[-2*Pd(1),0;0,-2*Pd(2)],[2*Sd(1),0;0,2*Sd(2)],[-2*Qd(1),0;0,-2*Qd(2)];1/Rdc,-1/Rdc,zeros(1,7)];%%为使雅可比矩阵不奇异，将JRX对角线元素为0设为1

JPX=zeros(2,9);JPX(1,4)=-1;JPX(2,5)=1;
JQX=zeros(2,9);JQX(1,8)=-1;JQX(2,9)=-1;
JRV_1=[-3*sqrt(2)/pi*W(1)*U(DCn(1)),0;0,-3*sqrt(2)/pi*W(2)*U(DCn(2))];JRV_2=[-3*sqrt(2)/pi*Id*U(DCn(1)),0;0,-3*sqrt(2)/pi*Id*U(DCn(2))];JRV=[JRV_1;zeros(2);JRV_2;zeros(3,2)];
J=[H,N,[zeros(6,9);JPX];M,L,[zeros(4,9);JQX];zeros(9,8),zeros(9,4),JRV,JRX];DD=-inv(J)*[DPa;DPt;DQa;DQt;Dd1;Dd2;Dd3;Dd4;Dd5]; %将增量加到U th上
   dth(Pn)=DD(1:6);    dth(DCn)=DD(7:8);    dU(Qn)=DD(9:12).*U(Qn);    dU(DCn)=DD(13:14).*U(DCn);    dVd=DD(15:16);    dW=DD(17);    dPd=DD(18:19);    dSd=DD(20:21);    dQd=DD(22:23);            th=th dth;    U=U dU;    Vd=Vd dVd;    W(1)=W(1) dW;    Pd=Pd dPd;    Sd=Sd dSd;    Qd=Qd dQd;        case 3JRX_1=eye(2);

JRX_2=[-3*sqrt(2)/pi*Kt(1)*U(DCn(1)),0;0,-3*sqrt(2)/pi*Kt(2)*U(DCn(2))];JRX_3=zeros(2);

JRX=[[1;0],JRX_2,zeros(2,1),JRX_3,JRX_3,[3/pi*Xc(1);3/pi*Xc(2)];[-Id;0],JRX_3,[0;1],JRX_3,JRX_3,[-Vd(1);-Vd(2)];[0;0],JRX_3,[0;0],JRX_1,JRX_3,[-3*sqrt(2)/pi*Kt(1)*U(DCn(1));-3*sqrt(2)/pi*Kt(2)*U(DCn(2))];...    [0;0],JRX_3,[0;-2*Pd(2)],[2*Sd(1),0;0,2*Sd(2)],[-2*Qd(1),0;0,-2*Qd(2)],[0;0];1/Rdc,zeros(1,7),-1];%%为使雅可比矩阵不奇异，将JRX对角线元素为0设为1

JPX=zeros(2,9);JPX(2,4)=1;
JQX=zeros(2,9);JQX(1,7)=-1;JQX(2,8)=-1;
JRV_1=[-3*sqrt(2)/pi*W(1)*U(DCn(1)),0;0,-3*sqrt(2)/pi*W(2)*U(DCn(2))];JRV_2=[-3*sqrt(2)/pi*Id*U(DCn(1)),0;0,-3*sqrt(2)/pi*Id*U(DCn(2))];JRV=[JRV_1;zeros(2);JRV_2;zeros(3,2)];


%% 雅可比矩阵J=[H,N,[zeros(6,9);JPX];M,L,[zeros(4,9);JQX];zeros(9,8),zeros(9,4),JRV,JRX];DD=-inv(J)*[DPa;DPt;DQa;DQt;Dd1;Dd2;Dd3;Dd4;Dd5]; %将增量加到U th上
   dth(Pn)=DD(1:6);    dth(DCn)=DD(7:8);    dU(Qn)=DD(9:12).*U(Qn);    dU(DCn)=DD(13:14).*U(DCn);    dVd=DD(15);    dW=DD(16:17);    dPd=DD(18);    dSd=DD(19:20);    dQd=DD(21:22);    dId=DD(23);        th=th dth;    U=U dU;    Vd(1)=Vd(1) dVd;    W=W dW;    Pd(2)=Pd(2) dPd;    Sd=Sd dSd;    Qd=Qd dQd;    Id=Id dId;     case 4JRX_1=eye(2);

JRX_2=[-3*sqrt(2)/pi*Kt(1)*U(DCn(1));0];JRX_3=zeros(2);

JRX=[JRX_1,JRX_2,zeros(2,1),JRX_3,JRX_3,[3/pi*Xc(1);3/pi*Xc(2)];[-Id,0;0,-Id],[0;0],[0;1],JRX_3,JRX_3,[-Vd(1);-Vd(2)];[0;0],JRX_3,[0;0],JRX_1,JRX_3,[-3*sqrt(2)/pi*Kt(1)*U(DCn(1));-3*sqrt(2)/pi*Kt(2)*U(DCn(2))];...    [0;0],JRX_3,[0;-2*Pd(2)],[2*Sd(1),0;0,2*Sd(2)],[-2*Qd(1),0;0,-2*Qd(2)],[0;0];1/Rdc,-1/Rdc,zeros(1,6),-1];%%为使雅可比矩阵不奇异，将JRX对角线元素为0设为1

JPX=zeros(2,9);JPX(2,4)=1;
JQX=zeros(2,9);JQX(1,7)=-1;JQX(2,8)=-1;
JRV_1=[-3*sqrt(2)/pi*W(1)*U(DCn(1)),0;0,-3*sqrt(2)/pi*W(2)*U(DCn(2))];JRV_2=[-3*sqrt(2)/pi*Id*U(DCn(1)),0;0,-3*sqrt(2)/pi*Id*U(DCn(2))];JRV=[JRV_1;zeros(2);JRV_2;zeros(3,2)];


%% 雅可比矩阵J=[H,N,[zeros(6,9);JPX];M,L,[zeros(4,9);JQX];zeros(9,8),zeros(9,4),JRV,JRX];DD=-inv(J)*[DPa;DPt;DQa;DQt;Dd1;Dd2;Dd3;Dd4;Dd5]; %将增量加到U th上
   dth(Pn)=DD(1:6);    dth(DCn)=DD(7:8);    dU(Qn)=DD(9:12).*U(Qn);    dU(DCn)=DD(13:14).*U(DCn);    dVd=DD(15:16);    dW=DD(17);    dPd=DD(18);    dSd=DD(19:20);    dQd=DD(21:22);    dId=DD(23);        th=th dth;    U=U dU;    Vd=Vd dVd;    W(1)=W(1) dW;    Pd(2)=Pd(2) dPd;    Sd=Sd dSd;    Qd=Qd dQd;    Id=Id dId;            case 5JRX_1=eye(2);

JRX_2=[0;-3*sqrt(2)/pi*Kt(2)*U(DCn(2))];JRX_3=zeros(2);

JRX=[JRX_1,JRX_2,zeros(2,1),JRX_3,JRX_3,[3/pi*Xc(1);3/pi*Xc(2)];[-Id,0;0,-Id],[0;0],[1,0;0,1],JRX_3,JRX_3;[0;0],JRX_3,JRX_3,JRX_1,JRX_3;...    [0;0],JRX_3,[-2*Pd(1),0;0,-2*Pd(2)],[2*Sd(1),0;0,2*Sd(2)],[-2*Qd(1),0;0,-2*Qd(2)];1/Rdc,-1/Rdc,zeros(1,6),0];%%为使雅可比矩阵不奇异，将JRX对角线元素为0设为1

JPX=zeros(2,9);JPX(2,4)=1;
JQX=zeros(2,9);JQX(1,7)=-1;JQX(2,8)=-1;
JRV_1=[-3*sqrt(2)/pi*W(1)*U(DCn(1)),0;0,-3*sqrt(2)/pi*W(2)*U(DCn(2))];JRV_2=[-3*sqrt(2)/pi*Id*U(DCn(1)),0;0,-3*sqrt(2)/pi*Id*U(DCn(2))];JRV=[JRV_1;zeros(2);JRV_2;zeros(3,2)];


%% 雅可比矩阵J=[H,N,[zeros(6,9);JPX];M,L,[zeros(4,9);JQX];zeros(9,8),zeros(9,4),JRV,JRX];DD=-inv(J)*[DPa;DPt;DQa;DQt;Dd1;Dd2;Dd3;Dd4;Dd5]; %将增量加到U th上
   dth(Pn)=DD(1:6);    dth(DCn)=DD(7:8);    dU(Qn)=DD(9:12).*U(Qn);    dU(DCn)=DD(13:14).*U(DCn);    dVd=DD(15:16);    dW=DD(17);    dPd=DD(18:19);    dSd=DD(20:21);    dQd=DD(22:23);            th=th dth;    U=U dU;    Vd=Vd dVd;    W(2)=W(2) dW;    Pd=Pd dPd;    Sd=Sd dSd;    Qd=Qd dQd;  end           AP=abs(DD);%取模值



%判断是否收敛
if(max(AP)<1e-5)>    break;end
kk=kk 1;     end if mm==200     disp('注意：结果不收敛！');end    disp(' 迭代次数：')    kk    disp('电压幅值：')    U      toc                       


2）子函数


Bus.con = [ ...   1  16.5  1  0  4  1;  2  18  1  0  5  1;  3  13.8  1  0  3  1;  4  230  1  0  2  1;  5  230  1  0  2  1;  6  230  1  0  2  1;  7  230  1  0  2  1;  8  230  1  0  2  1;  9  230  1  0  2  1; ];
Line.con = [ ...   9  8  100  230  60  0  0  0.0119  0.1008  0.209  0  0  0  0  0  1;  9  6  100  230  60  0  0  0.039  0.17  0.358  0  0  0  0  0  1;  7  5  100  230  60  0  0  0.032  0.161  0.306  0  0  0  0  0  1;  5  4  100  230  60  0  0  0.01  0.085  0.176  0  0  0  0  0  1;  6  4  100  230  60  0  0  0.017  0.092  0.158  0  0  0  0  0  1;  2  7  100  18    60  0  0.07826087  0  0.0625  0  0  0  0  0  0  1;  3  9  100  13.8  60  0  0.06  0  0.0586  0  0  0  0  0  0  1;  1  4  100  16.5  60  0  0.07173913  0  0.0576  0  0  0  0  0  0  1; ];
SW.con = [ ...   1  100  16.5  1.04  0  99  -99  1.1  0.9  0.8  1  1  1; ];
PV.con = [ ...   2  100  18  1.63  1.025  99  -99  1.1  0.9  1  1;  3  100  13.8  0.85  1.025  99  -99  1.1  0.9  1  1; ];
PQ.con = [ ...   6  100  230  0.9  0.3  1.2  0.8  0  1;  8  100  230  1  0.35  1.2  0.8  0  1;  5  100  230  1.25  0.5  1.2  0.8  0  1; ];
Hvdc.con = [ ...   7  8  100  230  230  45  300  0.3  0.1345  0.1257  1  1  20  25  0.0625  0.2  120  5  17  17  1.2  0.1  1.1  0  1  0.5  0.5  1  1; ];
Bus.names = {...   'Bus 1'; 'Bus 2'; 'Bus 3'; 'Bus 4'; 'Bus 5';   'Bus 6'; 'Bus 7'; 'Bus 8'; 'Bus 9'};


登录后查看更多
交直流系统等值模型
0
评论 0
收藏
侵权举报
声明：本文内容及配图由入驻作者撰写或者入驻合作网站授权转载。文章观点仅代表作者本人，不代表凡亿课堂立场。文章及其配图仅供工程师学习之用，如有内容图片侵权或者其他问题，请联系本站作侵删。
上一篇 物联网的安全面临着哪些危险？

热门评论0

相关文章

基于改进粒子群算法的微电网优化调度

基于改进粒子群算法的微电网优化调度
基于改进粒子群算法的微电网优化调度
基于改进粒子群算法的微电网优化调度

参考文献：摘 要：当今全球普遍面临着能源危机和环境污染的加重，污染严重的化石能源将逐渐被无污染的清洁能源代替，随着经济的发展和社会的进步，以传统的发电方式已经不能满足当今电力用户对电能高可靠性、稳定性的要求需求。以风、光、储、微型燃气轮机、燃料电池等分布式电源构成的微电网能有效的解决了能源危机和环境

电磁MATLAB
2024-08-15 17:30
浏览数325
基于matpower的电力系统潮流计算

基于matpower的电力系统潮流计算

MATPOWER是一个用MATLAB的M文件编写，用来解决电力潮流和优化潮流的问题的软件包。它是由美国康奈尔大学电力系统工程研究中心（PSERC of Cornell University）的RAY D. Zimmerman、Carlos E. Murillo-Sánchez和甘德强在Robert

电磁MATLAB
2024-10-18 14:56
浏览数161

基于两相静止坐标系上的状态方程的感应电动机起动动态计算matlab程序

基于两相静止坐标系上的状态方程的感应电动机起动动态计算matlab程序
基于两相静止坐标系上的状态方程的感应电动机起动动态计算matlab程序
基于两相静止坐标系上的状态方程的感应电动机起动动态计算matlab程序

1 静止两相正交坐标系中的动态数学模型1.1 定子绕组和转子绕组的3/2变换对静止的定子三相绕组和旋转的转子三相绕组进行相同的3/2变换，变换后的定子两相正交坐标系静止，而转子两相正交坐标系以角速度逆时针旋转。1.2 电压方程1.3 磁链方程1.4 转矩方程2 静止两相正交坐标系中的状态方程三相感应

电磁MATLAB
2024-08-15 17:23
浏览数289

Matlab的文件管理命令大全

Matlab的文件管理命令大全
Matlab的文件管理命令大全
Matlab的文件管理命令大全

Matlab是一款专业的科学计算和数据分析软件，广泛应用在工程、科学、金融等多种领域，它可以帮助用户进行高效数据可视化、数值计算、模拟建模等操作，为用户提供了一个良好完整的科学计算环节，今天我将向大家减少这款工具的文件命令，希望对小伙伴们有

电子芯期天
2023-03-24 11:39
浏览数1335

直流潮流计算matlab程序

直流潮流计算matlab程序
直流潮流计算matlab程序
直流潮流计算matlab程序

一、直流潮流计算原理直流潮流的特点是用电力系统的交流潮流（有功功率和无功功率）等值的直流电流来代替。甚至只用直流电路的解析法来分析电力系统的有功潮流，而不考虑无功分布对有功的影响。这样一来计算速度加快，但计算的准确度有所降低，本方法适用于对潮流计算准确度要求不高的计算场景。下面先对直流潮流法的原理进

电磁MATLAB
2024-08-15 17:22
浏览数512

电磁MATLAB

助力科研

硬件设计 Altium 单片机 全部

精彩课程

开班信息
设计指南
更多>>
PCB设计指南

2023-08-21 16:30:02

Altium Designer

2023-08-26 13:52:07

Cadence Allegro

2023-08-26 13:52:41

Mentor PADS

2023-08-26 13:52:59

PCB封装指南

2023-08-27 22:52:59

嵌入式单片机指南

2023-08-30 15:55:42

热门面试题
更多>>
中间继电器的工作原理及作用详解

2023-09-09 14:36:37

请解释PCB（印刷电路板）是什么，并简述其工作原理

2023-09-09 15:31:00

不同等级的高速PCB工程师，薪资差距有多大？

2023-09-09 15:41:33

电话面试常见的8个错误

2023-09-09 16:37:29

技术岗位面试如何与HR谈薪

2023-09-09 16:38:54

热门推荐
更多>>
简谈8051微控制器及汇编指令 1 时钟发生器在计算机到底什么作用？ 2 双T滤波器和其他滤波器有什么不同？ 3 π型滤波器和其他滤波器有什么不同？ 4 穿心电容如何滤波？三分钟说清！ 5 LC滤波器和其他滤波器有什么不同？ 6 C语言：如何分辨指针变量*p和取值运算*p？ 7 电池模拟器的常见参数有哪些？ 8 Allegro半自动出换pin表SKILL小程序 9 allegro如何自动绕线 10 【Allegro封装库设计50问解析】第15问 Allegro软件中的常规表贴焊盘应该如何创建呢？
  • 社交平台
关于凡亿 发展历程 联系我们 企业新闻 网站广告 APP下载

关注我们

下载APP

版权所有© 2022 湖南凡亿智邦电子科技有限公司湘ICP备18024493号-1

网站地图 xml地图 营业执照 资质证书