第一篇：电力系统谐波成因分析及谐波潮流计算 外文翻译
中文译文：
电力系统谐波成因分析及谐波潮流计算
首先一个理想的电力系统是由单一恒定频率与规定幅值的稳定电压供电的。但实际上，由于近年来随着科学技术的不断发展，在电力系统中大功率换流设备和调压装置的使用、高压直流输电的应用、大量非线性负荷的出现以及供电系统本身存在的非线性元件等引起系统中的电压波形畸变是越来越严重，对电力系统造成了很大的危害，例如：导致供电系统中的元件损耗增大、降低用电设备的使用寿命、干扰通信系统等。严重时甚至还能使设备损坏，自动控制失灵，继电保护误动作，从而造成停电事故等及其它问题。所“知己知彼，百战不殆”，因此,要实现对电网谐波的综合治理,就必须明确谐波的来源及电网在各种不同运行方式下谐波潮流的分布情况，以采取相应的措施限制和消除谐波，从而改善供电系统供电质量和确保系统的安全经济运行。
其次，电力系统中谐波源是多种多样的。主要有以下几种：
1、系统中的各种非线性用电设备如：换流设备、调压装置、电气化铁道、电弧炉、荧光灯、家用电器以及各种电子节能控制设备等是电力系统谐波的主要来源。即使这些设备提供理想的正弦波电压，它取用的电流也是非线性的，即有谐波电流存在。并且这些设备产生的谐波电流也将注入电力系统，导致系统各处电压产生谐波分量。这些设备的谐波含量决定于它本身的特征和工作条件，基本上与电力系统参数无关，可视为谐波恒流源。
2、供电系统本身存在的非线性元件是谐波的又一来源。这些非线性元件主要有变压器激磁支路、交直流换流站的可控硅控制元件、可控硅控制的电容器、电抗器组等。
3、用户电器设备产生的谐波分量也会对系统产生影响。如荧光灯、家用电器等的单容量不大，但数量很大且散布于各处，电力部门又难以管理的用电设备。如果这些设备的电流谐波含量过大，则会对电力系统造成严重影响，对该类设备的电流谐波含量，在制造时即应限制在一定的数量范围之内。
4、发电机发出的谐波电势。发电机发出额定电势的同时也会有谐波电势产生，其谐波电势取决于发电机本身的结构和工作条件，基本上与外接阻抗无关。故可视为谐波恒压源，但其值很小。
再次，电力系统谐波潮流计算，所谓电力系统谐波潮流计算，就是通过求解网络方程In=YnUn(n=3,5,7…...n：谐波次数。In为谐波源负荷注入电网的n次谐波电流列向量。Yn为电网的n次谐波导纳阵。Un为电网中各节点母线的n次谐波电压列向量)。求得电网中各节点（母线）的谐波电压，进而求得各支路中的谐波电流。当电力系统中存在有谐波源时，此时系统中各接点电压和支路电流均会有高次谐波。为了确定谐波电压和谐波电流在供电系统中的分布，需要对谐波阻抗构成的等效电路进行潮流计算，同时当整流装置供电系统中有容性元件存在时，还要根据各支路谐波阻抗的性质和大小，来检验有无谐振的情况。进行谐波潮流计算，首先必须确定电网元件的谐波阻抗。
3．1 电网各类元件的谐波阻抗：
1.同步发电机的谐波阻抗
合格的发电机的电势是纯正弦的，不含有高次谐波，其发电机电势只存在于基波网络。在高次谐波网络里，由于发电机谐波电势很小，此时可视发电机谐波电势为零。故其等值电路为连接机端与中性点的谐波电抗。其中 XGn=nXG1-------------（1）
式中 XG1为基波时发电机的零序、正序或负序电抗，有该次谐波的序特性决定 如果需要计及网络损耗，对于发电机，可将其阻抗角按85度估计，对于输电线，变压器和负荷等元件的等值发电机，可将其阻抗角按75度估计。
2.变压器的谐波阻抗
电力系统谐波的幅值常是随着频率的升高而衰减，故在基波潮流计算尤其是高压电网中，常忽略变压器的激磁支路和匝间电容。在计算谐波电流时，只考虑变压器的漏抗，且认为与谐波次数所认定的频率成正比。在一般情况下，变压器的等值电路就简化为一连接原副边节点其中的谐波电抗为变压器基波漏电抗。
在高次谐波的作用下，绕组内部的集肤效应和临近效应增大，这时变压器的电阻大致与谐波次数的平方成正比，此时的变压器谐波阻抗为： Zn=sqrt(n)RT1+jnXT1（3）其中RT1为基波时变压器的电阻。
对于三相绕组变压器，可采用星型等值电路，其谐波阻抗的计算方法通上。当谐波源注入的高次谐波电流三相不对称时，则要根据变压器的接线方式和各序阻抗计算出三相谐波阻抗。
3.电抗器的谐波阻抗
当只计及电抗器感抗时，对n次谐波频率为：XLn=Nxl*UN/sqrt(3)IN 4.输电线路的谐波阻抗
输电线路是具有均匀分布参数的电路，经过完全换位的输电线路可看作是三相对称的。
在潮流计算中，通常以集中参数的PI型等值电路表示。在计及分布特性的情况下，则： ZLn=Znsh(rnl)Yln/2=(chrnl-1)/(Znshrnl)ZN和RN分别为对于于该次谐波时线路的波阻抗和传播常数。其中 Zn=sqrt(Z0n/Y0n)Rn=sqrt(Z0nYon)Z0N和Y0N 分别为该次谐波时输电线路单位长度的阻抗和导纳
5.负荷的谐波阻抗
在谐波潮流计算时，基波部分可按节点注入功率看待，而在谐波网络中将它看作是恒定阻抗，近似地可认为综合负荷为一等值电动机。其综合负荷的谐波等值阻抗值为：
ZN=SQRT（N）R1+JNX1 其中 R1，X1 为基波等值电动机的负序电阻、电抗、其值可由该节点的基波电压、功率
值经换算求得。
零序电流一般不会进入负荷，因而在零序性的高次谐波网络里，可忽略负荷支路。
当确定了电路中各电气元件的谐波阻抗后，可以构成一个谐波作用的等效电路，以便进行计算，绘制谐波作用下的等效电路时应注意以下几个特点：
（1）、谐波作用的等效电路，均应以整流装置为中心，按照实际接线构成，于是整流装置视为谐波源，而电力系统的发电机不是以能源出现，而是作为谐波源的负载阻抗的一部分。
（2）、电路元件阻抗可以用有名值进行计算，也可以用标幺值进行计算。当采用有名值进行计算时，全部电路应折算到某一基准电压，便于分析和应用。
（3）一般计算中，元件的所有电阻均可忽略，但是当系统某一部分发生或接近并联或串联谐振时，此时的电阻影响却不能忽略。
（4）、在谐波电流近似计算中，所确定的是整流装置侧的总谐波电流，根据谐波作用等效电路，才能确定各支路谐波电流和电压的分布。3.2 谐波潮流计算
3.2.1 无容性元件网络的谐波潮流计算（1）、对称系统的谐波潮流计算
对称系统中三相情况相同，因此可以按一相情况来计算。
当确定了整流装置任一侧总谐波电流后，结合谐波等效电路，就可以确定系统网络中任一支路的谐波电流分布。然后再根据节点谐波电压和节点注入谐波电流的关系I=YU（其中，Y为谐波导纳阵），就可以确定各处的节点谐波电压了。进而可求出潮流功率。其计算步骤如下：
<1>、根据所给运行条件，以通常的潮流计算方法求解基波潮流。<2>、按谐波源工作条件，确定其它有关参数及需要计算的谐波次数。
<3>、计算各元件谐波参数，形成各次谐波网络节点导纳矩阵，并计算相应谐波网的注入电流。
<4>、由式IN=YNUN确定各节点的谐波电压，并计算各支路谐波功率。
其中，应注意有谐波仪测出的谐波注入电流，其相角是相对于基波电流的相角。故求出基波电流后，需将谐波注入电流相角进行修正。同样，系统节点的功率是基波功率与谐波功率之和，故基波注入功率也应进行修正。但线性负荷处的基波注入功率不必修正。
（2）、不对称系统谐波潮流计算
在不对称系统中，三相情况各不相同，而且相互影响，因此必须同时进行三相系统的计算。
不对称网络潮流的计算可将网络分为各次谐波网络，先计算基波网络，求得各节点基波电压后，按它计算各谐波潮流的各次注入电流，再按此谐波注入电流解算各次谐波的网络方程，求出各节点的各次谐波电压。四．总结
电力系统中谐波的出现，对于电力系统运行是一种“污染”。它们减少了系统电压正玄波形的质量，不仅严重地影响了电力系统自身，而且也损害了用户和周围的通信系统。因此对电力系统谐波的研究对于改善电能质量，抑制和消除谐波具有十分重要的意义。
Power system harmonic analysis and harmonic tidal calculation Firstly, the ideal electrical power system is supplied by a regulated voltage source that has a single and constant frequency and specified amplitude.But in fact, with the development of science and technology in recent years, the use of high-power commutation devices and regulating devices, the applications of HVDC，the emergence of a large number of non-linear load as well as the power supply system itself and so on that have cause a distortion which is more and more serious to the Voltage waveform, which has caused very great harm to the electrical power system.For example: it can Causes the Component Losses in the power supply system more serious, reduce the life of electrical equipment, disturb the communication system and so on.It can also damage the equipment, result in automatic control malfunction and cause the misoperation of relay protection.Thus cause the power outage accident and other questions.So-called “Know me and know thee, lose none and win all ”, therefore, if you want to Comprehensive manage the power system harmonic, you must make clear the harmonic sources and distribution grid of harmonic power flow in a variety of different operating modes , then you can take appropriate measures to limit and eliminate harmonics to improve power quality of the power supply system and ensure the safety of system with economic operation.Secondly, the source of electrical power system overtone has a varied forms.It mainly has the following several kinds:
1, system's each kind of nonlinear electric equipment.For example: Converter equipment, regulator device, electrified railway, electric-arc furnace, fluorescence the lamp, the domestic electric appliances and kinds of electronic energy conservation control device and so on are the important source of electrical power system overtone.Even if these equipment supplies the ideal sine wave voltage, the electric current it takes is also non-linear, namely the harmonic current exists.And the harmonic current the equipment produces will also pour into the electrical power system, which will cause existence of harmonic component in the system voltage of each place.The overtone that has nothing to do with the electrical power system parameter in these equipment decided its characteristic and working condition, which maybe regarded as source of the overtone with constant flow.2, the overtone generated by the nonlinear element in the power supply system itself is another
origin.There are many kinds of nonlinear element.For example: transformer excitation branch , Ac/dc converter thyristor control element, thyristor controlled capacitor, reactor, etc.3, The harmonic component produced by user electric equipment will also impact on the system.Such as fluorescent lamps, household appliances with a small capacity has a large quantity and spread in everywhere.So the power sector has difficulty in managing such equipment.If harmonic content current of these devices is too large, it will cause serious impact to power system,and the manufacturing number of the current harmonic content of this kind of equipment should be limited in a certain range.4, Harmonic voltage from generator is another source.Generator is sued to produce the rated voltage but at the same time also engender harmonic voltage which is determined by the structure and the working conditions of the generator.And harmonic voltage has nothing to do with external impedance that can be treated as the harmonic constant voltage source, but its value is very small.Thirdly, The calculation of power system harmonic power flow is completed by solve the network equations In = YnUn(n = 3, 5, 7...n: number of harmonics.In: Harmonic current’s n times column injected into power grid by load that is considered as the harmonic source.Yn: n times harmonic admittance matrix of power grid.Un: n times harmonic voltage of each bus nodes in the grid column).We can get the power grid harmonic voltage of each node(bus)then the harmonic current of each current branch.When there are harmonic sources exist in the power system, the system of contact voltages and branch currents will have higher harmonic.In order to determine the distribution of harmonic voltage and harmonic current in power supply system , we need do calculation of power flow for equivalent circuit consisted by harmonic impedance.At the same time, when capacitive component exist in rectifying device in the power supply system, also according to the nature and size of each branch of harmonic impedance to verify presence of resonance condition.To do harmonic power flow calculation, we must determine the harmonic impedance of power grid element first.3.1 Power grid harmonic impedance of all kinds of components: 1.Harmonic impedance of synchronous generator
Qualified potential is pure sine generator, do not contain high order harmonic, the generator voltage
is only the base network.In higher harmonic network, due to the harmonic voltage generator is very small, we can regard harmonic voltage generator as zero at this time.Therefore, its equivalent circuit is harmonic reactance that connected to the machine side and neutral.The XGn = nXG1--------------------------(1)
Type of XG1 is fundamental wave generator, positive sequence and negative sequence based on the zero sequence reactance decided by the characteristics of harmonic sequence.if needs to take lose of the network into consideration, estimates that the generator’s impedance angle according to 85 degrees, regarding the transmission line, part's and so on transformer and load equivalent generators, may its impedance angle according to 75 degrees.2.Harmonic impedance of transformer
Power system harmonic with the increase of frequency and amplitude attenuation is often change, so often ignore the excitation transformer branch and turn-to-turn capacitive when calculating the base wave method especially in high voltage grid.When calculating the harmonic current, only consider the transformer leakage reactance, and think that is proportional to the harmonic frequency determined by the frequency.In the general case, the transformer equivalent circuit is simplified as the former vice edge connection node and one of the fundamental harmonic reactance is transformer leakage reactance.Under the action of high order harmonic, the skin effect and the adjacent effect inside the winding increases, then the of the transformer’s resistance is proportional to the square of the number of times for the harmonic, the harmonic impedance for the transformer: Zn= SQRT(n)+ jnXT1 RT1----------------------(3)The RT1 is fundamental wave resistance of a transformer.For three-phase winding transformer, the radial equivalent circuit can be used, the calculation method of the harmonic impedance is same above.When the three-phase of higher harmonic current injected in harmonic sources are asymmetric, three-phase harmonic impedance can be calculated according to the connection mode of the transformer and each sequence impedance.3.Harmonic impedance of reactor.When only consider reactor inductance, the NTH harmonic frequency is:
XLn = Nxl * UN/SQRT(3)IN
4.Harmonic impedance of the transmission line
The distributed parameter of transmission lines are uniform, the transmission line is considered to be a three-phase symmetrical after complete transposition.In the calculation of power flow, it usually be expressed as PI type equivalent circuit with a lumped parameter.Under the consideration of the meter and the distribution characteristics, then: ZLn=Znsh(rnl)Yln/2=(chrnl-1)/(Znshrnl)Zn and RN: the line impedance and the harmonic propagation constant.The Zn = SQRT(Z0n/Y0n)
Rn = SQRT(Z0nYon)Z0N and Y0N: the harmonic impedance and admittance of per unit length of transmission lines 5.Harmonic impedance of load
In the calculation of harmonic power flow, the fundamental part can be regarded as the power injected by node, but it often be seen as a constant impedance in harmonic network, approximately integrated load can be thought of as an equivalent motor.Equivalent impedance values of the comprehensive load harmonic: ZN=SQRT（N）R1+JNX1 R1, the X1 : the fundamental negative sequence resistance of equivalent motor, reactance, its value can be calculated by base wave’s voltage and power of the node.Not to load in the zero sequence current, so negligible load branch can be ignored in the higher harmonic zero sequence sex network because there is no zero-sequence current in load Under normal circumstances.When the harmonic impedance of each electric component are determined in the circuit, a harmonic equivalent circuit can be formed and calculated, in order to draw equivalent circuit under the action of harmonic , should pay attention to the following features:
1.Equivalent circuit under the action of harmonic and shall be centered on rectifying device, according to the actual wiring, so regard the rectifier device as a source of harmonic, and power system energy generator is not to appear, but as part of the harmonic source load impedance.2.Circuit components impedance can be calculated by using a name value and also using the MAO values to calculate.When using a name values to calculate, all should be converted to a voltage reference circuits, in order to facilitate analysis and application.3.The resistance of all component can be ignored in general calculation, but when one part of the system or nearly parallel or series resonance occurs, the influence of the resistance cannot be ignored.4.Harmonic current is determined by total harmonic current of the rectifying device side in the approximate calculation, the selection of harmonic current and voltage distribution can be determined according to equivalent circuit of the harmonic.3.2 The calculation of harmonic power flow
3.2.1 The calculation of harmonic power flow without capacitive element network(1).Calculation of harmonic power flow in the symmetric system Three-phase system is the same in the Symmetric system, so you can make a calculation according to consideration of one phase.When the total harmonic current of the rectifying device either side is determined, the harmonic current distribution of the system at any branch in the network can be determined according to harmonic equivalent circuit.Then according to the node harmonic voltage and harmonic current injected by the node I = YU(among them, Y : harmonic admittance matrix), the harmonic voltage of various nodes can be determined.And then the power of the tidal current can be worked out.The calculation steps are as follows: < 1 >.The fundamental wave tide can be calculated by using the usual flow method according to the given operating conditions.< 2 >.According to the working conditions of harmonic source, other relevant parameters and harmonic frequency will be computed.< 3 >.Calculate the component of harmonic parameters, form every harmonic network node admittance matrix, and calculate the current injected in the corresponding harmonic net.< 4 >.Harmonic voltage of each node is determined by the type = YNUN IN and harmonic power of each branch is calculated.11
Among them, attention should be paid to harmonic current measured by the harmonic meter injected in current, the phase Angle is observed in current phase Angle.So after the calculation of the fundamental wave current, the phase Angle of the injected harmonic current must be amended.Also, the system power of a node is the sum of power of fundamental wave and harmonic power, so the injected fundamental wave power should also be revised.But fundamental wave power that injected of linear load don't have to be corrected.(2),Calculation of asymmetric system harmonic power flow
In the asymmetric three-phase system, the t conditions of each phase are not identical, and influence each other, so the calculation of three-phase system must be conducted at the same time.Calculation of the tide of asymmetric network can be conducted by dividing the network into every harmonic network, computing the fundamental network first, after get the voltage of base wave at each node, harmonic wave of all injected current is calculated according to it, then calculate the network equations of every harmonic according to harmonic current injected, work out the every harmonic voltage of each node.4.Summary The presence of harmonics in power system is a kind of “contamination” in the power system operation.They reduce the quality of the waveform of system voltage, not only affect the power system itself seriously, but also harm the user and the surrounding communication system.So the study of power system harmonic is of great significance in improving the power quality, restraining and eliminating harmonics.
第二篇：电力系统潮流计算
南 京 理 工 大 学
《电力系统稳态分析》
课程报告
姓名
XX
学 号： 5*** 自动化学院 电气工程
基于牛顿-拉夫逊法的潮流计算例题编程报学院(系): 专
业: 题
目: 任课教师 硕士导师 告
杨伟 XX
2015年6月10号
基于牛顿-拉夫逊法的潮流计算例题编程报告
摘要：电力系统潮流计算的目的在于：确定电力系统的运行方式、检查系统中各元件是否过压或者过载、为电力系统继电保护的整定提供依据、为电力系统的稳定计算提供初值、为电力系统规划和经济运行提供分析的基础。潮流计算的计算机算法包含高斯—赛德尔迭代法、牛顿-拉夫逊法和P—Q分解法等，其中牛拉法计算原理较简单、计算过程也不复杂，而且由于人们引入泰勒级数和非线性代数方程等在算法里从而进一步提高了算法的收敛性和计算速度。同时基于MATLAB的计算机算法以双精度类型进行数据的存储和运算, 数据精确度高，能进行潮流计算中的各种矩阵运算，使得传统潮流计算方法更加优化。
一 研究内容
通过一道例题来认真分析牛顿-拉夫逊法的原理和方法（采用极坐标形式的牛拉法），同时掌握潮流计算计算机算法的相关知识，能看懂并初步使用MATLAB软件进行编程，培养自己电力系统潮流计算机算法编程能力。
例题如下：用牛顿-拉夫逊法计算下图所示系统的潮流分布，其中系统中5为平衡节点，节点5电压保持U=1.05为定值，其他四个节点分别为PQ节点，给定的注入功率如图所示。计算精度要求各节点电压修正量不大于10-6。
二 牛顿-拉夫逊法潮流计算 1 基本原理
牛顿法是取近似解x(k)之后，在这个基础上，找到比x(k)更接近的方程的根，一步步地迭代，找到尽可能接近方程根的近似根。牛顿迭代法其最大优点是在方程f(x)=0的单根附近时误差将呈平方减少，而且该法还可以用来求方程的重根、复根。电力系统潮流计算，一般来说，各个母线所供负荷的功率是已知的，各个节点的电压是未知的（平衡节点外）可以根据网络结构形成节点导纳矩阵，然后由节点导纳矩阵列写功率方程，由于功率方程里功率是已知的，电压的幅值和相角是未知的，这样潮流计算的问题就转化为求解非线性方程组的问题了。为了便于用迭代法解方程组，需要将上述功率方程改写成功率平衡方程，并对功率平衡方程求偏导，得出对应的雅可比矩阵，给未知节点赋电压初值，将初值带入功率平衡方程，得到功率不平衡量，这样由功率不平衡量、雅可比矩阵、节点电压不平衡量（未知的）构成了误差方程，解误差方程，得到节点电压不平衡量，节点电压加上节点电压不平衡量构成节点电压新的初值，将新的初值带入原来的功率平衡方程，并重新形成雅可比矩阵，然后计算新的电压不平衡量，这样不断迭代，不断修正，一般迭代三到五次就能收敛。2 基本步骤和设计流程图
形成了雅克比矩阵并建立了修正方程式，运用牛顿-拉夫逊法计算潮流的核心问题已经解决，已有可能列出基本计算步骤并编制流程图。由课本总结基本步骤如下：
1)形成节点导纳矩阵Y；
2)设各节点电压的初值，如果是直角坐标的话设电压的实部e和虚部f；如果是极坐标的话则设电压的幅值U和相角a；
3)将各个节点电压的初值代入公式求修正方程中的不平衡量以及修正方程的系数矩阵的雅克比矩阵；
4)解修正方程式，求各节点电压的变化量，即修正量； 5)计算各个节点电压的新值，即修正后的值；
6)利用新值从第（3）步开始进入下一次迭代，直至达到精度退出循环； 7)计算平衡节点的功率和线路功率，输出最后计算结果； ① 公式推导
② 流程图
三
matlab编程代码
clear;
% 如图所示1，2，3，4为PQ节点,5为平衡节点
y=0;
% 输入原始数据，求节点导纳矩阵
y(1,2)=1/(0.07+0.21j);
y(4,5)=0;y(1,3)=1/(0.06+0.18j);
y(1,4)=1/(0.05+0.10j);
y(1,5)=1/(0.04+0.12j);
y(2,3)=1/(0.05+0.10j);
y(2,5)=1/(0.08+0.24j);
y(3,4)=1/(0.06+0.18j);
for i=1:5
for j=i:5
y(j,i)=y(i,j);
end
end
Y=0;
% 求节点导纳矩阵中互导纳
for i=1:5
for j=1:5
if i~=j
Y(i,j)=-y(i,j);
end
end
end
% 求节点导纳矩阵中自导纳
for i=1:5
Y(i,i)=sum(y(i,:));
end
Y
% Y为导纳矩阵
G=real(Y);
B=imag(Y);% 输入原始节点的给定注入功率
S(1)=0.3+0.3j;
S(2)=-0.5-0.15j;
S(3)=-0.6-0.25j;
S(4)=-0.7-0.2j;
S(5)=0;
P=real(S);
Q=imag(S);
% 赋初值，U为节点电压的幅值，a为节点电压的相位角
U=ones(1,5);
U(5)=1.05;
a=zeros(1,5);
x1=ones(8,1);
x2=ones(8,1);
k=0;
while max(x2)>1e-6
for i=1:4
for j=1:4
H(i,j)=0;
N(i,j)=0;
M(i,j)=0;
L(i,j)=0;
oP(i)=0;
oQ(i)=0;
end
end
% 求有功、无功功率不平衡量
for i=1:4
for j=1:5
oP(i)=oP(i)-U(i)*U(j)*(G(i,j)*cos(a(i)-a(j))+B(i,j)*sin(a(i)-a(j)));
oQ(i)=oQ(i)-U(i)*U(j)*(G(i,j)*sin(a(i)-a(j))-B(i,j)*cos(a(i)-a(j)));
end
oP(i)=oP(i)+P(i);
oQ(i)=oQ(i)+Q(i);
end
x2=[oP,oQ]';
% x2为不平衡量列向量
% 求雅克比矩阵
% 当i~=j时，求H,N,M,L
for i=1:4
for j=1:4
if i~=j
H(i,j)=-U(i)*U(j)*(G(i,j)*sin(a(i)-a(j))-B(i,j)*cos(a(i)-a(j)));
N(i,j)=-U(i)*U(j)*(G(i,j)*cos(a(i)-a(j))+B(i,j)*sin(a(i)-a(j)));
L(i,j)=H(i,j);
M(i,j)=-N(i,j);
end
end
end
% 当i=j时，求H,N,M,L
for i=1:4
for j=1:5
if i~=j H(i,i)=H(i,i)+U(i)*U(j)*(G(i,j)*sin(a(i)-a(j))-B(i,j)*cos(a(i)-a(j)));N(i,i)=N(i,i)-U(i)*U(j)*(G(i,j)*cos(a(i)-a(j))+B(i,j)*sin(a(i)-a(j)));
M(i,i)=M(i,i)-U(i)*U(j)*(G(i,j)*cos(a(i)-a(j))+B(i,j)*sin(a(i)-a(j)));
L(i,i)=L(i,i)-U(i)*U(j)*(G(i,j)*sin(a(i)-a(j))-B(i,j)*cos(a(i)-a(j)))
end
end
N(i,i)=N(i,i)-2*(U(i))^2*G(i,i);
L(i,i)=L(i,i)+2*(U(i))^2*B(i,i);
end
J=[H,N;M,L]
% J为雅克比矩阵
x1=-((inv(J))*x2);
% x1为所求△x的列向量
% 求节点电压新值，准备下一次迭代
for i=1:4
oa(i)=x1(i);
oU(i)=x1(i+4)*U(i);
end
for i=1:4
a(i)=a(i)+oa(i);
U(i)=U(i)+oU(i);
end
k=k+1;
end
k,U,a
% 求节点注入功率
i=5;
for j=1:5
P(i)=U(i)*U(j)*(G(i,j)*cos(a(i)-a(j))+B(i,j)*sin(a(i)-a(j)))+P(i);
Q(i)=U(i)*U(j)*(G(i,j)*sin(a(i)-a(j))-B(i,j)*cos(a(i)-a(j)))+Q(i);
end
S(5)=P(5)+Q(5)*sqrt(-1);
S
% 求节点注入电流
I=Y*U'
四
运行结果
节点导纳矩阵
经过五次迭代后的雅克比矩阵
迭代次数以及节点电压的幅值和相角（弧度数）
节点注入功率和电流
五 结果分析
在这次学习和实际操作过程里：首先，对电力系统分析中潮流计算的部分特别是潮流计算的计算机算法中的牛顿-拉夫逊法进行深入的研读，弄明白了其原理、计算过程、公式推导以及设计流程。牛顿-拉夫逊法是求解非线性方程的迭代过程，其计算公式为FJX，式中J为所求函数的雅可比矩阵；X为需要求的修正值；F为不平衡的列向量。利用x(*)=x(k+1)+X(k+1)进行多次迭代，通过迭代判据得到所需要的精度值即准确值x(*)。六 结论
通过这个任务，自己在matlab编程，潮流计算，word文档的编辑功能等方面均有提高，但也暴漏出一些问题：理论知识储备不足，对matlab的性能和特点还不能有一个全面的把握，对word软件也不是很熟练，相信通过以后的学习能弥补这些不足，达到一个新的层次。
第三篇：SD 126-1984 电力系统谐波管理暂行规定
电力系统谐波管理暂行规定
SD126～84 第一章 总则
第一条 电力系统中的谐波主要是治金、化工、电气化铁路等换流设备及其他非线性用电设备产生的。随着硅整流及可控硅换流设备的广泛使用和各种非线性负荷的增加，大量的谐波电流注入电网，造成电压正弦波形畸变，使电能质量下降，给发供电设备及用户用电设备带来严重危害。为保证向国民经济各部门提供质量合格的50赫兹电能，必须对各种非线性用电设备注入电网的谐波电流加以限制，以保证电网和用户用电设备的安全经济运行，特制订本规定。
第二条 本规定适于电力系统以及由电网供电的所有电力用户。
第三条 电网原有的谐波超过本规定的电压正弦波形畸变率极限值时，应查明谐波源并采取措施，把电压正弦波形畸变率限制在规定的极限值以内。在本规定颁发前已接入电网的非线性用电设备注入电网的谐波电流超过本规定的谐波电流允许值时，应制订改造计划并限期把谐波电流限制在允许范围以内。所需投资和设备由非线性用电设备的所属单位负责。
第四条 新建或扩建的非线性用电设备接入电网，必须按本规定执行。如用户的非线性用电设备接入电网，增加或改变了电网的谐波值及其分布，特别是使与电网连接点的谐波电压、电流升高，用户必须采取措施，把谐波电流限制在允许的范围内，方能接入电网运行。
第五条 进口设备和技术合作项目亦应执行本规定。但如对方的国家标准或企业标准的全部或部分规定比本规定严格，则应按对方较严格的规定执行。
第六条 谐波对通讯等的影响应按国内有关规定执行。
第七条 用户用电设备对谐波电压的要求较本规定的电压正弦波畴变率极限更严格时，由用户自行采取限制谐波电压的措施。
第二章 电压正弦波形畸变率极限值和谐波电流允许值
第八条 电网中任何一点的电压正弦波形畴变率均不得超过表1规定的极限值。
表1 电网电压正弦畸形畸变率极限值(相电压)如110千伏电网电压正弦波形畸变率已接近或超过1.5%，但经过测量和计算，下一级电网的电压正弦形畸变率未超过表1规定的极限值时，110千伏电网电压正弦波形畸变率可限制在2%以内。
第九条 限制用户非线性用电设备注入电网的谐波电流是控制电网电压正弦波形畴变的关键。任一用户注入电网连接点的各次谐波电流均不得超过表2规定的允许值。
第十条 用户与电网连接点原有的总电压正弦波形畸变率未超过表1规定的极限值的75%，且新增非线性用电设备向电网注入的谐波电流不超过表2规定的允许值，则允许该用电设备接入电网。
第十一条 用户与电网连接点原有的总电压正弦波形畸变率已超过表1规定的极限值的75%，或新增非线性用电设备向电网注入致谐波电流超过表2规定的允许值时，应由电力部门经过核算，确定该用电设备接入电网的技术措施和供电方案。
第十二条 单台换流设备或交流调压装置的容量不超过表3的规定值时，供电部门可不进行谐波核算即允许接入电网。但是，在公共连接点接入多个换流设备或交流调压装置，总容量超过表3的规定值时，应按第十条、第十一条和附录中(二)的有关规定办理。
第十三条 电流冲击持续时间不超过2秒钟，并且两次冲击之间的间隙时间不小于30秒钟，这种短时间的冲击电注所包含的谐波分量称为短时间谐波电流(或暂态谐波电流)。注入电网的短时间谐波电流，不属于本规定第九条和第十二条的限制范围。
表3 单台三相换流设备和交流调压装置接入电网的允许容量
注：因三相三脉冲换流器向低压电网注入直流，因此必须经过适当的隔离变压器才允许接入0.38千伏电网。第十四条 单相或三相不对称非线性用电设备接入电网，按第十条和第十一条的规定执行。但考核其注入电网的谐波值，应以谐波电流最大的一相作为依据。
用户注入电网的谐波电流允许值 表2 见附表中(2)之3。第十五条 对由110千伏及以上电压供电的用户在接入电网前，应计算非线性用电设备注入的谐波电流在电网中的分布和电压正弦波形畸变率，以便采取措施使电网各部分电压正弦波形畸变率均不超过表1规定的极限值。还应检验并采取措施防止发生谐振和谐波放大。
第十六条 计算电网各部分电压正弦波形畸变率时，应采用电网最小运行方式时该用户接入点的最小短路容量。
第三章 谐波管理
第十七条 各级电力部门对电网的谐波情况，应定期进行测量分析。当发现电网电压正弦波形畸变率超过表1的规定时，应查明谐波源并按第三条的规定，协助非线性用电设备所属单位采取措施，把注入电网的谐波电流限制在表2规定的允许值以下。在新的非线性用电设备接入电网前后，均要进行现场测量，检查谐波电流、电压正弦波形畸变率是否符合本规定。
第十八条 电力部门和用户均应校核接入电网的电力电容组是否会发生有害的并联谐振、串联谐振和谐波放大，防止电力设备因谐波过电流或过电压而损坏。为此，电力部门和用户所安装的电力电容器组，应根据实际存在的谐波情况，采取加装串联电抗器等措施，保证电力设备安全运行。
第十九条 应概括据谐波源的分布，在电网中谐波量较高的地点逐步设置谐波监测点。在该点测量谐波电压，并在向用户供电的线路的送电端测量谐波电流。
测量或计算谐波的次数应不少于19次。即需测量或计算2、3、4、5、6、7、8、9、10、11、12、13、14、15、16、17、18、19次谐波电压和谐波电流。
第二十条 在正常情况下，谐波测量应选拔在电网最小运行方式和非线性用电设备的运行周期中谐波发生量最大的时间内进行。谐波电压和电流应选取5次测量接近数值的算术平均值。
第二十一条 电力系统的运行方式和谐波值都是经常变化的，当谐波量已接近最大允许值时，应加强对电网发供电设备运行工况的监视，避免电器设备受谐波的影响而发生故障。在电网谐波量较高的地点，要逐步安装谐波警报指示器，以便进一步分析谐波情况，并采取措施，保证电力设备安全运行。
第四篇：电力系统潮流计算程序设计
电力系统潮流计算程序设计
姓名：韦应顺
学号：2011021052 电力工程学院
牛顿—拉夫逊潮流计算方法具有能够将非线性方程线性化的特点，而使用MATLAB语言是由于MATLAB语言的数学逻辑强，易编译。
【】【】1.MATLAB程序12
Function tisco %这是一个电力系统潮流计算的程序 n=input（‘n请输入节点数：n=’）； m=input(‘请输入支路数：m=’);ph=input(‘n请输入平衡母线的节点号：ph=’)； B1=input（‘n请输入支路信号：B1=’）;%它以矩阵形式存贮支路的情况，每行存贮一条支路 %第一列存贮支路的一个端点 %第二列存贮支路的另一个端点 %第三列存贮支路阻抗
%第四列存贮支路的对地导纳
%第五列存贮变压器的变比，注意支路为1 %第六列存贮支路的序号
B2=input(‘n请输入节点信息：B2=’)； %第一列为电源侧的功率 %第二列为负荷侧的功率 %第三列为该点的电压值
%第四列为该点的类型：1为PQ，2为PV节点，3为平衡节点 A=input（‘n请输入节点号及对地阻抗：A=’）； ip=input（‘n请输入修正值：ip=’）； %ip为修正值);Y=zeros（n）；
Y(p,q)=Y(p,q)-1./(B1(i3)*B1(i5);e=zeros（1，n）；
Y(p,q)=Y(p,q);f=zeros（1，n）；
no=2*ph=1； Y(q,q)=Y(q,q)+1./B1(i3)+B1(i4)/2；
End for i=1：n
G=real(Y);if A（i2）=0
B=imag(Y);p=A(i1)；
Y(p p)=1./A(i2);for i=1：n End e(i)=real(B2(i3));End f(i)=imag(B2(i3));For i=1:m S(i)=B2(i1)-B2(i2);p=B1(i1);V(i)=B2(i3);p=B1(i2);end Y(p,p)=Y(p,p)+1./(B1(i3)*B1(i5)^2+B1(i4)./2P=real(S);Q=imag(S);[C,D,DF]=xxf(G,B,e,f,P,Q,n,B2,ph,V,no);J=jacci(Y,G,B,P,Q,e,f,V,C,D,B2,n,ph,no);[De,Di]=hxf(J,D,F,ph,n,no);t=0;while
max(abs(De))>ip&max(abs(Dfi)>ip
t=t+1;
e=e+De；
f=f+Df;
[C,D,DF]=xxf(G,B,e,f,P,Q,n,B2,ph,V,no)；
J=jacci(Y,G,B,P,Q,e,f,V,C,D,B2,n,ph,no)；
[De,Df]=hxf(J,Df,ph,n,no);end v=e+f*j;for i=1:n hh(i)=conj(Y(ph,i)*v(i));end S(ph)=sum(hh)*v(ph);B2(ph,1)=S(ph);V=abs(v);
jd=angle(v)*180/p;resulte1=[A(:,1),real(v),imag(v),V,jd,real(S’)，imag(S’),real(B2(:1)),imag(B2(:1)),real(B2(:2)),imag(B2(:,2))];for i=1:m
a(i)=conj((v(B1(i1))/B1(i5)-v(B1(i2))/B1(i3));
b(i)=v(B1(i1))*a(i)-j*B1(i4)*v(B1(i))^2/2;
c(i)=-v(B1(i2))*a(i)-j*B1(i4)*v(B1(i2))^2/2;end result2=[B1(:,6),B1(:,1),B1(:,2),real(b’),imag(b’),real(c’),imag(c’), real(b’+c’),imag(b’+c’)];printcut(result1,S,b,c,result2);type resultm function [C,D,Df]=xxf(G,B,e,f,P,Q,n,B2,ph,V,no)%该子程序是用来求取Df for i=1:n
If
i=ph
C(i)=0;
D(i)=0;
For j=i:n
C(i)=C(i)+G(i,j)*e(j)-B(i,j)*f(j);D(i)=D(i)+G(i,j)*f(j)+B(i,j)*e(j);end
P1=C(i)*e(i)+D(i)*f(i);Q1=C(i)*f(i)-D(i)*e(i);V1=e(i)^2+f(i)^2;If
B2(i4)=2 p=2*i-1;
Df(p)=P(i)-P1;p=p+1;else p=2*i-1;
Df(p)=P(i)-P1;p=p+1;
Df(p)=Q(i)-Q1;end end end Df=Df’;If ph=n Df(no=[];end
function [De,Df]=hxf(J,Df,ph,n,no)%该子函数是为求取De Df DX=JDf;DX1=DX;
x1=length(DX1);if ph=n DX(no)=0;DX(no+1)=0;
For i=(no+2):(x1+2)DX(i)=DX1(i-2);End Else
DX=[DX1,0,0];End k=0;
[x,y]=size(DX);For i=1:2:x K=k+1;
Df(k)=DX(i);De(k)=DX(i+1);End End case 2 Function for j=1:n J=jacci(Y,G,B,PQ,e,f,V,C,D,B2,n,ph,no)X1=G(i,j)*f(i)-B(i,j)*e(i);
X2=G(i,j)*e(i)+B(i,j)*f(i);%该子程序是用来求取jacci矩阵
for i=1:n X3=0;switch B2(i4)X4=0;case 3 P=2*i-1;continue q=2*j-1;case 1 J(p,q)=X1;for j=1:n m=p+1;if
J=&J=ph J(m,q)=X3;X1=G(i)*f(i)-B(i,j)*e(i);q=q+1;X2=G(i,j)*e(i)+B(i,j)*f(i);J(p,q)=X2;X3=-X2;J(m,q)=X4;X4=X1;X1=D(i)+G(i,j)*f(i)-B(i,j)*e(i);p=2*i-1;X2=C(i)+G(i,j)*e(i)+B(i,j)*f(i);q=2*j-1;X3=0;J(p,q)=X1;X4=0;m=p+1;P=2*i-1;J(p,q)=X2;q=2*j-1;J(m,q)=X4;J(p,q)=X1;Else if j=&j=jph m=p+1;X1=D(i)+G(i,j)*f(i)-B(i,j)*e(i);J(m,q)=X3;X2=C(i)+G(i,j)*e(i)+B(i,j)*f(i);q=q+1;X3= C(i)+G(i,j)*e(i)-B(i,j)*f(i);J(p,q)=X2;X4= C(i)+G(i,j)*f(i)-B(i,j)*e(i);J(m,q)=X4;P=2*i-1;end q=2*j-1;end J(p,q)=X1;end m=p+1;end J(m,q)=X3;if ph=n q=q+1;J(no:)=[];J(p,q)=X2;J(no:)=[];J(m,q)=X4;J(:,no)=[];End J(:,no)=[];End
2实例验证 【例题】设有一系统网络结线见图1，各支路阻抗和各节点功率均已以标幺值标示于图1中，其中节点2连接的是发电厂，设节点1电压保持U1＝1.06定值，试计算其中的潮流分布，请输入节点数：n=5 请输入支路数：m=7 请输入平衡母线的节点号：ph=l 请输入支路信息：
BI=[ l 2 0．02+0．06i O l 1；1 3 0．08+0．24i 0 1 2；2 3 0．06+0．18i 0 l 3： 2 4 0．06+0．18i O l 4： 2 5 0．04+0．12i 0 l 5： 3 4 0．01+0．03i 0 l 6： 4 5 0．08+0．24i O 1 7] 请输入节点信息：
B2=[ 0 0 1．06 3；0．2+0．20i 0 1 1；一O．45一O．15i 0 l l；一0．4-0．05i 0 l 1；一0．6—0．1i 0 1 l] 请输入节点号及对地阻抗： A=[l 0；2 0；3 0；4 0；5 O ] 请输入修正值：ip=0．000 0l
参考文献
[1]陈珩．电力系统稳定分析[M]．北京：中国电力出版社，2002：139—187．
[2]郑阿奇．MATLAB实用教程[M]．北京：电子工业出版社，2005：1-243．
[3] 束洪春，孙士云，等．云电送粤交商流混联系统全过 程动态电压研究[J】．中国电力，2008，4l(10)：l-4． SHU Hong—ch吼，SUN Shi-yun，et a1．Research on fun prc'cess dyn锄ic Voltage stabil时of hybrid AC／DC poWer tmnsmission System舶m Yu衄an proVince to G啪gdong province【J】．Electric Power，2008，4l(10)： l-4．
[4] 朱新立，汤涌，等．大电网安全分析的全过程动态仿 真技术[J】．电网技术，2008，32(22)：23—28． SONG Xin—Ii，TANG Yof唱，et a1． Full dyn锄ic simulation for the stabilhy a眦lysis of large power system【J】．Power System融IlrIolo影，2008，32(22)： 23．28．
[5]Roytelm锄I，Shallidehpour S M．A comprehcnsivc long teml dynaIIlic simulation for powcr system recoVery【J】． IEEE Transactions 0n Power Systems，1994，9(3)． [6] 石雩梅，汪志宏，等．发电机励磁系统数学模型及参 数对电网动态稳定性分析结果影响的研究[J】．继电 器，2007，35(21)：22-27．
SHI Xue．mei，WANG Zlli-hon舀et a1．Iksearch on the innuence of g锄e翰to璐baScd ∞de诅iled excitation system models柚d parameterS t0 power铲id dyn锄ic stabil时【J】．Relay，2007，35(2 1)：22-27．
[7] 方思立，朱方．快速励磁系统对系统稳定的影响[J】．中 国电机工程学报，1986，6(1)：20．28．
FANG Si．1i，ZHU Fang．The effbct of f弧t．respon∞
excitation system on the stability of power netwofk【J】． Proceedings ofthe CSEE，1986，6(1)：20-28．
[8] 刘取．电力系统稳定性及发电机励磁控制[M】．北京： 中国电力出版社，2007．
LIU Qu．Power system S诅bility锄d generator excitation control【M】．BeUing：ChiIla Electric Powef Press，2007． [9] Dallachy J L，Anderson T．EXperience with rcplacing ro诅ting exciters wim static exciters【J】．1k InStitution of Electrical Engineers，1 996．
[10] 陈利芳，陈天禄．浅谈自并励励磁系统在大容量机组 中的应用【J】．继电器，2007，35(1)：8l培4． CHEN Li-f抽岛CHEN Tian—lIL Application of 辩l仁exci组tion mode in large capacity髫memtor unit【J】． ReIay'2007，35(1)：81-84．
[11] 方思立，刘增煌，孟庆和．大型汽轮发电机自并励励 磁系统的应用条件【J]．中国电力，1994，27(12)：61．63． FANG Si．Ii，LIU Zeng-hu锄g，MENG Qin争hc．m application conditions of large turbine generator self-excitation system【J】．Electric Powef，1994，27(12)： 61．63．
[12]梁小冰，黄方能．利用EMTDC进行长持续时间过程 的仿真研究【J】．电网技术，2002，26(9)：55．57． LIANG Xiao-bing，HUANG Fan争眦ng．How to cany out simulalion of long dul‘ation processes by use of EMTDC【J】．Power System 11echnology，2002，26(9)： 55-57．
[13]王卉，陈楷，彭哲，等．数字仿真技术在电力系统中 的应用及常用的几种数字仿真工具【J】．继电器，2004，32(21)：7l一75．
wANG Hui，CHEN Kai，PENG zhe，et a1．Application of digital simulation眦hniques棚d severaJ simulation tools in power system[J】．Relay，2004，32(21)：71·75．
[14]IEEE Power Engmeering Socie哆．IEEE std 421．5．2005 IEEE玎ccOmmended practice for excitation system models for power system stabiI时studies【s】．
第五篇：电力系统潮流计算程序
电力系统潮流计算c语言程序，两行，大家可以看看，仔细研究，然后在这个基础上修改。谢谢
#include “stdafx.h” #include #include“Complex.h” #include“wanjing.h” #include“gauss.h” using namespace std; int _tmain(int argc, _TCHAR* argv[]){ int i; //i作为整个程序的循环变量 int N=Bus::ScanfBusNo();//输入节点个数 int L=Line::ScanflineNo();//输入支路个数 if((L&&N)==0){return 0;} //如果找不到两个文件中的任意一个，退出 Line *line=new Line[L];//动态分配支路结构体 Line::ScanfLineData(line);//输入支路参数 Line::PrintfLineData(line,L);//输出支路参数 Bus *bus=new Bus[N];//动态分配结点结构体 for(int i=0;i bus[i].Sdelta.real=0; bus[i].Sdelta.image=0;} Bus::ScanfBusData(bus);//输入节点参数 Bus::PrintfBusData(bus,N);//输出结点参数 Complex **X;X=new Complex *[N];for(i=0;i Bus::JisuanNodeDnz(X,line,bus,L,N);//计算节点导纳矩阵 Bus::PrintfNodeDnz(X,N);//输出节点导纳矩阵 int NN=(N-1)*2;double **JacAug;JacAug=new double *[NN];for(i=0;i double *x;x=new double[NN];int count=1; LOOP: Bus::JisuanNodeI(X,bus,N);//计算节点注入电流 Bus::JisuanNodeScal(X,bus,N);//计算节点功率 Bus::JisuanNodeScal(X,bus,N);//计算节点功率 Bus::JisuanNodeSdelta(bus,N);//计算节点功率差值 Bus::PrintfNodeScal(X,bus,N);//输出节点功率差值 int icon=wehcon1(bus,N);//whether converbence看迭代是否结束 if(icon==1){ cout<<“icon=”< Bus::JisuanJacAug(JacAug,X,bus,N);//计算雅可比增广矩阵 // Bus::PrintfJacAug(JacAug,N); gauss::gauss_slove(JacAug,x,NN);//解方程组求出电压差值 Bus::ReviseNodeV(bus,x,N);//修正节点电压 // Bus::PrintfNodeV(bus,N); count++; goto LOOP;} else { for(i=0;i { int statemp,endtemp; Complex aa,bb,cc,dd,B; B.real=0; B.image=-line[i].B; statemp=line[i].start; endtemp=line[i].end; aa=Complex::productComplex(Complex::getconj(bus[statemp-1].V), B); bb=Complex::subComplex (Complex::getconj(bus[statemp-1].V), Complex::getconj(bus[endtemp-1].V)); cc=Complex::productComplex(bb , Complex::getconj(line[i].Y)); dd=Complex::CaddC(aa,cc); line[i].stoe=Complex::productComplex(bus[statemp-1].V,dd); aa=Complex::productComplex(Complex::getconj(bus[endtemp-1].V), B); bb=Complex::subComplex (Complex::getconj(bus[endtemp-1].V), Complex::getconj(bus[statemp-1].V)); cc=Complex::productComplex(bb , Complex::getconj(line[i].Y)); dd=Complex::CaddC(aa,cc); line[i].etos=Complex::productComplex(bus[endtemp-1].V,dd); } cout<<“icon=”< Bus::JisuanNodeScal(X,bus,N);//计算节点功率 for(i=0;i { bus[i].Scal.real = bus[i].Scal.real + bus[i].Load.real;//发电机功率=注入功率+负荷功率 bus[i].Scal.image= bus[i].Scal.image+ bus[i].Load.image; bus[i].V=Complex::Rec2Polar(bus[i].V); } cout<<“====节点电压===============发电机发出功率======”< for(i=0;i { cout<<“节点”<<(i+1)<<'t'; Complex::PrintfComplex(bus[i].V); coutt(bus[i].Scal.real); coutt(bus[i].Scal.image); cout< } cout<<“======线路传输功率==========”< for(i=0;i { int statemp,endtemp; statemp=line[i].start; endtemp=line[i].end; cout< Complex::PrintfComplex(Complex::ComDivRea(line[i].stoe,0.01)); Complex::PrintfComplex(Complex::ComDivRea(line[i].etos,0.01)); cout< } } return 0;} #include “stdafx.h” #include #include“Complex.h” #include“wanjing.h” #include“gauss.h” using namespace std; int _tmain(int argc, _TCHAR* argv[]){ int i; //i作为整个程序的循环变量 int N=Bus::ScanfBusNo();//输入节点个数 int L=Line::ScanflineNo();//输入支路个数 if((L&&N)==0){return 0;} //如果找不到两个文件中的任意一个，退出 Line *line=new Line[L];//动态分配支路结构体 Line::ScanfLineData(line);//输入支路参数 Line::PrintfLineData(line,L);//输出支路参数 Bus *bus=new Bus[N];//动态分配结点结构体 for(int i=0;i bus[i].Sdelta.real=0; bus[i].Sdelta.image=0;} Bus::ScanfBusData(bus);//输入节点参数 Bus::PrintfBusData(bus,N);//输出结点参数 Complex **X;X=new Complex *[N];for(i=0;i Bus::JisuanNodeDnz(X,line,bus,L,N);//计算节点导纳矩阵 Bus::PrintfNodeDnz(X,N);//输出节点导纳矩阵 int NN=(N-1)*2;double **JacAug;JacAug=new double *[NN];for(i=0;i double *x;x=new double[NN];int count=1; LOOP: Bus::JisuanNodeI(X,bus,N);//计算节点注入电流 Bus::JisuanNodeScal(X,bus,N);//计算节点功率 Bus::JisuanNodeScal(X,bus,N);//计算节点功率 Bus::JisuanNodeSdelta(bus,N);//计算节点功率差值 Bus::PrintfNodeScal(X,bus,N);//输出节点功率差值 int icon=wehcon1(bus,N);//whether converbence看迭代是否结束 if(icon==1){ cout<<“icon=”< Bus::JisuanJacAug(JacAug,X,bus,N);//计算雅可比增广矩阵 // Bus::PrintfJacAug(JacAug,N); gauss::gauss_slove(JacAug,x,NN);//解方程组求出电压差值 Bus::ReviseNodeV(bus,x,N);//修正节点电压 // Bus::PrintfNodeV(bus,N); count++; goto LOOP;} else { for(i=0;i { int statemp,endtemp; Complex aa,bb,cc,dd,B; B.real=0; B.image=-line[i].B; statemp=line[i].start; endtemp=line[i].end; aa=Complex::productComplex(Complex::getconj(bus[statemp-1].V), B); bb=Complex::subComplex (Complex::getconj(bus[statemp-1].V), Complex::getconj(bus[endtemp-1].V)); cc=Complex::productComplex(bb , Complex::getconj(line[i].Y)); dd=Complex::CaddC(aa,cc); line[i].stoe=Complex::productComplex(bus[statemp-1].V,dd); aa=Complex::productComplex(Complex::getconj(bus[endtemp-1].V), B); bb=Complex::subComplex (Complex::getconj(bus[endtemp-1].V), Complex::getconj(bus[statemp-1].V)); cc=Complex::productComplex(bb , Complex::getconj(line[i].Y)); dd=Complex::CaddC(aa,cc); line[i].etos=Complex::productComplex(bus[endtemp-1].V,dd); } cout<<“icon=”< Bus::JisuanNodeScal(X,bus,N);//计算节点功率 for(i=0;i { bus[i].Scal.real = bus[i].Scal.real + bus[i].Load.real;//发电机功率=注入功率+负荷功率 bus[i].Scal.image= bus[i].Scal.image+ bus[i].Load.image; bus[i].V=Complex::Rec2Polar(bus[i].V); } cout<<“====节点电压===============发电机发出功率======”< for(i=0;i { cout<<“节点”<<(i+1)<<'t'; Complex::PrintfComplex(bus[i].V); coutt(bus[i].Scal.real); coutt(bus[i].Scal.image); cout< } cout<<“======线路传输功率==========”< for(i=0;i { int statemp,endtemp; statemp=line[i].start; endtemp=line[i].end; cout< Complex::PrintfComplex(Complex::ComDivRea(line[i].stoe,0.01)); Complex::PrintfComplex(Complex::ComDivRea(line[i].etos,0.01)); cout< } } return 0;} #include class Complex//定义复数类 { public: double real;double image;int RecPolar;//0表示直角坐标，1表示极坐标 static Complex CaddC(Complex c1,Complex c2);//求两个复数和 static Complex subComplex(Complex c1,Complex c2);//求两个复数差 static Complex productComplex(Complex c1,Complex c2);//求两个复数积 static Complex divideComplex(Complex c1,Complex c2);//求两个复数商 static Complex ComDivRea(Complex c1,double r2);//除数 static Complex getconj(Complex c1);//求一个复数共轭 static Complex getinverse(Complex c1);//取倒数 static double getComplexReal(Complex c1);//求一个复数实部 static double getCompleximage(Complex c1);//求一个复数虚部 static void PrintfComplex(Complex c1);//显示一个复数 static void PrintfmultiComplex(Complex C,int N);//显示多个复数 static void zeroComplex(Complex c1);//将复数复零 static Complex Rec2Polar(Complex c1);//取极坐标 Complex(){ RecPolar=0;} }; Complex Complex::Rec2Polar(Complex c1)//极坐标表示 { Complex Node;Node.real=sqrt(c1.real*c1.real+c1.image*c1.image);Node.image=atan2(c1.image,c1.real)*180/3.1415926;Node.RecPolar=1;return Node;} Complex Complex::CaddC(Complex c1,Complex c2)//复数加法 { Complex Node; Node.real=c1.real+c2.real; Node.image=c1.image+c2.image; return Node;} Complex Complex::subComplex(Complex c1,Complex c2)//复数减法 { Complex Node; Node.real=c1.real-c2.real; Node.image=c1.image-c2.image; return Node;} Complex Complex::productComplex(Complex c1,Complex c2)//复数乘法 { Complex Node; Node.real=c1.real*c2.real-c1.image*c2.image; Node.image=c1.image*c2.real+c2.image*c1.real; return Node;} Complex Complex::divideComplex(Complex c1,Complex c2)//复数除法 { Complex Node; Node.real=(c1.real*c2.real+c1.image*c2.image)/(pow(c2.real,2)+pow(c2.image,2));Node.image=(c1.image*c2.real-c1.real*c2.image)/(pow(c2.real,2)+pow(c2.image,2));return Node;} Complex Complex::ComDivRea(Complex c1,double r1)//复数除数 { Complex Node;Node.real=c1.real/(r1);Node.image=c1.image/(r1);return Node;} Complex Complex::getconj(Complex c1)//取共轭 { Complex Node; Node.real=c1.real;Node.image=-c1.image; return Node;} Complex Complex::getinverse(Complex c1)//取倒数 { Complex Node;Node.real=1;Node.image=0;Node=(Complex::divideComplex(Node,c1));return Node;} double Complex::getComplexReal(Complex c1)//取实部 { return c1.real;} double Complex::getCompleximage(Complex c1)//取虚部 { return c1.image;} void Complex::PrintfComplex(Complex c1)//按直角坐标输出 { if(c1.RecPolar==0){ cout.precision(6); cout.width(8); cout.setf(ios::right); cout< ”; cout.precision(6); cout.width(8); cout.setf(ios::left); cout< ”;} else { cout< Complex::zeroComplex(Complex c1)//清零 { c1.real=0;c1.image=0;} class gauss { public: static void gauss_slove(double **a,double *x,int NN);static void gauss_output();}; void gauss::gauss_slove(double **a,double *x,int NN){ int n,i,j,k,*pivrow,**pivcol;double eps,pivot,sum,aik,al; n=NN;pivrow=new int[n];pivcol=new int *[n]; for(i=0;i pivot= fabs(a[k][k]); pivrow[k]=k;//行 pivcol[k][0]=k;pivcol[k][1]=k;//列n*2矩阵 for(i=k;i { for(j=k;j { if(pivot { pivot=fabs(a[i][j]); pivrow[k]=i;//行 pivcol[k][1]=j;//列 } } } if(pivot { cout<<“error”< getchar(); exit(0); } if(pivrow[k]!=k)//行变换 { for(j=k;j<(n+1);j++) { al=a[pivrow[k]][j]; a[pivrow[k]][j]=a[k][j]; a[k][j]=al; } } if(pivcol[k][1]!=k)//列变换 { for(i=0;i { al=a[i][pivcol[k][1]]; a[i][pivcol[k][1]]=a[i][k]; a[i][k]=al; } } if(k!=(n-1))//将矩阵化为上三角形 式 { for(i=(k+1);i { aik=a[i][k]; for(j=k;j<(n+1);j++) { a[i][j]-=aik*a[k][j]/a[k][k]; } } } } x[n-1]=a[n-1][n]/a[n-1][n-1];//解方程 for(i=(n-2);i>=0;i--){ sum=0; for(j=(i+1);j { sum +=a[i][j]*x[j];0.182709 0.016894-0.0310701 -0.0402051 0.156702 -0.0355909-0.0668055 -0.00703229-0.0886481 -0.0129814-0.0390805 -0.0135062-0.1023 -0.0460568 -0.0342827 -0.00382402-0.102896 -0.0184062 } x[i]=(a[i][n]-sum)/a[i][i];} for(k=(n-2);k>=0;k--){ al=x[pivcol[k][1]]; x[pivcol[k][1]]=x[pivcol[k][0]]; x[pivcol[k][0]]=al;} cout<<“节点电压修正量”< cout< } ====节点功率计算值==== 0.935261 -0.159048 0.573909 0.0789973-0.00289889 -0.00796623-0.0791247 -0.0168362-0.436255 -0.0580392 0.0359139 -0.0106592-0.229118 -0.0885419-0.136179 -0.148207 0.0446243 0.0111298-0.0223764 -0.00695775-0.0237482 -0.198318 -5.24266e-015 -0.0354071 -0.0925078 -1.05629e-015 -0.0391348 0.014529 0.00158644 -0.0258771 -0.109514 icon=1进行第2次迭代 节点电压修正量 =================-0.00164889-0.000540034-0.00261067-0.00532027-0.00235315-0.00600971-0.00189677-0.00643874-0.0023631-0.00650659-0.00170949-0.0074907-0.00164545-0.00485415-0.00493977-0.0119042-0.00331285-0.0175611-0.00207908 -0.00347744-0.0869347-9.48909e-015-0.0110778-0.0538236-7.53784e-016-0.0168097 7.049e-005-0.00146487-0.00458276 0.00251645 -0.00336375-0.00530645-0.0147816-0.000326161-0.00640487-0.00251701-0.0169829-0.00175286-0.0174333-0.0239063 -0.0119192-0.076014 -0.0160104-0.441997 -0.0750285 0.000250012 3.72542e-005-0.228052 -0.108844-0.100078 -0.105634 0.000410707 0.000378067-0.057497 -0.0195879 0.200039 0.0582563-0.00307326-0.0163809-0.00232773-0.0175806 8.74293e-005-0.0192018 0.000558996-0.0197776-0.000247851-0.0193784-0.00115346-0.0185848-0.00127275-0.0186244-0.00010108-0.0188966 0.000553585-0.0200901-3.76315e-005-0.0208303 0.00308341-0.0219386-0.00195916-0.0205356-0.00184757-0.0076401 0.00197593-0.0245534 0.00434657-0.027534 ====节点功率计算值==== 0.98623 -0.134163 0.583136 0.166278-0.111173 0.199792 -0.0621041 -0.0821379 -0.0350785 -0.0902383 -0.0320461 -0.0951562 -0.0220362 -0.175458 4.72557e-015 -0.0320661 -0.0871134 -7.03489e-017 -0.0350769 0.000273455 1.51804e-005 -0.0240417 -0.10604 icon=1进行第3次迭代 节点电压修正量 =================-2.67079e-005-2.30128e-006-2.20543e-005-6.00686e-005-2.33043e-005-6.85601e-005-3.22294e-005-2.61107e-005-2.80198e-005-6.6167e-005-2.34528e-005 -0.0739846 0.0227868-0.0158709-0.0248173-0.0179447-0.0578368-0.00890719-0.0337091-0.00693706-0.111601 1.21429e-014-0.0159145-0.0667319 9.24355e-016-0.0228592 7.10354e-005-6.6188e-006-0.00889343-0.0184098 -5.66132e-005-4.4646e-005-1.74668e-005-4.50947e-005-0.000181763-3.81763e-006-0.000286581-6.68993e-005-1.28441e-005-5.17172e-005-0.000223284-4.54717e-005-2.47586e-005 4.32335e-007-0.000258494 1.82635e-005-0.000272051-6.95195e-006-0.000251969 1.11318e-005-0.000279418 5.74737e-005-0.000307368 6.86998e-005-0.000320274 5.38112e-005-0.00031447 3.59531e-005-0.00030494 3.37607e-005-0.000307449 5.26532e-005-0.000310721 6.92761e-005-0.000350373 5.60942e-005-0.00040977 0.000123641-0.000440259 1.36149e-005-0.000426973-1.70227e-005-9.37794e-005 0.000113675-0.000544011 0.000176034-0.000636202 ====节点功率计算值==== 0.986878 -0.133979 0.583 0.167193-0.024 -0.012-0.076 -0.016-0.442 -0.0748606 1.43501e-008 1.07366e-008-0.228 -0.109 -0.0999999 -0.104049 4.51318e-008 8.98835e-008-0.0579999 -0.0199999 0.2 0.0591018-0.112 -0.0749997 0.2 0.0242519-0.062 -0.016-0.082 -0.025-0.035 -0.018 -0.0900001 -0.058-0.032 -0.00899997-0.095 -0.0339999-0.022 -0.00699998-0.175 -0.112 -6.07156e-015 -1.19217e-014-0.032 -0.016-0.087 -0.0669999 7.03078e-017 -9.23979e-016-0.035 -0.0229999 1.09492e-007 4.45699e-008 1.54958e-009 -2.01531e-010-0.024 -0.00899994-0.106 -0.0189996 icon=0,迭代结束。 ====节点电压===============发电机发出功率====== 节点1 1.05 0。 98.6878-13.3979 节点2 1.045 -1.846。 29.4193 节点3 1.02384-3.83352。 0 节 点25 1.01216-9.68486。 0 0 0 节点4 1.01637-4.55698。 0 节 点26 0.994393 -10.1089。 0 0 0 节点5 1.01 -6.48617。 节 点27 1.02012-9.42025。 0 11.5139 0 节点6 1.01332-5.38073。 0 节 点28 1.00992-5.86244。 0 0 0 节点7 1.00489-6.38368。 0 节 点29 1.00022-10.6579。 0 0 节点8 19.5951 节点9 0 节点10 0 节点11 5.91018 节点12 0 节点13 2.42519 节点14 0 节点15 0 节点16 0 节点17 0 节点18 0 节点19 0 节点20 0 节点21 0 节点22 0 节点23 0 节点24 0 1.01 -5.62974。 1.03905-6.78143。 1.03595-8.69362。 -4.5962。 1.04711-7.80323。 1.05 -6.34392。 1.03242-8.7401。 1.02788-8.86784。 1.03458-8.45044。 1.03051-8.83678。 1.01845-9.5141。 1.01604-9.70326。 1.02022-9.50938。 1.0237-9.17478。 1.02432-9.17024。 1.01802-9.36719。 1.01339-9.68362。 0 20 节 点30 0.988705 -11.5464。 0 0 0 ====== 线路传输功率========== 2to1 -57.7373 5.41674i 58.3454 0 -15.1827i 3to1 -39.659 -7.75964i 40.3424 1.78481i 4to2 -30.87 -9.74186i 31.4153 0 3.58352i 4to3 -37.0772 -7.78596i 37.259 6.55964i 5to2 -44.3717 -9.78456i 45.2968 0 4.84242i 6to2 -38.4766 -8.22625i 39.3252 0 2.87667i 6to4 -34.946 1.92384i 35.0885 0 -3.28202i 7to5 -0.16304 -6.41767i 0.171702 0 2.2985i 7to6 -22.637 -4.48233i 22.7745 0 1.44238i 8to6 -11.8939 -5.48098i 11.913 0 3.70557i 6to9 12.3737 -12.3826i -12.3737 0 13.0033i 6to10 10.9107 -3.80907i -10.9107 0 4.53223i 11to9 5.91018i 0 -5.08963i 10to9 -32.652 -2.3712i 32.652 0 3.46974i 4to12 23.5411 -11.5375i -23.5411 0 13.2407i 13to12 2.42519i 1.05 -1.90978i 1.66484i 14to12 -7.9019 -2.06732i 7.97894 30to29 -3.6702 -0.542564i 3.70398 2.22749i 0.606393i 15to12 -18.254 -5.74885i 18.4835 28to8 -1.89152 -3.79982i 1.89395 6.20089i-4.9239i 16to12-7.53872 -2.90237i 7.59633 28to6 -14.7868 -2.82565i 14.8234 3.02352i 0.294601i 15to14-1.69544 -0.461488i 1.70189 请按任意键继续...0.467323i 17to16-4.03014 1.10238i 18to15-6.08074 1.46028i 19to18-2.87549 0.478389i 20to19 6.6418-2.93222i 20to10 -8.8418 3.85077i 17to10-4.96987 4.76656i 21to10-16.1562 9.42843i 22to10-7.87782 4.21401i 22to21 1.34443-2.01837i 23to15-5.59369 2.25006i 24to22-6.48186 2.08163i 24to23-2.38596 0.579814i 25to24-0.167617 0.281364i 26to25 -3.5 2.3674i 27to25 3.39433-2.08638i 28to27 16.1446 3.13006i 29to27-6.10398 1.67047i 30to27-6.92979-1.07089i-1.37839i-0.467767i 2.96679i-3.66679i-4.72911i-9.18162i-4.10132i 2.01969i-2.17981i-2.00141i-0.56401i -0.28102i-2.29999i 2.11848i-2.10093i-1.50639i -1.3574i 4.03872 6.12096 2.88074 -6.62452 8.9242 4.98423 16.2709 7.93248 -1.34378 5.62846 6.53339 2.39369 0.167814 3.54513 -3.37751 -16.1446 6.19083 7.09313 高等电力系统分析 IEEE30节点潮流程序 班级：电研114班 姓名：王大伟 学号：2201100151
点击咨询 