2 网络层次分析
对于辐射型网络,前推回代法的基本原理是:① 假定节点电压不变,已知网络末端功率,由网络末端向首端计算支路功率损耗和支路功率,得到根节点注入功率;② 假定支路功率不变,已知根节点电压,由网络首端向末端计算支路电压损耗和节点电压。前推时,每条支路的功率都由该支路的下一层支路功率决定,回代时,节点电压都由上一层节点决定。这种特点一方面限制了不同层次间的功率前推和电压回代不能同时进行,另一方面也说明同一层次的支路功率之间没有前后关联,因此同一层次内完全可以实现功率或电压的并行计算。尤其对于大规模辐射型网络,由于分层数显著少于支路总数,所以分层后能够充分发挥并行计算的优势,提高计算速度。
以一个简单的11节点树状网为例,其节点和支路编号采用与网络结构无关的自然编号(即从1开始的自然数顺序编号),在这个网络中,支路1、2、5属于同一层次,当计算支路功率损耗和电压损耗时,彼此不相关,可以并行计算。同样,支路3、4、6、7、8、10也属同一层,其功率损耗和电压损耗也可以并行计算。这样,根据图1网络的拓扑结构,可以直观地看到网络支路共分为3层,且可以知道每一层的支路情况以及每一支路的送端节点和受端节点情况。显然,只要了解了这些信息,就能够分层实现功率前推和电压回代的并行计算,而且无需对节点和支路重新编号。
为了描述以上的网络层次信息,定义如下:
(1)网络层次矩阵L:
设网络分为Li层,每层包含的支路数最多为M,则网络层次矩阵L是1个(Li´M)矩阵,第i行的非零元素就是网络第i层包含的支路编号,非零元素的个数就是该层包含的支路数。从L1层到Li层代表了功率流动的方向,前推时从Li层到L1层,回代时从 L1层到 Li层。
(2)支路送端节点矩阵f和受端节点矩阵t:
由于原始数据中支路的首节点到末节点的方向不一定就是功率流向,因此必须根据功率方向来确定支路的送端节点和受端节点。每条支路上的功率都由该支路的送端流向受端,支路送端节点矩阵和受端节点矩阵都是一维矩阵,元素个数等于支路数,第i个元素就是支路i的送端(受端)节点编号。
(3)支路层次关联矩阵C:
设网络支路数为N,支路层次关联矩阵为1个(N´N)的矩阵。矩阵第i行j列元素为1,表示支路i与支路j为上下层关系,它们直接相连,且支路i的上层支路是支路j,支路j的下层支路是支路i。当支路间没有这种直接的上下层次关系时,对应的元素为0。
在上述几个矩阵中,以网络层次矩阵描述整个网络的支路分层情况。支路送端节点矩阵和受端节点矩阵反映每条支路与送端、受端节点的关联关系。支路层次关联矩阵反映的是支路之间的直接上下层次关系。
下面介绍一下如何分析网络结构,以形成这几个矩阵:
(1)形成网络层次矩阵L、支路送端节点矩阵f和受端节点矩阵t。
进行网络层次分析时,首先形成节点-支路关联矩阵。若节点数为N,则辐射型网络的支路数必定为N-1,节点-支路关联矩阵是1个N´N-1矩阵。当节点i与支路j相连时,则关联矩阵的i行j列元素为1,不相连时,则该元素为0。由此形成的节点-支路关联矩阵,每一列有两个非零元素,其对应的行号就是该列支路的两端节点编号;每一行的非零元素对应的列号就是与该行节点相连的支路编号式中 行表示节点1~11;列表示支路1~10。
从根节点7,即矩阵A的第7行出发,仅找到第9列的元素为1,即节点7仅与支路9相连。与根节点相连的所有支路都属于第1层支路,且根节点为送端节点,所以第1层支路为支路9,支路9的送端节点为7,受端节点为支路9的另一端节点,即矩阵第9列上另一个非零元素对应的节点1,这就是网络的第1层分析。
从网络第1层支路的所有受端节点出发,与它们相连的其他所有支路属于第2层支路,第2层支路的送端节点为第1层支路的受端节点。从节点1出发,查找矩阵A第1行的元素,找到第1、2、5列元素为1,故第2层支路为支路1、2、5,它们的送端节点为节点1,受端节点分别为相应列上另一个非零元素对应的节点。依次查找下去,沿着矩阵A中的轨迹可以整理出整个网络的层次结构和每条支路的送端、受端节点,其中实线表示第1层分析轨迹,虚线表示第2层分析轨迹,点划线表示第三层分析轨迹。
网络层次分析以后,形成的网络层次矩阵L(行表示L1~L3层)、支路送端节点矩阵f和受端节点矩阵t(2)形成支路层次关联矩阵C。
除了第1层支路没有上层支路外,任意1条支路只有1条直接相连的上层支路,而且始终遵循这样的原则:该支路的送端节点就是与其直接相连的上层支路的受端节点。通过支路送端节点和受端节点矩阵,可以很容易地找到任意一条支路的直接上层支路,比如由支路送端节点矩阵f找到任意支路i的送端节点bi,然后由支路受端节点矩阵t找到受端节点为bi的支路j,这就意味着支路i的上一层支路为支路j,即矩阵的i行j列元素为1。
查找每条支路的直接上层支路,可形成支路层次关联矩阵,图1网络的支路层次关联矩阵C为式中 行表示节点1~10;列表示支路1~10。
从支路层次关联矩阵可以查找任意支路的上层支路和下层支路。如果需要查找支路j的上一层支路,只需要知道矩阵C的第j行为1的元素所在的列就可以了;同样,如果需要查找支路j的下一层支路,只需要知道矩阵C的第j列为1的元素所在的行就可以了。例如,从矩阵C的第1行可知支路1的上一层支路为支路9,从矩阵第2列可知支路2的下一层支路为支路6、7,等等。
3 变压器模型
当网络中存在变压器时,通常采用图2(b)所示的P型等值电路来等效图2(a)的变压器支路。但是,通过多次计算分析发现,当网络中存在三绕组降压变压器时,由于中压侧等效绕组的阻抗普遍很小(常常是很小的负阻抗),所以当中压侧变比 时,将会产生很大的对地导纳,导致前推回代法不收敛。
下面以一简单的例子分析之,该例只有一个三绕组变压器的树状网络,高压端为根节点,中压和低压端接负荷,变压器型号为SFPZ9-180000/220 (220±8 ´1.25%/121/10.5, 180/180/90), 变压器参数为SB=100MVA;基准电压为220 kV / 110 kV /10kVY0 = 0.000748-j0.000799puVH = 1.05pu;SM = 0.09+ j0.03pu;SL =0.04-j0.04pu当三侧等效双绕组支路采用P型等值电路时,前推回代法不收敛。如果把中压侧的变比改为1,而阻抗不变, 或把中压侧的阻抗增大为低压侧或高压侧的阻抗,而变比不变,前推回代法都可以收敛。对同样的网络,我们又测试了其它型号的三绕组变压器,而且还改变了功率和电压,都得到同样的结论。所以可以确定不收敛的原因是中压侧的非标准变比和小阻抗联合作用产生的较大的对地导纳。
为了解决∏型等效模型产生的不收敛问题,本文根据理想变压器只改变电压、不改变传送功率的原理,提出了一种新的电压变换模型来处理变压器支路,并推导了在前推和回代时的公式,具体如下:
对于变压器支路,根据功率的流向,存在升压和降压两种方式,;4 分层前推回代法的主要步骤
同时考虑对地支路、线路支路、升压变压器和降压变压器支路的分层前推回代算法如下:
(1)功率前推
设支路受端计算电压V为式中 φ为与该支路相连的下层支路集合。
支路送端功率为
根据网络层次矩阵,从网络的第1层回代到第L层,逐层更新支路受端节点的电压,也即更新了下一层支路的送端节点电压。
式(5)~(10)中,Vf为支路送端节点电压;Vt为支路受端节点电压;由支路送端节点矩阵和受端节点矩阵可以容易得到;Y为支路受端节点对地导纳;S0为支路受端节点负荷;Z为支路阻抗;S¢是支路受端功率;S为支路送端功率;k为变压器支路变比;*表示共轭。
在MATLAB环境下,以上的功率前推和电压回代计算,都可以直接利用其快速的复数矩阵运算功能来实现。此时,式(5)~(10)中的变量都是复数的矩阵变量,它们可以直接进行相关的代数运算,其中,乘、除和乘方运算都使用点乘、点除和点乘方的方式,而取复数的共轭采用函数conj(·)。这样,只需简单的6句代码就可以实现式(5)~(10) 相应的潮流计算,代码量非常少,且相对单条支路功率前推和电压回代计算的循环实现方式,速度将会大幅度提高,且规模越大,速度提高的幅度越大。
5 算例
为了对比本文的分层前推回代法与文[3]的前推回代法,在MATLAB环境下进行了相应算法的程序编制,并分别以IEEE40节点树状网和一个实际的1338节点城市配电网络作为算例进行了计算,两种算法的计算结果完全一样,但分层前推回代法计算时间分别为0.03s和0.75s,文[3]的前推回代法计算时间分别为0.12s和91s。这显示出分层前推回代法在计算速度上的明显优势,并且网络规模越大,优势越显著。这是由于随着网络规模的增大,在供电半径的限制下网络层次不可能增大很多,因此相比较而言分层的效果更显著,例如IEEE40节点网络的39条支路分为10层,平均每层只有4条支路,而1338节点网络的1337条支路共分为38层,平均每层35条支路,最多的一层上有113条支路。
6 结论
利用辐射型网络同一层次之间的支路功率前推和电压回代相互独立的特点,本文提出了一种新颖的分层前推回代算法。该算法将网络支路按层次进行分类,并分层并行计算各层次的支路功率损耗和电压损耗,因而可大幅度提高配网潮流的计算速度。本文在MATLAB环境下,利用其快速的复数矩阵运算功能,实现了文中的分层前推回代算法,也取得了非常明显的速度效益。另外,本文还发现并讨论了当变压器支路阻抗过小时,利用∏型模型会产生数值巨大的对地导纳,由此会导致潮流不收敛。为此,本文根据理想变压器对功率和电压的变换原理,提出了一种有效的电压变换模型来处理变压器支路,从而改善了潮流算法的收敛特性。算例结果表明:该算法计算速度快、收敛性好,对于大规模辐射型网络,效果尤其明显。
在公式(18)中,和分别表示状态变量与其修正量组成的列向量;为方阵,一般叫作雅可比矩阵,第i行j列元素为 ,它的大小为第i个函数对第j个变量求偏导;k则表示阵元素都在处取;同时,F(X)是由n个函数组成的n维列向量;在极坐标下,节点电压可如下表示:在这里插入图片描述
(19)
若和为已知大小的功率,与从节点电压求得的有功和无功功率之差,为功率的不平衡量,则节点功率不平衡量可用如下公式计算:
在这里插入图片描述
(20)
节点功率可用各节点电压模值与相位表示,如下公式所示:
在这里插入图片描述
(21)
式(21)中,为节点i和j的相位差。
由以公式(18)-(21)推得牛顿法下,其潮流计算方程可写为:
在这里插入图片描述
(22)
公式(22)中,雅可比矩阵的各元素为
在这里插入图片描述
(23)
(24)
(25)
(26)
在这里插入图片描述
(27)
(28)
(29)
(30)
其中,节点导纳矩阵的元素由Gij 、Bij表示。
随着国内外配电系统自动化水平不断提高,电力行业人员也开始更加深入地研究配电网系统。配电网潮流计算作为DMS(配电管理系统)的重要基础,受到广大行业界人士的关注。因此,配电网潮流计算,已然成为配电网分析的重要内容。配电网与输电网相比,两者有明显不同,前者一般采用网格结构,线路参数R/X的值较大,三相负荷不对称程度明显。这些特点使得在输电网中计算有效,如牛顿法,不再适用于配电网。为此,有学者提出了适用于配电网的潮流算法,主要包括基于回路方程的潮流算法、前推回推法和改进的牛顿-拉夫逊法[17](简称改进的牛拉法)。其中,基于回路方程的方法具有较强的网格处理能力和良好的收敛性,但该方法的节点数和分支数处理非常复杂。前推回推法是针对配电网的树状特性,可以避免潮流计算中的病态条件,同时速度更快。然而,由于其公式和算法与牛顿潮流算法不同,其在其它方面(如潮流优化)的应用将受到限制。
改进牛顿法通过对传统法进行一定的近似,将J阵写成UDUT 的形式。U仅由网络拓扑决定,是一个上三角矩阵;D是一个对角矩阵。在牛拉法中,需要对J阵因子分解与前代回代,改进法则只有前推回代的计算过程。它很好地改善了传统法以及前推回推法。经过算例计算结果证明,改进法可以避免J阵病态,且拥有前推回代法的收敛速度、精度,又由于它属于牛顿型算法,所以该算法已经得到了广泛的运用[18]。
下面附带电力系统分析牛顿法算例及matlab程序:
网络结构如下:系统结构图
系统参数如下:
在上图所示的简单电力系统中,系统中节点1、2为PQ节点,节点3为PV节点,节点4为平衡节点,已给定P1s+jQ1s=-0.30-j0.18 P2s+jQ2s=-0.55-j0.13 P3s=0.5 V3s=1.10 V4s=1.05∠0°
容许误差ε=10-5
节点导纳矩阵:
导纳矩阵
各节点电压:
节点 e f v ζ
1.0.984637 -0.008596 0.984675 -0.500172
2.0.958690 -0.108387 0.964798 -6.450306
3.1.092415 0.128955 1.100000 6.732347
4.1.050000 0.000000 1.050000 0.000000
各节点功率:
节点 P Q
1-0.300000 -0.180000
2–0.550000 -0.130000
3 0.500000 -0.551305
4 0.367883 0.264698
matlab程序如下:
// 牛顿法潮流计算matlab程序
clc
Y=[1.042093-8.242876i -0.588235+2.352941i 3.666667i -0.453858+1.891074i
-0.588235+2.352941i 1.069005-4.727377i 0-0.480769+2.403846i
3.666667i 0 -3.333333i0
-0.453858+1.891074i -0.480769+2.403846i 00.934627-4.261590i]
%导纳矩阵
e=[1 1 1.1 1.05]%初始电压
f=zeros(4,1)
V=zeros(4,1)%节点电压
Ws=[-0.3 -0.18 -0.55 -0.13 0.5 1.1]%初始功率
W=zeros(6,1)
n=length(Y)%节点数
J=zeros(2*(n-1))%雅可比矩阵
delta_v=zeros(1,6)
delta_w=Ws
G=real(Y)
B=imag(Y)
S=zeros(4,2)
c=0%循环次数
m=input('请输入PQ节点数:')
while max(abs(delta_w))>10^-5
for i=1:(n-1)%以下为求取雅可比矩阵
for j=1:(n-1)
if (i~=j)
J(2*i-1,2*j-1)=-(G(i,j)*e(i)+B(i,j)*f(i))
J(2*i,2*j)=-J(2*i-1,2*j-1)
J(2*i-1,2*j)=B(i,j)*e(i)-G(i,j)*f(i)
J(2*i,2*j-1)=J(2*i-1,2*j)
end
end
end
for j=1:(n-2)
J(6,2*j-1)=0
J(6,2*j)=0
end%以上为非对角线元素
s1=0
s2=0
for i=1:(n-1)
for j=1:n
s1=s1+(G(i,j).*e(j)-B(i,j).*f(j))
s2=s2+(G(i,j).*f(j)+B(i,j).*e(j))
end
J(2*i-1,2*i-1)=-s1-G(i,i) *e(i)-B(i,i)*f(i)
J(2*i-1,2*i)=-s2+B(i,i) *e(i)-G(i,i)*f(i)
s1=0
s2=0
end
for i=1:m
for j=1:n
s1=s1+G(i,j).*f(j)+B(i,j).*e(j)
s2=s2+(G(i,j).*e(j)-B(i,j).*f(j))
end
J(2*i,2*i-1)=s1+B(i,i) *e(i)-G(i,i)*f(i)
J(2*i,2*i)=-s2+G(i,i) *e(i)+B(i,i)*f(i)
s1=0
s2=0
end
J(6,5)=-2*e(3)
J(6,6)=-2*f(3)%对角线元素求解
for i=1:m
for j=1:n
s1=s1+e(i)*(G(i,j).*e(j)-B(i,j).*f(j))+f(i)*(G(i,j).*f(j)+B(i,j).*e(j))
s2=s2+f(i)*(G(i,j).*e(j)-B(i,j).*f(j))-e(i)*(G(i,j).*f(j)+B(i,j).*e(j))
end
delta_w(2*i-1)=Ws(2*i-1)-s1
delta_w(2*i)=Ws(2*i)-s2
W(2*i-1)=s1
W(2*i)=s2
s1=0
s2=0
end
for j=1:n
s1=s1+e(3)*(G(3,j).*e(j)-B(3,j).*f(j))+f(3)*(G(3,j).*f(j)+B(3,j).*e(j))
end
delta_w(5)=Ws(5)-s1
delta_w(6)=(Ws(6)^2-(e(3)^2+f(3)^2))
W(5)=s1
W(6)=sqrt(e(3)^2+f(3)^2)%以上求功率差值
delta_v=-inv(J)*delta_w
for i=1:(n-1)
e(i)=e(i)+delta_v(2*i-1)
f(i)=f(i)+delta_v(2*i)
end%求电压差值
c=c+1
end
for x=1:4
V(x)=e(x)+f(x)*1i
end%节点电压
s1=0
for x=3:4
for j=1:4
s1=s1+conj(Y(x,j))*conj(V(j))
end
S(x,1)=real(V(x)*s1)
S(x,2)=imag(V(x)*s1)
s1=0
end%PV与平衡节点功率
for x=1:2
S(x,1)=W(2*x-1)
S(x,2)=W(2*x)
end%节点功率
c
J
V
S
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
运行结果如下:
潮流计算结果
文章知识点与官方知识档案匹配
算法技能树首页概览
34356 人正在系统学习中
电磁MATLAB
微信公众号
程序从vx公众号获取,助力科研!
打开CSDN APP,看更多技术内容
matlab潮流程序,潮流计算的MATLAB源程序_weixin_39760433的博客-CSDN...
1、简单潮流计算的牛顿拉夫逊程序,相关的原始数据数据数据输入格式如下:%B1是支路参数矩阵,第一列和第二列是节点编号。节点编号由小到大编写%对于含有变压器的支路,第一列为低压侧节点编号,第二列为高压侧节点%编号,将变压器的串联阻抗...
继续访问
...编程实现_某只旺仔的博客_电力系统潮流计算程序matlab
潮流计算是电力系统中最基本,应用最广泛的一种计算,是电力系统稳定计算和故障分析的基础。这个代码通过matlab实现了用极坐标表示的牛顿法和P-Q分解法进行潮流计算,使用IEEE14节点系统进行测试,计算结果和应用matpower的潮流计算完全一致。代码...
继续访问
<em>matlab</em>导数计算
<p>matlab导数计算,通过编程实战掌握具体应用。包括matlab导数计算前、matlab导数计算中、matlab导数计算后。</p>
继续访问
(完整版)基于MATLAB牛顿拉夫逊法进行潮流计算.doc
>%本程序的功能是用牛顿拉夫逊法进行潮流计算 n=input' 请输入节点数 :n=nl=input' 请输入支路数 :nl=isb=input' 请输入平衡母线节点号 :isb=pr=input' 请输入误差精度 :pr=B1=input' 请输入由各支路参数形成的矩阵 :B1=B2=input' 请输入各节点参数形成的矩阵 :B2=Y=zeros(ne=zeros
牛顿拉夫逊潮流计算matlab程序
基于牛顿拉夫逊进行潮流计算,求得各节点电压,各支路功率流动,内附程序输入说明以及案例。
珍藏多年的matlab潮流计算程序源代码集合,包含多个潮流计算程序
【达摩老生出品,必属精品,亲测校正,质量保证】 资源名:珍藏多年的matlab潮流计算程序源代码集合,包含多个潮流计算程序 资源类型:matlab项目全套源码 源码说明: 全部项目源码都是经过测试校正后百分百成功运行的,如果您下载后不能运行可联系我进行指导或者更换。 适合人群:新手及有一定经验的开发人员
牛顿-拉夫逊法潮流计算——采用极坐标形式
采用极坐标形式的Newton-Raphson法进行潮流计算,并提供IEEE 14、57节点及New England 39节点的计算数据。使用方法:运行/修改PowerFlow_Newton_Polar.m文件。
电力系统牛拉法潮流计算MATLAB程序
实现潮流计算的MATLAB代码,使用牛顿拉夫逊法,可更改结点数,支路数和导纳阵,通用性强
【课设/毕业设计】电力系统潮流计算(Matlab代码实现)
潮流计算是电力系统分析中的一种最基本的计算,对给定系统进行潮流计算可以得到各母线上的电压、网络中的功率分布及功率损耗等。本文介绍了潮流计算在电力系统分析中的作用及潮流计算的发展状况。通过对潮流计算所用的数学模型进行分析,建立潮流计算的基本方程。牛顿-拉夫逊法是目前广泛应用的一种潮流计算方法,本文阐述了牛顿-拉夫逊潮流计算的基本原理。详细展示了利用matlab平台编写潮流计算程序的具体过程,通过实例证明基于牛顿-拉夫逊法的潮流计算程序具有收敛速度快、占用内存小的优点。潮流计算的程序实现手段。
继续访问
传统牛顿拉夫逊法潮流计算matlab程序.zip
采用何仰赞《电力系统分析》中介绍的牛顿拉夫逊法matlab编程,程序还能实现N-1校核和线路网损分析,适合学习电气工程领域的学子。
基于matpower的电力系统潮流计算matlab程序
它是由美国康奈尔大学电力系统工程研究中心(PSERC of Cornell University)的RAY D. Zimmerman、Carlos E. Murillo-Sánchez和甘德强在Robert J. Thomas的指导下开发出来的,本章介绍的是MATPOWER4.0。每一个电网用变量名为“mpc”的结构体(structures)来定义,结构体mpc的不同字段用baseMVA、bus、branch、gen等来定义和返回电网的具体参数。列的数据类似于标准的IEEE 和PTI 列的数据格式。
继续访问
直流潮流计算matlab程序
直流潮流发的特点是用电力系统的交流潮流(有功功率和无功功率)等值的直流电流来代替。甚至只用直流电路的解析法来分析电力系统的有功潮流,而不考虑无功分布对有功的影响。这样一来计算速度加快,但计算的准确度有所降低,本方法适用于对潮流计算准确度要求不高的计算场景。θ为网络中各节点的电压相位角的向量;P为节点注入的有功功率向量。这就相当于线路两端的直流电位分别为θi和θj。(2)按照标幺值计算时,节点电压与其额定电压相差不大,故有:Ui≈Uj≈1.0;以IEEE9节点系统为算例,系统参数如下。
继续访问
牛顿-拉夫逊法潮流计算matlab程序,牛顿—拉夫逊法潮流计算MATLAB程序.doc
!!!!!!!!!!!!########################牛顿—拉夫逊法潮流计算程序By Yuluo%牛顿--拉夫逊法进行潮流计算n=input('请输入节点数:n=')n1=input('请输入支路数:n1=')isb=input('请输入平衡母线节点号:isb=')pr=input('请输入误差精度:pr=')B1=input('请输入由支路参数形成的矩阵:B1=')B...
继续访问
matlab程序 潮流计算,潮流计算matlab程序
《潮流计算matlab程序》由会员分享,可在线阅读,更多相关《潮流计算matlab程序(3页珍藏版)》请在人人文库网上搜索。1、clear%各节点参数:节点编号,类型,电压幅值,电压相位,注入有功,注入无功 %类型:1PQ节点,2PV节点,3平衡节点%本程序中将最后一个节点设为平衡节点R_1=1 1 1.0 0 0.2 0.2j2 1 1.0 0 -0.45 -0.15j3 1 1.0 0 ...
继续访问
matlab潮流计算函数,基于牛顿拉夫逊法潮流计算的matlab实验报告(含源程序和结果)...
《基于牛顿拉夫逊法潮流计算的matlab实验报告(含源程序和结果)》由会员分享,可在线阅读,更多相关《基于牛顿拉夫逊法潮流计算的matlab实验报告(含源程序和结果)(24页珍藏版)》请在人人文库网上搜索。1、基于牛顿拉夫逊法潮流计算的matlab实验报告一、 实验目的和要求1.学习掌握matlab的基本用法2.应用MATLAB语言编写具有一定通用性的牛顿-拉夫逊法潮流计算程序。要求:(1)潮流计...
继续访问
matlab的做潮流计算,Matlab实现潮流计算程序
关于Matlab趋势计算program.doc文件的实现,爱文共享的信息具有丰富的相关文档,每天都有成千上万的行业名人在该站点上共享最新信息.程序代码如下: 读取数据clcclearfilename ='txt'a = textread(filename)n = a(,)pinghengjd = a(,)phjddianya = a(,)jingdu = a(,)b = zeros (,)j = ...
继续访问
天然气潮流计算matlab程序
此外,由于天然气在管道内传输时会产生压力损耗,所以每隔一段都要设置压缩机来提升节点压力,保证天然气系统的正常运行。(3)如果节点中存在压缩机,则可以求得压缩机消耗的等效天然气流量,将其作为一个负荷加入到潮流计算中;管道或者节点的约束,例如流量、节点压力、压缩机消耗电能和压缩比等,如下面四式所示。(2)根据天然气节点压力与节点间流量的关系式求出前一个节点的节点压力;(1)由用户负荷求得微型燃气轮机的功率