现代电力系统分析.doc

(完整word)现代电力系统分析(完整word)现代电力系统分析(完整word)现代电力系统分析选用牛顿—拉佛森方法，利用matlab软件计算基于PQ节点情况下的潮流计算。一．所用公式其中二、程序流程图开始形成节点导纳矩阵输入原始数据设节点电压，i=1,2…,n,is置迭代次数置节点号i=1计算雅克比矩阵元素计算节点的，，节点的，求解修正方程式，得，雅克比矩阵是否已全部形成？计算平衡节点及PV节点功率求，迭代次数k=k+1i=i+1？潮流计算完成计算各节点电压的新值：牛顿—拉佛森例题中的简单模型系统系统单线图图1简单模型系统二、系统参数节点矩阵％(bus#）(volt)（ang）（p）(q）（bustype)bus=［11.000。00-1。60—0。801；21。000.00—2。00—1。001;31.000。00—3.70—1。301；41.050。005。000。002;51。050。000.000。003];线路矩阵％b#1b#2(R）（X)（G）（B)(K)line=[120.040.2500.250；130.100。3500.000;230。080.3000。250；530.000.0300.001。05;420。000.01500.001.05];三、计算结果：牛顿－拉夫逊法潮流计算结果节点计算结果：n节点节点电压节点相角（角度）节点注入功率10。862150-4.778511—1。600000+j—0.80000021。07791617.853530-2。000000+j-1。00000031.036411—4。281930-3.700000+j-1。30000041.05000021.8433195。000000+j1.81308451.0500000.0000002。579427+j2.299402n线路计算结果：n节点I节点J线路功率S（I,J）线路功率S（J，I）线路损耗dS（I，J）12-1。466181+j-0.4090761。584546+j0.6725560.118365+j0。26348013-0.133819+j—0.3909240.156788+j0。4713150。022969+j0。080391231。415454+j—0。244333—1。277360+j0。2031700.138093+j-0。041163532.579427+j2。299402-2.579427+j-1。9744850。000000+j0。324917425.000000+j1。813084-5.000000+j-1。4282230.000000+j0.384861导纳矩阵：Y=［1.3787—6。2917i-0.6240+3.9002i-0.7547+2。6415i00—0。6240+3。9002i1。4539-66。9808i-0。8299+3。1120i0+63.4921i0—0。7547+2。6415i-0.8299+3。1120i1。5846-35。7379i00+31。7460i00+63.4921i00—66.6667i0000+31。7460i00-33.3333i］图2。Matlab运行结果结果分析：此程序的运行结果和试验程序给出的结果是一致的。说明程序无误，但在精确度上有微小差异，这主要是和导纳矩阵的精确度以及显示精度有关。心得：本程序分模块进行,先是排序，再是求导纳阵，然后求雅阁比，再进行迭代运算，程序本身很简洁明了，运行的时候只需要在matlab里输入main就行了，然后打开BUS和line所在的。m文件，结果就会自动存在result文件中了,通过编写牛顿拉夫逊法matlab潮流计算程序复习了潮流计算的知识，也实现了计算机算法附录：实验源程序:Main函数:clear［dfile,pathname］=uigetfile（’*.m’，’SelectDataFile')；ifpathname==0error（’youmustselectavaliddatafile’)elselfile=length（dfile）；％stripoff.meval（dfile(1：lfile—2））；end;globalbus；globalline；［nb,mb］=size(bus)；［nl，ml]=size（line）;％计算bus和line矩阵的行数和列数[bus,line,nPQ,nPV,nodenum］=Num（bus,line);％对节点重新排序的子程序Y=y（bus,line)%计算节点导纳矩阵的子程序myf=fopen('Result。m',’w’)；fprintf（myf,’计算结果’）；fclose(myf)；％在当前目录下生成“Result.m”文件，写入节点导纳矩阵formatlongEP