主权项 |
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>Σ</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>Σ</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>·</mo><mi>ω</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>·</mo><mi>ω</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></mrow></math>]]></maths>选用高通数字FIR滤波器,选取凯塞窗为<maths num="0006"><![CDATA[<math><mrow><mi>ω</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>Σ</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为通带边缘频率。 |