资源描述：
《电力系统分析潮流实验报告.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、南昌大学实验报告学生姓名：学号：专业班级：实验类型：□验证□综合■设计□创新实验日期：实验成绩：电力系统潮流计算实验一、实验目的：本实验通过对电力系统潮流计算的计算机程序的编制与调试，获得对复杂电力系统进行潮流计算的计算机程序，使系统潮流计算能够由计算机自行完成，即根据已知的电力网的数学模型（节点导纳矩阵）及各节点参数，由计算程序运行完成该电力系统的潮流计算。通过实验教学加深学生对复杂电力系统潮流计算计算方法的理解，学会运用电力系统的数学模型，掌握潮流计算的过程及其特点，熟悉各种常用应用软件，熟悉硬件设备的使用方法，加强编制调试计算机程序的能力，提高工程计算的能力，学习如何将理论知
2、识和实际工程问题结合起来。二、实验内容：编制调试电力系统潮流计算的计算机程序。程序要求根据已知的电力网的数学模型（节点导纳矩阵）及各节点参数，完成该电力系统的潮流计算，要求计算出节点电压、功率等参数。1、在各种潮流计算的算法中选择一种，按照计算方法编制程序。2、将事先编制好的电力系统潮流计算的计算程序原代码由自备移动存储设备导入计算机。3、在相应的编程环境下对程序进行组织调试。4、应用计算例题验证程序的计算效果。三、实验程序：function[e,f,p,q]=flow_out(g,b,kind,e,f)%计算潮流后efpq的终值s=flow(g,b,kind,e,f);k=0;w
3、hilemax(abs(s))>10^-5J=J_out(g,b,kind,e,f);J_ni=inv(J);dv=J_ni*s;l=length(dv)/2;fori=1:le(i)=e(i)-dv(2*i-1);f(i)=f(i)-dv(2*i);ends=flow(g,b,kind,e,f);endl=length(e);fori=1:ls1=0;s2=0;forj=1:ls1=s1+g(i,j)*e(j)-b(i,j)*f(j);s2=s2+g(i,j)*f(j)+b(i,j)*e(j);endp(i)=e(i)*s1+f(i)*s2;q(i)=f(i)*s1-e(i)*s
4、2;endfunctions=flow(g,b,kind,e,f)%计算当前ef与规定的pqv的差值l=length(e);s=zeros(2*l-2,1);fori=1:(l-1)s1=0;s2=0;forj=1:ls1=s1+g(i,j)*e(j)-b(i,j)*f(j);s2=s2+g(i,j)*f(j)+b(i,j)*e(j);ends(2*i-1)=kind(2,i)-e(i)*s1-f(i)*s2;ifkind(1,i)==1s(2*i)=kind(3,i)-f(i)*s1+e(i)*s2;elses(2*i)=kind(3,i)^2-f(i)^2-e(i)^2;end
5、endfunctionJ=J_out(g,b,kind,e,f)%计算节点的雅克比矩阵l=length(e);J=zeros(2*l-2,2*l-2);fori=1:(l-1);ifkind(1,i)==1s=PQ_out(g,b,e,f,i);forj=1:(2*l-2)J(2*i-1,j)=s(1,j);J(2*i,j)=s(2,j);endelses=PV_out(g,b,e,f,i);forj=1:(2*l-2)J(2*i-1,j)=s(1,j);J(2*i,j)=s(2,j);endendendfunctionpq=PQ_out(g,b,e,f,i)%计算pq节点的雅克比
6、矩阵l=length(e);pq=zeros(2,2*l-2);forj=1:(l-1)ifj==is=0;fork=1:ls=s-(g(i,k)*e(k)-b(i,k)*f(k));endpq(1,2*i-1)=s-g(i,i)*e(i)-b(i,i)*f(i);s=0;fork=1:ls=s-(g(i,k)*f(k)+b(i,k)*e(k));endpq(1,2*i)=s+b(i,i)*e(i)-g(i,i)*f(i);s=0;fork=1:ls=s+(g(i,k)*f(k)+b(i,k)*e(k));endpq(2,2*i-1)=s+b(i,i)*e(i)-g(i,i)*f(
7、i);s=0;fork=1:ls=s-(g(i,k)*e(k)-b(i,k)*f(k));endpq(2,2*i)=s+g(i,i)*e(i)+b(i,i)*f(i);elsepq(1,2*j-1)=-(g(i,j)*e(i)+b(i,j)*f(i));pq(1,2*j)=b(i,j)*e(i)-g(i,j)*f(i);pq(2,2*j)=-pq(1,2*j-1);pq(2,2*j-1)=pq(1,2*j);endendfunctionpv=PV_out(g,b,e,
