1、电 气 工 程 及 其 自 动 化 专 业 课 程 设 计复 杂 网 络 N-R 法 潮 流 分 析 与 计 算 的 设 计学 生 学 号 ： 学 生 姓 名 ： 班 级 ： 指 导 教 师 ： 起 止 日 期 ： 哈 尔 滨 工 程 大 学 自 动 化 学 院课 程 设 计 报 告 撰 写 内 容一 、 设 计 要 求 （ 宋 体 ， 小 四 号 字 ， 加 黑 ）用 matlab 编 程 ， N_R 法 计 算 潮 流 分 布具 体 要 求 为 ：（ 1） 给 出 程 序 ， 并 给 出 注 释（ 2） 输 出 迭 代 次 数 ， 各 节 点 电 压 ， 各 支 路 电 流（ 3） 在 图
2、 中 标 明 功 率 流 向123 456T 1T 2L 2L 5L 3L 4节 点 数 据 如 下 表 所 示 （ 标 幺 值 ）1 2 3 4 5 6P 3 1.8 0.6 3.5 5Q 1 0.5 0.8 1.3V 1 1.050支 路 及 变 压 器 数 据线 路 T1 T2 L2 L3 L4 L5阻 抗 j0.04 j0.02 0.06+j0.025 0.01+j0.2 0.06+j0.5 0.05+j0.3导 纳 /2 j0.25 j0.25 j0.25 j0.25变 比 1.05:1 1.05:1精 度 要 求 ： 0.0001- III -二 、 设 计 方 案 （ 要 求 给
3、 出 详 细 的 设 计 思 路 及 其 必 要 的 论 证 ）（1.）潮流计算的方法（1）高斯雅克比迭代法（2）高斯-塞得尔法（对初值要求底，迭代次数多）（3）牛顿-拉夫逊法（使用广泛）（4）PQ 快速分解法（提升运算速度）目前广泛应用的潮流计算方法都是基于节点电压法的，以节点导纳矩阵Y作为电力网络的数学模型。节点电压Ui 和节点注入电流Ii 由节点电压方程YV=I （1）根据S=VI(I为I 的共轭)可得非线性的节点方程YV=I=(S/V) （2）在实际的电力系统中，已知的运行条件不是节点的注入电流，而是负荷和发电机的功率，而且这些功率一般不随节点电压的变化而变化。由于各节点注入功率与注入
4、电流的关系为SiPijQi =ViIi，因此可将式（2）改写为Ii=Si/Vi=Pi+jQi/Vi (i= 1, 2,3 n) （3）式中，Pi 和Qi 分别为节点i 向网络注入的有功功率和无功功率，当i 为发电机节点时Pi0；当i 为负荷节点时Pi0；当i 为无源节点Pi0，Qi0；Vi 和Ii 分别为节点电压相量Vi 和节点注入电流相量Ii 的共轭。式（3）亦即潮流计算的基本方程式,它可以在直角坐标也可以在极坐标上建立2n 个实数形式功率方程式。发电机Pi、Qi 为正，负荷Pi、Qi 为负。展开YV=I 为Ii=YijVj=YiiVi+YijVi( i=1 2 3 n) (4)将式（4）代
5、入式（3），得n 维的非线性复数的电压方程组潮流计算的基本方程为(Pi-jQi)/ Vi= YiiVi+YijVi （i=1, 2, n） (5)（2.）变量的分类假设系统中有n 个节点，构成n 个复数方程，2n 个实数方程，变量总数为6n 个。a)不可控变量（2n 个）：负荷消耗的有功功率Li P 和无功功率Li Q .由于该类变量无法控制，取决于用户，而且出现事先没有预计的变动，使系统偏离原始运行状态，因此又称为不可控变量或扰动变量。b)控制变量（2n 个）:发电机发出的有功功率Gi P 和无功功率Gi Q ，因为该类变量可控。也称独立变量。c)状态变量（2n 个）：母线电压或节点电压的幅
6、值大小i V 与相角大小i ，又称依从变量或因变量。并且i V 受Gi P 控制， i 受Gi Q 控制。其中2n 个扰动变量是给定的，2n 个控制变量和2n 个状态变量中给定两个，求另外两个。（3.）变量的约束条件a)扰动变量没有约束条件。b)控制变量约束条件：为满足发电机的技术经济特性指标。c)状态变量的i V 的约束条件：保证良好的电能质量。d) 状态变量的i 的约束条件：保证系统的稳定运行。(4.)系统节点的分类，根据给定的控制变量和状态变量进行分类如下：（1）PQ 节点（即负荷节点）：给定Pi、Qi，求Vi 和i （ i i e , f ）。通常变电所都是这一类型的节点，由于没有发电
7、设备，因而发电功率为零电力系统中的绝大多数节点属于这一节点。其包含变电站节点（即联络节点或浮游节点）。（2）PV 节点（即调节节点、电压控制节点）：给定Pi 和Vi，求Qi 和i （ i i e , f ）。这类节点必须有足够的可调无功容量，用以维持给定的电压幅值。一般时选择有一顶武功储备的发电厂和具有可调无功电源设备的变电所作为PV 节点。在电力系统中，这类节点数很少。（3）平衡节点（即松弛节点、参考节点、基准节点）：给定Vi 和i （ i =0），求Pi 和Qi。（只有一个）有功功率不能给定，这个节点承担了系统的有功功率平衡。同时其电压幅值也是给定的，相位为零。- V -(5. )P-Q
8、分解法是从改进和简化牛顿法潮流程序的基础上提出来的，它的基本思想是：把节点功率表示为电压向量的极坐标方程式，抓住主要矛盾，以有功功率误差作为修正电压向量角度的依据，以无功功率误差作为修正电压幅值的依据，把有功功率和无功功率迭代分开来进行。牛顿法潮流程序的核心是求解修正方程式，当节点功率方程式采取极坐标系统时，修正方程式展开为：P = H + NV/ VQ = J + LV /V以上方程式是从数学上推倒出来的，并没有考虑电力系统这个具体对象的特点。电力系统中有功功率主要与各节点电压向量的角度有关，无功功率则主要受各节点电压幅值的影响。大量运算经验也告诉我们，矩阵N 及J 中各元素的数值相对是很小
9、的，因此对牛顿法的第一步简化就是把有功功率和无功功率分开来进行迭代，即将式(4)化简为：P = HQ = LV /V (5)这样，由于我们把2n 阶的线性方程组变成了二个n 阶的线性方程组，对牛顿法的第二个化简，也是比较关键的一个化简，即把式(5)中的系数矩阵简化为在迭代过程中不变的对称矩阵。众所周知，一般线路两端电压的相角差是不大的(通常不超过1020 度)，因此可以认为：(6)此外，与系统各节点无功功率相应的导纳Li B 必定远远小于该节点自导纳的虚部，即：因此， (7)考虑到以上关系后，式(5)中系数矩阵中的元素表达式可以化简为：（8）这样，式(5)中系数矩阵可以表示为：（9）进一步可以
10、把它们表示为以下矩阵的乘积：（10）将它代入(5)中，并利用乘法结合率，我们可以把修正方程式变为：将以上两式的左右两侧用以下矩阵左乘就可得到- VII -以上两式就是P-Q 分解法达到修正方程式，其中系数矩阵只不过是系统导纳矩阵的虚部，因而是对称矩阵，而且在迭代过程中维持不变。它们与功率误差方程式构成了P-Q 分解法迭代过程中基本计算公式，其迭代步骤大致是：(1)给定各节点电压向量的电压初值V i (0) ,i (0);(2)根据(12)计算各节点有功功率误差Pi ，并求出;Pi/Vi(3)解修正方程式(11)，并进而计算各节点电压向量角度的修正量i (4)修正各节点电压向量角度i ；(5)根
11、据式(16)计算各节点无功功率误差i Q ，并求出 / ; i i Q V(6)解修正方程式(11)，求出各节点电压幅值的修正量i V(7)修正各节点电压幅值i V (i ) (i 1) (i 1)i i i V =V V (18)(8)返回(2)进行迭代，直到各节点功率误差及电压误差都满足收敛条件。P-Q 分解法与牛顿法潮流程序的主要差别表现在它们的修正方程式上。P-Q分解法通过对电力系统具体特点的分析，对牛顿法修正方程式的雅可比矩阵进行了有效的简化和改进，有以下三个特点：（1）在提高计算速度和减少内存方面的作用是明显的，不再叙述。（2）使我们得到以下好处。首先，因为修正方程式的系数矩阵就是
12、导纳矩阵的虚部，因此在迭代过程中不必象牛顿法那样进行形成雅可比矩阵的计算，这样不仅是仅减少了运算量，而且也大大简化了程序。其次，由于系数矩阵在迭代过程中维持不变，因此在求解修正方程式时，可以迅速求得修正量，从而显著提高了迭代速度。（3）可以使我们减少形成因子表时的运算量，而且由于对称矩阵三角分解后，其上三角矩阵和下三角矩阵有非常简单的关系，所以在计算机中可以只储存上三角矩阵或下三角矩阵，从而也进一步节约了内存。三 、 设 计 内 容%本 程 序 的 功 能 是 用 牛 顿 拉 夫 逊 法 进 行 潮 流 计 算% B1 矩 阵 ： 1、 支 路 首 端 号 ； 2、 末 端 号 ； 3、 支
13、路 阻 抗 ； 4、 线 路 对 地 电 纳 (或 变 压 器 导纳 )；% 5、 支 路 的 变 比 ； 6、 支 路 首 端 处 于 K 侧 为 1， 1 侧 为 0；% 7、 线 路 /变 压 器 标 识 （ 0/1） 变 压 器 参 数 当 支 路 首 端 处 于 K 侧 标 识 为 1 时 归 算 至 末 端侧 ， 0 归 算 至 首 端 侧% B2 矩 阵 ： 1、 该 节 点 发 电 机 功 率 ； 2、 该 节 点 负 荷 功 率 ；% 3、 PQ 节 点 电 压 初 始 值 ， 或 平 衡 节 点 及 PV 节 点 电 压 的 给 定 值% 4、 节 点 所 接 无 功 补
14、偿 并 联 电 容 （ 感 ） 的 电 纳% 5、 节 点 分 类 标 号 ： 1 为 平 衡 节 点 （ 应 为 1 号 节 点 ） ； 2 为 PQ 节 点 ； 3 为 PV 节 点 ；clear;isb=1; %input(请 输 入 平 衡 母 线 节 点 号 ： isb=);pr=0.0001; %input(请 输 入 误 差 精 度 ： pr=);%-n=6;%input(请 输 入 节 点 数 ： n=);nl=6;%input(请 输 入 支 路 数 ： nl=);B1=1 2 0+0.04i 0 1.05 1 1;2 3 0.06+0.025i 0+0.5i 1 0 0;2
15、 5 0.01+0.2i 0+0.5i 1 0 0;3 4 0.06+0.50i 0+0.5i 1 0 0; 4 5 0.05+0.3i 0 1 0 0;6 5 0+0.02i 0 1.05 1 1B2=0 0 1 0 1;0 3+1i 1.00 0 2;0 1.8+0.50i 1.00 0 2;0 0.6+0.8i 1.00 0 2;0 3.5+1.3i 1.00 0 2;0 -5+0i 1.05 0 3 - IX -%input(请 输 入 各 节 点 参 数 形 成 的 矩 阵 ： B2=);%X=1 0;2 0;3 0;4 0;5 0;6 0%-%n=4;%input(请 输 入 节
16、点 数 ： n=);nl=4;%input(请 输 入 支 路 数 ： nl=);%B1=1 2 4+16i 0 1 0 0;1 3 4+16i 0 1 0 0;2 3 2+8i 0 1 0 0;2 4 1.49+48.02i 0 11/110 0 1 %input(请 输 入 由 支 路 参 数 形 成 的 矩 阵 ： B1=);%B2=0 0 115 0 1;0 0 110 0 2;0 20+4i 110 0 2;0 10+6i 10 0 2 %input(请 输 入 各 节 点 参 数形 成 的 矩 阵 ： B2=);%-Y=zeros(n);e=zeros(1,n);f=zeros(1
17、,n);V=zeros(1,n);sida=zeros(1,n);S1=zeros(nl);% % %-求 导 纳 矩 阵 -%for i=1:n% if X(i,2)=0;% p=X(i,1);% Y(p,p)=1/X(i,2);%end%endfor i=1:nl %从 1 到 n1（ 总 支 路 数 ）if B1(i,7)=1 %-如 果 是 变 压 器 支 路 -if B1(i,6)=0 %左 节 点 (首 端 )处 于 1 侧p=B1(i,1);q=B1(i,2);else %左 节 点 (首 端 )处 于 K 侧p=B1(i,2);q=B1(i,1);endY(p,q)=Y(p,q
18、)-1./(B1(i,3)*B1(i,5); %非 对 角 元Y(q,p)=Y(p,q); %非 对 角 元Y(q,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)2); %对 角 元 K 侧Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4); %对 角 元 1 侧 +励 磁 导 纳 else %-否 则 为 线 路 支 路 -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(q,q)=Y(q,q)+1./B1(i,3)+B1(i,4)./2.0000; %对 角
19、 元 j 侧 +线 路 电 纳 的 一 半Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)./2.0000; %对 角 元 i 侧 +线 路 电 纳 的 一 半endenddisp(导 纳 矩 阵 Y=);disp(Y);%-给 定 各 节 点 初 始 电 压 及 给 定 各 节 点 注 入 功 率 -G=real(Y);B=imag(Y); %分 解 出 导 纳 阵 的 实 部 和 虚 部 for i=1:n %给 定 各 节 点 初 始 电 压 的 实 部 和 虚 部 e(i)=real(B2(i,3);f(i)=imag(B2(i,3);V(i)=abs(B2(i,3);
20、 %PV、 平 衡 节 点 及 PQ 节 点 电 压 模 值 endfor i=1:n %给 定 各 节 点 注 入 功 率 S(i)=B2(i,1)-B2(i,2); %i 节 点 注 入 功 率 SG-SL B(i,i)=B(i,i)+B2(i,4); %i 节 点 无 功 补 偿 量 (电 纳 值 )end%=用 牛 顿 -拉 夫 逊 法 迭 代 求 解 非 线 性 代 数 方 程 （ 功 率 方 程 ）=P=real(S);Q=imag(S); %分 解 出 各 节 点 注 入 的 有 功 和 无 功 功 率ICT1=0;IT2=1;N0=2*n;N1=N0+1;a=0; %迭 代 次
21、 数 ICT1、 a； 不 满 足 收 敛 要 求 的 节 点 数 IT2while IT2=0 % N0=2*n 雅 可 比 矩 阵 的 阶 数 ； N1=N0+1 扩 展 列IT2=0;a=a+1;JZ=Jacobi 矩 阵 第 (,num2str(a),)次 消 去 运 算 ;JZ1=Jacobi 矩 阵 第 (,num2str(a),)次 回 代 运 算 ;JZ0=功 率 方 程 第 (,num2str(a),)次 差 值 ： ;%-求 取 各 个 节 点 的 功 率 及 功 率 偏 差 及 PV 节 点 的 电 压 偏 差 -for i=1:n %n 个 节 点 2n 行 (每 节
22、点 两 个 方 程 P 和 Q 或 U) p=2*i-1;m=p+1;C(i)=0;D(i)=0;for j1=1:n %第 i 行 共 n 列 (n 个 节 点 间 互 导 纳 及 节 点 电 压 相 乘 即 电 流 )C(i)=C(i)+G(i,j1)*e(j1)-B(i,j1)*f(j1);% (Gij*ej-Bij*fj)D(i)=D(i)+G(i,j1)*f(j1)+B(i,j1)*e(j1);% (Gij*fj+Bij*ej)end%求 i 节 点 有 功 和 无 功 功 率 P,Q的 计 算 值 P1=C(i)*e(i)+f(i)*D(i);%节 点 功 率 P 计 算 ei (Gij*ej-Bij*fj)+fi (Gij*fj+Bij*ej)Q1=C(i)*f(i)-e(i)*D(i);%节 点 功 率 Q 计 算 fi (Gij*ej-Bij*fj)-ei (Gij*fj+Bij*ej)V2=e(i)2+f(i)2; %电 压 模 平 方
