电力系统在线暂态分析计算

电力系统在线暂态分析计算

电力系统暂态稳定是指电力系统受到大干扰后，各同步发电机保持同步运行并过

渡到新的或恢复到原来稳态运行方式的能力。通过暂态稳定计算可以判别系统的暂态

功角稳定性、暂态电压稳定性和暂态频率稳定性。

1.1暂态稳定计算问题的数学描述

暂态稳定计算首先是以某一稳态的潮流计算结果作为初始状态。此外，暂态稳定

计算的数学模型包括一次电网的数学描述(网络方程)和发电机，负荷，无功补偿，直

流输电，发电机的调压器、调速器、电力系统稳定器、继电保护等一次设备和二次装

置动态特性的数学描述(微分方程)以及各种可能发牛的扰动方式和稳定措施的模拟

等。因此PSASP暂态稳定计算(ST)的数学模型可归为以下三个部分：

①电网的数学模型，即网络方程：

X=F(X,y)(5-1)

其中：

F=5/,……/)'

X=(X»2，……为网络方程求解的变量

②发电机、负荷等一次设备和二次自动装置的数学模型，即微分方程：

y=G(X,y)(5-2)

其中：

G=(gi'g?，...

Y=(乂，为，……‘,为微分方程求解的变量

③扰动方式和稳定措施的模拟，如电网的简单故障或复杂故障及冲击负荷、

快关汽门、切机、切奂荷、切线路等。这些因素的作用结果是改变X,Yo

1.2暂态稳定计算的主要功能和特点

PSASP暂态稳定计算的主要功能可概括为以下几个方面:

1）一般模型的计算功能

＞可计算交直流混合电力系统

＞可考虑变电站（主接线）内部的开关状态对系统网络结构的影响

＞程序提供的常用系统元件模型（固定模型）如下：

＞同步电机模型（7种）

＞励磁调节器模型

＞原动机调速器模型

＞电力系统稳定器（PSS）模型

＞感应电动机及综合动态负荷模型

＞静态负荷模型

＞静止无功补偿器模型

＞直流输电模型

2）复杂故障方式的计算功能

＞可同时考虑多处三相对称故障，包括三相短路，三相断线，串联电容保护

三相击穿等

＞可同时考虑多处不对称故障，包括单相短路，单相负荷投入，两相短路,

两相接地，单相断线，两相断线，串联电容器不对称击穿等

＞在不对称故障下，可考虑输电线零序互感的影响

＞既可做暂态稳定计算，也可做短路电流计算（在程序中做不解微分方程处

理），及动态过程中的复杂故障短路电流计算

＞可做输电线工频暂态过电压及潜供电流计算分析

＞能给出三相不平衡方式下序电压，序电流，相电压，相电流的分右

3）扰动方式和稳定措施模拟的计算功能

在PSASP暂态稳定中，除网络故障外，还设置了一些扰动方式，用以模拟电力系

统的某些冲击和稳定措施，其功能可概括如下：

＞冲击负荷对电力系统影响的动态仿真

＞负荷功率随机波动的模拟

＞励磁回路电压波动的模拟

＞发电机出力跟踪调节的模拟

＞调速器汽门快速调整的模拟

＞调压器励磁电压快速调整的模拟

＞按不同准则（时间、角度、高频率、频率差、过电流）减少母线上的部分和

全部机组的出力

＞按不同准则（时间、低电压、低频率、频率差、线路过电流、异步电动机自

身过电流）切除母线上的部分和全部负荷

＞按不同准则（时间、低电压、低频率、角度差、零序过电流、正序过电流）

开断线路开关

4）通过用户自定义模型的方法，用户可建立各种模型以实现所需的计算功能

＞电力系统各种一次设备模型，如不同型号的同步电机，异步电机，静止无

功补偿器等

＞电力系统各种自动装置的模型，如：各种类型的调压静、调速器、电力系

统稳定器（PSS）及各种继电保护和安全自动装置等

＞随不同工程而异的超高压直流输电线路及其控制系统的模型

＞灵活交流输电系统（FACTS）的元件的模型，如：可控硅串联补偿器（TCSC）,

统一潮流控制器（UPFC）等

1.3暂态稳定计算的元件数学模型

1.3.1励磁系统模型库

励磁系统模型是同步发电机励磁系统物理过程的数学描述。作为电力系统机电暂

态过程数学模型的重要组成部分，主要应用于电力系统稳定的计算分析。

我国电力系统稳定t-算分析，长期采用暂态电势（Eq，）恒定的模型，它不能准确地

反映各种不同励磁系统的物理过程。新的稳定分析导则建议在安全稳定分析中发电机

采用Eq“、Ed”电势变化的详细模型，同时也要求在暂态稳定分析中要考虑各种励磁系

统模型。PSASP用户自定义系统模型库包含很多种司步发电机励磁系统模型，可以用

来模拟我国电力系统中已有的绝大部分的励磁系统，满足电力系统稳定计算的需要。

类别励磁系统特点模型号

具有复式励磁加电压校正器的直流励磁机励磁系1104

统，带并联校正

直

流具有快速相复励调节器的直流励磁机励磁系统，带1105

励并联校正

磁带由恒定电源供电的电子型励磁调节器的直流励1106

机磁机励磁系统，带串联校正

励带由机端供电的电子型励磁调节器的直流励磁机1116

磁励磁系统，带串联校正

系带由恒定电源供电的电子型励磁调节器的直流励1107

统磁机励磁系统，带并联校正

带由机端供电的电子型励磁调节器的直流励磁机1117

励磁系统，带并联校正

有副励磁机的交流励磁机不可控整流器励磁系统，1108

带并联校正

无副励磁机的交流励磁机不可控整流器励磁系统，1118

交带并联校正

流有副励磁机的交流励磁机静止整流器（有刷）励磁1109

励系统，昔串联校正

磁无副励磁机的交流励磁机静止整流器（有刷）励磁1119

机系统，带串联校正

励

有副励磁机的交流励磁机不可控整流器（无刷）励1110

磁磁系统，带并联校正

系

无副励磁机的交流励磁机不可控整流器（无刷）励1120

统

磁系统，带并联校正

交流励磁机可控整流器励磁系统（它励可控整流器1111

励磁系统），带串联校正

静电压源可控整流器励磁系统（自并励励磁系统），1112

止带串联校正

励

交流侧串联相加的复励源可控整流器励磁系统（自1113

磁复励励磁系统），带串联校正

系

交流侧串联相加的复励源可控整流器励磁系统（自1114

统

复励励磁系统），带串联校正，较好地模拟了可控

整流器特性

表5-1给出了各模型所适用的励磁系统。

几点说明：

(1)直流励磁机饱和系数

对于直流励磁机励磁系统，励磁机饱和系数工的计算公式为:

=

Cl和C2的确定方法如下：

在励磁机特性曲线上取两点，相应的Efd为Efdmax和0.7Efdmax,求得Sei和

Se2,则有:

se\=Ge

S,？=C}e

联立求解可得到系数Cl和C2。

如图5-1所示，对每一点Se的计算公式为

Se=(A-B)/B

图5-1励磁机饱和特性

(2)不可控整流器的模型

不可控整流器的模型如图5-2所示。

1-0.58/

图5-2不可控整流器数学模型

1.2.1.2模型说明

(1)模型1104

1104励磁系统数学模型的框图如图5-3所示。励磁调节器的电源取自发电机定子

电流和机端电压，两者在直流侧相加，为了提高励磁控制系统的稳定性，采用了以发

电机磁场电压为信号的并联校正。

这种模型主要用来模拟具有复式励磁加电压校正器的直流励磁机励磁系统。

图5-3中各参数物理意义如下：

Ki电流复励分量系数

Ta校正器时间常数

Ka校正添增益

Ke励磁机增益，一般为1

Te励磁机时间常数

Kf并联校正增益

Tf并联校正时间常数

Tr电压测量回路时间常数

Vamax校正器最大输出

Vamin校正器最小输出

VrmaxAER最大输出

VrminAER最小输出

Xc.调差电抗，Xc.>0时为正调差，Xc<0时为负调差。

Ke自励系数，反映直流励磁机自励分量所占份额，完全它励时Ke=l,完全自

励时，Ke=0o

Se直流励磁机的饱和系数，不计饱和时Se=0,越是饱和，Sc越大。

Vt发电机端电压

It发电机定子电流

。发电机的运行功率因数角

Efd发电机的磁场电压

1104励磁系统自定义模型框图如图5-4所示。

图5-31104励磁系统数学模型

图5-41104励磁系统自定义模型框图

(2)模型1105

1105励磁系统数学模型的框图如图5-5所示。励磁调节器的电源取自发电机定子

电流和端电压，与1104模型不同的是，两个分量在交流侧相加，因而具有相位补偿作

用。另外，电压校正器所需功率少，其时间常数比1104模型的校正器要小的多。为了

提高励磁控制系统的稳定性，采用了并联校正环节，以发电机磁场电压为输入信号。

这种模型主要用来模拟具有快速相复励调节器的直流励磁机励磁系统c

图5-5中各参数物理意义如下：

Kv电压复励系数

Ki电流复励系数

Ta放大器时间常数

Ka放大器放大倍数

Tr测量回路时间常数

Vamax放大器最大输出

Vamin放大器最小输出

Vrmax调节器最大输出

Vrmin调节器最小输出

Ke励磁机增益，一般为1.0

Te励磁机时间常数

Se励磁机的饱和系数

Ke励磁机的自励系数

Vt发电机端电压

It发电机定子电流

*发电机功率因数角

Kf并联校正环节增益

Tf并联校正环节时间常数

Xc调差电抗

1105励磁系统自定义模型框图如图5-6所示。

图5-51105励磁系统数学模型

图5-61105励磁系统自定义模型框图

(3)模型1106

1106励磁系统数学模型的框图如图5-7所示。采用半导体电子型调节器，功率放

大器由电力电子元件组成，反应速度快。调节器的功率电源取自恒定电压。帝串联校

正(PID调节)和发电机磁场电压硬反馈，后者用以提高励磁系统的响应速度。

这种模型主要用来模拟带半导体电子型励磁调节器，并采用恒定电源供电的直流

励磁机励磁系统。

图5-7中各参数物理意义如下：

Tr电压测量回路时间常数

Ti、

T[

Tc>串联校正回路时间常数，一般有Tb<Tc<Td<Ti

Th-

Kv积分校正选择因子，Kv=0时为纯积分型校正，Kv=l时为比例积分型校正

Kj校正回路的直流增益

Ka功率放大单元的增益

Ta功率放大单元的时间常数

Vrmax调节器最大输出

Vrmin调节器最小输出

Te励磁机时间常数

Ke励磁机自励系数

Se励磁机饱和系数

Kg用来调节时间常数补偿度的比例反馈系数

Xc调差电抗

Vt发电机定子电压

It发电机定子电流

夕发电机功率因数角

Efd发电机的磁场电压

1106励磁系统自定义模型框图如图5-8所示。

图5-71106励磁系统数学模型

图5-81106励磁系统自定义模型框图

(4)模型1116

1116励磁系统数学模型的框图如图5-9所示。它是由1106模型演变而来的，与

1106模型的唯一区别是，1116励磁系统的励磁调节器的电源取自发电机端电压，因

此，调节器的最大输出和最小输出也和发电机端电压成比例。

1116模型用来模拟励磁调节器由机端变压器供电的直流励磁机励磁系统。

图5-9中的各参数物理意义如下：

Vrmax发电机机端电压为额定值时，励磁调节器最大输出电压

Vrmin发电机机端电压为额定值时，励磁调节器最小输出电压

其它参数的物理意义与1106模型完全相同,

1116励磁系统自定义模型框图如图570所示。

图5-91116励磁系统数学模型

图5T01116励磁系统自定义模型框图

(5)模型1107

和1106励磁系统数学模型相同，1107励磁系统数学模型主要用来模拟带电子型

励磁调节器的直流励磁机励磁系统。励磁调节器的功率元件用恒定电源供电，但提高

励磁系统稳定性的方法是采用了并联校正，信号取自发电机磁场电压。

1107励磁系统数学模型的框图如图5Tl所示。

图5-11中的参数物理意义如下：

Kf并联校正环节增益

Tf并联校正环节时间常数

其它参数的物理意义与1106模型完全相同°

1107励磁系统自定义模型框图如图5-12所示。

图5-111107励磁系统数学模型

图5-121107励磁系统自定义模型框图

(6)模型1117

1117励磁系统数学模型的框图如图5-13所示。其主要部分和1107励磁系统模型

是相同的，唯一的区别是，励磁调节器的功率元件由发电机端供电，它的最大输出和

最小输出是发电机端电压的函数，分别为Vt•Vrmax和Vt•Vrmin,Vrmax^Vrmin分

别是发电机端电压为额定值时的最大输出和最小输出。

1117励磁系统数学模型主要用来模拟有由机端供电的电子型励磁调节器的直流

励磁机励磁系统。

1117励磁系统自定义模型框图如图5-14所示。

图5T31117励磁系统数学模型

图5-141117励磁系统自定义模型框图

(7)模型1108

图5-15为1108励磁系统数学模型的框图。

本模型模拟的是具有不可控整流器的交流励磁机励磁系统，励磁调节器带并联校

正环节，校正信号取自Efd,功率电源取自同轴副励磁机，另外还有励磁机时间常数

的补偿反馈，信号取Efe。

图5-15中各参数物理意义如下：

Ve整流换相电抗后励磁机电压

Te交流励磁机空载时磁场绕组时间常数

Kd交流励磁机负载电流电枢反应的去磁系数，是交流励磁机电抗的函数

Kc与换相电抗相关的整流器负荷系数

Fex整流器的负荷因子

Kf并联校正环节增益

Tf并联校正环节时间常数

Tc和Tb励磁系统暂态增益衰减用时间常数，实际的调节器中，很少有这一环

节或只有增益而无时间常数，此时，可填Tb二Tc的任意数，但Tb不要为零

Kj暂态增益衰减环节的增益

Ifd发电机励磁电流，即转子电流，在标幺值系统中，它和同步电势Eq相等。

Tr、Xc、Kf、Kg、Ka、Ta与1106励磁系统中的物理意义相同。

1108励磁系统自定义模型框图如图5-16所示。

1+药

图5T51108励磁系统数学模型

18tfD

图5-161108励磁系统自定义模型框图

(8)模型1118

1118励磁系统数学模型用来模拟没有副励磁机的交流励磁机不可控整流器励磁

系统。此时，励磁调节器的功率电源由发电机机端电压供给，调节器的最大输出值、

最小输出值都与发电机机端电压有关。图5-17为其框图。Vrmax>Vrmin分别为发电

机电压为额定值时，调节器最大输出和最小输出，其它参数的物理意义与1108模型完

全相同。

1118励磁系统自定义模型框图如图5-18所示。

图5-171118励磁系统数学模型

图5-181118励磁系统自定义模型框图

（9）模型1109

图5-19为1109励磁系统数学模型的框图，励磁调节港中采用串联校正（PID调

节）。励磁机时间常数补偿信号取自发电机磁场电压Efd。

本模型主要用于模拟带串联校正的交流励磁机静止整流渊（有刷）励磁系统。

图5-19中各参数物理意义如下：

Tr电压测量回路时间常数

Ti、

Td

TeI串联校正回路时间常数，一般有Tb<Tc<Td<Ti

Tb•

Kv积分校正选择因子，Kv=0时为纯积分型校正，Kv二1时为比例积分型校正

KJ校正回路的直流增益

Ka功率放大单元的增益

Ta功率放大单元的时间常数

Te交流励磁机空载时磁场绕组时间常数

Kd交流励磁机负载电流电枢反应的去磁系数，是交流励磁机电抗的函数

Kg用来调节时间常数补偿度的比例反馈系数

Kc与换相电抗相关的整流器负荷系数

Xc调差电抗

Ve整流换相电抗后励磁机电压

Eex整流器的负荷因子

Vrmax调节器最大输出

Vrmin调节器最小输出

1109励磁系统自定义模型框图如图5-20所示。

图5-191109励磁系统数学模型

图5-201109励磁系统自定义模型框图

(10)模型1119

图5-21为1119励磁系统数学模型的框图，1119励磁系统用来模拟没有囿励磁机

的交流励磁机静止整流器(有刷)励磁系统，调节器采用串联校正。

图5-21中各参数的物理意义除Vrmax、Vrmin外,与1109模型完全相同。Yrmax、

Vrniin分别是发电机电压为额定值时调节器的最大输出和最小输出。

1119励磁系统自定义模型框图如图5-22所示。

图5-211119励磁系统数学模型

图5-221119励磁系统自定义模型框图

(11)模型1110

图5-23为1110励磁系统数学模型的框图。与1108模型一样，1110励磁系统模

型用来模拟由副励磁机向调节器供电的不可控整流交流励磁机励磁系统，不同之处是,

并联校正的信号、时间常数补偿信号均取自调节器的输出电流(即交流励磁机的励磁

电流)。它主要应用于无刷励磁系统。

图5-23中各参数的物理意义1108模型完全相同。

1110励磁系统自定义模型框图如图5-24所示。

图5-231110励磁系统数学模型

图5-241110励磁系统自定义模型框图

(12)模型1120

图5-25为1120励磁系统数学模型的框图。1120励磁系统模型用来模拟没有副励

磁机的交流励磁机不可挖整流器励磁系统，与1H0模型一样，采用并联校正，调节器

功率单元由机端变供电，并联校正和时间常数补偿信号均取自调节器的输出电流而不

是发电机的磁场电压，它主要也是应用于无刷励磁系统。除Vnnax、Vrmin外，图5-

25中各参数的物理意义，与1110模型完全相同。Vrmax、Vrmin分别是发电机电压为

额定值时调节器的最大输出和最小输出。

1120励磁系统自定义模型框图如图5-26所示。

图5-251120励磁系统数学模型

图5-261120励磁系统自定义模型框图

（13）模型1111

图5-27为1111励磁系统数学模型的框图。这种励磁系统，没有副励磁机，交流

励磁机采用自励恒压系统，其输出电压能保持不变，发电机励磁的调节是通过由交流

励磁机供电的可控整流器导通角的控制而实现的，乂称交流励磁机可控整流器励磁系

统（它励可控整流器励磁系统）。调节器采用串联校正（PID调节）。

图5-27中各参数物理意义如下：

Tr电压测量回路时间常数

7n

T,

Tc|串联校正回路的时间常数，一般有Ti>Td>Tc>Tb

ThJ

Kv积分校正选择因子，Kv=0时为纯积分型校正，Kv=l时为比例积分型校正。

KJ串联校正环节增益

Kc与换弧电抗相关的整流器负荷系数

Xc调差电抗

Vrmax调节器最大输出

Vrmin调节器最小输出

1111励磁系统自定义模型框图如图5-28所示。

图5-271111励磁系统数学模型

图5-281111励磁系统自定义模型框图

(14)模型1112

图5-29为1112励磁系统数学模型的框图。这种模型模拟的励磁系统称为电压源

可控整流器励磁系统，又称自并励励磁系统。与1111励磁系统模型不同的是，可控整

流器是由发电机端的整流变压器(又称机端变)供电，没有旋转部件，是静止励磁系

统的一种。它的最大输出电压是发电机端电压的函数，成正比例的关系。

图5-29中，Vrmax.Vrmin分别为发电机电压为额定值时调节器的最大输出和最

小输出。

Kc的意义与1111模型相同，由整流变压器漏抗决定。由于整流变压器的漏抗较

小，工作在调节特性的第1段，换相压降用Ifd・Kc来模拟。其它各参数的物理意义,

与1111模型完全相同。

1112励磁系统自定义模型框图如图5-30所示。

v

图5-291112励磁系统数学模型

图5-301112励磁系统自定义模型框图

(15)模型1113

图5-31为1113励磁系统数学模型的框图。一种由发电机端电压及定子电流供电

的可控整流器励磁系统称为复励源可控整流器励磁系统，它也是静止励磁系统的一种。

可控整流器的供电电压由两部分组成，一部分是与发电机端电压成比例的分量

Vl=Ki.Vt,由并联变压器提供：另一部分为与发电机定子电流成比例的分量

V2=jKiTtXu,由串联变压器提供。向量相加后向整流器供电。由于串联变压器有较大

的漏抗，换弧压降由整流器负荷因子Fex来模拟。

图5-31中，

Fox整流器的负荷因子

Ve电压分量和电流分量相加后的有效值

Xu串联变压器的互感抗

其它各参数的物理意义，与1111模型完全相同。

1113励磁系统自定义模型框图如图5-32所示。

图5-311113型励磁系统数学模型

图5-321113型励磁系统自定义模型框图

（16）模型1114

图5-33为1114励磁系统数学模型的框图。它与1113励磁系统模型一样，用来

模拟交流侧串联相加的复励源可控整流器励磁系统（自复励励磁系统）。不同之处是

较好地模拟了可控整流器的特性，图中Vcmax表示考虑换相压降后整流可能输出的最

大电压（整流时控制角a=0°,逆变时a=180°）,只用作限幅用。图5-33中各参数

的物理意义与1113模型完全相同。

1114励磁系统自定义模型框图如图5-34所示。

图5-331114型励磁系统数学模型

图5-341114型励磁系统自定义模型框图

132低频率、低电压减负荷继电器模型

(1)De_Loadbyf(编号:6121)

功能：低频率减负荷继电器

模型原理图:

模型图:

(2)De_LoadbyV(编号:6122)

功能：低电压减负荷继电器

模型原理图：

10111213

133直流输电控制器模型

HVDC_Mode13（编号：6203）

功能：模拟直流输电线路动态，两侧换流器准稳态和控制器（整流侧功率调节器,

逆变侧熄弧角调节器）动态的模型，适用于电力系统机电暂态稳定分析和小干扰稳定

分析。

模型示意图:

VDIOIDOGOAOXCIXCJTNITNJ

(U(L)(L)(L)(L)(L)(L)(L)

VTIR(BI)

》1TR(B1)

VTH(BI)9

VTIR(B2)

直流线路及其控»rri(Ri)

VT1KB2)

制系统暂态稳定

VT(Bl)>ITR(B2)

自定义模型

VT(B2)

>ITKB2)

TM7(BI)

RCDLILsQCIQCJ

(D(L)(L)

图中：L直流输电线路

Bl直流线整流侧交流母线

B2直流线逆变侧交流母线

VDIO整流侧直流电压潮流稳态值

ID0直流电流潮流稳态值

G0逆变侧熄弧角潮流稳态值

A0整流侧点燃角潮流稳态值

QCI,QCJ两侧换流站无功补偿容量

XCI,XCJ两侧换相电抗

TNLTNJ两侧变压器变比

RCD直流线路电阻

LI,Ls直流线路和平波电抗器电感

VT1R,VT1I两侧交流母线电压

1TR,ITI两侧母线注入电流

TM7外接系统稳定信号的入口

［说明］

①可以加电力系统稳定器的UD模型参与计算，ID模型号为6xxx,其输出定

义为TM7(B1),Bl为整流侧交流母线。

②模型参数

直流线数据

线路和换流站参数直接取直流线数据；

线路和换流站运行参数由UD模型输入信息得到。

调节器数据

调节器数据取公用参数库直流线调节器参数：

整流侧母线(B1)按整流侧调节器参数组号取值；

逆变侧母线(B2)按逆变侧熄弧角调节器参数组号取值。

③该模型不具备暂态稳定直流故障的能力及其他调节功能。

1.4暂态稳定的计算方法

如上所述,暂态稳定的数学模型可归结为网络方程和微分方程联立求解,即：

X=F(X,Y)

［y=G(x,y)(5-3)

其具体的处理方法是:采用梯形稳积分的迭代法,求解微分方程；采用直接三角分

解和迭代相结合的方法求解网络方程；微分方程和网络方程两者交替迭代,直至收敛,

以完成一个时段t的求解。

(1)微分方程的梯形隐积分迭代法

微分方程

y=G(x,y)(5-4)

的求解方法原理，与下面的单变量微分方程式的求解方法是一致的。

设微分方程式

dt(5-5)

当'〃处函数值Yn已知时，可按下式求出。+产/，，+A处的函数值Yn+1:

%=%+0'/()⑺力

(5-6)

上式中的定积分相当于下图中阴影部分的面积。

图5-35梯形积分法的几何解释

当步长At足够小时，函数f(y,t)在tn到tn+1之间的曲线可以近似地用直线代

替，如图中虚线所示。这样,阴影部分的面积就可以用梯形ABCD的面积来代替，因此，

式(5-6)可以改写为：

心=匕+白[/(匕乙)+/(匕”小)]停7、

乙in-/；

即是梯形积分法的差分方程。也就是把微分方程转换成代数方程求解。由于式(5-

7)等号的右侧也含有待求量YnH,这种隐式形式很难直接求解,通常采用如下的迭代

方法：

针)=%+等/50+/(燔4)

(5-8)

其中K为迭代次数，并设窗=儿这样，按式(4-9),由比;求染，再力瞪求

叶之依此类推，直至

,屋+1)V{K}

|Zn+I-Z«+l

时，即求得n+1时段的值

Y一y(K+i)

/n+l—1n+l

式(5-8),即是梯形隐积分的迭代方程式,可以根据函数f具体表达式对式(5-8)

进行整理,使之更有利于收敛。

为了简化叙述,现设暂态稳定的梯形隐积分方程如下：

yd'+D=G(x,yg)(5-