本发明涉及一种电路网络状态方程的迭代求解方法，具体涉及一种加速电路网络状态方程迭代求解的方法。
背景技术：
目前，在电力系统分析软件中一般采用lu分解等直接解法对电路网络方程进行求解。直接解法一般要求生成满阵，不能充分利用电路网络方程中系数矩阵的稀疏特性，会占用较大的内存空间。直接解法可以有效求解一千阶以下的电路方程，而对于拥有上万个节点的大规模电路网络则不具有任何优势。
共轭梯度法(以及带有预调件的共轭梯度法)是一种成熟的迭代求解稀疏矩阵方程的算法，但其要求待解矩阵方程中的系数矩阵具有正定性和对称性。使用状态方程法对电路网络进行分析时，其迭代求解式中的系数矩阵并不满足正定性和对称性，而无法直接使用共轭梯度算法。
技术实现要素：
为解决上述现有技术中的不足，本发明提供一种加速电路网络状态方程迭代求解的方法，该方法通过变换矩阵方程的形式，使待解方程的系数矩阵满足正定性、稀疏性和对称性，从而可以调用共轭梯度等成熟的快速迭代算法，以达到求解加速的目的。
本发明的目的是采用下述技术方案实现的：
本发明提供一种加速电路网络状态方程迭代求解的方法，其改进之处在于，所述方法包括下述步骤：
(1)列写电路的系统状态方程，并通过变量代换，将其变换一阶形式的状态方程；
(2)利用直接积分法对一阶状态方程进行迭代求解；
(3)在直接积分法迭代过程中调用共轭梯度法或带预条件的共轭梯度算法实现加速求解。
进一步地，所述步骤(1)中，电路的系统状态方程表达式如下：
其中：
通过变量代换，将系统状态方程①变换为一阶状态方程，如下：
其中：
k1、k2均表示系数矩阵；x为状态变量；e是单位对角矩阵。
进一步地，所述步骤(2)中，对于标准形式的一阶动态方程②，用直接积分法进行迭代求解，直接积分法的迭代格式如下式④所示：
其中：δt是第n个时刻与n+1时刻之间的时间步长；β是一个可供选择的系数，当β取值为0时，称为前差法(或前向欧拉法)，稳定性为条件稳定，数值精度为一阶精度；当β取值为0.5时，称为crank-nicolson法(或梯形法)，稳定性为无条件稳定，数值精度为二阶精度；当β取值为1时，称为后差法(或后向欧拉法)，稳定性为无条件稳定，数值精度为一阶精度；xn代表了第n步，即t＝nδt时刻的系统状态变量；xn+1代表了第n+1步，即t＝(n+1)δt时刻的系统状态变量；rn和rn+1分别代表了第n步，即t＝nδt时刻和第n+1步，即t＝(n+1)δt时刻的系统输入向量。
进一步地，所述步骤(3)包括下述步骤：
1)引入中间变量a和b，将迭代式④拆解为式⑤和⑥；
2)利用式⑤计算中间变量a和b；
3)式⑥拆分成式⑦和⑧，再将式⑧代入⑦中消去其中的变量ψn+1，得到式⑨；
4)调用预调件的共轭梯度算法(共轭梯度法是一种常见的矩阵方程求解方法)对式⑨中的矩阵逆乘向量运算进行求解，得到状态变量
5)利用式⑧获得状态变量ψn+1，最终得到状态变量
进一步地，所述步骤1)中，将式③代入式④中，并按照计算顺序将迭代式④在等号左右两边分解成两个部分，如式⑤和⑥，其中的a和b是推导过程引入的中间变量，均为n阶列向量；
其中：kc是电容系数矩阵；kr是电阻系数矩阵；kl是电感系数矩阵；e是单位对角矩阵；式⑤是利用前一时刻的状态量
进一步地，所述步骤3)包括：
其中：kc是电容系数矩阵；kr是电阻系数矩阵；kl是电感系数矩阵；e是单位对角矩阵；a和b均为中间变量；由于kc、kl和kr都满足正定性和对称性，且β和δt均大于零，所以
与最接近的现有技术相比，本发明提供的技术方案具有的优异效果是：
1.本发明通过变换直接积分法迭代格式的求解步骤，使求解过程中的矩阵逆乘向量运算满足了共轭梯度迭代求解算法所要求的矩阵特性，使求解过程可以调用现已成熟的、高速的 共轭梯度算法来完成求解加速。
2.本发明在求解过程中保持了状态方程稀疏性的特点，可以使用成熟的稀疏存储和求解技术，可以在有限的内存空间内实现大规模节点数的电路方程仿真。
附图说明
图1是本发明提供的加速电路网络状态方程迭代求解的方法的流程图。
具体实施方式
下面结合附图对本发明的具体实施方式作进一步的详细说明。
以下描述和附图充分地示出本发明的具体实施方案，以使本领域的技术人员能够实践它们。其他实施方案可以包括结构的、逻辑的、电气的、过程的以及其他的改变。实施例仅代表可能的变化。除非明确要求，否则单独的组件和功能是可选的，并且操作的顺序可以变化。一些实施方案的部分和特征可以被包括在或替换其他实施方案的部分和特征。本发明的实施方案的范围包括权利要求书的整个范围，以及权利要求书的所有可获得的等同物。在本文中，本发明的这些实施方案可以被单独地或总地用术语“发明”来表示，这仅仅是为了方便，并且如果事实上公开了超过一个的发明，不是要自动地限制该应用的范围为任何单个发明或发明构思。
本发明提供一种加速电路网络状态方程迭代求解的方法，其流程图如图1所示，包括下述步骤：
(1)一般线性电路(由电阻、电容、电感以及电源等基本元件组成的，具有n+1个节点的电路，其中0号节点定义为参考地电位节点)列写电路的系统状态方程，并通过变量代换，将其变为一阶形式的状态方程；
电路的系统状态方程表达式如下：
其中：
一阶状态方程表达式如下：
其中：
x为状态变量；e是单位对角矩阵。
(2)利用直接积分法对一阶状态方程进行迭代求解；对于标准形式的一阶动态方程②，用直接积分法进行迭代求解，直接积分法的迭代格式如下式④所示：
其中：δt是第n个时刻与n+1时刻之间的时间步长；β是一个可供选择的系数，当β取不同的值时，相关的名称如下表所述：
表1β的选择与数值稳定性
xn和xn+1分别是t＝nδt和t＝(n+1)δt时刻的状态变量；rn和rn+1分别是t＝nδt和t＝(n+1)δt时刻的系统输入向量。
将式③代入式④中，并按照计算顺序将迭代式④在等号左右两边分解成两个部分，如式⑤和⑥，其中的a和b是推导过程引入的中间变量，均为n阶列向量；
其中：kc是电容系数矩阵；kr是电阻系数矩阵；kl是电感系数矩阵；e是单位对角矩阵；式⑤是利用前一时刻的状态量
如果能够加速式⑥的求解速度，也就加速了整个迭代的求解过程。本专利即给出一种加速式④迭代运算的方法。
(3)在直接积分法迭代过程中调用共轭梯度法或带预条件的共轭梯度算法实现加速求解，对于满足稀疏性的高阶正定对称系数矩阵方程，一般采用已经成熟的共轭梯度等算法来加速计算。迭代式④中的各系数矩阵满足稀疏性，却不满足正定性和对称性。为了加速迭代式④的计算，采用如下步骤：
1)引入中间变量a和b，将迭代式④拆解为式⑤和⑥；
2)利用式⑤计算中间变量a和b；
3)式⑥拆分成式⑦和⑧，再将式⑧代入⑦中消去其中的变量ψn+1，得到式⑨；
其中：kc是电容系数矩阵；kr是电阻系数矩阵；kl是电感系数矩阵；e是单位对角矩阵；a和b均为中间变量；由于kc、kl和kr都满足正定性和对称性，且β和δt均大于零，所以
4)调用预调件的共轭梯度算法对式⑨中的矩阵逆乘向量运算进行求解，得到状态变量
5)利用式⑧获得状态变量ψn+1，最终得到状态变量
本发明提供的加速电路网络状态方程迭代求解的方法，采用直接积分法迭代式分步求解的方法，即：通过方程变换获得满足正定对称特性的系数矩阵的方法；在直接积分法迭代过 程中调用共轭梯度法(或带预条件的共轭梯度算法)等快速稀疏矩阵算法实现求解加速的方法。
以上实施例仅用以说明本发明的技术方案而非对其限制，尽管参照上述实施例对本发明进行了详细的说明，所属领域的普通技术人员依然可以对本发明的具体实施方式进行修改或者等同替换，这些未脱离本发明精神和范围的任何修改或者等同替换，均在申请待批的本发明的权利要求保护范围之内。