电力系统分析matlab仿真_电力系统稳定性分析

电力系统分析matlab仿真_电力系统稳定性分析基于Wirtinger不等式的时滞电力系统稳定性判定方法【专利摘要】本发明公开了一种基于Wirtinger不等式的时滞电力系统稳定性判定方法,用于分析电力系统所能承受的最大时滞稳定裕度。该方法的具体步骤如下:首先,建立考虑时滞影响的电力系统模型。然后,针对所建模型构建Lyapunov泛函,在泛函的求导过程中通过采用Wirtinger不等式进行放缩,以减少判据的保守性。最后将所得判据用一组线性矩阵不…

大家好,又见面了,我是你们的朋友全栈君。如果您正在找激活码,请点击查看最新教程,关注关注公众号 “全栈程序员社区” 获取激活教程,可能之前旧版本教程已经失效.最新Idea2022.1教程亲测有效,一键激活。

Jetbrains全系列IDE稳定放心使用

基于Wirtinger不等式的时滞电力系统稳定性判定方法

【专利摘要】本发明公开了一种基于Wirtinger不等式的时滞电力系统稳定性判定方法,用于分析电力系统所能承受的最大时滞稳定裕度。该方法的具体步骤如下:首先,建立考虑时滞影响的电力系统模型。然后,针对所建模型构建Lyapunov泛函,在泛函的求导过程中通过采用Wirtinger不等式进行放缩,以减少判据的保守性。最后将所得判据用一组线性矩阵不等式(LMI)表示。

【专利说明】

基于W i rt i nger不等式的时滞电力系统稳定性判定方法

技术领域

[0001 ]本发明涉及一种基于Wirtinger不等式的时滞电力系统稳定性判定方法,适用于 解决互联电力系统广域控制策略中的延时问题,属于电力系统技术领域。

【背景技术】

[0002] 随着电网规模的不断扩大,单纯的依靠本地控制信息已经无法满足电网性能的要 求。因此电网稳定控制策略逐步从本地控制趋向于全局控制,特别是近年来广域测量系统 (WAMS)的快速发展,使电力系统全局控制策略的实现成为可能。

[0003] 在此基础上,现代电力系统控制器将从全局角度出发,利用本地信号和远端信号 作为反馈信号设计广域控制器。由于远端信号的引入,信号的延时将变得不可避免,已有研 究表明,即使很小的时滞都可能对电力系统稳定性产生影响。因此考虑电力系统所能承受 的最大时滞,对于电网的安全稳定运行具有十分重要的意义。

[0004] 目前关于时滞系统的研究,具有多种求解方法,常采用的求解方法是构造 Lyapunov泛函,基于Lyapunov稳定性理论,得到系统稳定判据,最后借助线性矩阵不等式 (LMI)来求解时滞稳定裕度。但是由上述方法得到的判据,实际是对求导后的Lyapunov泛函 放缩所得到的充分条件,因此造成所得判据具有一定的保守性。所以降低保守性问题成为 该方法研究的重点,也是难点之一。

【发明内容】

[0005] 发明目的:针对现有技术中存在的问题,本发明提出一种基于Wirtinger不等式的 时滞电力系统稳定性判定方法,首先构造全新的Lyapunov泛函,将时滞下限不为零考虑进 判据中,然后利用Wirtinger不等式放缩技巧,降低结果的保守性。

[0006] 技术方案:一种基于Wirtinger不等式的时滞电力系统稳定性判定方法,包括如下 步骤: j x(t) = Ax(t) + A,)X{t – ci(t))

[0007] (1)建立包含广域控制回路的时滞电力系统模型],、,、’ \式中:

[x(j) = (p(j) /eH/{0,0] hi彡d(tXh2,办),y彡1。其中0彡hi

[0008] (2)给定稳定判定条件:

[0009] 若存在正定矩阵P G R4nX4n;正定矩阵Qi G RnXn,i = 1,2,3;正定矩阵Zj G RnXn,j = 1, 2;矩阵XkGR2nX2n,k=l,2使下列矩阵不等式成立,则时滞电力系统是渐近稳定的:

[0012] 其中:(1)。= //e((7/ /)(72) + Q + 4 + 4 + 匀允2

[0013] r:!-^ Gf Gl G[]r

[0014] G3 =ei_e2,G4=ei+e2_2es,Gs= e2_e4,G6= e2+e4_2e7

[0015] r2=[Gf Ggr 0i

[0016] Gy = e3_e2,G8 = e3+e2_2e6,G9 = e2_e4,Gio = e2+e4_2e7

[0017] Qj -diag^O.-Q^J> Q: = ^^g(〇2^°2,;’Q: °.i/,);

[0018] Q,,=diag{0,,-{\-d)0,,Q5tt) _ Z^h2 xdiagiZ^Q.J + Uh-l\)^diag{Z^0,:)

[0019] 之=c//tzg(Z丨,3Z|),之2 =r”ag(Z2,3Z2)

[0020] G^\e[ d(!)el {d{l)-h{)el {h.-dU^J 一 T

[0021 ] (?2 = A^C-, ~ (1 – cl)~ ~ cl)c!j -+ (1- d)C^

[0022] ei=[0i-i I O7-1], 1 = 1, , 7

[0023] 对于矩阵A,He(A)=A+AT。其中I代表单位矩阵。

[0024] (3)利用Matlab中的线性矩阵(LMI)工具箱判断给定时滞d(t)是否满足步骤(2)给 出的判定条件,若满足,则可判定在延时d(t)条件下的时滞电力系统是渐近稳定的。

[0025] 时滞电力系统模型乂(〇 =也的+ .^^_

【附图说明】

[0026] 图1为本发明所提的时滞电力系统稳定性判定方法流程图;

[0027]图2为包含WADC的时滞电力系统;

[0028]图3为四机两区域电力系统;

[0029]图4为d(t)=0ms情况下四机两区域电力系统响应图;

[0030]图5为d(t) = 110ms情况下四机两区域电力系统响应图。

【具体实施方式】

[0031]下面结合具体实施例,进一步阐明本发明,应理解这些实施例仅用于说明本发明 而不用于限制本发明的范围,在阅读了本发明之后,本领域技术人员对本发明的各种等价 形式的修改均落于本申请所附权利要求所限定的范围。

[0032] (1)电力系统时滞模型的建立

[0033]通常情况下,电力系统可由一组微分代数方程描述,在系统运行点附近对其线性 化,最终系统可表示为:

[i, (/) = A,x, (/) + _ KU ⑴ | j(,) = ChV|(0

[0035] 式中:xi为系统状态变量;Ai为系统状态矩阵;Bi为系统输入矩阵;Ci为系统输出矩 阵;u为系统控制输入,y为系统控制输出。

[0036] 系统的控制策略采用基于WAMS的广域阻尼控制,其中广域阻尼控制器(WADC)的控 制输入信号含有系统远端信号,其状态方程可表示为: f.v {t) = Ax U) + Bu (t)

[0037] ^ f ( (2) UV⑴=(、弋⑴+ Z). W?

[0038] 式中:x。为控制器的状态变量;A。为控制系统的状态矩阵;B。为控制系统的输入矩 阵;C。为控制系统的输出矩阵;u。为控制系统输入,y。为系控制系统输出。D。为标量,反映了 输出y。与输入u。之间的直接关联。

[0039] 图2给出了包含广域控制回路的时滞电力系统结构关系图,其中的变延时d(t)是 由于远端信号传输所引起的传输延时。根据图2所示的关系图,可以得到:

[0040] (3) |w(/) = V(.(〇

[0041 ]其中:hi彡d(tXh2,

[0042]进一步,时滞电力系统模型可表示为:

[0043] i(/) = Ax{() + Adx{t-d{t)) (4). T 「4 B.C 1 「5,1) C, 0]

[0044] 式中,4= 1 e 1 L 丨 c」 L〇 A J L sc\

[0045] (2)电力系统延时依赖稳定性判据方法

[0046]求解时滞稳定判据过程中,常用的放缩技巧是利用Jensen不等式。这种方法虽然 可行但增加了结果的保守性。本发明所提方法采用一种全新的不等式-Wirtinger不等式进 行放缩,可大大降低所得结果的保守性。首先,给出本发明所提方法用到的两个重要的引 理。

[0047]引理1:对于给定的正定矩阵M>0,以下不等式对于在区间[a,b]上连续可微函数x 都成立:

[0051 ]引理2:对于给定的正定矩阵R>0,矩阵I,W2和标量a G (〇,1),定义对于所有的I, 函数 0 (a,R):

[0053] 如果存在矩阵乂,使A X >0,那么以下不等式成立: -* ‘R XTW£~

[0054] min 0(a., R)> ” (仙 Lmi」L*及」L:i^i

[0055] 构造如下Lyapunov泛函: V ~ if Pq i- [ x1 (s)Q,x(s)cls + I x! (s)()2x(s)cis + [ x1 (s)Q,x(s)ds r n h-h h–h2

[0056] 〇 , _h t (5) + J [ i:1 (a0Z + | ^ j” x1 (s)Z lx(.s)dsci0

[0057] 式中:PGR4nX4n;QiGRnXn,i = l,2,3;ZjGRnXn,j = l,2,

[0058] = x1 (t) | ^xr(s)d$ I . xJ (syds

[0059] SeiGRM’eftOi-丨I 〇7-」,1 = 1,…,7为分块坐标矩阵,可得如下稳定性判据:

[0060] 判据:若存在正定矩阵P G R4nX4n;正定矩阵Qi G RnXn,i = 1,2,3 ;正定矩阵Zj G RnXn, j = l,2;矩阵XkGR2nX2n,k=l,2使下列矩阵不等式成立,则时滞电力系统是渐近稳定的:

(6)

[0063] 其中:c丨)(1 = /”,((,丫/%) + Q + 4 + 4 + (7:丨f

[0064] ri: = [Gf GTA Gl Gl^

[0065] 〇3 =ei_e2,G4=ei+e2_2es,Gs= e2_e4,G6= e2+e4_2e7

[0066] r2: = [(3f Gj G^J

[0067] 〇7 = e3_e2,G8 = e3+e2_2e6,G9 = e2_e4,Gio = e2+e4_2e7

[0068] Q = , 4 =t/Vag(a,02;,,

[0069] (), -^K2′,^°5″) 7 Z ^h? y.diag(Zi,0,ii) + (h2-h^xdiagiZ,,%^

[0070] 氧-―g(Z丨,耳),之二-g(Z2,3Z2)

[0071] 0;=[< (/(0< (d ⑴-(/z?-

[0072] Gj – – (1.-一(1 一 -4+:(1 –

[0073] ei=[0i-i I 07-i] ,i = l ,7 [0074]那么系统(4)渐近稳定。

[0075] 对于矩阵A,He(A)=A+AT。其中I代表单位矩阵。

[0076]证明:

[0077] 对于判据中的Lyapunov泛函进行求导可得:

[0081 ] 然后对求导后的Lyapunov泛函(7)的最后两项利用引用1,2处理,以

[0083]针对式(9)分别使用引理1可得: 处理方法相同,具体操作如下: (9)

[0085] 其中:

[0086] 111 = 63^1 ^12=64^1 ^21 = G5Cl^22 = G6ll [0087] 针对上式,使用引理2可得:

[0090]另一项-「了 办采用同样的处理方法:

[0092] 将处理后的两项代入到求导后的Lyapunov泛函中去可得:

[0096]系统渐近稳定,需要使所求Lyapunov泛函导数小于0,判据得证。

[0097]注:判据中的不等式(6)依赖于d (t)和,无法直接使用LMI工具箱求解。但 是不等式(6)是关于d(t)和的凸函数,所以只要使上述不等式在cKUzln, d(t) = -jii : d(t) – hx . d(t) = j.i \ d{l) -h2 , d{t) = -pi ; d{l)二 h2、?'(〇 = “上都成立即 可。

[0098]下面介绍本发明的一个实施例:

[0099]四机两区域电力系统如图3所示,1号发电机上安装有广域阻尼控制器,选择《13作 为控制器反馈信号。广域阻尼控制器常规超前滞后WADC,如下式所示:

(10)

[0101] 其中:Tw=10s,Ti = 0.324s,T2 = 0.212s

[0102] 在111 = 0^ = 0情况下,利用本发明所提方法得到四机两区域电力系统在1=10与 Ka = 22时系统的稳定裕度分别为346.4ms与101.0ms。

[0103]设置系统在母线3处发生三相短路故障,持续200ms。通过图4中无时滞条件下的系 统响应,可以看出广域控制器WADC优化了系统的性能,消除了内部振荡。图5中给出了系统 中存在110ms延时情况下的系统响应,可以看出延时110ms情况下,K a= 10系统是稳定的,而 Ka = 22系统是不稳定的,呈现发散趋势。这也符合本发明所提方法求得Ka = 22条件下的时滞 稳定裕度101 .〇ms。

【主权项】

1. 一种基于Wirtinger不等式的时滞电力系统稳定性判定方法,其特征在于,包括如下 步骤: (1) 建立包含广域控制回路的时滞电力系统模型?中:In < d (t)^h2, d{t) < μ -μ^Ι ;(2) 给定稳定判定条件: 若存在正定矩阵P e R4nX4n;正定矩阵& e RnXn,i = 1,2,3;正定矩阵^ e RnXn,j = 1,2;矩 阵XkeR2nX2n,k=l,2使下列矩阵不等式成立,则时滞电力系统是渐近稳定的:ei=[Oi-1 I 〇7-1],1 = 1,···,7 对于矩阵A,He(A)=A+AT。其中I代表单位矩阵。 (3) 利用Matlab中的线性矩阵(LMI)工具箱判断给定时滞d(t)是否满足步骤(2)给出的 判定条件,若满足,则可判定在延时d(t)条件下的时滞电力系统是渐近稳定的。2. 如权利要求1所述的基于Wirtinger不等式的时滞电力系统稳定性判定方法,其特征 在于,时滞电力系统模型式中:χ=[χι X。]1″, XI为系统状态变量;X。为控制器的状态变量;Αι为系统 状态矩阵;Βι为系统输入矩阵;A。为控制系统的状态矩阵;C。为控制系统的输出矩阵;Ci为系 统输出矩阵;B。为控制系统的输入矩阵。

【文档编号】H02J3/00GK105958476SQ201610299335

【公开日】2016年9月21日

【申请日】2016年5月6日

【发明人】孙永辉, 李宁, 卫志农, 孙国强, 张世达, 郭敏, 秦晨

【申请人】河海大学

版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 举报,一经查实,本站将立刻删除。

发布者:全栈程序员-用户IM,转载请注明出处:https://javaforall.cn/187007.html原文链接:https://javaforall.cn

【正版授权,激活自己账号】: Jetbrains全家桶Ide使用,1年售后保障,每天仅需1毛

【官方授权 正版激活】: 官方授权 正版激活 支持Jetbrains家族下所有IDE 使用个人JB账号...

(0)


相关推荐

  • python chr()和ord()_Python函数ord

    python chr()和ord()_Python函数ordord()函数主要用来返回对应字符的ascii码,chr()主要用来表示ascii码对应的字符他的输入时数字,可以用十进制,也可以用十六进制。例如:printord('a)#97pri

  • 对Spring的IOC和AOP的理解,可做面试题用「建议收藏」

    对Spring的IOC和AOP的理解,可做面试题用「建议收藏」IOC:控制反转,也叫DI依赖注入,它并不是一种技术实现,而是一种设计思想。在实际项目开发中,我们往往是通过类与类之间的相互协作来完成特定的业务逻辑,这个时候,每个类都要管理与自己有交互的类的引用和依赖,这就使得代码的维护异常困难并且耦合度过高,而IOC的出现正是为了解决这个问题,IOC将类与类的依赖关系写在配置文件中,程序在运行时根据配置文件动态加载依赖的类,降低的类与类之间的耦合度。它的原…

  • 大数据应用实例[通俗易懂]

    大数据应用实例[通俗易懂]通过对全国各地不同时间段的购买华硕笔记本电脑的销售量进行大数据处理,将其通过汇总分析,即可制成图表的格式,使用户能对华硕笔记本销售量进行研究,可以查看出华硕笔记本最近的销售同比(与上一年的销售比较)与环比(与上月的销售比较),购买的人数地域分布量,及购买的人群年龄分布,性别分布等等首先就是对华硕笔记本最近的销售同比(与上一年的销售比较)与环比(与上月的销售比较),整体同比及移…

  • Aria2使用详细教程「建议收藏」

    Aria2使用详细教程「建议收藏」随着百度网盘限速,迅雷限速,逼着用户花钱充会员,越来越多的激活成功教程方法和替代工具出现,Aria2就是这样一种工具,Aria2的简介就不赘述了,网上多的是。

  • Hexo-Matery主题细致美化

    Hexo-Matery主题细致美化Hexo-Matery主题美化在一番瞎改js代码后,终于无法忍受next主题,于是愤然投入Matery大家庭,结果证明,香!下面是我记录的配置Matery主题的流程,仅供后来的师傅们参考。大家可以来我Hexo博客主页看看具体效果我的小站下载安装Matery主题首先在hexo官网主题里面搜索Matery主题点击进入作者大大的github,然后下载master分支的最新稳定版的代码,解压缩后将hexo-theme-matery的文件夹重命名为matery,复制到Hexo的theme

  • Oracle中关于bitmap index的使用问题

    Oracle中关于bitmap index的使用问题

发表回复

您的电子邮箱地址不会被公开。

关注全栈程序员社区公众号