电力系统分析潮流计算实验报告.docVIP

电力系统分析潮流计算实验报告 姓名：XXXXXX 学号：XXXXXXXXXX 班级：XXXXXXXX 一、实验目的 掌握潮流计算计算机算法的方法，熟悉MATLAB的程序调试方法。 二、实验准备 根据课程内容，熟悉MATLAB软件的使用方法，自行学习MATLAB程序的基础语法，并根据所学知识编写潮流计算牛顿拉夫逊法(或PQ分解法) 的计算程序，用相应的算例在MATLAB上进行计算、调试和验证。 三、实验要求 每人一组，在实验课时内，调试和修改运行程序，用算例计算输出潮流结果。 四、程序流程 五、实验程序 %本程序的功能是用牛拉法进行潮流计算 %原理介绍详见鞠平著《电气工程》 %默认数据为鞠平著《电气工程》例8.4所示数据 %B1是支路参数矩阵 %第一列和第二列是节点编号。节点编号由小到大编写 %对于含有变压器的支路，第一列为低压侧节点编号，第二列为高压侧节点编号 %第三列为支路的串列阻抗参数，含变压器支路此值为变压器短路电抗 %第四列为支路的对地导纳参数，含变压器支路此值不代入计算 %第五烈为含变压器支路的变压器的变比，变压器非标准电压比 %第六列为变压器是否是否含有变压器的参数，其中“1”为含有变压器，“0”为不含有变压器 %B2为节点参数矩阵 %第一列为节点注入发电功率参数 %第二列为节点负荷功率参数 %第三列为节点电压参数 %第四列 %第五列 %第六列为节点类型参数，“1”为平衡节点，“2”为PQ节点，“3”为PV节点参数 %X为节点号和对地参数矩阵 %第一列为节点编号 %第二列为节点对地参数 %默认算例 % n=4; % n1=4; % isb=4; % pr=0.00001; % B1=[1 2 0.1667i 0 0.8864 1;1 3 0.1302+0.2479i 0.0258i 1 0;1 4 0.1736+0.3306i 0.0344i 1 0;3 4 0.2603+0.4959i 0.0518i 1 0]; % B2=[0 0 1 0 0 2;0 -0.5-0.3i 1 0 0 2;0.2 0 1.05 0 0 3;0 -0.15-0.1i 1.05 0 0 1]; % X=[1 0;2 0.05i;3 0;4 0]; clear; clc; num=input(是否采用默认数据？(1-默认数据;2-手动输入)); if num==1 n=4; n1=4; isb=4; pr=0.00001; B1=[1 2 0.1667i 0 0.8864 1;1 3 0.1302+0.2479i 0.0258i 1 0;1 4 0.1736+0.3306i 0.0344i 1 0;3 4 0.2603+0.4959i 0.0518i 1 0]; B2=[0 0 1 0 0 2;0 -0.5-0.3i 1 0 0 2;0.2 0 1.05 0 0 3;0 -0.15-0.1i 1.05 0 0 1]; X=[1 0;2 0.05i;3 0;4 0]; else n=input(请输入节点数:n=); n1=input(请输入支路数:n1=); isb=input(请输入平衡节点号:isb=); pr=input(请输入误差精度:pr=); B1=input(请输入支路参数:B1=); B2=input(请输入节点参数:B2=); X=input(节点号和对地参数:X=); end Times=1; %迭代次数 %创建节点导纳矩阵 Y=zeros(n); for i=1:n1 if B1(i,6)==0 %不含变压器的支路 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(p,p)=Y(p,p)+1/B1(i,3)+0.5*B1(i,4); Y(q,q)=Y(q,q)+1/B1(i,3)+0.5*B1(i,4); else %含有变压器的支路 p=B1(i,1); q=B1(i,2); Y(p,q)=Y(p,q)-B1(i,5)/B1(i,3); Y(q,p)=Y(p,q); Y(p,p)=Y(p,p)+B1(i,5)/B1(i,3)+(1-B1(i,5))/B1(i,3); Y(q,q)=Y(q,q)+B1(i,5)/B1(i,3)+(B1(i,5)*(B1(i,5)-1))