MATLAB牛顿拉夫逊法算潮流分析.doc

MATLAB牛顿拉夫逊法算潮流分析%本程序的功能是用牛顿——拉夫逊法进行潮流计算%B1矩阵:1、支路首端号;2、末端号;3、支路阻抗;4、线路对地电纳(或变压器导纳);%5、支路的变比;6、支路首端处于K侧为1,1侧为0;%7、线路/变压器标识(0/1)变压器参数当支路首端处于K侧标识为1时归算至末端侧,0归算至首端侧%B2矩阵:1、该节点发电机功率;2、该节点负荷功率;%3、PQ节点电压初始值,或平衡节点及PV节点电压的给定值%4、节点所接无功补偿并联电容(感)的电纳%5、节点分类标号:1为平衡节点(应为1号节点);2为PQ节点;3为PV节点;clear;isb=1;%input('请输入平衡母线节点号:isb=');pr=1e-5;%input('请输入误差精度:pr=');%---------------------------------------------------n=10;%input('请输入节点数:n=');nl=10;%input('请输入支路数:nl=');B1=[+-4i100;+-4i100;+-4i100;+16i7e-4i100;+-4i100;+-4i100;+;+;+;+]%input('请输入由支路参数形成的矩阵:B1=');B2=[;0022002;0022002;0022002;0022002;120023103;+;+;+;+]%input('请输入各节点参数形成的矩阵:B2=');%-------------------------------------------------------------%n=4;%input('请输入节点数:n=');%nl=4;%input('请输入支路数:nl=');%B1=[124+16i0100;%134+16i0100;%232+8i0100;%+]%input('请输入由支路参数形成的矩阵:B1=');%B2=[0011501;%0011002;%020+4i11002;%010+6i1002]%input('请输入各节点参数形成的矩阵:B2=');%-------------------------------------------------------------Y=zeros(n);e=zeros(1,n);f=zeros(1,n);V=zeros(1,n);sida=zeros(1,n);S1=zeros(nl);%%%-----------求导纳矩阵------------------------fori=1:nl%从1到n1(总支路数)ifB1(i,7)==1%-----------如果是变压器支路--------ifB1(i,6)==0%左节点(首端)处于1侧p=B1(i,1);q=B1(i,2);else%左节点(首端)处于K侧p=B1(i,2);q=B1(i,1);endY(p,q)=Y(p,q)-1./(B1(i,3)*B1(i,5));%非对角元Y(q,p)=Y(p,q);%非对角元Y(q,q)=Y(q,q‎‎)+1./(B1(i,3)*B1(i,5)^2);%对角元K侧Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4);%对角元1侧+励磁导纳else%------------否则为线路支路--------------------p=B1(i,1);q=B1(i,2);Y(p,q)=Y(p,q)-1./B1(i,3);%非对角元Y(q,p)=Y(p,q);%非对角元Y(q,q)=Y(q,q)+1./B1(i,3)+B1(i,4)./;%对角元j侧+线路电纳的一半Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)./;%对角元i侧+线路电纳的一半endenddisp('导纳矩阵Y=');disp(Y);%-----------给定各节点初始电压及给定各节点注入功率--------------------------G=real(Y);B=imag(Y);%分解出导纳阵的实部和虚部fori=1:n%给定各节点初始电压的实部和虚部e(i)=real(B2(i,3))