极坐标法潮流计算.docx
课程设计论文任务书学生姓名年级专业题目名称采用极坐标下的牛顿-拉夫逊计算设计时间课程名称潮流计算课程设计课程编号121202306设计地点综合仿真实验室一、课程设计（论文）目的1 .掌握电力系统极坐标下的牛顿-拉夫逊计算的根本原理；2 .掌握并能熟练运用一门计算机语言（MATLAB语言或FoRTRAN或C语言或C+语言）；3 .采用计算机语言对极坐标下的牛顿-拉夫逊计算进行计算机编程计算。通过课程设计，使学生稳固电力系统潮流计算的根本原理与方法，掌握潮流计算的数值求解方法（节点导纳矩阵，修正方程），开发系统潮流计算的计算程序。让学生掌握用计算机仿真分析电力系统的方法。同时，通过软件开发，也有助于计算机操作能力和软件开发能力的提高。二、技术参数和条件在图所示的简单电力系统中，系统中节点1、2为PQ节点，节点3为PV节点，节点4为平衡节点，已给定Sl=-0.4-/0.3,$2=-。3-/（）.2,6=0.4,V3=1.02,V4=1.05,6>4=O0,网络各元件参数的标幺值如表2所示，给定电压的初始值如表2所示，收敛系数£=0.00001。试求：采用极坐标下的牛顿-拉夫逊计算图1网络的潮流分布。三、任务和要求任务：熟练掌握计算机语言，并采用计算机编程进行以下计算：根据电力系统网络推导电力网络数学模型,写出节点导纳矩阵；掌握潮流计算的数值求解方法（节点导纳矩阵，修正方程），开发系统潮流计算的计算程序。要求：1 .手工计算，手写，采用A4纸，得出计算结果。2.编写程序：它包括程序源代码；程序说明；局部程序的流程图；程序运行结果，电子版。注：1.此表由指导教师填写，经系、教研室审批，指导教师、学生签字后生效：2.此表1式3份，学生、指导教师、教研室各1份。四、参考资料和现有根底条件（包括实验室、主要仪器设备等）1何仰赞等.电力系统分析M.武汉：华中理工大学出版社，2002.32西安交通大学等.电力系统计算M.北京：水利电力出版社，1993.12五、进度安排2010年12月20日：下达课程设计的方案书，任务书，设计题目及分组情况。2010年12月21日-23日：学生完成潮流计算的手工计算。2010年12月24日：讲述课程设计编程的思路、要求。举例：用MATLAB软件编写的局部程序。2010年12月25日-30日：学生编写程序。2011年1月1日-3日：上机调试程序，得出正确结果。2011年1月4日-5日：整理课程设计报告。2011年1月6日：学生辩论。六、教研室审批意见教研室主任（签字）：年月日七I、主管教学主任意见主管主任（签字）：年月日八、备注指导教师（签字）：学生（签字）:设计主题题目一：在以下图所示的简单电力系统中，系统中节点1、2为PQ节点，节点3为PV节点，节点4为平衡节点，已给定4=-0.4-川.3,52=-0.3-0.2,尸3=04，V3=1.02,匕=1.05,=(T,网络各元件参数的标幺值如表1所示，给定电压的初始值如表2所示,收敛系数£=0.00001。试求：采用极坐标下的牛顿-拉夫逊计算图示网络的潮流计算。表1网络各元件参数的标幺值支路电阻电抗输电线路;先变压器变比k120.020.060.01130.010.030.01一230.030.07一240.00.050.9625340.020.05表2各节点电压初值标幺值参数3节点i1234潮力W)=暧)+族/JlL1.00+j0.01.0+j0.01.0+j0.01.05+j0.0计算流程图本次课程设计采用极坐标下的牛顿-拉夫逊计算网络的潮流计算。其牛顿-拉夫逊潮流计算程序框图如下所示。n2=inputC请输入PQ结点数：n2=t);isb=input('请输入平衡结点：isb=');pr=input('请输入精确度：pr=t);K=input，请输入变比矩阵:K=);C=input('请输入支路阻抗矩阵：C=);y=input('请输入支路导纳矩阵：y-);U=inputC请输入结点电压矩阵：U-);S=inputC请输入各结点的功率：S=);Z=zeros(1,n);N=zeros(n1+n2,n2);L=ZeroS(n2,n2);QT1=zeros(l,n1+n2);form=l:nforR=l:nC(m,m)=C(m,m)+y(m,R);ifK(m,R)-=OC(m,m)=C(m,m)+l(C(m,R)/(K(m,R)*(K(m,R)-l);C(R,R)=C(R,R)+l(C(m,R)(l-K(m,R);C(m,R)=C(m,R)K(m,R);C(R,m)=C(m,R);endendendform=l:nforR=l:nifm二RZ(m)=Z(m)+1C(m,R);endendform=l:nforR=l:nifm=RY(m,m)=C(m,m)+Z(m);elseY(m,R)=-lC(m,R);endendenddispC结点导纳矩阵:);disp(Y);dispC迭代中的雅克比矩阵厂);G=real(Y);B=imag(Y);O=angle(U);Ul=abs(U);k=0;PR=I;P=real(三);Q=imag(三);whilePR>prform=l:n2UD(m)=Ul(m);endform=knl+n2forR=l:nPT(R)=U1(m)*U1(R)*(G(m,R)*cos(O(m)-O(R)÷B(m,R)*sin(O(m)-O(R);endPT1(m)=sum(PT);PP(m)=P(m)-PTl(m);PPl(k+l,m)=PP(m);endform=l:n2forR=l:nQT(R)=Ul(m)*Ul(R)*(G(m,R)*sin(O(m)-O(R)-B(m,R)*cos(O(m)-O(R);endQTl(m)=sum(QT);QQ(m)=Q(m)-QTl(m);QQl(k+l,m)=QQ(m);endPR1=max(abs(PP);PR2=max(abs(QQ);PR=max(PRl,PR2);form=knl+n2forR=knl+n2ifm=RH(m,m)=U1(m)2*B(m,m)+QTl(m);elseH(m,R)=-U1(m)*U1(R)*(G(m,R)*sin(O(m)-O(R)-B(m,R)*cos(O(m)-O(R);endendendform=knl+n2forR=I:n2ifm=RN(m,m)=-U1(m)A2*G(m,m)-PT1(m);elseN(m,R)=-Ul(m)*Ul(R)*(G(m,R)*cos(O(m)-O(R)+B(m,R)*sin(O(m)-O(R);endendendform=l:n2forR=knl+n2ifm=RJ(m,m)=U1(m)2*G(m,m)-PTl(m);elseJ(m,R)=Ul(m)*Ul(R)*(G(m,R)*cos(O(m)-O(R)+B(m,R)*sin(O(m)-O(R);endendendform=l:n2forR=I:n2ifm=R1.(m,m)=U1(m)2*B(m,m)-QT1(m);else1.(m,R)=-U1(m)*U1(R)*(G(m,R)*sin(O(m)-O(R)-B(m,R)*cos(O(m)-O(R);endendendJJ=HN;JL;disp(JJ);PQ=PP;QQ;DA=-inv(JJ)*PQ;DA1=DA,;form=knl+n2OO(m)=DAl(m);endform=n:nl+n2+n2UUl(m-nl-n2)=DAl(m);endUD2=diag(UD);UU=UUPUD2;form=knl+n2O(m)=O(m)+OO(m);endform=l:n2Ul(m)=Ul(m)+UU(m);endform=knl+n2o(k+l,m)=180pi*O(m);endform=l:n2u(k+l,m)=Ul(m);endk=k+l;endform=l:nb(m)=Ul(m)*cos(O(m);c(m)=U1(m)*sin(O(m);endU=b+i*c;forR=l:nPHl(R)=U(isb)*conj(Y(isb,R)*conj(U(R);endPH=Sum(PHl);form=l:nforR=l:nifm二RCl(m,R)=lC(m,R);elseCl(m,m)=C(m,m);endendendform=l:nforR=l:nif(C(m,R)=inf)&(m=R)SS(m,R)=Ul(m)A2*conj(Cl(m,m)+U(m)*(conj(U(m)-conj(U(R)*conj(Cl(m,R)；endendenddisp('迭代中的P：)disp(PPl);dispC迭代中的4Q:)disp(QQl);dispC迭代中相角:);d中P(O);dispC迭代中电压的模:);disp(u);disp('平衡结点的功率:');disp(PH);dispC全部线路功率分布)disp(SS);注意：matlab默认输出结果保存4位小数，可在显示屏上输入>>fomatlong此时小数点后面保存14位小数5.2程序结果请输入结点数：n=4请输入PV结点数：nl=l请输入PQ结点数：112=2请输入平衡结点：isb=4请输入精确度：pr=0.00001请输入变比矩阵:K=0000;0000.9625;0000;0000请输入支路阻抗矩阵：C=00.02+0.06i0.01+0.03iinf;0.02÷0.06i00.03+0.07i0.0+0.05i;0.01+0.03i0.03+0.07i00.02+0.05i;inf0.0+0.05i0.02+0.05i01注：inf表示两者未连接请输入支路导纳矩阵：y=00.01i0.0Ii0;0.01i000;0.01i000;0000请输入结点电压矩阵：U=l+0il+0i1.02+0i1.05+0i请输入各结点的功率：S=-0.4-0.3i-0.3-0.2i0.40结点导纳矩阵：15.0000-44.9800i-5.0000+15.0000i-10.0000+30.0000i0-5.0000+15.0000i10.1724-45.587Ii-5.1724+12.0690i0+19.2500i-10.0000+30.0000i-5.1724+12.0690i22.0690-59.3003i-6.8966+17.2414i00+19.2500i-6.8966+17.2414i6.8966-37.241
