CN116045897A - 一种全方位太阳角测量方法及其系统 - Google Patents

一种全方位太阳角测量方法及其系统 Download PDF

Info

Publication number
CN116045897A
CN116045897A CN202211664741.0A CN202211664741A CN116045897A CN 116045897 A CN116045897 A CN 116045897A CN 202211664741 A CN202211664741 A CN 202211664741A CN 116045897 A CN116045897 A CN 116045897A
Authority
CN
China
Prior art keywords
solar angle
sampling
solar
sampling module
fyass
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202211664741.0A
Other languages
English (en)
Inventor
钱方亮
聂章海
程卫强
程朝阳
张磊
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai Aerospace Control Technology Institute
Original Assignee
Shanghai Aerospace Control Technology Institute
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shanghai Aerospace Control Technology Institute filed Critical Shanghai Aerospace Control Technology Institute
Priority to CN202211664741.0A priority Critical patent/CN116045897A/zh
Publication of CN116045897A publication Critical patent/CN116045897A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C1/00Measuring angles

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Photometry And Measurement Of Optical Pulse Characteristics (AREA)

Abstract

本发明公开了一种全方位太阳角测量方法及其系统,该系统通过根据卫星对日定向需求情况以及对日方向配置多台模拟太阳角计,并设计其安装方位,根据多台模拟太阳角计测量信息,实现全方位太阳角的计算测量。其优点是:本发明利用多台模拟太阳角计可以实现全方位的太阳角计算,不需要陀螺信息,综合实现了0‑1式太阳敏感器和太阳角计功能,简化了对日定向所需单机配置,降低了系统成本。

Description

一种全方位太阳角测量方法及其系统
技术领域
本发明涉及太阳角测量领域,具体涉及一种全方位太阳角测量方法及其系统。
背景技术
卫星捕获太阳及对日定向一般使用太阳敏感器,0-1式太阳敏感器信息用于搜索太阳,测量太阳方位角的太阳敏感器(也称太阳角计)用于捕获太阳并对日定向。由于需根据太阳敏感器信号控制卫星按一定逻辑搜索太阳,还需要陀螺测量卫星角速度,以控制执行机构(飞轮或喷气)实现卫星搜索太阳并对日定向。然而上述测量方式所需测量工具及数据众多,特别是需要陀螺额外测量卫星角速度,在实际应用中操作比较繁琐。
发明内容
本发明的目的在于提供一种全方位太阳角测量方法及其系统,该方法利用多台模拟太阳角计可以实现全方位的太阳角计算,不需要陀螺测量信息,综合实现了0-1式太阳敏感器和太阳角计的功能,简化了对日定向所需单机配置,降低了系统成本。
为了达到上述目的,本发明通过以下技术方案实现:
一种全方位太阳角测量系统,包含:
多台模拟太阳角计,所述模拟太阳角计包含多个采样模块,多个采样模块的采样角度覆盖全方位三维角度,随着采样模块的太阳入射角发生变化,所述采样模块输出的信号强度随之变化,基于各个采样模块输出的信号强度可实现全方位太阳角的计算测量。
可选的,包含:第一模拟太阳角计,其安装于卫星本体对日方向中远离太阳的第一面的边缘,所述第一模拟太阳角计沿周向设置有多个采样模块;
第二模拟太阳角计,其安装于卫星本体对日方向中靠近太阳的第二面的边缘,所述第二模拟太阳角计沿周向设置有多个采样模块;
多个采样模块的采样角度覆盖全方位三维角度,随着采样模块的太阳入射角发生变化,所述采样模块输出的信号强度随之变化,基于各个采样模块输出的信号强度可实现全方位太阳角的计算测量。
可选的,卫星本体的对日方向平行于卫星本体的Z轴,所述第一模拟太阳角计安装于卫星本体的-Z面,第二模拟太阳角计安装于卫星本体的+Z面;
所述第一模拟太阳角计和第二模拟太阳角计的水平面坐标相同;
可选的,所述第一模拟太阳角计包含第一采样模块、第二采样模块、第三采样模块和第四采样模块,所述第四采样模块与第二采样模块相对设置,所述第四采样模块位于外侧,所述第二采样模块位于内侧,所述第一采样模块、第二采样模块、第三采样模块和第四采样模块的采样电压信号分别为V11、V12、V13、V14
所述第二模拟太阳角计包含第五采样模块、第六采样模块、第七采样模块和第八采样模块,所述第六采样模块与第八采样模块相对设置,所述第六采样模块位于外侧,所述第八采样模块位于内侧,所述第五采样模块、第六采样模块、第七采样模块和第八采样模块的采样电压信号分别为V21、V22、V23、V24
可选的,当所有采样模块的电池电压和大于0.1V时,为光照区;
当所有采样模块的电池电压和小于0.1V时,为阴影区;
当所有采样模块的电池电压和等于0.1V时,光照阴影标志不变更。
可选的,当所有采样模块的电池电压和小于0.1V时,太阳角置0;
当所有采样模块的电池电压和大于或等于0.1V时,全方位太阳角包含滚动太阳角Gdass、滚动伪速率Gdass_ω、俯仰太阳角Fyass和俯仰伪速率Fyass_ω,其中,
(1)滚动太阳角和滚动伪速率:
当V14≥V24且V14+V12≥V24+V22且V12≥V22按公式(1)计算滚动太阳角Gdass
Figure BDA0004014324920000021
当V24≥V14且V12≥V22且V12+V24≥V22+V14按公式(2)计算滚动太阳角Gdass
Figure BDA0004014324920000031
当V22+V24≥V14+V12且V22≥V12且V24≥V14且V22≥V24按公式(3)计算滚动太阳角Gdass
Figure BDA0004014324920000032
当V22+V24≥V14+V12且V22≥V12且V24≥V14且V22<V24按公式(4)计算滚动太阳角Gdass
Figure BDA0004014324920000033
当V22≥V12且V22+V14≥V12+V24且V14≥V24按公式(5)计算滚动太阳角Gdass
Figure BDA0004014324920000034
滚动太阳角Gdass进行滤波,即
Figure BDA0004014324920000035
其中,k-1表示上一节拍的滤波值,N为滤波系数,T为控制周期,模拟太阳角计的滚动伪速率Gdass_ω为
Gdass_ω=(Gdass-Gdss_f(k-1))/N   (7);
(2)俯仰太阳角Fyass和俯仰伪速率Fyass_ω:
当V13≥V21且V13+V11≥V23+V21且V11≥V23按公式(8)计算俯仰太阳角Fyass
Figure BDA0004014324920000036
当V21≥V13且V11≥V23且V11+V21≥V23+V13按公式(9)计算俯仰太阳角Fyass
Figure BDA0004014324920000037
当V23+V21≥V13+V11且V23≥V11且V21≥V13且V23≥V21按公式(10)计算俯仰太阳角Fyass
Figure BDA0004014324920000038
当V23+V21≥V13+V11且V23≥V11且V21≥V13且V23<V21按公式(11)计算俯仰太阳角Fyass
Figure BDA0004014324920000041
当V23≥V11且V23+V13≥V11+V21且V13≥V21按公式(12)计算俯仰太阳角Fyass
Figure BDA0004014324920000042
俯仰太阳角Fyass进行滤波,即
Figure BDA0004014324920000043
其中,k-1表示上一节拍的滤波值,N为滤波系数,T为控制周期;
模拟太阳角计的俯仰伪速率Fyass_ω为
Fyass_ω=(Fyass-Fyss_f(k-1))/N   (14)。
可选的,一种前述的全方位太阳角测量系统的太阳角测量方法,包含:
多台模拟太阳角计的各个采样模块进行采样并分别输出各自的信号强度;
基于各个采样模块输出的信号强度数据,进行全方位太阳角的计算测量。
可选的,当所有采样模块的电池电压和大于0.1V时,为光照区;
当所有采样模块的电池电压和小于0.1V时,为阴影区;
当所有采样模块的电池电压和等于0.1V时,光照阴影标志不变更。
可选的,所述全方位太阳角测量系统包含第一模拟太阳角计和第二模拟太阳角计;所述第一模拟太阳角计安装于卫星本体对日方向中远离太阳的第一面的边缘,所述第一模拟太阳角计沿周向设置有多个采样模块;所述第二模拟太阳角计安装于卫星本体对日方向中靠近太阳的第二面的边缘,所述第二模拟太阳角计沿周向设置有多个采样模块;
多个采样模块的采样角度覆盖全方位三维角度,随着采样模块的太阳入射角发生变化,所述采样模块输出的信号强度随之变化,基于各个采样模块输出的信号强度可实现全方位太阳角的计算测量。
可选的,当所有采样模块的电池电压和小于0.1V时,太阳角置0;
当所有采样模块的电池电压和大于或等于0.1V时,太阳角测量包含滚动太阳角Gdass、滚动伪速率Gdass_ω、俯仰太阳角Fyass和俯仰伪速率Fyass_ω,V11、V12、V13、V14分别为第一模拟太阳角计的第一采样模块、第二采样模块、第三采样模块和第四采样模块的采样电压信号,V21、V22、V23、V24分别为第二模拟太阳角计的第五采样模块、第六采样模块、第七采样模块和第八采样模块的采样电压信号,其中,
(1)滚动太阳角和滚动伪速率:
当V14≥V24且V14+V12≥V24+V22且V12≥V22按公式(1)计算滚动太阳角Gdass
Figure BDA0004014324920000051
当V24≥V14且V12≥V22且V12+V24≥V22+V14按公式(2)计算滚动太阳角Gdass
Figure BDA0004014324920000052
当V22+V24≥V14+V12且V22≥V12且V24≥V14且V22≥V24按公式(3)计算滚动太阳角Gdass
Figure BDA0004014324920000053
当V22+V24≥V14+V12且V22≥V12且V24≥V14且V22<V24按公式(4)计算滚动太阳角Gdass
Figure BDA0004014324920000054
当V22≥V12且V22+V14≥V12+V24且V14≥V24按公式(5)计算滚动太阳角Gdass
Figure BDA0004014324920000055
滚动太阳角Gdass进行滤波,即
Figure BDA0004014324920000056
其中,k-1表示上一节拍的滤波值,N为滤波系数,T为控制周期,
模拟太阳角计的滚动伪速率Gdass_ω为
Gdass_ω=(Gdass-Gdss_f(k-1))/N   (7);
(2)俯仰太阳角Fyass和俯仰伪速率Fyass_ω:
当V13≥V21且V13+V11≥V23+V21且V11≥V23按公式(8)计算俯仰太阳角Fyass
Figure BDA0004014324920000061
当V21≥V13且V11≥V23且V11+V21≥V23+V13按公式(9)计算俯仰太阳角Fyass
Figure BDA0004014324920000062
当V23+V21≥V13+V11且V23≥V11且V21≥V13且V23≥V21按公式(10)计算俯仰
太阳角Fyass
Figure BDA0004014324920000063
当V23+V21≥V13+V11且V23≥V11且V21≥V13且V23<V21按公式(11)计算俯仰
太阳角Fyass
Figure BDA0004014324920000064
当V23≥V11且V23+V13≥V11+V21且V13≥V21按公式(12)计算俯仰太阳角Fyass
Figure BDA0004014324920000065
俯仰太阳角Fyass进行滤波,即
Figure BDA0004014324920000066
其中,k-1表示上一节拍的滤波值,N为滤波系数,T为控制周期;
模拟太阳角计的俯仰伪速率Fyass_ω为
Fyass_ω=(Fyass-Fyss_f(k-1))/N   (14)。
本发明与现有技术相比具有以下优点:
本发明的一种全方位太阳角测量方法及其系统中,根据卫星对日定向需求情况以及对日方向配置多台模拟太阳角计(太阳敏感器),并设计其安装方位,根据各台模拟太阳角计测量信息,计算滚动轴和俯仰轴的太阳角及伪速率。本发明利用多台模拟太阳角计可以实现全方位的太阳角计算,不需要陀螺信息,综合实现了0-1式太阳敏感器和太阳角计功能,简化了对日定向所需单机配置,降低了系统成本。
附图说明
图1为本发明的一种卫星本体-Z面示意图;
图2为本发明的一种卫星本体+Z面示意图。
具体实施方式
以下结合附图,通过详细说明一个较佳的具体实施例,对本发明做进一步阐述。
基于前述问题,本发明提供了一种全方位太阳角测量系统,该系统包含多台模拟太阳角计,所述模拟太阳角计包含多个采样模块,多个采样模块的采样角度覆盖全方位三维角度,随着采样模块的太阳入射角发生变化,所述采样模块输出的信号强度随之变化,基于各个采样模块输出的信号强度可实现全方位太阳角的计算测量。
如图1和图2结合所示,为本实施例的一种全方位太阳角测量系统所在的卫星本体示意图,该系统具体包含第一模拟太阳角计(模拟太阳角计1)和第二模拟太阳角计(模拟太阳角计2),所述第一模拟太阳角计和第二模拟太阳角计均为金字塔型太阳敏感器,以便全方位捕获太阳角相关信息。所述第一模拟太阳角计安装于卫星本体对日方向中远离太阳的第一面的边缘,所述第一模拟太阳角计沿周向设置有多个采样模块。所述第二模拟太阳角计安装于卫星本体对日方向中靠近太阳的第二面的边缘,所述第二模拟太阳角计沿周向设置有多个采样模块。多个采样模块的采样角度覆盖全方位三维角度,随着采样模块的太阳入射角发生变化,所述采样模块输出的信号强度随之变化,各个采样模块反应不同分区的太阳角的变化,基于各个采样模块输出的信号强度可实现全方位太阳角的计算测量。
需要说明的是,所述模拟太阳角计并不仅限为上述金字塔型太阳敏感器,在其他实施例中,其也可为其他类型的结构,例如平面型太阳敏感器或其他立体结构型太阳敏感器。进一步的,该系统包含的模拟太阳角计的个数也不仅限于上述,本发明对其个数不做限制,只要可实现多个采样模块的采样角度覆盖全方位三维角度均可。示例地,在另一实施例中,该系统包含四个模拟太阳角计,所述模拟太阳角计为平面型太阳敏感器,各个模拟太阳角计的各个采样模块的采样角度覆盖全方位三维角度。
进一步的,星体坐标系的三维方向如图1和图2所示,在本实施例中,卫星本体的对日方向平行于卫星本体的Z轴,所述第一模拟太阳角计安装于卫星本体的-Z面即第一面,第二模拟太阳角计安装于卫星本体的+Z面即第二面。
在本实施例中,所述采样模块为电池模块,根据各个采样模块的采样电压值确定卫星光照阴影状态。进一步的,所述第一模拟太阳角计和第二模拟太阳角计的水平面坐标相同,即两者对应设置。所述第一模拟太阳角计包含第一采样模块、第二采样模块、第三采样模块和第四采样模块(分别为图1中模拟太阳角计1的1234),所述第四采样模块与第二采样模块相对设置,所述第四采样模块位于外侧,所述第二采样模块位于内侧,所述第一采样模块、第二采样模块、第三采样模块和第四采样模块的采样电压信号分别为V11、V12、V13、V14。所述第二模拟太阳角计包含第五采样模块、第六采样模块、第七采样模块和第八采样模块(分别为图2中模拟太阳角计2的1234),所述第六采样模块与第八采样模块相对设置,所述第六采样模块位于外侧,所述第八采样模块位于内侧,所述第五采样模块、第六采样模块、第七采样模块和第八采样模块的采样电压信号分别为V21、V22、V23、V24。可以理解的是,所述第一模拟太阳角计或第二模拟太阳角计不仅限于包含上述四个采样模块,在其他实施例中,还可以包含其他个数的采样模块,本发明对此不加以限制,只要可实现对周向不同方位的太阳入射角的测量即可。
当所有采样模块的电池电压和大于0.1V时,为光照区。当所有采样模块的电池电压和小于0.1V时,为阴影区。当所有采样模块的电池电压和等于0.1V时,光照阴影标志不变更。
在本实施例中,当所有采样模块的电池电压和小于0.1V时,太阳角置0。
进一步的,当所有采样模块的电池电压和大于或等于0.1V时,全方位太阳角包含滚动太阳角Gdass、滚动伪速率Gdass_ω、俯仰太阳角Fyass和俯仰伪速率Fyass_ω。
其中,(1)滚动太阳角和滚动伪速率的计算如下:
当V14≥V24且V14+V12≥V24+V22且V12≥V22按公式(1)计算滚动太阳角Gdass
Figure BDA0004014324920000091
当V24≥V14且V12≥V22且V12+V24≥V22+V14按公式(2)计算滚动太阳角Gdass
Figure BDA0004014324920000092
当V22+V24≥V14+V12且V22≥V12且V24≥V14且V22≥V24按公式(3)计算滚动太阳角Gdass
Figure BDA0004014324920000093
当V22+V24≥V14+V12且V22≥V12且V24≥V14且V22<V24按公式(4)计算滚动太阳角Gdass
Figure BDA0004014324920000094
当V22≥V12且V22+V14≥V12+V24且V14≥V24按公式(5)计算滚动太阳角Gdass
Figure BDA0004014324920000095
滚动太阳角Gdass进行滤波,即
Figure BDA0004014324920000096
其中,k-1表示上一节拍即上一采样周期的滤波值,N为滤波系数,T为控制周期,
模拟太阳角计(第一模拟太阳角计或第二模拟太阳角计)的滚动伪速率Gdass_ω为
Gdass_ω=(Gdass-Gdss_f(k-1))/N   (7)。
(2)俯仰太阳角Fyass和俯仰伪速率Fyass_ω的计算如下:
当V13≥V21且V13+V11≥V23+V21且V11≥V23按公式(8)计算俯仰太阳角Fyass
Figure BDA0004014324920000097
当V21≥V13且V11≥V23且V11+V21≥V23+V13按公式(9)计算俯仰太阳角Fyass
Figure BDA0004014324920000101
当V23+V21≥V13+V11且V23≥V11且V21≥V13且V23≥V21按公式(10)计算俯仰
太阳角Fyass
Figure BDA0004014324920000102
当V23+V21≥V13+V11且V23≥V11且V21≥V13且V23<V21按公式(11)计算俯仰
太阳角Fyass
Figure BDA0004014324920000103
当V23≥V11且V23+V13≥V11+V21且V13≥V21按公式(12)计算俯仰太阳角Fyass
Figure BDA0004014324920000104
俯仰太阳角Fyass进行滤波,即
Figure BDA0004014324920000105
其中,k-1表示上一节拍即上一周期的滤波值,N为滤波系数,T为控制周期;
模拟太阳角计(第一模拟太阳角计或第二模拟太阳角计)的俯仰伪速率Fyass_ω为
Fyass_ω=(Fyass-Fyss_f(k-1))/N   (14)。
基于同一发明构思,本发明还提供了一种前述全方位太阳角测量系统的太阳角测量方法,该方法包含:多台模拟太阳角计的各个采样模块进行采样并分别输出各自的信号强度;基于各个采样模块输出的信号强度数据,进行全方位太阳角的计算测量。
如前所述,当所有采样模块的电池电压和大于0.1V时,为光照区;当所有采样模块的电池电压和小于0.1V时,为阴影区;当所有采样模块的电池电压和等于0.1V时,光照阴影标志不变更。
在本实施例中,当所有采样模块的电池电压和小于0.1V时,太阳角置0。
进一步的,当所有采样模块的电池电压和大于或等于0.1V时,太阳角测量包含滚动太阳角Gdass、滚动伪速率Gdass_ω、俯仰太阳角Fyass和俯仰伪速率Fyass_ω,V11、V12、V13、V14分别为第一采样模块、第二采样模块、第三采样模块和第四采样模块的采样电压信号,V21、V22、V23、V24分别为第五采样模块、第六采样模块、第七采样模块和第八采样模块的采样电压信号。
其中,(1)滚动太阳角和滚动伪速率的计算如下:
当V14≥V24且V14+V12≥V24+V22且V12≥V22按公式(1)计算滚动太阳角Gdass
Figure BDA0004014324920000111
当V24≥V14且V12≥V22且V12+V24≥V22+V14按公式(2)计算滚动太阳角Gdass
Figure BDA0004014324920000112
当V22+V24≥V14+V12且V22≥V12且V24≥V14且V22≥V24按公式(3)计算滚动太阳角Gdass
Figure BDA0004014324920000113
当V22+V24≥V14+V12且V22≥V12且V24≥V14且V22<V24按公式(4)计算滚动太阳角Gdass
Figure BDA0004014324920000114
当V22≥V12且V22+V14≥V12+V24且V14≥V24按公式(5)计算滚动太阳角Gdass
Figure BDA0004014324920000115
滚动太阳角Gdass进行滤波,即
Figure BDA0004014324920000116
其中,k-1表示上一节拍即上一采样周期的滤波值,N为滤波系数,T为控制周期,
模拟太阳角计(第一模拟太阳角计或第二模拟太阳角计)的滚动伪速率Gdass_ω为
Gdass_ω=(Gdass-Gdss_f(k-1))/N   (7)。
(2)俯仰太阳角Fyass和俯仰伪速率Fyass_ω的计算如下:
当V13≥V21且V13+V11≥V23+V21且V11≥V23按公式(8)计算俯仰太阳角Fyass
Figure BDA0004014324920000121
当V21≥V13且V11≥V23且V11+V21≥V23+V13按公式(9)计算俯仰太阳角Fyass
Figure BDA0004014324920000122
当V23+V21≥V13+V11且V23≥V11且V21≥V13且V23≥V21按公式(10)计算俯仰太阳角Fyass
Figure BDA0004014324920000123
当V23+V21≥V13+V11且V23≥V11且V21≥V13且V23<V21按公式(11)计算俯仰太阳角Fyass
Figure BDA0004014324920000124
当V23≥V11且V23+V13≥V11+V21且V13≥V21按公式(12)计算俯仰太阳角Fyass
Figure BDA0004014324920000125
俯仰太阳角Fyass进行滤波,即
Figure BDA0004014324920000126
其中,k-1表示上一节拍的滤波值,N为滤波系数,T为控制周期;
模拟太阳角计的俯仰伪速率Fyass_ω为
Fyass_ω=(Fyass-Fyss_f(k-1))/N   (14)。
综上所述,本发明的一种全方位太阳角测量方法及其系统中,根据卫星对日定向需求情况以及对日方向配置多台模拟太阳角计(太阳敏感器),并设计其安装方位,根据多台模拟太阳角计测量信息,计算滚动轴和俯仰轴的太阳角及伪速率。本发明利用多台模拟太阳角计可以实现全方位的太阳角计算,不需要陀螺信息,综合实现了0-1式太阳敏感器和太阳角计功能,简化了对日定向所需单机配置,降低了系统成本。
尽管本发明的内容已经通过上述优选实施例作了详细介绍,但应当认识到上述的描述不应被认为是对本发明的限制。在本领域技术人员阅读了上述内容后,对于本发明的多种修改和替代都将是显而易见的。因此,本发明的保护范围应由所附的权利要求来限定。

Claims (10)

1.一种全方位太阳角测量系统,其特征在于,包含:
多台模拟太阳角计,所述模拟太阳角计包含多个采样模块,多个采样模块的采样角度覆盖全方位三维角度,随着采样模块的太阳入射角发生变化,所述采样模块输出的信号强度随之变化,基于各个采样模块输出的信号强度可实现全方位太阳角的计算测量。
2.如权利要求1所述的全方位太阳角测量系统,其特征在于,包含:
第一模拟太阳角计,其安装于卫星本体对日方向中远离太阳的第一面的边缘,所述第一模拟太阳角计沿周向设置有多个采样模块;
第二模拟太阳角计,其安装于卫星本体对日方向中靠近太阳的第二面的边缘,所述第二模拟太阳角计沿周向设置有多个采样模块;
多个采样模块的采样角度覆盖全方位三维角度,随着采样模块的太阳入射角发生变化,所述采样模块输出的信号强度随之变化,基于各个采样模块输出的信号强度可实现全方位太阳角的计算测量。
3.如权利要求2所述的全方位太阳角测量系统,其特征在于,
卫星本体的对日方向平行于卫星本体的Z轴,所述第一模拟太阳角计安装于卫星本体的-Z面,第二模拟太阳角计安装于卫星本体的+Z面;
所述第一模拟太阳角计和第二模拟太阳角计的水平面坐标相同。
4.如权利要求2所述的全方位太阳角测量系统,其特征在于,
所述第一模拟太阳角计包含第一采样模块、第二采样模块、第三采样模块和第四采样模块,所述第四采样模块与第二采样模块相对设置,所述第四采样模块位于外侧,所述第二采样模块位于内侧,所述第一采样模块、第二采样模块、第三采样模块和第四采样模块的采样电压信号分别为V11、V12、V13、V14
所述第二模拟太阳角计包含第五采样模块、第六采样模块、第七采样模块和第八采样模块,所述第六采样模块与第八采样模块相对设置,所述第六采样模块位于外侧,所述第八采样模块位于内侧,所述第五采样模块、第六采样模块、第七采样模块和第八采样模块的采样电压信号分别为V21、V22、V23、V24
5.如权利要求2或4所述的全方位太阳角测量系统,其特征在于,
当所有采样模块的电池电压和大于0.1V时,为光照区;
当所有采样模块的电池电压和小于0.1V时,为阴影区;
当所有采样模块的电池电压和等于0.1V时,光照阴影标志不变更。
6.如权利要求4所述的全方位太阳角测量系统,其特征在于,
当所有采样模块的电池电压和小于0.1V时,太阳角置0;
当所有采样模块的电池电压和大于或等于0.1V时,全方位太阳角包含滚动太阳角Gdass、滚动伪速率Gdass_ω、俯仰太阳角Fyass和俯仰伪速率Fyass_ω,其中,
(1)滚动太阳角和滚动伪速率:
当V14≥V24且V14+V12≥V24+V22且V12≥V22按公式(1)计算滚动太阳角Gdass
Figure FDA0004014324910000021
当V24≥V14且V12≥V22且V12+V24≥V22+V14按公式(2)计算滚动太阳角Gdass
Figure FDA0004014324910000022
当V22+V24≥V14+V12且V22≥V12且V24≥V14且V22≥V24按公式(3)计算滚动太阳角Gdass
Figure FDA0004014324910000023
当V22+V24≥V14+V12且V22≥V12且V24≥V14且V22<V24按公式(4)计算滚动太阳角Gdass
Figure FDA0004014324910000024
当V22≥V12且V22+V14≥V12+V24且V14≥V24按公式(5)计算滚动太阳角Gdass
Figure FDA0004014324910000025
滚动太阳角Gdass进行滤波,即
Figure FDA0004014324910000026
其中,k-1表示上一节拍的滤波值,N为滤波系数,T为控制周期,模拟太阳角计的滚动伪速率Gdass_ω为
Gdass_ω=(Gdass-Gdss_f(k-1))/N    (7);
(2)俯仰太阳角Fyass和俯仰伪速率Fyass_ω:
当V13≥V21且V13+V11≥V23+V21且V11≥V23按公式(8)计算俯仰太阳角Fyass
Figure FDA0004014324910000031
当V21≥V13且V11≥V23且V11+V21≥V23+V13按公式(9)计算俯仰太阳角Fyass
Figure FDA0004014324910000032
当V23+V21≥V13+V11且V23≥V11且V21≥V13且V23≥V21按公式(10)计算俯仰太阳角Fyass
Figure FDA0004014324910000033
当V23+V21≥V13+V11且V23≥V11且V21≥V13且V23<V21按公式(11)计算俯仰太阳角Fyass
Figure FDA0004014324910000034
当V23≥V11且V23+V13≥V11+V21且V13≥V21按公式(12)计算俯仰太阳角Fyass
Figure FDA0004014324910000035
俯仰太阳角Fyass进行滤波,即
Figure FDA0004014324910000036
其中,k-1表示上一节拍的滤波值,N为滤波系数,T为控制周期;
模拟太阳角计的俯仰伪速率Fyass_ω为
Fyass_ω=(Fyass-Fyss_f(k-1))/N    (14)。
7.一种如权利要求1~6任一项所述的全方位太阳角测量系统的太阳角测量方法,其特征在于,包含:
多台模拟太阳角计的各个采样模块进行采样并分别输出各自的信号强度;
基于各个采样模块输出的信号强度数据,进行全方位太阳角的计算测量。
8.如权利要求7所述的全方位太阳角测量系统的太阳角测量方法,其特征在于,
当所有采样模块的电池电压和大于0.1V时,为光照区;
当所有采样模块的电池电压和小于0.1V时,为阴影区;
当所有采样模块的电池电压和等于0.1V时,光照阴影标志不变更。
9.如权利要求7所述的全方位太阳角测量系统的太阳角测量方法,其特征在于,
所述全方位太阳角测量系统包含第一模拟太阳角计和第二模拟太阳角计;所述第一模拟太阳角计安装于卫星本体对日方向中远离太阳的第一面的边缘,所述第一模拟太阳角计沿周向设置有多个采样模块;所述第二模拟太阳角计安装于卫星本体对日方向中靠近太阳的第二面的边缘,所述第二模拟太阳角计沿周向设置有多个采样模块;
多个采样模块的采样角度覆盖全方位三维角度,随着采样模块的太阳入射角发生变化,所述采样模块输出的信号强度随之变化,基于各个采样模块输出的信号强度可实现全方位太阳角的计算测量。
10.如权利要求9所述的全方位太阳角测量系统的太阳角测量方法,其特征在于,
当所有采样模块的电池电压和小于0.1V时,太阳角置0;
当所有采样模块的电池电压和大于或等于0.1V时,太阳角测量包含滚动太阳角Gdass、滚动伪速率Gdass_ω、俯仰太阳角Fyass和俯仰伪速率Fyass_ω,V11、V12、V13、V14分别为第一模拟太阳角计的第一采样模块、第二采样模块、第三采样模块和第四采样模块的采样电压信号,V21、V22、V23、V24分别为第二模拟太阳角计的第五采样模块、第六采样模块、第七采样模块和第八采样模块的采样电压信号,其中,
(1)滚动太阳角和滚动伪速率:
当V14≥V24且V14+V12≥V24+V22且V12≥V22按公式(1)计算滚动太阳角Gdass
Figure FDA0004014324910000051
当V24≥V14且V12≥V22且V12+V24≥V22+V14按公式(2)计算滚动太阳角Gdass
Figure FDA0004014324910000052
当V22+V24≥V14+V12且V22≥V12且V24≥V14且V22≥V24按公式(3)计算滚动太阳角Gdass
Figure FDA0004014324910000053
当V22+V24≥V14+V12且V22≥V12且V24≥V14且V22<V24按公式(4)计算滚动太阳角Gdass
Figure FDA0004014324910000054
当V22≥V12且V22+V14≥V12+V24且V14≥V24按公式(5)计算滚动太阳角Gdass
Figure FDA0004014324910000055
滚动太阳角Gdass进行滤波,即
Figure FDA0004014324910000056
其中,k-1表示上一节拍的滤波值,N为滤波系数,T为控制周期,
模拟太阳角计的滚动伪速率Gdass_ω为
Gdass_ω=(Gdass-Gdss_f(k-1))/N    (7);
(2)俯仰太阳角Fyass和俯仰伪速率Fyass_ω:
当V13≥V21且V13+V11≥V23+V21且V11≥V23按公式(8)计算俯仰太阳角Fyass
Figure FDA0004014324910000057
当V21≥V13且V11≥V23且V11+V21≥V23+V13按公式(9)计算俯仰太阳角Fyass
Figure FDA0004014324910000058
当V23+V21≥V13+V11且V23≥V11且V21≥V13且V23≥V21按公式(10)计算俯仰太阳角Fyass
Figure FDA0004014324910000061
当V23+V21≥V13+V11且V23≥V11且V21≥V13且V23<V21按公式(11)计算俯仰太阳角Fyass
Figure FDA0004014324910000062
当V23≥V11且V23+V13≥V11+V21且V13≥V21按公式(12)计算俯仰太阳角Fyass
Figure FDA0004014324910000063
俯仰太阳角Fyass进行滤波,即
Figure FDA0004014324910000064
其中,k-1表示上一节拍的滤波值,N为滤波系数,T为控制周期;
模拟太阳角计的俯仰伪速率Fyass_ω为
Fyass_ω=(Fyass-Fyss_f(k-1))/N    (14)。
CN202211664741.0A 2022-12-23 2022-12-23 一种全方位太阳角测量方法及其系统 Pending CN116045897A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211664741.0A CN116045897A (zh) 2022-12-23 2022-12-23 一种全方位太阳角测量方法及其系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211664741.0A CN116045897A (zh) 2022-12-23 2022-12-23 一种全方位太阳角测量方法及其系统

Publications (1)

Publication Number Publication Date
CN116045897A true CN116045897A (zh) 2023-05-02

Family

ID=86112473

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211664741.0A Pending CN116045897A (zh) 2022-12-23 2022-12-23 一种全方位太阳角测量方法及其系统

Country Status (1)

Country Link
CN (1) CN116045897A (zh)

Similar Documents

Publication Publication Date Title
US9612157B2 (en) Apparatus and method for navigation
CN101556155B (zh) 小卫星姿态确定系统及方法
ITRM20120664A1 (it) Bussola solare elettronica ad alta precisione.
CN102914306B (zh) 一种双探头星敏感器及其设计方法
CN101865692B (zh) 偏振光栅导航传感器
KR20100119007A (ko) 태양의 위치를 계산하는 연산부를 구비한 태양광 추적장치
CN102927982A (zh) 双光谱的自主导航敏感器及其设计方法
CN102706363A (zh) 一种高精度星敏感器的精度测量方法
US5844232A (en) Slit-type sun sensor having a plurality of slits
CN108072372A (zh) 一种用于微小卫星的集成姿态测量系统
CN105486315A (zh) 遥感卫星对月绝对定标姿态调整方法
JPH05240655A (ja) 3軸安定宇宙船および太陽捕捉方法
CN116045897A (zh) 一种全方位太阳角测量方法及其系统
CN102749089B (zh) 一种三探头星敏感器姿态的确定方法
Liebe Solar compass chip
CN111912402A (zh) 基于地磁信息辅助gps的高旋转载体的测姿方法及装置
Heywood The computation of solar radiation intensities Part 2—Solar radiation on inclined surfaces
CN107631666B (zh) 一种基于地磁和太阳光角的弹体滚转角检测系统及方法
JP2014058176A (ja) 地球センサ及び地球センサを用いた地球エッジ検出方法
RU2354591C2 (ru) Система определения максимальной выходной мощности солнечных батарей космического аппарата
US9783323B2 (en) Ram angle and magnetic field sensor (RAMS)
He et al. A mimu/polarized camera/gnss integrated navigation algorithm for uav application
CN104715143A (zh) 基于卫星图像数据的临边像素点的地球切线高度估计方法
CN106595657B (zh) 轴对称飞行器姿态测量装置及其测量方法
JP2000234925A (ja) 姿勢検出装置

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination