[发明专利]应用于保密通信的受控Shimizu-Morioka系统与T系统的广义混沌同步方法在审
申请号: | 201811137686.3 | 申请日: | 2018-09-28 |
公开(公告)号: | CN109474427A | 公开(公告)日: | 2019-03-15 |
发明(设计)人: | 张端;孙莹 | 申请(专利权)人: | 浙江工业大学 |
主分类号: | H04L9/12 | 分类号: | H04L9/12;H04L9/00 |
代理公司: | 杭州斯可睿专利事务所有限公司 33241 | 代理人: | 王利强 |
地址: | 310014 浙江省*** | 国省代码: | 浙江;33 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | 一种应用于保密通信的受控Shimizu‑Morioka系统与T系统的广义混沌同步方法,所述方法包括以下步骤:1)广义混沌同步问题描述;2)响应系统的状态变换和反馈;3)驱动系统的状态转换;4)广义同步。本发明提供一种应用于保密通信的受控Shimizu‑Morioka系统与T系统的广义混沌同步方法,以T系统为驱动系统,以单输入的受控Shimizu‑Morioka系统为响应系统,采用微分几何中向量场的李导数方法设计一种混沌同步算法,实现广义同步,控制品质较高。 | ||
搜索关键词: | 混沌同步 受控 保密通信 驱动系统 响应系统 应用 问题描述 状态变换 状态转换 单输入 向量场 导数 算法 反馈 | ||
【主权项】:
1.一种应用于保密通信的受控Shimizu‑Morioka系统与T系统的广义混沌同步方法,其特征在于,包括以下步骤:1)广义混沌同步问题描述驱动系统为T系统,形式如下:
其中x=(x1,x2,x3)T是状态变量,a、b和c是已知的正实数参数,2a‑b≠0;以受控Shimizu‑Morioka系统为响应系统,形式如下:
其中ξ=(ξ1,ξ2,ξ3)T是状态变量,u是标量输入,α和β为系统中已知的正实数参数,此β与式(1)中的b等值,即β=b;广义混沌同步要实现的目标是:在驱动系统(1)与响应系统(2)初值分别为x(t0)和ξ(t0),响应系统轨迹经过状态反馈u=u(x,ξ,t) (3)其中t表示时间,和相空间之间的状态变换ξ=T(x) (4)后趋向于驱动系统的轨迹,即
这里||·||代表
空间中向量的2‑范数;2)响应系统的状态变换和反馈对响应系统(2)作如下状态变换η=S(ξ)其中η=(η1,η2,η3)T
所以,这是一个线性变换,MS为3阶方阵,此线性变换的逆变换为
以η为状态,系统表示为
作反馈u=‑(1‑η1)η2+αη3+u0 (9)再考虑到β=b,系统简化为
该系统属于受控的下三角系统,三阶下三角系统的一般形式为
其中w为输入控制量;另一方面观察系统(10)的后两个等式实际形成了线性系统形式,所以系统(10)为实现部分线性化;3)驱动系统的状态转换为了找寻驱动系统(1)的状态变换以简化系统,先为此系统加上控制量成为
其中v为加入的输入控制量;系统(12)作反馈v=(a‑c)x1+ax1x3+v1 (13)系统可简化为
考虑将系统(14)通过状态变换和进一步的反馈转换为更为简单与系统(10)更为相似的形式,以便于设计广义同步控制方法;由于系统(10)为下三角形式的受控常微分方程,希望系统(14)能转换为同样后者相似形式,为此,记系统(14)的漂移向量场为
以及输入向量场为
令向量场
计算如下向量场李括号
注意
在全局范围内秩为2,并且
说明此分布对合;令
计算如下向量场李括号
在x1=0或者2a‑b=0时
是秩仍为2,这也说明系统(14)不可能实现状态反馈线性化;但是,当2a‑b≠0时,
仅在一个零测度集内秩为2,此集合之外秩均为3,所以系统(14)可以通过状态变换等价转换为下三角系统(11),然而,仍需探究系统(14)究竟能转换为何种下三角系统的问题,并且希望得到形式上较为简单的下三角系统,为此,注意到
此时全局范围内分布
的秩为3并且对合,令
取如下分布Δ0=span{X0};Δ1=span{X0,X1};Δ2=span{X0,X1,X2}, (23)分布Δ0,Δ1,Δ2及X0,X1,X2具有如下性质:①可验证[X0,X1]=0,[X1,X2]=0以及[X0,X2]=0②由①,Δ0,Δ1,Δ2均为对合分布;③由①,存在状态变换h=(h1(x),h2(x),h3(x))T=H(x)=H(x1,x2,x3)满足
④由于
说明性质③中的状态变换h下,系统必定仍然具有下三角形式;上述性质③也意味着满足以下3个偏微分方程组,第1组为
其中h1(x)为光滑函数,符号”L”表示做李导数,第2组为
其中h2(x)为光滑函数,第3组为
其中h3(x)为光滑函数,上述3组偏微分方程的可行解分别为
在h=(h1,h2,h3)T状态下系统成为
该系统实际上已具有下三角系统如系统(11)的形式,但从系统(29)的第二个式子看,上述系统通过如下状态变换还可进一步简化
如用y=(y1,y2,y3)T状态表示x状态则为
在y状态下写出系统
对比系统(12)和系统(32)可知经过状态变化(30)系统(1)成为
上述系统的前2个方程在形式于驱动系统的等价形式(34)已实现一致;4)广义同步现在考虑系统(33)与系统(35)的同步问题,令二者状态差为e=η‑y=(e1,e2,e3)T,则
设计反馈
系统表示为
对于上述系统的子系统
根据线性系统的经典方法设计如下控制器:
该控制器下系统(39)将在的有限时间控制内,即t1时刻实现e2(t1)=e3(t1)=0,设计一种控制器从t0时刻开始,经有限时间实现e2(t1)=e3(t1)=0,并保证此过程中控制量有连续的一阶导数并过渡到0;首先,设计预想的e2(t)为
其中p(t)为一元多项式,由于要求t1时刻到达系统(39)的原点以及u1在t>t0范围内有连续的一阶导数,这意味着e2(t)在t1时有连续的三阶导数,实际上e2(t)和其一、二、三阶导数再t1时刻为保证连续均只能为0,即
再考虑系统(39)的t0时刻应满足
由于式(42)和式(43)共给出6个条件,所以p(t0)应为5次多项式,再利用式(42)得
其中C0和C1为待定系数,利用式(43)的第1个式子得到
再由式(43)的第2个式子
整理得到
该e2(t)满足式(42)和式(43)的各项要求,那么
以及
明显e2(t1)=e3(t1)=u1(t1)=0;在时间t1之后,系统(38)的第一个方程成为
此方程明显是大范围渐进稳定的,从而系统(38)大范围渐进稳定,说明系统(50)与系统(33)在此控制律下实现同步;回到系统(1)与系统(2)的广义同步问题,系统(2)经过状态变换(6)成为系统(51),系统(1)作了反馈和状态变换后成为系统(33),其间的状态变换需综合式(28)以及式(30)
并命名此状态变换为y=(y1,y2,y3)T=Y((x1,x2,x3)T)=Y(x),而控制律可见式(37),其中u1的表达式见式(49)。
下载完整专利技术内容需要扣除积分,VIP会员可以免费下载。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于浙江工业大学,未经浙江工业大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/patent/201811137686.3/,转载请声明来源钻瓜专利网。