[发明专利]一种无拖曳卫星姿态通道的抗干扰滤波方法有效

专利信息
申请号: 201610195609.8 申请日: 2016-03-31
公开(公告)号: CN105629986B 公开(公告)日: 2017-06-06
发明(设计)人: 郭雷;迟伟杰;乔建忠;朱玉凯;李文硕 申请(专利权)人: 北京航空航天大学
主分类号: G05D1/08 分类号: G05D1/08
代理公司: 北京科迪生专利代理有限责任公司11251 代理人: 杨学明,顾炜
地址: 100191*** 国省代码: 北京;11
权利要求书: 查看更多 说明书: 查看更多
摘要:
搜索关键词: 一种 拖曳 卫星 姿态 通道 抗干扰 滤波 方法
【权利要求书】:

1.一种无拖曳卫星姿态通道的抗干扰滤波方法,其特征在于包括以下步骤:首先,建立含有多源干扰的无拖曳卫星姿态通道的数学模型;其次,建立外部环境干扰力矩的数学模型;再次,设计针对无拖曳卫星姿态通道的抗干扰滤波器;最后,设计最优控制器,进行反馈控制;具体步骤如下:

(1)建立含有多源干扰的无拖曳卫星姿态通道的数学模型:

<mrow><mover><mi>x</mi><mo>&CenterDot;</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mi>M</mi><mi>x</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><mi>N</mi><mo>&lsqb;</mo><mi>u</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><mi>w</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><mi>d</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>d</mi><mi>J</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>&rsqb;</mo></mrow>

y(t)=Hx(t)+v(t)

其中,x(t)为含有多源干扰的无拖曳卫星姿态通道的数学模型的状态变量,y(t)为含有多源干扰的无拖曳卫星姿态通道的数学模型的输出变量,系数矩阵H=[I3×3 03×3],无拖曳卫星的标称转动惯量矩阵J0=diag{J0x,J0y,J0z},ω0为轨道角速度,u(t)为无拖曳卫星所受到的控制力矩,d(t)为无拖曳卫星所受到的外部环境干扰力矩,dJ(t)为惯量不确定性带来的等价干扰,w(t)和v(t)分别表示无拖曳卫星姿态通道的执行机构噪声和量测噪声;

(2)建立外部环境干扰力矩的数学模型:

<mfenced open = "{" close = ""><mtable><mtr><mtd><mover><mi>&psi;</mi><mo>&CenterDot;</mo></mover><mo>(</mo><mi>t</mi><mo>)</mo><mo>=</mo><mi>P</mi><mi>&psi;</mi><mo>(</mo><mi>t</mi><mo>)</mo><mo>+</mo><mi>&xi;</mi><mo>(</mo><mi>t</mi><mo>)</mo></mtd></mtr><mtr><mtd><mi>d</mi><mo>(</mo><mi>t</mi><mo>)</mo><mo>=</mo><mi>Q</mi><mi>&psi;</mi><mo>(</mo><mi>t</mi><mo>)</mo></mtd></mtr></mtable></mfenced>

其中,ψ(t)为外部环境干扰力矩的数学模型的状态变量,P和Q为系数矩阵,ξ(t)为模型不确定性项;

(3)建立增广的含有多源干扰的无拖曳卫星姿态通道的数学模型:

<mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mover><mover><mi>x</mi><mo>&OverBar;</mo></mover><mo>&CenterDot;</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mover><mi>M</mi><mo>&OverBar;</mo></mover><mover><mi>x</mi><mo>&OverBar;</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><mover><mi>N</mi><mo>&OverBar;</mo></mover><mo>&lsqb;</mo><mi>u</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><mi>w</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><mi>d</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>d</mi><mi>J</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>&rsqb;</mo><mo>+</mo><mover><mi>G</mi><mo>&OverBar;</mo></mover><mi>&xi;</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><mover><mi>y</mi><mo>&OverBar;</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mover><mi>H</mi><mo>&OverBar;</mo></mover><mover><mi>x</mi><mo>&OverBar;</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><mi>v</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mfenced>

其中,系数矩阵

(4)设计针对无拖曳卫星姿态通道的抗干扰滤波器:

<mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mover><mover><mover><mi>x</mi><mo>&OverBar;</mo></mover><mo>^</mo></mover><mo>&CenterDot;</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mover><mi>M</mi><mo>&OverBar;</mo></mover><mover><mover><mi>x</mi><mo>&OverBar;</mo></mover><mo>^</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><mover><mi>N</mi><mo>&OverBar;</mo></mover><mi>u</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><mi>K</mi><mo>&lsqb;</mo><mover><mi>y</mi><mo>&OverBar;</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><mover><mover><mi>y</mi><mo>&OverBar;</mo></mover><mo>^</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>&rsqb;</mo></mrow></mtd></mtr><mtr><mtd><mrow><mover><mover><mi>y</mi><mo>&OverBar;</mo></mover><mo>^</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mover><mi>H</mi><mo>&OverBar;</mo></mover><mover><mover><mi>x</mi><mo>&OverBar;</mo></mover><mo>^</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mfenced>

其中,为的估计值,为的估计值,K为针对无拖曳卫星姿态通道的抗干扰滤波器的增益矩阵;

(5)设计最优控制器,进行反馈控制:

<mrow><mi>u</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mi>L</mi><mover><mi>x</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><mi>Q</mi><mover><mi>&psi;</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow>

其中,L为最优控制器的增益矩阵,为含有多源干扰的无拖曳卫星姿态通道的数学模型状态变量x(t)的估计值,为外部环境干扰力矩的数学模型状态变量ψ(t)的估计值,V为对称正定矩阵,即V=VT>0,未知矩阵S为下列Riccati方程的解:

<mrow><mi>S</mi><mover><mi>M</mi><mo>&OverBar;</mo></mover><mo>+</mo><msup><mover><mi>M</mi><mo>&OverBar;</mo></mover><mi>T</mi></msup><mi>S</mi><mo>-</mo><mi>S</mi><mover><mi>N</mi><mo>&OverBar;</mo></mover><msup><mi>V</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mover><mi>N</mi><mo>&OverBar;</mo></mover><mi>S</mi><mo>+</mo><mi>W</mi><mo>=</mo><mn>0</mn></mrow>

其中,W为对称半正定矩阵,即W=WT≥0。

下载完整专利技术内容需要扣除积分,VIP会员可以免费下载。

该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于北京航空航天大学,未经北京航空航天大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服

本文链接:http://www.vipzhuanli.com/pat/books/201610195609.8/1.html,转载请声明来源钻瓜专利网。

×

专利文献下载

说明:

1、专利原文基于中国国家知识产权局专利说明书;

2、支持发明专利 、实用新型专利、外观设计专利(升级中);

3、专利数据每周两次同步更新,支持Adobe PDF格式;

4、内容包括专利技术的结构示意图流程工艺图技术构造图

5、已全新升级为极速版,下载速度显著提升!欢迎使用!

请您登陆后,进行下载,点击【登陆】 【注册】

关于我们 寻求报道 投稿须知 广告合作 版权声明 网站地图 友情链接 企业标识 联系我们

钻瓜专利网在线咨询

周一至周五 9:00-18:00

咨询在线客服咨询在线客服
tel code back_top