CN106125715A - 一种基于阻抗回比矩阵的mimo级联系统稳定性分析方法 - Google Patents
一种基于阻抗回比矩阵的mimo级联系统稳定性分析方法 Download PDFInfo
- Publication number
- CN106125715A CN106125715A CN201610454774.0A CN201610454774A CN106125715A CN 106125715 A CN106125715 A CN 106125715A CN 201610454774 A CN201610454774 A CN 201610454774A CN 106125715 A CN106125715 A CN 106125715A
- Authority
- CN
- China
- Prior art keywords
- mimo
- return ratio
- impedance
- ratio matrix
- stability
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B23/00—Testing or monitoring of control systems or parts thereof
- G05B23/02—Electric testing or monitoring
- G05B23/0205—Electric testing or monitoring by means of a monitoring system capable of detecting and responding to faults
- G05B23/0218—Electric testing or monitoring by means of a monitoring system capable of detecting and responding to faults characterised by the fault detection method dealing with either existing or incipient faults
- G05B23/0243—Electric testing or monitoring by means of a monitoring system capable of detecting and responding to faults characterised by the fault detection method dealing with either existing or incipient faults model based detection method, e.g. first-principles knowledge model
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B2219/00—Program-control systems
- G05B2219/20—Pc systems
- G05B2219/24—Pc safety
- G05B2219/24065—Real time diagnostics
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Automation & Control Theory (AREA)
- Measurement Of Resistance Or Impedance (AREA)
Abstract
本发明公开了一种基于阻抗回比矩阵的MIMO级联系统稳定性分析方法,以高铁车网级联系统为对象,包括推导计算系统的等效前级输出阻抗,推导计算系统的等效后级输入导纳,计算系统整体传递函数并获得回比矩阵,估计回比矩阵特征值并限制其分布区域,特征值分布区域设置禁区,获得降低保守性的MIMO系统稳定性判据。推导计算系统的等效前级输出阻抗和后级输入导纳是获得系统整体传递函数的前提,是进一步获得级联系统回比矩阵的必要条件;盖尔圆定理是设置禁区限制回避矩阵特征值范围,进一步降低判据保守性的基础。该判据方法适用于MIMO级联系统的稳定性分析,可以简单有效地分析MIMO级联系统稳定性,并且保守性相较于现有的奇异值判据、范数判据等都更小。
Description
技术领域
本发明涉及系统稳定性分析是自控控制领域,尤其是基于阻抗回比矩阵的多输入多输出(MIMO)级联系统稳定性分析方法。
背景技术
实际生产应用中存在各种类型的自动控制系统,但其只有在满足稳定性条件下才能正常工作,因此系统稳定性分析是自控控制领域一个至关重要的问题。
经典控制理论中的方法,如Nyquist判据、Bode图分析等,可以对单输入单输出(SISO)系统稳定性进行分析。Middlebrook对于级联SISO系统提出了基于阻抗回比函数的稳定性判据,随后不少学者基于禁区的概念对Middlebrook阻抗判据进行了推广,减小了判据的保守性。而实际生产应用中的系统大多为多输入多输出(MIMO)系统,MIMO系统稳定性分析往往通过线性系统理论的方法,需要对系统状态空间模型进行精确建模。Belkhayat等人将SISO级联系统阻抗稳定性分析方法推广于MIMO级联系统,提出奇异值判据、范数判据等方法,相较于线性系统理论的方法更为简单,但是这类判据保守性较大。
本发明针对上述问题,提出了基于阻抗回比矩阵的MIMO级联系统稳定性分析方法,缩小现有判据的禁区范围。该判据可以简单有效地分析MIMO级联系统稳定性,并且保守性相较于现有的奇异值判据、范数判据等都更小。
发明内容
本发明的目的是提出一种基于阻抗回比矩阵的MIMO级联系统稳定性分析方法。该方法通过分析MIMO级联系统的前级输出阻抗及后级输入导纳,计算阻抗回比矩阵,基于广义Nyquist判据进行车网级联系统稳定性分析。该方法将传统单输入单输出系统阻抗稳定性分析方法拓展于多输入多输出系统,得出一种简单的稳定性判据,具有一定保守性。
本发明的目的是通过如下的手段实现的。
一种基于阻抗回比矩阵的MIMO级联系统稳定性分析方法,包括阻抗回比矩阵计算和系统稳定性分析两部分;通过分析MIMO级联系统的前级输出阻抗及后级输入导纳,计算阻抗回比矩阵,基于广义Nyquist判据进行车网级联系统稳定性分析;包含如下的主要步骤:
1)、获取计算前级系统的等效输出阻抗:
在两相旋转坐标系(dq系)下,MIMO级联系统的前级系统复频域下输出阻抗矩阵:
其中,Rs和Ls为牵引网等效电阻值和电感值,ω为基波角频率;
2)、计算后级系统的等效输入导纳:
对于CRH5型车牵引变流器,动车组在dq坐标下的导纳矩阵可以表示为
YLdq=2maYin=nYin
其中m为动车组台数,a为单台动车组牵引动力单元个数,Yin为PWM变流器的输入导纳矩阵;
3)、计算系统整体传递函数并获得回比矩阵:
式中:Ldq为系统阻抗回比矩阵,其元素分别为Ldd、Ldq、Lqd和Lqq。ZSdq为前级系统等效输出阻抗,其元素分别为Zdd、Zdq、Zqd和Zqq。YLdq为后级系统等效输入导纳,其元素分别为Ydd、Ydq、Yqd和Yqq,s表示复频域;
4)、估计阻抗回比矩阵特征值并限制其分布区域:
利用盖尔圆定理限制特征值的分布区域,特征值满足以下两式
|λi-Ldd|<|Ldq|或|λi-Lqq|<|Lqd|
|λi-Ldd|<|Lqd|或|λi-Lqq|<|Ldq|
式中:λi为系统阻抗回比矩阵Ldq的第i个特征值,Ldd、Ldq、Lqd和Lqq分别为Ldq的元素;
5)、特征值分布区域设置禁区,获得低保守性的MIMO级联系统稳定性判据:
相应的四个子判据为
Re{Ldd}-|Ldq|+1>0
Re{Lqq}-|Lqd|+1>0
Re{Ldd}-|Lqd|+1>0
Re{Lqq}-|Ldq|+1>0
式中,Ldd、Ldq、Lqd和Lqq分别为系统阻抗回比矩阵Ldq的四个元素,Re表示元素的实部。
与现有技术相比,本发明的有益效果是:
1、本发明将传统SISO系统阻抗稳定性分析方法拓展于MIMO系统,克服了传统线性系统理论的方法应用于MIMO系统时需要对系统状态空间模型进行精确建模的问题,实现过程更加简便。
2、本发明进一步缩小判据的禁区,相对于现有的用于MIMO系统的奇异值判据、范数判据等方法保守性更小。
附图说明
图1为本发明的全并联复线AT牵引网等效电路。
图2为本发明的CRH5型车单个牵引动力单元示意图。
图3为本发明的等效电路图。
图4为本发明的CRH5型车牵引变流器dq电流控制框图。
图5为本发明的车网级联系统拓扑图。
图6为本发明的基于禁区的MIMO稳定性判据示意图。
图7为本发明7台车接入牵引网的稳定性分析幅频图。
图8为本发明8台车接入牵引网的稳定性分析幅频图。
图9为本发明的多车接入牵引网低频振荡的仿真再现图。
具体实施方式
下面结合附图和实例对本发明的实施例做详细说明。
本实施例在以本发明技术方案为前提下进行实施,给出了详细的实施过程,但本发明的保护范围不限于下述的实施例。
本实施例以高速铁路中CRH5型动车组接入的车网级联系统为例。
包括前级牵引网侧系统等效源阻抗的计算,后级动车组侧系统等效负荷导纳计算,车网级联系统的传递函数的计算并得出回比矩阵,估计回比矩阵特征值并限制其分布区域,设置禁区并获得低保守性的系统稳定性分析结果。
本实施例的MIMO级联系统稳定性分析由以下各步组成。
1)、计算前级系统的等效输出阻抗;
前级系统的等效输出阻抗通过推导牵引网等效阻抗获得,此处对中国高铁较常用的全并联复线AT供电系统进行建模和计算。图1为全并联复线AT牵引网等效电路。基于广义对称分量法,全并联复线AT牵引网的单位阻抗变换可以推导为:
式中,CW和AF分别表示接触线和AT馈线,U和D分别表示上行和下行,A4为
根据图1的等效电路,全并联复线AT牵引网的等效阻抗可以推导为
其中,Z’AT为AT漏抗的两倍,x是列车与前方相邻AT站的距离,l是列车与牵引变电站之间的距离,DAT是相邻AT站的距离。
将前级系统输出阻抗Zs分解到两相旋转坐标系(dq系)下,得到MIMO级联系统的前级系统复频域下输出阻抗矩阵:
其中,Rs和Ls为牵引网等效电阻值和电感值,ω为基波角频率。
2)、计算后级系统的等效输入导纳;
后级系统的输入导纳可以通过推导列车的数学模型获得,CRH5型车单个牵引动力单元如图2所示。逆变器和电机部分对稳定性分析影响微弱,故可等效为一个电阻以简化建模,且由于两重化整流器结构一致,对单个变流器的直流侧等效电容和负荷分别为原来的1/2和2倍,因此简化后CRH5型车牵引变流器等效电路如图3所示。
假设da和db分别表示a桥和b桥的开关状态,则电压和电流的关系可以表示为
因为采用了dq电流控制,输入电压un和输入电流in经dq解耦后得到ud,uq和id,iq。un和in被视为α分量,通过下式坐标变换
vα=vd cosωt-vq sinωt
可将电压电流分解到dq系下:
在静态工作点,上式可以进一步表示为
由于Eq和Iq均为零,Id,Dd和Dq稳态值可以表示为下式
为了推导整流器的小信号模型,对每个变量的稳态值加上一个小信号值
整理可得到矩阵形式的整流器小信号模型
CRH5型车网侧整流器常规控制方法为dq电流控制,如图4,可以得到下式
其中,iq*=0,HI为电流环PI控制器传递函数。id*可由电压环PI控制器输出得到,为
其中HU为电压环PI控制器传递函数。
图4的SPWM能够映射到dq结构下,假设载波幅值为UC,则[dd,dq]T和[ud,uq]T的关系为
进一步可得到
将上式代入小信号模型,并通过Laplace变换可得
则PWM变流器的输入导纳矩阵被表示为
动车组在dq坐标下的导纳矩阵可以表示为
YLdq=2maYin=nYin
其中m为动车组台数,a为单台动车组牵引动力单元个数。
3)、计算系统整体传递函数并获得回比矩阵;
车网系统可以被考虑为一个级联系统,整个系统的传递函数能够通过网侧输出阻抗和车侧输入导纳推导得到。根据图5由KVL电压定律可得
通过Laplace变换可得到车网级联系统的传递函数矩阵为
Gcl=YLdq(I2+ZSdqYLdq)-1=YLdq(I2+Ldq)-1
则回比矩阵定义为
4)、估计阻抗回比矩阵特征值并限制其分布区域;
利用盖尔圆定理限制特征值的分布区域,特征值位于以对角元素为圆心,以同一行的非对角元素和为半径的区域,因此Ldq(s)的特征值λi位于图6所示的盖尔圆G1或者G2区域中,满足下式
|λi-Ldd|<|Ldq|或|λi-Lqq|<|Lqd|
由于Ldq的所有特征值在Ldq T的n个盖尔圆内,即Ldq的列盖尔圆,因此特征值同样满足下式
|λi-Ldd|<|Lqd|或|λi-Lqq|<|Ldq|
5)、特征值分布区域设置禁区,获得低保守性的MIMO级联系统稳定性判据;
为了降低保守性,一个基于禁区的判据被推广到MIMO系统,通过限制Ldq(s)的特征值在(-1+j0)右侧,使得绕过(-1+j0)的Nyquist围线仍旧为0,保持系统稳定,如图6所示。
对于Ldq盖尔圆对特征值的限制,响应的系统稳定性充分条件为
对于Ldq列盖尔圆对特征值的限制,相应的系统稳定性充分条件为
结合以上两式,可以得到低保守性的MIMO级联系统稳定性判据为
相应的四个子判据为
Re{Ldd}-|Ldq|+1>0
Re{Lqq}-|Lqd|+1>0
Re{Ldd}-|Lqd|+1>0
Re{Lqq}-|Ldq|+1>0
计算参数如表I所示
表I
车网系统参数
提出的判据用幅频图表示,四条曲线分别表示四个子判据,当两条红线和蓝线均大于零时系统稳定,由图7和图8可得,7台车接入牵引网时系统稳定,当接入8台时系统不稳定,结合图9低频振荡现象发生的情况来看,判据能够准确确定车网级联系统的稳定性情况。
Claims (4)
1.一种基于阻抗回比矩阵的MIMO级联系统稳定性分析方法,包括阻抗回比矩阵计算和系统稳定性分析两部分;通过分析MIMO级联系统的前级输出阻抗及后级输入导纳,计算阻抗回比矩阵,基于广义Nyquist判据进行车网级联系统稳定性分析;包含如下的主要步骤:
1)、获取计算前级系统的等效输出阻抗:
在两相旋转坐标系(dq系)下,MIMO级联系统的前级系统复频域下输出阻抗矩阵:
其中,Rs和Ls为牵引网等效电阻值和电感值,ω为基波角频率;
2)、计算后级系统的等效输入导纳:
对于CRH5型车牵引变流器,动车组在dq坐标下的导纳矩阵可以表示为
YLdq=2maYin=nYin
其中m为动车组台数,a为单台动车组牵引动力单元个数,Yin为PWM变流器的输入导纳矩阵;
3)、计算系统整体传递函数并获得回比矩阵:
式中:Ldq为系统阻抗回比矩阵,其元素分别为Ldd、Ldq、Lqd和Lqq。ZSdq为前级系统等效输出阻抗,其元素分别为Zdd、Zdq、Zqd和Zqq。YLdq为后级系统等效输入导纳,其元素分别为Ydd、Ydq、Yqd和Yqq,s表示复频域;
4)、估计阻抗回比矩阵特征值并限制其分布区域:
利用盖尔圆定理限制特征值的分布区域,特征值满足以下两式
|λi-Ldd|<|Ldq|或|λi-Lqq|<|Lqd|
|λi-Ldd|<|Lqd|或|λi-Lqq|<|Ldq|
式中:λi为系统阻抗回比矩阵Ldq的第i个特征值,Ldd、Ldq、Lqd和Lqq分别为Ldq的元素;
5)、特征值分布区域设置禁区,获得低保守性的MIMO级联系统稳定性判据:
相应的四个子判据为
Re{Ldd}-|Ldq|+1>0
Re{Lqq}-|Lqd|+1>0
Re{Ldd}-|Lqd|+1>0
Re{Lqq}-|Ldq|+1>0
式中,Ldd、Ldq、Lqd和Lqq分别为系统阻抗回比矩阵Ldq的四个元素,Re表示元素的实部。
2.根据权利要求1所述的基于阻抗回比矩阵的MIMO级联系统稳定性分析方法,其特征在于,所述步骤5)中,该稳定性判据应用于MIMO级联系统稳定性分析。
3.根据权利要求1所述的基于阻抗回比矩阵的MIMO级联系统稳定性分析方法,其特征在于,在所述步骤3)、4)、5)中,该稳定性判据可以通过阻抗回比矩阵进行稳定性分析。
4.根据权利要求1所述的基于阻抗回比矩阵的MIMO级联系统稳定性分析方法,其特征在于,所述步骤5)中通过设置禁区减少了该稳定性判据的保守性。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610454774.0A CN106125715B (zh) | 2016-06-21 | 2016-06-21 | 一种基于阻抗回比矩阵的mimo级联系统稳定性分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610454774.0A CN106125715B (zh) | 2016-06-21 | 2016-06-21 | 一种基于阻抗回比矩阵的mimo级联系统稳定性分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106125715A true CN106125715A (zh) | 2016-11-16 |
CN106125715B CN106125715B (zh) | 2018-06-19 |
Family
ID=57471337
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610454774.0A Expired - Fee Related CN106125715B (zh) | 2016-06-21 | 2016-06-21 | 一种基于阻抗回比矩阵的mimo级联系统稳定性分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106125715B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107565545A (zh) * | 2017-07-28 | 2018-01-09 | 西南交通大学 | 一种基于阻抗回比矩阵的牵引车‑网耦合系统稳定性判据计算方法 |
CN108847670A (zh) * | 2018-07-30 | 2018-11-20 | 西南交通大学 | 一种双馈风机网侧变换器的谐波不稳定分析方法 |
CN108983635A (zh) * | 2018-06-19 | 2018-12-11 | 兰州交通大学 | 一种crh5型动车组投入多车网电气耦合系统稳定性分析方法 |
CN112636380A (zh) * | 2020-12-21 | 2021-04-09 | 中国科学院电工研究所 | 一种基于盖尔圆理论的交直流配用电系统特征值分析方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060271232A1 (en) * | 2005-05-31 | 2006-11-30 | Hitachi Via Mechanics, Ltd. | Digital servo control unit and laser machining apparatus |
US20100292813A1 (en) * | 2009-05-18 | 2010-11-18 | Imb Controls Inc. | Method and Apparatus for Tuning a PID Controller |
CN102176118A (zh) * | 2011-01-27 | 2011-09-07 | 西北工业大学 | 多输入-多输出飞行器等价稳定裕度鲁棒确定方法 |
CN105223810A (zh) * | 2015-09-06 | 2016-01-06 | 西南交通大学 | 一种基于多变量控制的高铁低频振荡过电压阻尼方法 |
-
2016
- 2016-06-21 CN CN201610454774.0A patent/CN106125715B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060271232A1 (en) * | 2005-05-31 | 2006-11-30 | Hitachi Via Mechanics, Ltd. | Digital servo control unit and laser machining apparatus |
US20100292813A1 (en) * | 2009-05-18 | 2010-11-18 | Imb Controls Inc. | Method and Apparatus for Tuning a PID Controller |
CN102176118A (zh) * | 2011-01-27 | 2011-09-07 | 西北工业大学 | 多输入-多输出飞行器等价稳定裕度鲁棒确定方法 |
CN105223810A (zh) * | 2015-09-06 | 2016-01-06 | 西南交通大学 | 一种基于多变量控制的高铁低频振荡过电压阻尼方法 |
Non-Patent Citations (3)
Title |
---|
孙青林等: "基于频域的多变量广义预测控制(MIMO-GPC)稳定性分析", 《中国工程科学》 * |
曹玉腾等: "再入飞行器多变量稳定裕度研究", 《航天控制》 * |
顾树生等: "一种基于奇异值分析的多变量系统设计方法", 《控制与决策》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107565545A (zh) * | 2017-07-28 | 2018-01-09 | 西南交通大学 | 一种基于阻抗回比矩阵的牵引车‑网耦合系统稳定性判据计算方法 |
CN108983635A (zh) * | 2018-06-19 | 2018-12-11 | 兰州交通大学 | 一种crh5型动车组投入多车网电气耦合系统稳定性分析方法 |
CN108983635B (zh) * | 2018-06-19 | 2022-05-13 | 兰州交通大学 | 一种crh5型动车组投入多车网电气耦合系统稳定性分析方法 |
CN108847670A (zh) * | 2018-07-30 | 2018-11-20 | 西南交通大学 | 一种双馈风机网侧变换器的谐波不稳定分析方法 |
CN108847670B (zh) * | 2018-07-30 | 2021-04-09 | 西南交通大学 | 一种双馈风机网侧变换器的谐波不稳定分析方法 |
CN112636380A (zh) * | 2020-12-21 | 2021-04-09 | 中国科学院电工研究所 | 一种基于盖尔圆理论的交直流配用电系统特征值分析方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106125715B (zh) | 2018-06-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106125715A (zh) | 一种基于阻抗回比矩阵的mimo级联系统稳定性分析方法 | |
CN107565545B (zh) | 一种基于阻抗回比矩阵的牵引车-网耦合系统稳定性判据计算方法 | |
CN101719182B (zh) | 一种交直流电力系统分割并行电磁暂态数字仿真方法 | |
CN106410848A (zh) | 一种电力电子多馈入电力系统小干扰稳定性评估方法 | |
CN104393590B (zh) | 电网暂态失稳模式智能辨识方法 | |
CN105095566A (zh) | 一种基于小波分析和svm的逆变器故障诊断方法 | |
CN112149280B (zh) | 含svg的新能源的多馈入系统电网强度获得方法 | |
Nagarajan et al. | Application of minimum spanning tree algorithm for network reduction of distribution systems | |
CN103226636B (zh) | 基于定性推理的放大电路故障传播分析方法 | |
CN105071447A (zh) | 一种交直流混联多馈入直流受端电网强度获得方法 | |
CN104578153A (zh) | 一种电网拓扑结构识别方法 | |
CN102104252A (zh) | 一种适用于电磁暂态分析的电力系统动态等值方法 | |
CN106682407A (zh) | 基于戴维南等值和支路传输功率极限的电压稳定评估方法 | |
CN102163845B (zh) | 基于功率矩法的分布式电源优化配置方法 | |
CN106208099A (zh) | 一种基于二层规划的电力系统无功优化方法及其应用 | |
CN106097156A (zh) | 配电网线损计算方法和系统 | |
Lin et al. | Analysis of influencing factors of short circuit ratio of multi-infeed AC/DC power systems | |
CN106096765A (zh) | 分散式风电群输电规划方案的评价方法 | |
CN106877309B (zh) | 主从控制模式下直流微电网的稳定性判定方法 | |
CN102914718A (zh) | 基于响应成分和振荡特征辨识的低频振荡类型判别方法 | |
CN105469316A (zh) | 配电网任意两点间理论线损计算方法及系统 | |
CN103427433B (zh) | 一种混合双极直流输电系统直流回路阻抗的计算方法 | |
CN105119279A (zh) | 一种分布式电源规划方法及其系统 | |
CN102184296B (zh) | 基于实测数据的电气化铁路冲击负荷建模方法 | |
CN105044517B (zh) | 基于潮流特征值的智能变电站保护向量诊断方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180619 Termination date: 20210621 |