电力系统分析潮流计算最终完整版.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=[120、1667i00、88641;130、1302+0、2479i0、0258i10;140、1736+0、3306i0、0344i10;340、2603+0、4959i0、0518i10];

%B2=[001002;0-0、5-0、3i1002;0、201、05003;0-0、15-0、1i1、05001];

%X=[10;20、05i;30;40];

clear;

clc;

num=input(就是否采用默认数据？(1-默认数据;2-手动输入));

ifnum==1

n=4;

n1=4;

isb=4;

pr=0、00001;

B1=[120、1667i00、88641;130、1302+0、2479i0、0258i10;140、1736+0、3306i0、0344i10;340、2603+0、4959i0、0518i10];

B2=[001002;0-0、5-0、3i1002;0、201、05003;0-0、15-0、1i1、05001];

X=[10;20、05i;30;40];

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);

fori=1:n1

ifB1(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,