电力系统潮流上机计算报告书写要求

时间：2024.3.15

电力系统认识实习报告

系别：电力工程系

班级：电力0903

姓名：陈跃

学号：200904000603

20xx年1月7日

一、程序说明

包括：程序设计思想、程序流程图、程序使用说明。

二、给定题目的手算过程（迭代两次）

包括：原题目、节点导纳矩阵、雅克比矩阵、第一次和第二次迭代

结果。

三、给定题目的程序计算结果

包括：原题目、节点导纳矩阵、雅克比矩阵、程序输入和输出文件

（误差0.0001）。

四、编程特色与创新

包括：程序能够完成的基本功能；程序能够完成的高级功能（如：

是否包括非标准变比变压器支路，是否采用了稀疏矩阵技术，

是否增加了人机对话界面，程序的通用性和实用性如何）。

五、总结

包括：手算结果与程序计算结果的分析比较；本次上机体会，如：

独立编程体会、跟踪调试技能的掌握情况、C语言中结构体、

指针、文件输入输出的掌握情况等。

报告要求：

1．报告中除上面的第“三”项外，其他部分必须手写（最好使用

黑色水笔）。

2．报告统一采用A4打印纸书写（留出页边距: 1.5~2厘米）。不使

用实验报告纸。

3．封面按上述格式书写。

4．装订统一在左侧1厘米，二个钉。

5．上述五部分内容必须齐全，各部分内容可以扩充。

6．报告书写要求字迹清楚，不得潦草，流程图中的框要用尺子画。

7．报告必须与本人提交程序吻合，否则取消成绩。

8．报告不得有雷同，否则全部取消成绩。

第二篇：电力系统潮流计算仿真报告

Beijing Jiaotong University

电力系统潮流计算

仿真报告

姓名：TYP

班级：电气0906

学号：09291183

指导老师：吴俊勇

完成日期：2012.6.24

一、实验内容

电力系统潮流计算是研究电力系统稳态运行情况的一种基本电气计算。它的任务是根据给定的运行条件和网路结构确定整个系统的运行状态，如各母线上的电压（幅值及相角）、网络中的功率分布以及功率损耗等。电力系统潮流计算的结果是电力系统稳定计算和故障分析的基础。

对于简单系统，可以将其分为开式网络和闭式网络手工计算。对于复杂电力系统，根据定解条件，应用牛顿—拉夫逊法进行计算，在手工计算中，由于涉及大量变量、微分方程、矩阵计算，求解很烦琐，而且容易出错，计算不同系统时需要重新计算。故而我们可以借助计算机来进行潮流计算，方便快捷且准确率高。

二、计算机潮流计算方法

我们常用牛顿—拉夫逊法来进行潮流计算。

牛顿—拉夫逊法(简称牛顿法)在数学上是求解非线性代数方程式的有效方法，其要点是把非线性方程式的求解过程变成反复地对相应的线性方程式进行求解的过程，即通常所称的逐次线性化过程。

1、基本原理

从几何意义上，牛顿—拉夫逊法实质上就是切线法，是一种逐步线性化的方法。

2、牛顿—拉夫逊法潮流求解过程

以下讨论的是用直角坐标形式的牛顿—拉夫逊法潮流的求解过程。当采用直角坐标时，潮流问题的待求量为各节点电压的实部和虚部两个分量，由于平衡节点的电压向量是给定的，因此待求量共2(n-1)需要2(n-1)个方程式。事实上，除了平衡节点的功率方程式在迭代过程中没有约束作用以外，其余每个节点都可以列出两个方程式。

求解过程大致可以分为以下步骤：

(1)形成节点导纳矩阵；

(2)将各节点电压设初值；

(3)将节点初值代入相关求式，求出修正方程式的常数项向量；

(4)将节点电压初值代入求式，求出雅可比矩阵元素；

(5)求解修正方程，求修正向量；

(6)求取节点电压的新值；

(7)检查是否收敛，如不收敛，则以各节点电压的新值作为初值自第3步重新开始进行狭义次迭代，否则转入下一步；

(8)计算支路功率分布，PV节点无功功率和平衡节点功率。

其用程序仿真的过程可以由以下流程图来简单表示出来：

三、Matlab程序

1、原代码

%潮流计算

fprintf('开始潮流计算\n');

fprintf('请输入待求网络的相应参数\n');

%参数输入部分

n=input('网络中的节点数:n=');

L=input('网络中的支路数:L=');

ss=input('平衡节点ss=');

pr=input('误差精度:pr=');

X1=input('支路参数:X1=');

X2=input('节点参数:X2=');

X=input('节点号和对地参数:X=');

fprintf('参数输入部分结束\n\n');

Y=zeros(n);

%置迭代次数

mm=1;

%创建节点导纳矩阵

for i=1:L

if X1(i,6)==0 %不含变压器的支路

p=X1(i,1);

q=X1(i,2);

Y(p,q)=Y(p,q)-1/X1(i,3);

Y(q,p)=Y(p,q);

Y(p,p)=Y(p,p)+1/X1(i,3)+0.5*X1(i,4);

Y(q,q)=Y(q,q)+1/X1(i,3)+0.5*X1(i,4);

else %含有变压器的支路

p=X1(i,1);

q=X1(i,2);

Y(p,q)=Y(p,q)-1/(X1(i,3)*X1(i,5));

Y(q,p)=Y(p,q);

Y(p,p)=Y(p,p)+1/X1(i,3);

Y(q,q)=Y(q,q)+1/(X1(i,5)^2*X1(i,3));

end

OrgS=zeros(2*n-2,1);

DetaS=zeros(2*n-2,1); %将OrgS、DetaS初始化

%创建OrgS，用于存储初始功率参数

h=0;

j=0;

for i=1:n %对PQ节点的处理

if i~=ss&&X2(i,6)==2

h=h+1;

for j=1:n

OrgS(2*h-1,1)=OrgS(2*h-1,1)+real(X2(i,3))*(real(Y(i,j))*real(X2(j,3))-imag(Y(i,j))*imag(X2(j,3)))+imag(X2(i,3))*(real(Y(i,j))*imag(X2(j,3))+imag(Y(i,j))*real(X2(j,3)));

OrgS(2*h,1)=OrgS(2*h,1)+imag(X2(i,3))*(real(Y(i,j))*real(X2(j,3))-imag(Y(i,j))*imag(X2(j,3)))-real(X2(i,3))*(real(Y(i,j))*imag(X2(j,3))+imag(Y(i,j))*real(X2(j,3)));

end

for i=1:n %对PV节点的处理，注意这时不可再将h初始化为0

if i~=ss&&X2(i,6)==0

h=h+1;

for j=1:n

OrgS(2*h-1,1)=OrgS(2*h-1,1)+real(X2(i,3))*(real(Y(i,j))*real(X2(j,3))-imag(Y(i,j))*imag(X2(j,3)))+imag(X2(i,3))*(real(Y(i,j))*imag(X2(j,3))+imag(Y(i,j))*real(X2(j,3)));

OrgS(2*h,1)=OrgS(2*h,1)+imag(X2(i,3))*(real(Y(i,j))*real(X2(j,3))-imag(Y(i,j))*imag(X2(j,3)))-real(X2(i,3))*(real(Y(i,j))*imag(X2(j,3))+imag(Y(i,j))*real(X2(j,3)));

end

OrgS;

%创建PVU 用于存储PV节点的初始电压

PVU=zeros(n-h-1,1);

t=0;

for i=1:n

if X2(i,6)==0

t=t+1;

PVU(t,1)=X2(i,3);

end

PVU;

%创建DetaS，用于存储有功功率、无功功率和电压幅值的不平衡量

h=0;

for i=1:n %对PQ节点的处理

if i~=ss&&X2(i,6)==2

h=h+1;

DetaS(2*h-1,1)=real(X2(i,2))-OrgS(2*h-1,1);

DetaS(2*h,1)=imag(X2(i,2))-OrgS(2*h,1);

end

t=0;

for i=1:n %对PV节点的处理，注意这时不可再将h初始化为0

if i~=ss&&X2(i,6)==0

h=h+1;

t=t+1;

DetaS(2*h-1,1)=real(X2(i,2))-OrgS(2*h-1,1);

DetaS(2*h,1)=real(PVU(t,1))^2+imag(PVU(t,1))^2-real(X2(i,3))^2-imag(X2(i,3))^2;

end

DetaS;

%创建I，用于存储节点电流参数

i=zeros(n-1,1);

h=0;

for i=1:n

if i~=ss

h=h+1;

I(h,1)=(OrgS(2*h-1,1)-OrgS(2*h,1)*sqrt(-1))/conj(X2(i,3));

end

%创建Jacbi(雅可比矩阵)

Jacbi=zeros(2*n-2);

h=0;

k=0;

for i=1:n %对PQ节点的处理

if X2(i,6)==2

h=h+1;

for j=1:n

if j~=ss

k=k+1;

if i==j %对角元素的处理

Jacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(X2(i,3))+real(Y(i,j))*imag(X2(i,3))+imag(I(h,1));

Jacbi(2*h-1,2*k)=real(Y(i,j))*real(X2(i,3))+imag(Y(i,j))*imag(X2(i,3))+real(I(h,1));

Jacbi(2*h,2*k-1)=-Jacbi(2*h-1,2*k)+2*real(I(h,1));

Jacbi(2*h,2*k)=Jacbi(2*h-1,2*k-1)-2*imag(I(h,1));

else %非对角元素的处理

Jacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(X2(i,3))+real(Y(i,j))*imag(X2(i,3));

Jacbi(2*h-1,2*k)=real(Y(i,j))*real(X2(i,3))+imag(Y(i,j))*imag(X2(i,3));

Jacbi(2*h,2*k-1)=-Jacbi(2*h-1,2*k);

Jacbi(2*h,2*k)=Jacbi(2*h-1,2*k-1);

end

if k==(n-1) %将用于内循环的指针置于初始值，以确保雅可比矩阵换行

k=0;

end

k=0;

for i=1:n %对PV节点的处理

if X2(i,6)==0

h=h+1;

for j=1:n

if j~=ss

k=k+1;

if i==j %对角元素的处理

Jacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(X2(i,3))+real(Y(i,j))*imag(X2(i,3))+imag(I(h,1));

Jacbi(2*h-1,2*k)=real(Y(i,j))*real(X2(i,3))+imag(Y(i,j))*imag(X2(i,3))+real(I(h,1));

Jacbi(2*h,2*k-1)=2*imag(X2(i,3));

Jacbi(2*h,2*k)=2*real(X2(i,3));

else %非对角元素的处理

Jacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(X2(i,3))+real(Y(i,j))*imag(X2(i,3));

Jacbi(2*h-1,2*k)=real(Y(i,j))*real(X2(i,3))+imag(Y(i,j))*imag(X2(i,3));

Jacbi(2*h,2*k-1)=0;

Jacbi(2*h,2*k)=0;

end

if k==(n-1) %将用于内循环的指针置于初始值，以确保雅可比矩阵换行

k=0;

end

Jacbi;

%求解修正方程，获取节点电压的不平衡量

DetaU=zeros(2*n-2,1);

DetaU=inv(Jacbi)*DetaS;

DetaU;

%修正节点电压

j=0;

for i=1:n %对PQ节点处理

if X2(i,6)==2

j=j+1;

X2(i,3)=X2(i,3)+DetaU(2*j,1)+DetaU(2*j-1,1)*sqrt(-1);

end

for i=1:n %对PV节点的处理

if X2(i,6)==0

j=j+1;

X2(i,3)=X2(i,3)+DetaU(2*j,1)+DetaU(2*j-1,1)*sqrt(-1);

end

X2;

%开始循环**********************************************************************

while abs(max(DetaU))>pr

OrgS=zeros(2*n-2,1); %!!!初始功率参数在迭代过程中是不累加的，所以在这里必须将其初始化为零矩阵

h=0;

j=0;

for i=1:n

if i~=ss&&X2(i,6)==2

h=h+1;

for j=1:n

OrgS(2*h-1,1)=OrgS(2*h-1,1)+real(X2(i,3))*(real(Y(i,j))*real(X2(j,3))-imag(Y(i,j))*imag(X2(j,3)))+imag(X2(i,3))*(real(Y(i,j))*imag(X2(j,3))+imag(Y(i,j))*real(X2(j,3)));

OrgS(2*h,1)=OrgS(2*h,1)+imag(X2(i,3))*(real(Y(i,j))*real(X2(j,3))-imag(Y(i,j))*imag(X2(j,3)))-real(X2(i,3))*(real(Y(i,j))*imag(X2(j,3))+imag(Y(i,j))*real(X2(j,3)));

end

for i=1:n

if i~=ss&&X2(i,6)==0

h=h+1;

for j=1:n

OrgS(2*h-1,1)=OrgS(2*h-1,1)+real(X2(i,3))*(real(Y(i,j))*real(X2(j,3))-imag(Y(i,j))*imag(X2(j,3)))+imag(X2(i,3))*(real(Y(i,j))*imag(X2(j,3))+imag(Y(i,j))*real(X2(j,3)));

OrgS(2*h,1)=OrgS(2*h,1)+imag(X2(i,3))*(real(Y(i,j))*real(X2(j,3))-imag(Y(i,j))*imag(X2(j,3)))-real(X2(i,3))*(real(Y(i,j))*imag(X2(j,3))+imag(Y(i,j))*real(X2(j,3)));

end

OrgS;

%创建DetaS

h=0;

for i=1:n

if i~=ss&&X2(i,6)==2

h=h+1;

DetaS(2*h-1,1)=real(X2(i,2))-OrgS(2*h-1,1);

DetaS(2*h,1)=imag(X2(i,2))-OrgS(2*h,1);

end

t=0;

for i=1:n

if i~=ss&&X2(i,6)==0

h=h+1;

t=t+1;

DetaS(2*h-1,1)=real(X2(i,2))-OrgS(2*h-1,1);

DetaS(2*h,1)=real(PVU(t,1))^2+imag(PVU(t,1))^2-real(X2(i,3))^2-imag(X2(i,3))^2;

end

DetaS;

%创建I

i=zeros(n-1,1);

h=0;

for i=1:n

if i~=ss

h=h+1;

I(h,1)=(OrgS(2*h-1,1)-OrgS(2*h,1)*sqrt(-1))/conj(X2(i,3));

end

%创建Jacbi

Jacbi=zeros(2*n-2);

h=0;

k=0;

for i=1:n

if X2(i,6)==2

h=h+1;

for j=1:n

if j~=ss

k=k+1;

if i==j

Jacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(X2(i,3))+real(Y(i,j))*imag(X2(i,3))+imag(I(h,1));

Jacbi(2*h-1,2*k)=real(Y(i,j))*real(X2(i,3))+imag(Y(i,j))*imag(X2(i,3))+real(I(h,1));

Jacbi(2*h,2*k-1)=-Jacbi(2*h-1,2*k)+2*real(I(h,1));

Jacbi(2*h,2*k)=Jacbi(2*h-1,2*k-1)-2*imag(I(h,1));

else

Jacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(X2(i,3))+real(Y(i,j))*imag(X2(i,3));

Jacbi(2*h-1,2*k)=real(Y(i,j))*real(X2(i,3))+imag(Y(i,j))*imag(X2(i,3));

Jacbi(2*h,2*k-1)=-Jacbi(2*h-1,2*k);

Jacbi(2*h,2*k)=Jacbi(2*h-1,2*k-1);

end

if k==(n-1)

k=0;

end

k=0;

for i=1:n

if X2(i,6)==0

h=h+1;

for j=1:n

if j~=ss

k=k+1;

if i==j

Jacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(X2(i,3))+real(Y(i,j))*imag(X2(i,3))+imag(I(h,1));

Jacbi(2*h-1,2*k)=real(Y(i,j))*real(X2(i,3))+imag(Y(i,j))*imag(X2(i,3))+real(I(h,1));

Jacbi(2*h,2*k-1)=2*imag(X2(i,3));

Jacbi(2*h,2*k)=2*real(X2(i,3));

else

Jacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(X2(i,3))+real(Y(i,j))*imag(X2(i,3));

Jacbi(2*h-1,2*k)=real(Y(i,j))*real(X2(i,3))+imag(Y(i,j))*imag(X2(i,3));

Jacbi(2*h,2*k-1)=0;

Jacbi(2*h,2*k)=0;

end

if k==(n-1)

k=0;

end

Jacbi

DetaU=zeros(2*n-2,1);

DetaU=inv(Jacbi)*DetaS;

DetaU

%修正节点电压

j=0;

for i=1:n

if X2(i,6)==2

j=j+1;

X2(i,3)=X2(i,3)+DetaU(2*j,1)+DetaU(2*j-1,1)*sqrt(-1);

end

for i=1:n

if X2(i,6)==0

j=j+1;

X2(i,3)=X2(i,3)+DetaU(2*j,1)+DetaU(2*j-1,1)*sqrt(-1);

end

mm=mm+1; %迭代次数加1

end

2、数据输入

开始潮流计算

请输入待求网络的相应参数

网络中的节点数:n=22

网络中的支路数:L=40

平衡节点ss=1

误差精度:pr=0.0001

支路参数:X1=[1,7,0.0+0.015i,0.0+0.0i,1,1.05;

2,9.,0.0+0.0217i,0.0+0.0i,1,1.075;

3,2.,0.0+0.0124i,0.0+0.0i,1,1.1;

4,19,0.0+0.064i,0.0+0.0i,1,1.025;

5,18,0.0+0.0375i,0.0+0.0i,1,1.05;

6,17,0.0+0.0337i,0.0+0.0i,1,1;

7,8,0.0106+0.074i,0.0+0.0i,1,1;

7,9,0.0147+0.104i,0.0+0.0i,1,1;

8,9,0.0034+0.0131i,0.0+0.0i,1,1;

8,22,0.0537+0.19i,0.0+0.0i,1,1;

8,22,0.0+0.1653i,0.0+0.0i,1,0;

9,10,0.0+0.002i,0.0+0.0i,1,1;

9,22,0.0559+0.218i,0.0+0.0i,1,1;

9,22,0.0+0.1954i,0.0+0.0i,1,0;

10,11,0.0+0.018i,0.0+0.0i,1,1;

11,11,0.0-1.3675i,0.0+0.0i,1,0;

11,12,0.0033+0.0343i,0.0+0.0i,1,1;

11,12,0.0+1.8797i,0.0+0.0i,1,0;

12,12,0.0-1.3675i,0.0+0.0i,1,0;

12,13,0.00245+0.0255i,0.0+0.0i,1,1;

12,13,0.0+1.395i,0.0+0.0i,1,0;

12,15,0.0+0.018i,0.0+0.0i,1,1;

13,17,0.0+0.01i,0.0+0.0i,1,1;

14,15,0.0-0.002i,0.0+0.0i,1,1;

14,19,0.0034+0.02i,0.0+0.0i,1,1;

16,16,0.0+0.5018i,0.0+0.0i,1,0;

16,17,0.0+0.001i,0.0+0.0i,1,0.975;

16,18,0.0033+0.0333i,0.0+0.0i,1,1;

16,19,0.0578+0.218i,0.0+0.0i,1,1;

16,19,0.0+0.1887i,0.0+0.0i,1,0;

16,20,0.0165+0.0662i,0.0+0.0i,1,1;

16,20,0.0+0.2353i,0.0+0.0i,1,0;

16,21,0.0374+0.178i,0.0+0.0i,1,1;

16,21,0.0+0.164i,0.0+0.0i,1,0;

19,21,0.0114+0.037i,0.0+0.0i,1,1;

20,22,0.0214+0.0859i,0.0+0.0i,1,1;

20,22,0.0+0.3008i,0.0+0.0i,1,0;

21,22,0.015+0.0607i,0.0+0.0i,1,1;

21,22,0.0+0.2198i,0.0+0.0i,1,0]

节点参数:X2=[0, 0.0+0.0i ,1.0,0.0,0,1;

0,6.0+3.2i,1.0,0.0,0,2;

0,3.1+0.0i,1.0,1.0,0,0;

0,1.6+0.7i,1.0,0.0,0,2;

0,4.3+3.34i,1.0,0.0,0,2;

0,-0.01+0.0i,1.0,0.0,0,0;

0,0.0+0.0i,1.0,0.0,0,2;

0,-2.87-1.44i,1.0,0.0,0,2;

0,-3.176-2.21i,1.0,0.0,0,2;

0,0.183+0.0i,1.0,0.0,0,2;

0,0.0+0.0i,1.0,0.0,0,2;

0,0.0+0.0i,1.0,0.0,0, 2;

0,-5-2.9i,1.0,0.0,0, 2;

0,0.0+0.0i,1.0,0.0,0, 2;

0,-4.3-2.6i,1.0,0.0,0, 2;

0,-0.864-0.662i,1.0,0.0,0, 2;

0,-0.719-0.474i,1.0,0.0,0, 2;

0,-0.7-0.5i,1.0,0.0,0, 2;

0,-2.265-1.69i,1.0,0.0,0, 2]

节点号和对地参数:

X=[1,0;2,0;3,0;4,0;5,0;6,0;7,0;8,0;9,0;10,0;11,0;12,0;13,0;14,0;15,0;16,0;17,0;18,0;19,0;20,0;21,0;22,0]

参数输入部分结束

四、仿真结果

由于结果太长，粘入后一部分。

DetaU =

1.0e+152 *

0.0000

-0.0000

0.8278

1.2904

-0.6757

-0.5324

0.0138

0.0600

-0.3561

0.7466

-0.0086

-0.5668

-0.0710

0.7725

-0.1085

-4.6600

-8.8086

0.7060

-1.1456

0.3951

-0.0042

0.0181

-0.4711

-0.5285

4.4101

1.0732

5.4030

1.8329

0.1459

1.2651

0.0512

1.2975

-0.1952

0.0423

-0.5872

0.4488

0.1269

0.6914

-0.1404

0.4546

0.1120

0.0489

X2 =

1.0e+152 *

Columns 1 through 3

0 0 0.0000

0 0.0000 + 0.0000i -0.0000 + 0.0000i

0 0.0000 -0.6061 - 2.4840i

0 0.0000 + 0.0000i 1.2324 + 0.8115i

0 0.0000 + 0.0000i 2.1905 + 0.4803i

0 -0.0000 0.0700 + 0.1603i

0 0 0.2226 + 0.0678i

0 -0.0000 - 0.0000i 0.4741 + 0.1168i

0 -0.0000 - 0.0000i -0.7678 + 0.3106i

0 0.0000 0.7060 - 0.1106i

0 0 -3.7987 + 3.0737i

0 0 1.4448 - 5.2542i

0 0 2.1504 + 0.6668i

0 0 -0.7980 - 0.1097i

0 0 -1.6845 - 1.2253i

0 -0.0000 - 0.0000i 1.0346 + 4.5579i

0 0 1.7606 + 5.3485i

0 -0.0000 - 0.0000i 0.3183 + 0.0898i

0 -0.0000 - 0.0000i 0.4975 + 0.0130i

0 -0.0000 - 0.0000i 0.0975 - 0.2849i

0 -0.0000 - 0.0000i 0.1583 - 0.8069i

0 -0.0000 - 0.0000i 0.5199 - 0.1895i

Columns 4 through 6

0 0 0.0000

0.0000 0 0

0 0 0.0000

0 0 0

0 0 0.0000

Jacbi =

1.0e+155 *

Columns 1 through 6

-0.0843 0.1860 0.0000 0.0000 0 0

0.1860 0.0843 -0.0000 0.0000 0 0

0 0 0 0 0.0192 -0.0225

0 0 0 0 0.0029 0.0193

0 0 0 0 0 0

0.0354 0.0143 0 0 0 0

-0.0143 0.0354 0 0 0 0

0 0 0 0 0 0

0 0 0 0 -0.0078 0.0002

0 0 0 0 -0.0002 -0.0078

0 0 0 0 0 0

0.0489 -0.2003 -0.1723 -0.4288 0 0

0 0 -0.0050 -0.0012 0 0

0 0 0 0 0 0

Columns 7 through 12

0 0 0 0 0 0

-0.0068 0.0292 0 0 0 0

0.0548 0.1236 0 0 0 0

0 0 0 0 0.0179 -0.0053

0 0 0 0 0.0054 0.0223

0 0 0 0 -0.0065 0.0006

0 0 0 0 -0.0006 -0.0065

0 0 0 0 0.0068 0.0040

0 0 0 0 -0.0040 0.0068

0 0 0 0 0 0

0 0 -0.0522 0.1587 0 0

0 0 -0.1587 -0.0522 0 0

-0.0085 0.0024 0 0 0 0

-0.0024 -0.0085 0 0 0 0

0 0 0 0 0 0

0 0 0.0516 0.0080 0 0

0 0 0.0003 0.0001 0 0

Columns 13 through 18

0 0 0.0000 0.0000 0 0

0 0 -0.0000 0.0000 0 0

0 0 0 0 0 0

-0.0031 0.0005 -0.0022 0.0003 0 0

-0.0005 -0.0031 -0.0003 -0.0022 0 0

-0.0827 0.2322 -0.0361 -0.0005 0 0

0.2339 0.1785 0.0005 -0.0361 0 0

0.0491 0.0365 -3.4291 1.1630 0.3839 0.1553

-0.0365 0.0491 1.5906 2.4649 -0.1553 0.3839

0 0 -0.3530 -0.0553 -1.1595 0.0968

0 0 0.0553 -0.3530 -0.0261 1.9439

0 0 0 0 0.2110 0.1708

0 0 0 0 -0.1708 0.2110

0 0 0 0 0 0

-0.0054 -0.0028 -0.0047 -0.0024 0 0

0.0028 -0.0054 0.0024 -0.0047 0 0

0 0 0 0 0 0

Columns 19 through 24

0 0 0 0 0 0

-0.0392 -0.0061 0 0 0 0

0.0061 -0.0392 0 0 0 0

2.7833 1.2855 0.1032 0.1010 0 0

1.8245 -3.4062 -0.1010 0.1032 0 0

-0.0279 -0.1586 -5.6838 0.8997 -0.0376 -0.2133

0.1586 -0.0279 -0.4125 5.9710 0.2133 -0.0376

0 0 -0.0876 0.0184 0.2256 -0.0086

0 0 -0.0184 -0.0876 0.1614 0.3797

0 0 0 0 0 0

0 0 0.0936 -0.0681 0 0

0 0 0.0681 0.0936 0 0

0 0 0 0 0 0

0 0 0 0 -0.1761 0.5348

0 0 0 0 -0.5348 -0.1761

0 0 0 0 0 0

Columns 25 through 30

0 0 0 0 0 0

0 0 -0.0803 -0.2919 0 0

0 0 0.2919 -0.0803 0 0

0 0 0 0 0 0

0.2653 -0.1190 -0.3990 0.0549 0 0

-0.0068 0.4533 -0.0549 -0.3990 0 0

-0.8422 0.6127 1.2177 -1.1488 0 0

-0.6127 -0.8422 -0.0596 0.2796 0 0

0 0 0 0 NaN NaN

0 0 0 0 -1.7606 5.3485

0 0 0 0 -5.3485 -1.7606

0 0 0 0 -0.0097 0.0017

0 0 0 0 -0.0017 -0.0097

-0.0243 -0.0035 0 0 -0.0048 -0.0004

0.0035 -0.0243 0 0 0.0004 -0.0048

0 0 0 0 -0.0008 -0.0056

0 0 0 0 0.0056 -0.0008

0 0 0 0 -0.0009 -0.0094

0 0 0 0 0.0094 -0.0009

0 0 0 0 0 0

Columns 31 through 36

0 0 0 0 0 0

0 0 0 0 -0.0193 0.0127

0 0 0 0 -0.0127 -0.0193

0 0 -0.0584 0.0128 0 0

0 0 -0.0128 -0.0584 0 0

0 0 0 0 0 0

-0.2150 0.0667 0 0 0 0

-0.0667 -0.2150 0 0 0 0

0 0 0 0 0.0397 0.0013

0 0 0 0 -0.0013 0.0397

0 0 0 0 0 0

-1.0346 4.5579 -0.0442 0.1325 -0.0151 0.0425

-4.5579 -1.0346 -0.1325 -0.0442 -0.0425 -0.0151

NaN NaN 0 0 0 0

0 0 0.0111 -0.0127 0 0

0 0 -0.0044 0.0253 0 0

0 0 0 0 0.0462 0.0110

0 0 0 0 -0.0033 0.0522

0 0 0 0 0 0

0 0 0 0 0.0022 -0.0211

0 0 0 0 0.0211 0.0022

0 0 0 0 0 0

-0.0021 0.0048 0 0 0 0

0 0 0 0 0 0

Columns 37 through 42

0 0 0 0 0 0

0 0 0 0 -0.0053 0.0006

0 0 0 0 -0.0006 -0.0053

0 0 0 0 0.0069 0.0038

0 0 0 0 -0.0038 0.0069

0 0 0 0 0 0

-0.0353 0.0805 -0.0170 0.0511 0 0

-0.0805 -0.0353 -0.0511 -0.0170 0 0

0 0 0 0 0 0

0 0 -0.0124 -0.0035 0 0

0 0 0.0035 -0.0124 0 0

-0.0849 -0.0002 0 0 -0.0006 -0.0043

-0.0201 0.0878 0 0 0.0043 -0.0006

0 0 -0.1598 0.0021 -0.0001 -0.0168

0 0 -0.0926 0.1573 0.0168 -0.0001

-0.0069 -0.0041 -0.0097 -0.0058 -0.0152 -0.0084

0.0041 -0.0069 0.0058 -0.0097 -0.0385 0.0686

0 0 0 0 0 0

X2 =

Columns 1 through 3

0 0 1.0000

0 6.0000 + 3.2000i NaN + NaNi

0 3.1000 NaN + NaNi

0 1.6000 + 0.7000i NaN + NaNi

0 4.3000 + 3.3400i NaN + NaNi

0 -0.0100 NaN + NaNi

0 0 NaN + NaNi

0 -2.8700 - 1.4400i NaN + NaNi

0 -3.1760 - 2.2100i NaN + NaNi

0 0.1830 NaN + NaNi

0 0 NaN + NaNi

0 -5.0000 - 2.9000i NaN + NaNi

0 0 NaN + NaNi

0 -4.3000 - 2.6000i NaN + NaNi

0 -0.8640 - 0.6620i NaN + NaNi

0 -0.7190 - 0.4740i NaN + NaNi

0 -0.7000 - 0.5000i NaN + NaNi

0 -2.2650 - 1.6900i NaN + NaNi

Columns 4 through 6

0 0 1.0000

0 0 2.0000

1.0000 0 0

0 0 2.0000

0 0 0

0 0 2.0000

mm =

五、实验心得

经过这次潮流计算编程的相关学习，我对于利用Matlab计算电力系统潮流的方法和思想都有了一定的了解，同时也知道了Matlab编程的基本思想。在使用Matlab仿真的过程中我遇到了不少问题，比如矩阵输入的疏忽使得行列不对等就让我重新检查了好几遍源程序和输入。还有在参数设定的时候前后没有统一也造成了很大的麻烦。由此可见，不管做什么都要认真细致，不可忽略细节。还有可能因为误差不同而造成的结果差异较大而不断进行改进和检查，反复研究，在进步的同时也加深对Matlab软件的了解，使用的也会更为顺手。

总体说来，感觉收获很大，不仅加深了对电力系统潮流计算步骤的理解，也使对Matlab工具的应用更加熟练，相信以后的学习过程中这次的经历也会带来很大的帮助，加深了学习的乐趣和自己学习的热情。

“路漫漫其修远兮，吾将上下而求索”!这次大作业的Matlab仿真并不是学习的结束，而是新的开始！认真对待每一份作业，并且花时间去琢磨和钻研，总会有收获的。

更多相关推荐：

计算机实验二实验报告: 深圳大学实验报告课程名称：计算机导论项目名称：操作系统与工具软件学院：专业：指导教师：报告人学号：实验时间：20**年10月15日提交时间：20**年10月18日教务处制一、实验目的与要求1．掌握Windows…
计算机实验报告模板: 实验报告实验序号实验项目名称
计算机实验报告: 成绩计算机原理实验室实验报告课程汇编语言与微机原理姓名马路伟专业测控技术与仪器学号1220xx111日期20xx年11月太原工业学院电子工程系实验一输入输出实验123计算机原理实验室实验报告实验二分支程序实验4...
计算机实验报告格式: 电子信息工程学系实验报告适用于计算机课程课程名称面向对象程序设计实验项目名称Visualstudioc60集开发环境的使用实验时间班级计教101姓名蔡静学号实验目的1熟悉并学习使用C程序编译平台VC602掌握如...
计算机网络上机实验报告: 计算机网络课程设计论文设计论文题目计算机网络综合实习课程设计学院名称信息科学与技术学院专业名称通信工程学生姓名杜立华学生学号20xx13070112任课教师王权海陈红娟设计论文成绩教务处制20xx年4月18日通...
电大计算机实验报告: 实验一WindowsXP基本操作11实验目的掌握使用WindowsXP的基本方法12实验要求1观察WindowsXP的基本方法2掌握WindowsXP的启动和关闭方法3掌握鼠标器的操作方法4掌握窗口图标菜单工具...
大一计算机实验报告: 大学计算机基础实验报告20xx20xx学年第1学期实验成绩

计算机实验报告: 四川大学网络教育学院实验报告实验名称计算机综合实践班级姓名学号DH1091420xx四川大爱网络教育学院计算机综合实践实验报告实验内容1熟悉使用DEBUG调试工具借助DEBUG调试工具来发现汇编语言程序的错误所...
计算机模拟手工实验报告: 计算机模拟手工实验学生实验报告学院商学院课程名称会计模拟手工会计专业班级财务管理2班姓名肖鹏学号20xx1222238郑州大学西亚斯国际学院商学院学生实验报告第一部分实验概况与内容一实验的目的及要求实验目的通过...
计算机模拟手工实验报告: 计算机模拟手工实验学生实验报告学院商学院课程名称计算机模拟手工实验专业班级姓名学号学生实验报告第一部分实验概况与内容一实验的目的及要求1实验目的本实验以模拟企业的实际会计工作为基础按照企业会计制度和企业会计准则...

大学计算机基础实验报告模版: 武汉商学院信息工程系大学计算机基础课程实验报告信息工程系20xx年编制武汉商学院信息工程系实验报告编写说明1实验名称填写以所用实验指导书内容为准2实验目的和要求实验项目摘要填写以所用实验指导书内容为准3实验预习...
大学计算机基础实验报告模板: 计算机基础实践报告学院计算机与信息学院专业班级电子信息工程二班学号20xx2661学生姓名陈雪莹实验的目的和要求通过上机实验能够熟练掌握计算机的基本操作技能该实验安排内容包括Windows的基本操作使用Word...

热门关注

北理大学计算机实验基础实验二实验报告表答案

计算机网络编程实验报告

《财务核算计算机模拟》实验报告书模板

《大学计算机基础》实验报告模板

电力系统潮流上机计算报告书写要求

第二篇：电力系统 潮流计算仿真报告

第二篇：电力系统潮流计算仿真报告