电力系统分析潮流计算实验报告
姓名：XXXXXX 学号：XXXXXXXXXX 班级：XXXXXXXX
一、实验目的
掌握潮流计算计算机算法的方法，熟悉MATLAB的程序调试方法。
二、实验准备
根据课程内容，熟悉MATLAB软件的使用方法，自行学习MATLAB程序的基础语法，并根据所学知识编写潮流计算牛顿拉夫逊法(或PQ分解法) 的计算程序，用相应的算例在MATLAB上进行计算、调试和验证。
三、实验要求
每人一组，在实验课时内，调试和修改运行程序，用算例计算输出潮流结果。
四、程序流程
五、实验程序
%本程序的功能是用牛拉法进行潮流计算 %原理介绍详见鞠平著《电气工程》
%默认数据为鞠平著《电气工程》例8.4所示数据
±是支路参数矩阵
%第一列和第二列是节点编号。节点编号由小到大编写
%对于含有变压器的支路，第一列为低压侧节点编号，第二列为高压侧节点编号 %第三列为支路的串列阻抗参数，含变压器支路此值为变压器短路电抗 %第四列为支路的对地导纳参数，含变压器支路此值不代入计算 %第五烈为含变压器支路的变压器的变比，变压器非标准电压比
%第六列为变压器是否是否含有变压器的参数，其中“1”为含有变压器，“0”为不含有变压器
2为节点参数矩阵
%第一列为节点注入发电功率参数 %第二列为节点负荷功率参数 %第三列为节点电压参数 %第四列 %第五列
%第六列为节点类型参数，“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