极坐标法潮流计算.docx

上传人:王** 文档编号:1113355 上传时间:2024-03-31 格式:DOCX 页数:14 大小:81.25KB
下载 相关 举报
极坐标法潮流计算.docx_第1页
第1页 / 共14页
极坐标法潮流计算.docx_第2页
第2页 / 共14页
极坐标法潮流计算.docx_第3页
第3页 / 共14页
极坐标法潮流计算.docx_第4页
第4页 / 共14页
极坐标法潮流计算.docx_第5页
第5页 / 共14页
极坐标法潮流计算.docx_第6页
第6页 / 共14页
极坐标法潮流计算.docx_第7页
第7页 / 共14页
极坐标法潮流计算.docx_第8页
第8页 / 共14页
极坐标法潮流计算.docx_第9页
第9页 / 共14页
极坐标法潮流计算.docx_第10页
第10页 / 共14页
亲,该文档总共14页,到这儿已超出免费预览范围,如果喜欢就下载吧!
资源描述

《极坐标法潮流计算.docx》由会员分享,可在线阅读,更多相关《极坐标法潮流计算.docx(14页珍藏版)》请在优知文库上搜索。

1、课程设计论文任务书学生姓名年级专业题目名称采用极坐标下的牛顿-拉夫逊计算设计时间课程名称潮流计算课程设计课程编号121202306设计地点综合仿真实验室一、课程设计(论文)目的1 .掌握电力系统极坐标下的牛顿-拉夫逊计算的根本原理;2 .掌握并能熟练运用一门计算机语言(MATLAB语言或FoRTRAN或C语言或C+语言);3 .采用计算机语言对极坐标下的牛顿-拉夫逊计算进行计算机编程计算。通过课程设计,使学生稳固电力系统潮流计算的根本原理与方法,掌握潮流计算的数值求解方法(节点导纳矩阵,修正方程),开发系统潮流计算的计算程序。让学生掌握用计算机仿真分析电力系统的方法。同时,通过软件开发,也有助

2、于计算机操作能力和软件开发能力的提高。二、技术参数和条件在图所示的简单电力系统中,系统中节点1、2为PQ节点,节点3为PV节点,节点4为平衡节点,已给定Sl=-0.4-/0.3,$2=-。3-/().2,6=0.4,V3=1.02,V4=1.05,64=O0,网络各元件参数的标幺值如表2所示,给定电压的初始值如表2所示,收敛系数=0.00001。试求:采用极坐标下的牛顿-拉夫逊计算图1网络的潮流分布。三、任务和要求任务:熟练掌握计算机语言,并采用计算机编程进行以下计算:根据电力系统网络推导电力网络数学模型,写出节点导纳矩阵;掌握潮流计算的数值求解方法(节点导纳矩阵,修正方程),开发系统潮流计算

3、的计算程序。要求:1 .手工计算,手写,采用A4纸,得出计算结果。2.编写程序:它包括程序源代码;程序说明;局部程序的流程图;程序运行结果,电子版。注:1.此表由指导教师填写,经系、教研室审批,指导教师、学生签字后生效:2.此表1式3份,学生、指导教师、教研室各1份。四、参考资料和现有根底条件(包括实验室、主要仪器设备等)1何仰赞等.电力系统分析M.武汉:华中理工大学出版社,2002.32西安交通大学等.电力系统计算M.北京:水利电力出版社,1993.12五、进度安排2010年12月20日:下达课程设计的方案书,任务书,设计题目及分组情况。2010年12月21日-23日:学生完成潮流计算的手工

4、计算。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

5、.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计算流程图本次课程设计采用极坐标下的牛顿-拉夫逊计算网络的潮流计算。其牛顿-拉夫逊潮流计算程序框图如下所示。n

6、2=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

7、,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)

8、;k=0;PR=I;P=real(三);Q=imag(三);whilePRprform=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)=

9、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)*co

10、s(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);endendendJ

11、J=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)=

12、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:)

13、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.

14、020.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

展开阅读全文
相关资源
猜你喜欢
相关搜索

当前位置:首页 > 高等教育 > 微积分

copyright@ 2008-2023 yzwku网站版权所有

经营许可证编号:宁ICP备2022001189号-2

本站为文档C2C交易模式,即用户上传的文档直接被用户下载,本站只是中间服务平台,本站所有文档下载所得的收益归上传人(含作者)所有。装配图网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对上载内容本身不做任何修改或编辑。若文档所含内容侵犯了您的版权或隐私,请立即通知装配图网,我们立即给予删除!