发明名称 |
一种提高有效载荷观测效率的卫星编队构形保持控制方法 |
摘要 |
本发明公开了一种提高有效载荷观测效率的卫星编队构形保持控制方法,其包括如下步骤:据工程实际建立具体约束条件;通过相对导航信息获取编队参数调整量(Δa,Δe,Δu),计算编队保持控制量;根据工程约束条件确定γ可行的取值范围;根据评价方程优化γ,然后计算编队控制量,对编队构形进行保持控制,实现燃料消耗优化。本发明能在近圆非赤道轨道条件下的编队卫星构形保持控制,同时满足实际有效载荷的观测需求。 |
申请公布号 |
CN105353767A |
申请公布日期 |
2016.02.24 |
申请号 |
CN201510883179.4 |
申请日期 |
2015.12.03 |
申请人 |
上海新跃仪表厂 |
发明人 |
完备;杜耀珂;陈敏;沈阳;何煜斌;崔佳;汪礼成;郑科宇;吴敬玉 |
分类号 |
G05D1/10(2006.01)I;G05B13/04(2006.01)I |
主分类号 |
G05D1/10(2006.01)I |
代理机构 |
上海汉声知识产权代理有限公司 31236 |
代理人 |
郭国中;刘翠 |
主权项 |
一种提高有效载荷观测效率的卫星编队构形保持控制方法,其特征在于,包括如下步骤:步骤一、根据工程实际建立具体约束条件;步骤二、根据相对导航信息获取编队参数调整量(Δa,Δe,Δu),代入下式计算编队保持控制量:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mi>ΔV</mi><mn>1</mn></msub><mo>=</mo><mfrac><mn>1</mn><mrow><mo>(</mo><mrow><mn>4</mn><mo>+</mo><mn>4</mn><mi>γ</mi><mo>/</mo><mi>π</mi></mrow><mo>)</mo></mrow></mfrac><mrow><mo>[</mo><mrow><mfrac><mrow><mi>n</mi><mi>a</mi></mrow><mrow><mn>4</mn><mi>cos</mi><mi>γ</mi></mrow></mfrac><mrow><mo>(</mo><mrow><mn>1</mn><mo>+</mo><mn>2</mn><mi>γ</mi><mo>/</mo><mi>π</mi></mrow><mo>)</mo></mrow><mi>Δ</mi><mi>e</mi><mo>-</mo><mfrac><mi>n</mi><mn>4</mn></mfrac><mrow><mo>(</mo><mrow><mn>3</mn><mo>+</mo><mn>4</mn><mi>γ</mi><mo>/</mo><mi>π</mi></mrow><mo>)</mo></mrow><mi>Δ</mi><mi>a</mi><mo>-</mo><mfrac><mrow><mi>n</mi><mi>a</mi></mrow><mrow><mn>3</mn><mi>π</mi></mrow></mfrac><mrow><mo>(</mo><mrow><mi>Δ</mi><mi>u</mi><mo>+</mo><mi>d</mi><mi>u</mi></mrow><mo>)</mo></mrow></mrow><mo>]</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>ΔV</mi><mn>2</mn></msub><mo>=</mo><mfrac><mn>1</mn><mrow><mo>(</mo><mrow><mn>4</mn><mo>+</mo><mn>4</mn><mi>γ</mi><mo>/</mo><mi>π</mi></mrow><mo>)</mo></mrow></mfrac><mrow><mo>[</mo><mrow><mo>-</mo><mfrac><mrow><mi>n</mi><mi>a</mi></mrow><mrow><mn>4</mn><mi>cos</mi><mi>γ</mi></mrow></mfrac><mrow><mo>(</mo><mrow><mn>3</mn><mo>+</mo><mn>2</mn><mi>γ</mi><mo>/</mo><mi>π</mi></mrow><mo>)</mo></mrow><mi>Δ</mi><mi>e</mi><mo>-</mo><mfrac><mi>n</mi><mn>4</mn></mfrac><mrow><mo>(</mo><mrow><mn>3</mn><mo>+</mo><mn>4</mn><mi>γ</mi><mo>/</mo><mi>π</mi></mrow><mo>)</mo></mrow><mi>Δ</mi><mi>a</mi><mo>-</mo><mfrac><mrow><mi>n</mi><mi>a</mi></mrow><mrow><mn>3</mn><mi>π</mi></mrow></mfrac><mrow><mo>(</mo><mrow><mi>Δ</mi><mi>u</mi><mo>+</mo><mi>d</mi><mi>u</mi></mrow><mo>)</mo></mrow></mrow><mo>]</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>ΔV</mi><mn>3</mn></msub><mo>=</mo><mfrac><mn>1</mn><mrow><mo>(</mo><mrow><mn>4</mn><mo>+</mo><mn>4</mn><mi>γ</mi><mo>/</mo><mi>π</mi></mrow><mo>)</mo></mrow></mfrac><mrow><mo>[</mo><mrow><mfrac><mrow><mi>n</mi><mi>a</mi></mrow><mrow><mn>4</mn><mi>cos</mi><mi>γ</mi></mrow></mfrac><mrow><mo>(</mo><mrow><mn>3</mn><mo>+</mo><mn>2</mn><mi>γ</mi><mo>/</mo><mi>π</mi></mrow><mo>)</mo></mrow><mi>Δ</mi><mi>e</mi><mo>+</mo><mfrac><mi>n</mi><mn>4</mn></mfrac><mrow><mo>(</mo><mrow><mn>7</mn><mo>+</mo><mn>8</mn><mi>γ</mi><mo>/</mo><mi>π</mi></mrow><mo>)</mo></mrow><mi>Δ</mi><mi>a</mi><mo>+</mo><mfrac><mrow><mi>n</mi><mi>a</mi></mrow><mrow><mn>3</mn><mi>π</mi></mrow></mfrac><mrow><mo>(</mo><mrow><mi>Δ</mi><mi>u</mi><mo>+</mo><mi>d</mi><mi>u</mi></mrow><mo>)</mo></mrow></mrow><mo>]</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>ΔV</mi><mn>4</mn></msub><mo>=</mo><mfrac><mn>1</mn><mrow><mo>(</mo><mrow><mn>4</mn><mo>+</mo><mn>4</mn><mi>γ</mi><mo>/</mo><mi>π</mi></mrow><mo>)</mo></mrow></mfrac><mrow><mo>[</mo><mrow><mo>-</mo><mfrac><mrow><mi>n</mi><mi>a</mi></mrow><mrow><mn>4</mn><mi>cos</mi><mi>γ</mi></mrow></mfrac><mrow><mo>(</mo><mrow><mn>1</mn><mo>+</mo><mn>2</mn><mi>γ</mi><mo>/</mo><mi>π</mi></mrow><mo>)</mo></mrow><mi>Δ</mi><mi>e</mi><mo>+</mo><mfrac><mi>n</mi><mn>4</mn></mfrac><mrow><mo>(</mo><mrow><mn>7</mn><mo>+</mo><mn>8</mn><mi>γ</mi><mo>/</mo><mi>π</mi></mrow><mo>)</mo></mrow><mi>Δ</mi><mi>a</mi><mo>+</mo><mfrac><mrow><mi>n</mi><mi>a</mi></mrow><mrow><mn>3</mn><mi>π</mi></mrow></mfrac><mrow><mo>(</mo><mrow><mi>Δ</mi><mi>u</mi><mo>+</mo><mi>d</mi><mi>u</mi></mrow><mo>)</mo></mrow></mrow><mo>]</mo></mrow></mrow></mtd></mtr></mtable></mfenced><mo>;</mo></mrow>]]></math><img file="FDA0000866786700000011.GIF" wi="1542" he="583" /></maths>此时的控制量为γ的函数;其中,第一次喷气是在纬度幅角为<img file="FDA0000866786700000012.GIF" wi="399" he="135" />时刻,喷气量将产生的切向速度增量是ΔV<sub>1</sub>;第二次喷气是在纬度幅角为<img file="FDA0000866786700000013.GIF" wi="498" he="142" />时刻,喷气量将产生的切向速度增量是ΔV<sub>2</sub>;第三次喷气是在纬度幅角为<img file="FDA0000866786700000014.GIF" wi="523" he="142" />时刻,喷气量将产生的切向速度增量是ΔV<sub>3</sub>;第四次喷气是在纬度幅角为<img file="FDA0000866786700000015.GIF" wi="517" he="143" />时刻,喷气量将产生的切向速度增量是ΔV<sub>4</sub>。可以看出γ就是相对于u<sub>0</sub>(arctan(δΔe<sub>Y</sub>/δΔe<sub>X</sub>))偏置的角度,从而保证编队控制的纬度幅角不会在有效载荷的观测范围。同时考虑触发时间与起控时间存在时间差的影响,记触发时间对应的纬度幅角为u<sub>0</sub>,起控时间对应的纬度幅角为u<sub>1</sub>,则du=‑1.5(u<sub>1</sub>‑u<sub>0</sub>)Δa步骤三、根据步骤1所得的工程约束条件确定γ可行的取值范围;考虑极端情况,成像区域为[u<sub>0</sub>,u<sub>0</sub>+β°],则根据前述约束条件得到γ必须满足如下条件才能保证4次脉冲控制均不在有效载荷观测范围γ≥γ<sub>1</sub>°,其中,(90>γ<sub>1</sub>>0)步骤四、得到最优γ(即γ=γ<sub>1</sub>°),代入以下优化方程,实现燃料消耗优化;f<sub>min</sub>=|ΔV<sub>1</sub>|+|ΔV<sub>2</sub>|+|ΔV<sub>3</sub>|+|ΔV<sub>4</sub>|。 |
地址 |
200233 上海市徐汇区宜山路710号 |