潮流计算课程设计

《潮流计算课程设计》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、%本程序的功能是用牛顿——拉夫逊法进行潮流计算%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;

2、%input('请输入平衡母线节点号：isb=');pr=1e-5;%input('请输入误差精度：pr=');%---------------------------------------------------n=10;%input('请输入节点数：n=');nl=10;%input('请输入支路数：nl=');B1=[123.4+12.8i1.4e-4i100;145.1+19.2i2.1e-4i100;234.25+16i1.75e-4i100;454.25+16i7e-4i100;135.1+19.2i2.1e-4i100;

3、645.95+22.4i9.8e-4i100;271.78+53.89i038.5/23101;381.49+48.02i011/23101;491.49+48.02i011/23101;5102.46+70.17i038.5/23101]%input('请输入由支路参数形成的矩阵：B1=');B2=[00225.501;0022002;0022002;0022002;0022002;120023103;061.11+37.87i3502;047.53+29.46i1002;054.32+33.66i1002;040.74+25.25i

4、3502]%input('请输入各节点参数形成的矩阵：B2=');%-------------------------------------------------------------%n=4;%input('请输入节点数：n=');%nl=4;%input('请输入支路数：nl=');%B1=[124+16i0100;%134+16i0100;%232+8i0100;%241.49+48.02i011/11001]%input('请输入由支路参数形成的矩阵：B1=');%B2=[0011501;%0011002;%020+4i1

5、1002;%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%--------

6、---如果是变压器支路--------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%------------否则为线路支

7、路--------------------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)./2.0;%对角元j侧+线路电纳的一半Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)./2.0;%对角元i侧+线路电纳的一半endenddisp('导纳矩阵Y=');disp(Y);%-----------给定各节点初始电压及给定各节点注入功率-----------------

8、---------G=real(Y);B=imag(Y);%分解出导纳阵的实部和虚部fori=1:n%给定各节点初始电压的实部和虚部e(i)=real(B2(i,3));f(i)=imag(B2(i,3))