发明名称 基于光纤陀螺惯性测量系统的船舶横纵荡信息测量方法
摘要 本发明的目的在于提供基于光纤陀螺惯性测量系统的船舶横纵荡信息测量方法。实时采集光纤陀螺和加速度计的输出信号,得到船舶实时姿态信息,从而构成载体坐标系与地理坐标系的关系矩阵,由船舶的主航向角得到载体坐标系与半固定坐标系之间的转换矩阵,利用输出信号得到地理坐标系上第n个采样点时船舶实时速度信息,通过坐标转化得到半固定坐标系x轴和y轴上第n个采样点时的速度信息,将上述速度信息进行积分,得到半固定坐标系下第n个采样点时总的位移量以及半固定坐标系下x轴上总的位移量,对上述位移量进行滤波,得到船舶的纵荡信息和横荡信息。本发明不仅增加了原有惯性测量系统功能,同时可以提高系统的导航参数测量精度。
申请公布号 CN101881620A 申请公布日期 2010.11.10
申请号 CN201010209309.3 申请日期 2010.06.25
申请人 哈尔滨工程大学 发明人 奔粤阳;李倩;高伟;龚晶;于飞
分类号 G01C21/18(2006.01)I;G01C21/10(2006.01)I 主分类号 G01C21/18(2006.01)I
代理机构 代理人
主权项 1.基于光纤陀螺惯性测量系统的船舶横纵荡信息测量方法,其特征是:(1)对船用光纤陀螺捷联惯性系统充分预热,并实时采集光纤陀螺和加速度计的输出信号;(2)利用光纤陀螺和加速度计的输出,得到船舶实时姿态信息,包括纵摇角α、横摇角β、艏摇角γ,再由姿态信息构成载体坐标系b与地理坐标系t的关系矩阵<maths num="0001"><![CDATA[<math><mrow><msubsup><mi>C</mi><mi>b</mi><mi>t</mi></msubsup><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>C</mi><mn>11</mn></msub></mtd><mtd><msub><mi>C</mi><mn>12</mn></msub></mtd><mtd><msub><mi>C</mi><mn>13</mn></msub></mtd></mtr><mtr><mtd><msub><mi>C</mi><mn>21</mn></msub></mtd><mtd><msub><mi>C</mi><mn>22</mn></msub></mtd><mtd><msub><mi>C</mi><mn>23</mn></msub></mtd></mtr><mtr><mtd><msub><mi>C</mi><mn>31</mn></msub></mtd><mtd><msub><mi>C</mi><mn>32</mn></msub></mtd><mtd><msub><mi>C</mi><mn>33</mn></msub></mtd></mtr></mtable></mfenced><mo>,</mo></mrow></math>]]></maths>其中C<sub>11</sub>=cosβcosγ-sinβsinαsinγC<sub>12</sub>=-cosαsinγC<sub>13</sub>=sinβcosγ+cosβsinαsinγC<sub>21</sub>=cosβsinγ+sinβsinαcosγC<sub>22</sub>=cosαcosγ;C<sub>23</sub>=sinβsinγ-cosβsinαcosγC<sub>31</sub>=-sinβcosαC<sub>32</sub>=sinαC<sub>33</sub>=cosβcosα(3)由船舶航行操控时设定的船舶航迹向与地理北向的夹角,即船舶的主航向角<img file="FSA00000177100300012.GIF" wi="49" he="32" />得到地理坐标系t与半固定坐标系d的关系矩阵<img file="FSA00000177100300013.GIF" wi="475" he="188" />进而得到载体坐标系b与半固定坐标系d之间的转换矩阵<img file="FSA00000177100300014.GIF" wi="215" he="51" />(4)利用陀螺和加速度计的输出,得到地理坐标系t上第n个采样点时船舶实时速度信息,包括地理坐标系x轴上的速度<img file="FSA00000177100300015.GIF" wi="140" he="63" />地理坐标系y轴上速度<img file="FSA00000177100300016.GIF" wi="139" he="64" />地理坐标系z轴上速度<img file="FSA00000177100300017.GIF" wi="136" he="62" />(5)利用矩阵<img file="FSA00000177100300018.GIF" wi="73" he="51" />将地理坐标系t上船舶实时速度信息通过坐标转化,得到半固定坐标系x轴上第n个采样点时的速度信息<img file="FSA00000177100300021.GIF" wi="118" he="63" />和半固定坐标系y轴上第n个采样点时的速度信息<img file="FSA00000177100300022.GIF" wi="120" he="65" /><img file="FSA00000177100300023.GIF" wi="711" he="219" />(6)将半固定坐标系d下x轴上的速度信息<img file="FSA00000177100300024.GIF" wi="120" he="62" />和y轴上的速度信息<img file="FSA00000177100300025.GIF" wi="118" he="65" />进行一次积分,得到半固定坐标系d下第n个采样点时总的位移量,即半固定坐标系d下y轴上总的位移量<img file="FSA00000177100300026.GIF" wi="144" he="65" />以及半固定坐标系d下x轴上总的位移量<img file="FSA00000177100300027.GIF" wi="118" he="63" /><maths num="0002"><![CDATA[<math><mrow><msubsup><mi>s</mi><mi>y</mi><mi>d</mi></msubsup><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><mi>h</mi><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msubsup><mi>v</mi><mi>y</mi><mi>d</mi></msubsup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>,</mo></mrow></math>]]></maths><maths num="0003"><![CDATA[<math><mrow><msubsup><mi>s</mi><mi>x</mi><mi>d</mi></msubsup><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><mi>h</mi><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msubsup><mi>v</mi><mi>x</mi><mi>d</mi></msubsup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow></math>]]></maths>其中k=1,2,...n,<img file="FSA000001771003000210.GIF" wi="106" he="56" />表示半固定坐标系d下x轴上第k个采样点时的速度信息,<img file="FSA000001771003000211.GIF" wi="105" he="56" />表示半固定坐标系d下y轴上第k个采样点时的速度信息,<img file="FSA000001771003000212.GIF" wi="104" he="56" />与<img file="FSA000001771003000213.GIF" wi="105" he="56" />可以在第k个采样点时由步骤1到步骤5的过程测量得到并加以保存,h为光纤陀螺惯性测量系统的采样周期;(7)对步骤(6)中得到的半固定坐标系d下总的位移量<img file="FSA000001771003000214.GIF" wi="118" he="65" />与<img file="FSA000001771003000215.GIF" wi="117" he="63" />进行滤波,得到船舶沿艏向和右舷方向的即时高频位移量,即船舶的纵荡信息<img file="FSA000001771003000216.GIF" wi="91" he="53" />和横荡信息<img file="FSA000001771003000217.GIF" wi="115" he="47" /><maths num="0004"><![CDATA[<math><mrow><msub><mover><mi>s</mi><mo>^</mo></mover><mi>y</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><msubsup><mi>s</mi><mi>y</mi><mi>d</mi></msubsup><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>&CenterDot;</mo><mi>&omega;</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>,</mo></mrow></math>]]></maths><maths num="0005"><![CDATA[<math><mrow><msub><mover><mi>s</mi><mo>^</mo></mover><mi>x</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><msubsup><mi>s</mi><mi>x</mi><mi>d</mi></msubsup><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>&CenterDot;</mo><mi>&omega;</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></mrow></math>]]></maths>选用高通数字FIR滤波器,选取凯塞窗为<maths num="0006"><![CDATA[<math><mrow><mi>&omega;</mi><mo>[</mo><mi>n</mi><mo>]</mo><mo>=</mo><mfrac><mrow><msub><mi>I</mi><mn>0</mn></msub><mo>[</mo><mn>6</mn><msqrt><mn>1</mn><mo>-</mo><msup><mrow><mo>(</mo><mfrac><mrow><mn>2</mn><mi>n</mi></mrow><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></mfrac><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mn>2</mn></msup></msqrt><mo>]</mo></mrow><mn>70</mn></mfrac><mo>,</mo></mrow></math>]]></maths>其中函数定义为<maths num="0007"><![CDATA[<math><mrow><msub><mi>I</mi><mn>0</mn></msub><mo>[</mo><mi>x</mi><mo>]</mo><mo>=</mo><mn>1</mn><mo>+</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mn>20</mn></munderover><msup><mrow><mo>[</mo><mfrac><msup><mrow><mo>(</mo><mi>x</mi><mo>/</mo><mn>2</mn><mo>)</mo></mrow><mi>j</mi></msup><mrow><mi>j</mi><mo>!</mo></mrow></mfrac><mo>]</mo></mrow><mn>2</mn></msup><mo>,</mo></mrow></math>]]></maths>并且有<maths num="0008"><![CDATA[<math><mrow><mi>N</mi><mo>=</mo><mn>5.2</mn><mfrac><msub><mi>f</mi><mi>s</mi></msub><mi>f</mi></mfrac><mo>,</mo></mrow></math>]]></maths>其中f<sub>s</sub>为采样频率,f为通带边缘频率。
地址 150001 黑龙江省哈尔滨市南岗区南通大街145号哈尔滨工程大学科技处知识产权办公室