一种斜置冗余惯导系统标定误差模型建模方法与流程

文档序号:12444827阅读:来源:国知局

技术特征:

1.一种斜置冗余惯导系统标定误差模型建模方法,其特征在于,该方法包括以下步骤:

S1针对斜置冗余系统不同配置方案进行分析,得到系统的配置矩阵和变换矩阵:

在非正交配置的系统中,惯性原件敏感得到的信息需要先投影到惯导系统的载体系再进行解算,设冗余惯导系统的冗余度为n,即系统由n个陀螺和n个加速度计组成,n≥3,则n冗余度的斜置冗余惯导系统中,陀螺输出为:

<mrow> <msup> <mi>&omega;</mi> <mi>s</mi> </msup> <mo>=</mo> <msup> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msubsup> <mi>&omega;</mi> <mn>1</mn> <mi>s</mi> </msubsup> </mtd> <mtd> <msubsup> <mi>&omega;</mi> <mn>2</mn> <mi>s</mi> </msubsup> </mtd> <mtd> <mo>...</mo> </mtd> <mtd> <msubsup> <mi>&omega;</mi> <mi>n</mi> <mi>s</mi> </msubsup> </mtd> </mtr> </mtable> </mfenced> <mi>T</mi> </msup> </mrow>

ωs表示器件系(s系)中的角增量真值,表示第i个陀螺在器件系(s系)的角增量真值,i=1,2,…,n;

加速度计输出为:

<mrow> <msup> <mi>f</mi> <mi>s</mi> </msup> <mo>=</mo> <msup> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msubsup> <mi>f</mi> <mn>1</mn> <mi>s</mi> </msubsup> </mtd> <mtd> <msubsup> <mi>f</mi> <mn>2</mn> <mi>s</mi> </msubsup> </mtd> <mtd> <mo>...</mo> </mtd> <mtd> <msubsup> <mi>f</mi> <mi>n</mi> <mi>s</mi> </msubsup> </mtd> </mtr> </mtable> </mfenced> <mi>T</mi> </msup> </mrow>

fs表示s系中的速度增量真值,fis表示第i个加速度计在s系的速度增量真值,i=1,2,…,n;

在b系中,载体系右、前、上三个方向的角增量和速度增量分别为可得测量方程:

<mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <msup> <mi>&omega;</mi> <mi>s</mi> </msup> <mo>=</mo> <msub> <mi>H</mi> <mi>g</mi> </msub> <msup> <mi>&omega;</mi> <mi>b</mi> </msup> </mtd> </mtr> <mtr> <mtd> <mrow> <msup> <mi>f</mi> <mi>s</mi> </msup> <mo>=</mo> <msub> <mi>H</mi> <mi>f</mi> </msub> <msup> <mi>f</mi> <mi>b</mi> </msup> </mrow> </mtd> </mtr> </mtable> </mfenced>

其中陀螺和加速度计的配置矩阵Hg,Hf为n×3阶的矩阵,可得变换矩阵Mg,Mf,M=(HTH)-1HT);

S2对斜置冗余惯导系统的分立式标定误差模型进行分析:

惯性器件的测量误差主要由以下三部分构成:

a常值零偏

陀螺和加速度计的常值零偏ε,为:

ε=[ε1 ε2 … εn]T,εi表示第i个陀螺的常值零偏

表示第i个加速度计的常值零偏

b标度因数误差

由于标度因数不准确引入的误差即为器件的标度因数误差,陀螺和加速度计的标度因数误差ΔΚg,ΔΚf为:

表示第i个陀螺的标度因数误差

表示第i个加速度计的标度因数误差

c安装误差

斜置冗余惯导系统与传统正交系统相同,不可避免地存在安装误差,导致实际配置矩阵与理论值H产生误差:

<mrow> <mover> <mi>H</mi> <mo>~</mo> </mover> <mo>=</mo> <mi>H</mi> <mo>+</mo> <mi>&Delta;</mi> <mi>H</mi> <mo>=</mo> <mi>H</mi> <mo>+</mo> <mi>&delta;</mi> <mi>&alpha;</mi> <mo>&CenterDot;</mo> <mi>P</mi> <mo>+</mo> <mi>&delta;</mi> <mi>&beta;</mi> <mo>&CenterDot;</mo> <mi>Q</mi> </mrow>

其中:

<mfenced open = "" close = ""> <mtable> <mtr> <mtd> <mrow> <mi>P</mi> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <msub> <mi>cos&beta;</mi> <mn>1</mn> </msub> <msub> <mi>sin&alpha;</mi> <mn>1</mn> </msub> </mrow> </mtd> <mtd> <mrow> <mo>-</mo> <msub> <mi>cos&beta;</mi> <mn>1</mn> </msub> <msub> <mi>cos&alpha;</mi> <mn>1</mn> </msub> </mrow> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>cos&beta;</mi> <mn>2</mn> </msub> <msub> <mi>sin&alpha;</mi> <mn>2</mn> </msub> </mrow> </mtd> <mtd> <mrow> <mo>-</mo> <msub> <mi>cos&beta;</mi> <mn>2</mn> </msub> <msub> <mi>cos&alpha;</mi> <mn>2</mn> </msub> </mrow> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mo>&CenterDot;</mo> </mtd> <mtd> <mo>&CenterDot;</mo> </mtd> <mtd> <mo>&CenterDot;</mo> </mtd> </mtr> <mtr> <mtd> <mo>&CenterDot;</mo> </mtd> <mtd> <mo>&CenterDot;</mo> </mtd> <mtd> <mo>&CenterDot;</mo> </mtd> </mtr> <mtr> <mtd> <mo>&CenterDot;</mo> </mtd> <mtd> <mo>&CenterDot;</mo> </mtd> <mtd> <mo>&CenterDot;</mo> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>cos&beta;</mi> <mi>n</mi> </msub> <msub> <mi>sin&alpha;</mi> <mi>n</mi> </msub> </mrow> </mtd> <mtd> <mrow> <mo>-</mo> <msub> <mi>cos&beta;</mi> <mi>n</mi> </msub> <msub> <mi>cos&alpha;</mi> <mi>n</mi> </msub> </mrow> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> </mtable> </mfenced> </mrow> </mtd> <mtd> <mrow> <mi>Q</mi> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <mo>-</mo> <msub> <mi>sin&beta;</mi> <mn>1</mn> </msub> <msub> <mi>cos&alpha;</mi> <mn>1</mn> </msub> </mrow> </mtd> <mtd> <mrow> <mo>-</mo> <msub> <mi>sin&beta;</mi> <mn>1</mn> </msub> <msub> <mi>sin&alpha;</mi> <mn>1</mn> </msub> </mrow> </mtd> <mtd> <mrow> <msub> <mi>cos&beta;</mi> <mn>1</mn> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>-</mo> <msub> <mi>sin&beta;</mi> <mn>2</mn> </msub> <msub> <mi>cos&alpha;</mi> <mn>2</mn> </msub> </mrow> </mtd> <mtd> <mrow> <mo>-</mo> <msub> <mi>sin&beta;</mi> <mn>2</mn> </msub> <msub> <mi>sin&alpha;</mi> <mn>2</mn> </msub> </mrow> </mtd> <mtd> <mrow> <msub> <mi>cos&beta;</mi> <mn>2</mn> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mo>&CenterDot;</mo> </mtd> <mtd> <mo>&CenterDot;</mo> </mtd> <mtd> <mo>&CenterDot;</mo> </mtd> </mtr> <mtr> <mtd> <mo>&CenterDot;</mo> </mtd> <mtd> <mo>&CenterDot;</mo> </mtd> <mtd> <mo>&CenterDot;</mo> </mtd> </mtr> <mtr> <mtd> <mo>&CenterDot;</mo> </mtd> <mtd> <mo>&CenterDot;</mo> </mtd> <mtd> <mo>&CenterDot;</mo> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>-</mo> <msub> <mi>sin&beta;</mi> <mi>n</mi> </msub> <msub> <mi>cos&alpha;</mi> <mi>n</mi> </msub> </mrow> </mtd> <mtd> <mrow> <mo>-</mo> <msub> <mi>sin&beta;</mi> <mi>n</mi> </msub> <msub> <mi>sin&alpha;</mi> <mi>n</mi> </msub> </mrow> </mtd> <mtd> <mrow> <msub> <mi>cos&beta;</mi> <mi>n</mi> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow> </mtd> </mtr> </mtable> </mfenced>

βi为第i个惯性器件的安装轴hi与载体系xb-yb平面的夹角,αi为安装轴hi在b系xb-yb平面投影与xb轴的夹角,δαi和δβi为第i个惯性器件实际轴向和理论轴向的误差角,δα为δαi构成的对角矩阵,δβ为δβi构成的对角矩阵,P为αi角误差的变换矩阵,Q为βi角误差的变换矩阵;综合常值零偏、标度因数误差和安装误差,可得冗余系统的陀螺实际输出如下:

<mrow> <msup> <mover> <mi>&omega;</mi> <mo>~</mo> </mover> <mi>s</mi> </msup> <mo>=</mo> <msup> <mi>&omega;</mi> <mi>s</mi> </msup> <mo>+</mo> <mrow> <mo>(</mo> <msub> <mi>&Delta;K</mi> <mi>g</mi> </msub> <msub> <mi>H</mi> <mi>g</mi> </msub> <mo>+</mo> <msub> <mi>K</mi> <mi>g</mi> </msub> <msub> <mi>G</mi> <mi>g</mi> </msub> <mo>)</mo> </mrow> <msubsup> <mi>N</mi> <msub> <mi>g</mi> <mrow> <mi>i</mi> <mi>b</mi> </mrow> </msub> <mi>b</mi> </msubsup> <mo>+</mo> <msup> <mi>&epsiv;</mi> <mi>s</mi> </msup> </mrow>

上式中为陀螺s系中敏感的角增量,该角增量有误差,Kg为陀螺的标度因数,为陀螺在b系中的脉冲输出,εs为陀螺在s系的常值零偏,Hg为陀螺的配置矩阵,Gg为陀螺的安装误差矩阵;

冗余系统加速度计的实际输出如下:

<mrow> <msup> <mover> <mi>f</mi> <mo>~</mo> </mover> <mi>s</mi> </msup> <mo>=</mo> <msup> <mi>f</mi> <mi>s</mi> </msup> <mo>+</mo> <mrow> <mo>(</mo> <msub> <mi>&Delta;K</mi> <mi>f</mi> </msub> <msub> <mi>H</mi> <mi>f</mi> </msub> <mo>+</mo> <msub> <mi>K</mi> <mi>f</mi> </msub> <msub> <mi>G</mi> <mi>f</mi> </msub> <mo>)</mo> </mrow> <msubsup> <mi>N</mi> <msub> <mi>f</mi> <mrow> <mi>i</mi> <mi>b</mi> </mrow> </msub> <mi>b</mi> </msubsup> <mo>+</mo> <msup> <mo>&dtri;</mo> <mi>s</mi> </msup> </mrow>

上式中为加速度计在s系中敏感的角增量,该角增量有误差,Kf为加速度计的标度因数,为加速度计在b系中的脉冲输出,为加速度计在s系的常值零偏,Hf为加速度计的配置矩阵,Gf为加速度计的安装误差矩阵;

综上可得,在s系中陀螺和加速度计的分立式标定误差模型δωs,δfs

<mrow> <msup> <mi>&delta;&omega;</mi> <mi>s</mi> </msup> <mo>=</mo> <msup> <mover> <mi>&omega;</mi> <mo>~</mo> </mover> <mi>s</mi> </msup> <mo>-</mo> <msup> <mi>&omega;</mi> <mi>s</mi> </msup> <mo>=</mo> <mrow> <mo>(</mo> <msub> <mi>&Delta;K</mi> <mi>g</mi> </msub> <msub> <mi>H</mi> <mi>g</mi> </msub> <mo>+</mo> <msub> <mi>K</mi> <mi>g</mi> </msub> <msub> <mi>G</mi> <mi>g</mi> </msub> <mo>)</mo> </mrow> <msubsup> <mi>N</mi> <msub> <mi>g</mi> <mrow> <mi>i</mi> <mi>b</mi> </mrow> </msub> <mi>b</mi> </msubsup> <mo>+</mo> <msup> <mi>&epsiv;</mi> <mi>s</mi> </msup> </mrow>

<mrow> <msup> <mi>&delta;f</mi> <mi>s</mi> </msup> <mo>=</mo> <msup> <mover> <mi>f</mi> <mo>~</mo> </mover> <mi>s</mi> </msup> <mo>-</mo> <msup> <mi>f</mi> <mi>s</mi> </msup> <mo>=</mo> <mrow> <mo>(</mo> <msub> <mi>&Delta;K</mi> <mi>f</mi> </msub> <msub> <mi>H</mi> <mi>f</mi> </msub> <mo>+</mo> <msub> <mi>K</mi> <mi>f</mi> </msub> <msub> <mi>G</mi> <mi>f</mi> </msub> <mo>)</mo> </mrow> <msubsup> <mi>N</mi> <msub> <mi>f</mi> <mrow> <mi>i</mi> <mi>b</mi> </mrow> </msub> <mi>b</mi> </msubsup> <mo>+</mo> <msup> <mo>&dtri;</mo> <mi>s</mi> </msup> </mrow>

S3对斜置冗余惯导系统的系统级标定误差建模进行分析:

在分立式标定完成后,各器件的标度因数误差和安装误差均为小值,可以认为采用分立式标定后的配置矩阵转换后的系统输出在一个非常接近系统载体系的伪b系,因此在b系到s系的n×3阶的安装误差可以简化成伪b系和b系之间的3×3阶的安装误差角,器件的标度因数误差还是以每个器件单独的标度因数来分析,可得陀螺在b系中的系统级标定的误差模型δωb

<mrow> <msup> <mi>&delta;&omega;</mi> <mi>b</mi> </msup> <mo>=</mo> <msup> <mover> <mi>&omega;</mi> <mo>~</mo> </mover> <mi>b</mi> </msup> <mo>-</mo> <msup> <mi>&omega;</mi> <mi>b</mi> </msup> <mo>=</mo> <msub> <mi>&Delta;G</mi> <mi>g</mi> </msub> <msub> <mi>M</mi> <mi>g</mi> </msub> <mrow> <mo>(</mo> <mi>d</mi> <mi>i</mi> <mi>a</mi> <mi>g</mi> <mo>(</mo> <mrow> <msub> <mi>K</mi> <mi>g</mi> </msub> <mo>+</mo> <msub> <mi>&Delta;K</mi> <mi>g</mi> </msub> </mrow> <mo>)</mo> <msubsup> <mi>N</mi> <msub> <mi>g</mi> <mrow> <mi>i</mi> <mi>s</mi> </mrow> </msub> <mi>s</mi> </msubsup> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>M</mi> <mi>g</mi> </msub> <mrow> <mo>(</mo> <mi>d</mi> <mi>i</mi> <mi>a</mi> <mi>g</mi> <mo>(</mo> <mrow> <msub> <mi>&Delta;K</mi> <mi>g</mi> </msub> </mrow> <mo>)</mo> <msubsup> <mi>N</mi> <msub> <mi>g</mi> <mrow> <mi>i</mi> <mi>s</mi> </mrow> </msub> <mi>s</mi> </msubsup> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>M</mi> <mi>g</mi> </msub> <msup> <mi>&epsiv;</mi> <mi>s</mi> </msup> </mrow>

其中Mg为陀螺的变换矩阵,认为是无误差的,为陀螺在s系中的脉冲输出,为陀螺输出在伪b系和b系的安装误差矩阵;

同样可得加速度计在b系中的系统级标定的误差模型δfb

<mrow> <msup> <mi>&delta;f</mi> <mi>b</mi> </msup> <mo>=</mo> <msup> <mover> <mi>f</mi> <mo>~</mo> </mover> <mi>b</mi> </msup> <mo>-</mo> <msup> <mi>f</mi> <mi>b</mi> </msup> <mo>=</mo> <msub> <mi>&Delta;G</mi> <mi>f</mi> </msub> <msub> <mi>M</mi> <mi>f</mi> </msub> <mrow> <mo>(</mo> <mi>d</mi> <mi>i</mi> <mi>a</mi> <mi>g</mi> <mo>(</mo> <mrow> <msub> <mi>K</mi> <mi>f</mi> </msub> <mo>+</mo> <msub> <mi>&Delta;K</mi> <mi>f</mi> </msub> </mrow> <mo>)</mo> <msubsup> <mi>N</mi> <msub> <mi>f</mi> <mrow> <mi>i</mi> <mi>s</mi> </mrow> </msub> <mi>s</mi> </msubsup> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>M</mi> <mi>f</mi> </msub> <mrow> <mo>(</mo> <mi>d</mi> <mi>i</mi> <mi>a</mi> <mi>g</mi> <mo>(</mo> <mrow> <msub> <mi>&Delta;K</mi> <mi>f</mi> </msub> </mrow> <mo>)</mo> <msubsup> <mi>N</mi> <msub> <mi>f</mi> <mrow> <mi>i</mi> <mi>s</mi> </mrow> </msub> <mi>s</mi> </msubsup> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>M</mi> <mi>f</mi> </msub> <msup> <mo>&dtri;</mo> <mi>s</mi> </msup> </mrow>

其中Mf为加速度计的变换矩阵,认为是无误差的,为加速度计在s系中的脉冲输出,为加速度计输出在伪b系和b系的安装误差矩阵。

当前第2页1 2 3 
网友询问留言 已有0条留言
  • 还没有人留言评论。精彩留言会获得点赞!
1