• 四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 2019-05-13
  • 银保监会新规剑指大企业多头融资和过度融资 2019-05-12
  • 韩国再提4国联合申办世界杯 中国网友无视:我们自己来 2019-05-11
  • 中国人为什么一定要买房? 2019-05-11
  • 十九大精神进校园:风正扬帆当有为 勇做时代弄潮儿 2019-05-10
  • 粽叶飘香幸福邻里——廊坊市举办“我们的节日·端午”主题活动 2019-05-09
  • 太原设禁鸣路段 设备在测试中 2019-05-09
  • 拜耳医药保健有限公司获第十二届人民企业社会责任奖年度企业奖 2019-05-08
  • “港独”没出路!“梁天琦们”该醒醒了 2019-05-07
  • 陈卫平:中国文化内涵包含三方面 文化复兴表现在其中 2019-05-06
  • 人民日报客户端辟谣:“合成军装照”产品请放心使用 2019-05-05
  • 【十九大·理论新视野】为什么要“建设现代化经济体系”?   2019-05-04
  • 聚焦2017年乌鲁木齐市老城区改造提升工程 2019-05-04
  • 【专家谈】上合组织——构建区域命运共同体的有力实践者 2019-05-03
  • 【华商侃车NO.192】 亲!楼市火爆,别忘了买车位啊! 2019-05-03
    • / 11
    • 下载费用:30 金币  

    重庆时时彩是国外的吗: 一种折射星的模拟方法.pdf

    摘要
    申请专利号:

    重庆时时彩单双窍门 www.4mum.com.cn CN201410201313.3

    申请日:

    2014.05.14

    公开号:

    CN103968835A

    公开日:

    2014.08.06

    当前法律状态:

    授权

    有效性:

    有权

    法律详情: 授权|||实质审查的生效IPC(主分类):G01C 21/02申请日:20140514|||公开
    IPC分类号: G01C21/02 主分类号: G01C21/02
    申请人: 哈尔滨工程大学
    发明人: 钱华明; 孙龙; 钱林琛; 蔡佳楠; 黄蔚; 沈忱; 吴永慧; 李浩; 李通
    地址: 150001 黑龙江省哈尔滨市南岗区南通大街145号哈尔滨工程大学科技处知识产权办公室
    优先权:
    专利代理机构: 代理人:
    PDF完整版下载: PDF下载
    法律状态
    申请(专利)号:

    CN201410201313.3

    授权公告号:

    ||||||

    法律状态公告日:

    2017.02.15|||2014.09.03|||2014.08.06

    法律状态类型:

    授权|||实质审查的生效|||公开

    摘要

    本发明涉及一种折射星的模拟方法,其特征在于:步骤1:扫描星表,判断当前恒星是否为折射星;若判断为折射星,则进入步骤2,否则继续扫描;步骤2:计算星光折射角;步骤3:求解折射星视位置的赤经、赤纬;步骤4:判断折射星的折射星光是否被星敏感器捕捉到;若判断折射星的折射星光被星敏感器捕捉到,则进入步骤5,否则返回步骤1;步骤5:进行星点成像中心计算并生成模拟折射星图。

    权利要求书

    权利要求书
    1.  一种折射星的模拟方法,其特征在于:
    步骤1:扫描星表,判断当前恒星是否为折射星;若判断为折射星,则进入步骤2,否则继续扫描;
    步骤2:计算星光折射角;
    步骤3:求解折射星视位置的赤经、赤纬;
    步骤4:判断折射星的折射星光是否被星敏感器捕捉到;若判断折射星的折射星光被星敏感器捕捉到,则进入步骤5,否则返回步骤1;
    步骤5:进行星点成像中心计算并生成模拟折射星图。

    2.  根据权利要求1所述的折射星的模拟方法,其特征在于:步骤1中,符合下列方程的恒星判断为折射星,
    α-θR<arccos(r|r|&CenterDot;s)<β]]>
    式中,s为星表中记录的恒星星光的方向矢量;θR为折射高度为20km时对应的星光折射角,可通过大气折射模型计算得到;r=[x y z]为载体在地心赤道惯性坐标系的位置矢量,为其模长,代表载体到地心的距离;α和β可表示为,
    α=arcsin(ha+Rer);β=arcsin(hb+Rer);]]>
    式中,ha=20km,hb=50km分别为根据平流层厚度所选择的最小折射高度和最大折射高度;Re为地球半径。

    3.  根据权利要求2所述的折射星的模拟方法,其特征在于:步骤2中,星光折射角通过解如下方程组获得,
    ha=r2-u2+utanR-Reha=57.081+2.531e[0.981ln(R)-8.689]-6.441ln(R)]]>
    式中,ha为折射视高度;u=|r·u|=|xsx+ysy+zsz|,u=[sx sy sz]T为折射前的星光矢量,sx、sy、sz为恒星在天球坐标系的方向矢量,可通过查找星表得到;R为星光折射角。

    4.  根据权利要求3所述的折射星的模拟方法,其特征在于:步骤3中,恒星视位置的赤经、赤纬通过解如下方程组获得,
    P&CenterDot;u=cosRarccos(P&CenterDot;r|r|)=R+arccos(u&CenterDot;r|r|)]]>
    式中,P=[cosγ cosδ sinγ cosδ sinδ]T为折射星光在地心赤道坐标系的方向矢量,γ和δ分别为折射星视位置的赤经和赤纬。

    5.  根据权利要求4所述的折射星的模拟方法,其特征在于:步骤5中,具体包括如下步骤,
    步骤5.1:星等到灰度的模拟;
    将恒星在星敏感器中成像的亮度用灰度表示,在基本星表中星等的范围是-1至11之间的浮点型数,取0到6星等的恒星,令d为观测星的星等,g为所成像的灰度值,则模拟恒星在像平面的灰度为:
    g=100+20·(6-d)
    步骤5.2:星点成像中心模拟;
    通过如下公式计算折射星成像中心位置(xs,ys),
    xs=Nx2&CenterDot;1tan(θFOV/2)&CenterDot;cosδsin(γ-γ0)sinδsinδ0+cosδcosδ0cos(γ-γ0)ys=Ny2&CenterDot;1tan(θFOV/2)&CenterDot;sinδcosδ0-cosδsinδ0cos(γ-γ0)sinδsinδ0+cosδcosδ0cos(γ-γ0)]]>
    式中,Nx、Ny分别为星敏感器像平面两个方向的像素数目;θFOV为星敏感器的视场角;γ0、δ0分别星敏感器的视轴指向所对应的赤经、赤纬;γ和δ分别为折射星视位置的赤经和赤纬;
    步骤5.3:星象点模拟;
    星像能量分布公式如下,
    f(x,y)=g2πσ2exp[-(x-xs)2+(y-ys)22σ2]]]>
    式中,σ代表星点扩散半径表示星点弥散斑的大小,(x,y)为像平面的像元坐标,f(x,y)为像元(x,y)对应的灰度值;
    步骤5.4:星图生成;
    根据星敏感器的参数设定好星点扩散半径、背景灰度和背景噪声后生成的数字星图可表示为,
    P(x,y)=f(x,y)+B+NB(x,y)+N(x,y)
    式中,B为背景的灰度值,NB(x,y)为背景噪声,设置为均值为0方差为σB的高斯白噪声;N(x,y)为随机噪声,设置为均值为0,方差为σN的高斯白噪声。

    说明书

    说明书一种折射星的模拟方法
    技术领域
    本发明涉及一种折射星的模拟方法。
    背景技术
    随着载人航天、深空探测等空间技术的发展,各类航天器对自主导航能力的要求越来越高?;谛敲舾衅鞯奶煳牡己椒椒ㄓ捎诘己骄雀?、误差不随时间积累等特点成为了空间应用中的有效自主导航手段。在影响天文导航精度的因素中,地平的敏感精度是很重要的一点,根据星敏感器敏感地平方式的不同,可分为直接敏感地平和利用星光折射间接敏感地平两种方法。
    直接敏感地平法受地平仪精度的制约,导航精度较低。利用星光折射间接敏感地平的方法是一种低成本、高精度的自主导航方法。该方法的研究工作可以追溯到20世纪60年代,美国在实施Apollo计划的过程中,就对利用天体掩星、星光在大气中的折射、星光穿越大气时的衰减等实现自主导航的方案进行了研究.1975年由美国海军研究局和美国国防部高级研究计划局共同投资,麻省理工学院Draper实验室对星光折射/星光色散自主导航方案进行了研究和论证,结果表明在一个轨道周期可观测40颗折射星的理想条件下,导航精度可以达到100m。90年代初投入使用的MADAN(multi-mission attitude determination andautonomous navigation)导航系统(多任务姿态确定和自主导航系统)便利用了星光折射原理。二十世纪80年代初期,法国也进行了星光折射法自主导航的研究。1985年和1986年,CNES多次释放平流层气球对星光折射进行了实际测量,在此基础上,对大气折射的精确模型、测量方案、自然环境对系统观测的约束、误差分配和系统性能优化等方面进行了深入的分析和仿真试验,当时预计该系统导航精度为300m。
    星光折射间接敏感地平法精度高且成本低,是一种很有发展前途的导航方法,因此吸引了大量的学者对其进行研究。然而,由于航天实验费用巨大,往往不可 能进行实时的星空拍摄,因此目前的折射星导航的研究,大多数都是基础理论的研究,实验验证比较困难?;褂屑僖徊糠盅芯砍晒墙岷闲枪庹凵涞募负卧砗痛笃凵淠P椭苯蛹扑阈枪庹凵浣?,然后使用星光折射角来研究星光折射导航。这种方法直接使用了计算得到的星光折射角,没有对折射星图进行模拟,忽略了环境噪声、星敏感器噪声、星点位置噪声等众多噪声的影响,与真实情况差距较大,因此这些研究成果对于实际应用还有一定的距离。鉴于此,本发明提供了一种折射星的模拟方法,而且能够根据实际情况考虑星敏感器噪声、环境噪声的影响,为折射星导航的研究提供了一种更加符合实际情况且操作简单,成本低廉的方法。
    发明内容
    本发明目的在于提供一种不但能够模拟折射星在星敏感器像平面的位置,而且能够考虑星敏感器噪声,且操作简单,成本低廉的折射星的模拟方法。
    实现本发明目的技术方案:
    步骤1:扫描星表,判断当前恒星是否为折射星;若判断为折射星,则进入步骤2,否则继续扫描;
    步骤2:计算星光折射角;
    步骤3:求解折射星视位置的赤经、赤纬;
    步骤4:判断折射星的折射星光是否被星敏感器捕捉到;若判断折射星的折射星光被星敏感器捕捉到,则进入步骤5,否则返回步骤1;
    步骤5:进行星点成像中心计算并生成模拟折射星图。
    步骤1中,符合下列方程的恒星判断为折射星,
    α-θR<arccos(r|r|&CenterDot;s)<β]]>
    式中,s为星表中记录的恒星星光的方向矢量;θR为折射高度为20km时对应的星光折射角,可通过大气折射模型计算得到;r=[x y z]为载体在地心赤道惯性坐标系的位置矢量,为其模长,代表载体到地心的距离; α和β可表示为,
    α=arcsin(ha+Rer);β=arcsin(hb+Rer);]]>
    式中,ha=20km,hb=50km分别为根据平流层厚度所选择的最小折射高度和最大折射高度;Re为地球半径。
    步骤2中,星光折射角通过解如下方程组获得,
    ha=r2-u2+utanR-Reha=57.081+2.531e[0.981ln(R)-8.689]-6.441ln(R)]]>
    式中,ha为折射视高度;u=|r·u|=|xsx+ysy+zsz|,u=[sx sy sz]T为折射前的星光矢量,sx、sy、sz为恒星在天球坐标系的方向矢量,可通过查找星表得到;R为星光折射角。
    步骤3中,恒星视位置的赤经、赤纬通过解如下方程组获得,
    P&CenterDot;u=cosRarccos(P&CenterDot;r|r|)=R+arccos(u&CenterDot;r|r|)]]>
    式中,P=[cosγ cosδ sinγ cosδ sinδ]T为折射星光在地心赤道坐标系的方向矢量,γ和δ分别为折射星视位置的赤经和赤纬。
    步骤5中,具体包括如下步骤,
    步骤5.1:星等到灰度的模拟;
    将恒星在星敏感器中成像的亮度用灰度表示,在基本星表中星等的范围是-1至11之间的浮点型数,取0到6星等的恒星,令d为观测星的星等,g为所成像的灰度值,则模拟恒星在像平面的灰度为:
    g=100+20·(6-d)
    步骤5.2:星点成像中心模拟;
    通过如下公式计算折射星成像中心位置(xs,ys),
    xs=Nx2&CenterDot;1tan(θFOV/2)&CenterDot;cosδsin(γ-γ0)sinδsinδ0+cosδcosδ0cos(γ-γ0)ys=Ny2&CenterDot;1tan(θFOV/2)&CenterDot;sinδcosδ0-cosδsinδ0cos(γ-γ0)sinδsinδ0+cosδcosδ0cos(γ-γ0)]]>
    式中,Nx、Ny分别为星敏感器像平面两个方向的像素数目;θFOV为星敏 感器的视场角;γ0、δ0分别星敏感器的视轴指向所对应的赤经、赤纬;γ和δ分别为折射星视位置的赤经和赤纬;
    步骤5.3:星象点模拟;
    星像能量分布公式如下,
    f(x,y)=g2πσ2exp[-(x-xs)2+(y-ys)22σ2]]]>
    式中,σ代表星点扩散半径表示星点弥散斑的大小,(x,y)为像平面的像元坐标,f(x,y)为像元(x,y)对应的灰度值;
    步骤5.4:星图生成;
    根据星敏感器的参数设定好星点扩散半径、背景灰度和背景噪声后生成的数字星图可表示为,
    P(x,y)=f(x,y)+B+NB(x,y)+N(x,y)
    式中,B为背景的灰度值,NB(x,y)为背景噪声,设置为均值为0方差为σB的高斯白噪声;N(x,y)为随机噪声,设置为均值为0,方差为σN的高斯白噪声。
    本发明具有的有益效果:
    本发明首先扫描星表,判断当前恒星是否为折射星;若判断为折射星,则计算星光折射角;求解折射星视位置的赤经、赤纬;然后,判断折射星的折射星光是否被星敏感器捕捉到;若判断折射星的折射星光被星敏感器捕捉到,则进行星点成像中心计算并生成模拟折射星图。本发明不但能够模拟折射星在星敏感器像平面的位置,而且能够根据实际情况考虑星敏感器噪声、环境噪声的影响,使得折射星的模拟和折射角的计算更加符合实际,并且操作简单,成本低廉。
    附图说明
    图1是折射星的判断原理图;
    图2是星光折射原理图;
    图3是恒星矢量与主轴平面夹角示意图。
    具体实施方式
    步骤1:扫描星表,判断当前恒星是否为折射星;若判断为折射星,则进入步骤2,否则继续扫描;
    图1中:ha-最小折射高度20km,hb-最大折射高度50km,r-载体位置矢量,Re-地球半径s1、s2-恒星星光的方向矢量,θR-折射高度为20km时的星光折射角。
    在使用折射星导航时,根据平流层的厚度折射高度一般选取20km-50km,即图1中ha=20km,hb=50km;假设某恒星星光的方向矢量为s,则由图1得符合下列方程的恒星可以被选为折射星
    α-θR<arccos(r|r|&CenterDot;s)<β---(1)]]>
    其中,θR为折射高度为20km时对应的星光折射角,可通过大气折射模型计算得到;r=[x y z]为载体在地心赤道惯性坐标系的位置矢量,为其模长,代表载体到地心的距离;α和β如图1所示,可表示为
    α=arcsin(ha+Rer);β=arcsin(hb+Rer);]]>
    步骤2:计算星光折射角;
    图2中,R-星光折射角,ha-折射视高度,hg-折射高度,O-地心,u-星光折射前的方向矢量
    若某颗恒星利用式(1)判断为折射星之后,就计算该恒星的星光折射角R;如图2所示,根据折射高度与星光折射角之间的几何关系可得如下方程
    ha=r2-u2+utanR-Re---(2)]]>
    其中,u=|r·u|=|xsx+ysy+zsz|,u=[sx sy sz]T为折射前的星光矢量,sx、sy、sz为恒星在天球坐标系的方向矢量,可通过查找星表得到;
    同时根据大气折射模型也可以得到折射视高度与星光折射角之间的关系
    ha=57.081+2.531e[0.981ln(R)-8.689]-6.441ln(R) (3)
    对于公式(2),在进行折射星模拟时可以通过载体的真实轨迹得到载体的位置矢量r,同时从星表中可以得到u,地球半径Re已知,所以在公式(2)中只有折射视高度ha与星光折射角R未知;将其与公式(3)结合通过解方程组就可以得到星光折射角R;
    步骤3:求解折射星视位置的赤经、赤纬;
    假设折射星视位置的赤经、赤纬为γ和δ,则折射星光在地心赤道坐标系的方向矢量P可以表示如下
    P=[cosγ cosδ sinγ cosδ sinδ]T (4)
    星光折射角为R已经在上一步中得到,则有
    P·u=cosR (5)
    同时根据图2中的几何关系得
    arccos(R&CenterDot;r|r|)=R+arccos(u&CenterDot;r|r|)---(6)]]>
    将公式(5)、(6)结合通过解方程组就可以得到折射星视位置的赤经γ、赤纬δ;
    步骤4:判断折射星的折射星光是否被星敏感器捕捉到;若判断折射星的折射星光被星敏感器捕捉到,则进入步骤5,否则返回步骤1;
    图3中:P-恒星星光方向矢量,εx、εy、εz-分别为恒星星光方向矢量与主轴平面的夹角,Oxyz-地心赤道坐标系。
    如图3所示,折射星光矢量P与主轴平面的夹角分别为εx、εy、εz,且有
    &epsiv;x&epsiv;y&epsiv;z=arcsin(cosγcosδ)arcsin(sinγcosδ)δ---(7)]]>
    根据载体的真实轨??傻玫皆靥宓恼媸底颂?,通过载体与星敏感器的安装矩阵可以得到星敏感器的视轴指向所对应的赤经γ0、赤纬δ0,,根据公式(7)得到视轴指向与主平面的夹角εx0、εy0、εz0,则符合下述公式的折射星光可以被视场大小为θFOV的星敏感器捕获到
    |&epsiv;x-&epsiv;x0|θFOV/2|&epsiv;y-&epsiv;y0|θFOV/2|&epsiv;z-&epsiv;z0|θFOV/2---(8)]]>
    步骤5:进行星点成像中心计算并生成模拟折射星图。
    折射星图模拟的的具体步骤如下:
    步骤5.1:星等到灰度的模拟
    一般将恒星在星敏感器中成像的亮度用灰度表示,恒星的星等越小,对应的 亮度越强,相应的灰度值越大,在基本星表中星等的范围是-1至11之间的浮点型数,取0到6星等的恒星,令d为观测星的星等,g为所成像的灰度值,则模拟恒星在像平面的灰度为
    g=100+20·(6-d) (9)
    步骤5.2:星点成像中心模拟
    若判断折射星光被星敏感器捕获到后,使用如下公式计算折射星成像中心位置(xs,ys)
    xs=Nx2&CenterDot;1tan(θFOV/2)&CenterDot;cosδsin(γ-γ0)sinδsinδ0+cosδcosδ0cos(γ-γ0)ys=Ny2&CenterDot;1tan(θFOV/2)&CenterDot;sinδcosδ0-cosδsinδ0cos(γ-γ0)sinδsinδ0+cosδcosδ0cos(γ-γ0)---(10)]]>
    其中,Nx、Ny分别为星敏感器像平面两个方向的像素数目。
    步骤5.3:星象点模拟
    对于CCD成像装置,恒星为点目标,像点为CCD相机焦平面上的光学系统点扩散函数。
    一般而言,光学系统不是理想系统,星像能量分布近似为如下二维高斯分布:
    f(x,y)=g2πσ2exp[-(x-xs)2+(y-ys)22σ2]---(11)]]>
    其中,σ代表星点扩散半径表示星点弥散斑的大小,(x,y)为像平面的像元坐标,f(x,y)为像元(x,y)对应的灰度值;
    步骤5.4:星图生成
    根据星敏感器的参数设定好星点扩散半径、背景灰度和背景噪声后生成的数字星图可表示为;
    P(x,y)=f(x,y)+B+NB(x,y)+N(x,y) (12)
    其中,B为背景的灰度值,NB(x,y)为背景噪声一般设置为均值为0方差为σB的高斯白噪声。另外出背景噪声外,还有成像器件噪声和散粒噪声等,它们也可以用高斯分布的随机噪声来表示,即N(x,y)也可以设置为均值为0,方差为σN的 高斯白噪声。

    关 键 词:
    一种 折射 模拟 方法
      专利查询网所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
    关于本文
    本文标题:一种折射星的模拟方法.pdf
    链接地址://www.4mum.com.cn/p-6142718.html
    关于我们 - 网站声明 - 网站地图 - 资源地图 - 友情链接 - 网站客服客服 - 联系我们

    [email protected] 2017-2018 www.4mum.com.cn网站版权所有
    经营许可证编号:粤ICP备17046363号-1 
     


    收起
    展开
  • 四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 2019-05-13
  • 银保监会新规剑指大企业多头融资和过度融资 2019-05-12
  • 韩国再提4国联合申办世界杯 中国网友无视:我们自己来 2019-05-11
  • 中国人为什么一定要买房? 2019-05-11
  • 十九大精神进校园:风正扬帆当有为 勇做时代弄潮儿 2019-05-10
  • 粽叶飘香幸福邻里——廊坊市举办“我们的节日·端午”主题活动 2019-05-09
  • 太原设禁鸣路段 设备在测试中 2019-05-09
  • 拜耳医药保健有限公司获第十二届人民企业社会责任奖年度企业奖 2019-05-08
  • “港独”没出路!“梁天琦们”该醒醒了 2019-05-07
  • 陈卫平:中国文化内涵包含三方面 文化复兴表现在其中 2019-05-06
  • 人民日报客户端辟谣:“合成军装照”产品请放心使用 2019-05-05
  • 【十九大·理论新视野】为什么要“建设现代化经济体系”?   2019-05-04
  • 聚焦2017年乌鲁木齐市老城区改造提升工程 2019-05-04
  • 【专家谈】上合组织——构建区域命运共同体的有力实践者 2019-05-03
  • 【华商侃车NO.192】 亲!楼市火爆,别忘了买车位啊! 2019-05-03