CN107132155A - 基于实测大气密度的短时大气模型修正方法 - Google Patents
基于实测大气密度的短时大气模型修正方法 Download PDFInfo
- Publication number
- CN107132155A CN107132155A CN201710299136.0A CN201710299136A CN107132155A CN 107132155 A CN107132155 A CN 107132155A CN 201710299136 A CN201710299136 A CN 201710299136A CN 107132155 A CN107132155 A CN 107132155A
- Authority
- CN
- China
- Prior art keywords
- atmospheric
- density
- actual measurement
- models
- value
- 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
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N9/00—Investigating density or specific gravity of materials; Analysing materials by determining density or specific gravity
- G01N9/36—Analysing materials by measuring the density or specific gravity, e.g. determining quantity of moisture
Abstract
本发明提供了一种基于实测大气密度的短时大气模型修正方法,选定大气密度模型,采集计算设定时长的实测大气密度,分别对大气模型密度值和实测密度值进行线性拟合,得到实测大气密度拟合曲线和模型密度拟合曲线,计算实测大气密度拟合曲线和模型密度拟合曲线的夹角,然后对模型密度拟合曲线逐点进行旋转,得到旋转后的时间值t′和旋转后的选定大气模型密度值。本发明将实测大气密度值引入轨道计算,较好的提高了航天器轨道预报精度。
Description
技术领域
本发明属于航天测量与控制领域,涉及一种大气模型修正方法。
背景技术
近地轨道航天器运动受大气阻力摄动影响非常明显,大气阻力摄动的决定性因素就是大气密度值,大气密度的不均匀性导致了航天器在不同飞行高度受到不同的大气阻力。目前用作轨道计算的大气密度值都是通过模型计算得到,常用的大气模型有DTM、JACHIA、MSIS等系列。通常意义上,大气密度计算结果误差一般在10%到20%,在太阳活动活跃期和地磁暴期间,大气密度计算误差甚至会增加到一倍或更大,这严重影响了在轨运行航天器的轨道计算精度。通常,航天器上装载有大气密度探测器,探测器主要部件包括:取样室、输出电路、感测区等,这样就可以测得取样室内的温度和压力,进而计算出轨道空间的大气密度,经测量计算有助于对大气密度变化规律的直接把握。目前的大气密度探测结果主要用于研究空间大气成分或者环境演化分析,将探测结果用于航天器轨道计算及对现有大气模型进行修正的方法并不多。
发明内容
为了克服现有技术的不足,本发明提供一种基于实测大气密度的短时大气模型修正方法,减少理论大气密度模型计算带来的误差,提高短时低轨道航天器的预报精度(太阳活动和地磁活动均处于平静期)。
本发明解决其技术问题所采用的技术方案包括以下步骤:
(1)选定大气密度模型,采集计算设定时长的实测大气密度,分别对大气模型密度值和实测密度值进行线性拟合,令得到实测大气密度拟合曲线和模型密度拟合曲线,其中,ρmod_fit为模型密度线性拟合值,ρreal_fit为实测大气密度线性拟合值,a1和b1是MSIS2000大气模型线性拟合的常数项和一次项系数,a2和b2是实测大气密度线性拟合的常数项和一次项系数,t为时间;
(2)计算实测大气密度拟合曲线和模型密度拟合曲线的夹角然后对模型密度拟合曲线逐点进行旋转,得到旋转后的时间值t′和旋转后的选定大气模型密度值其中,ρmod为选定大气模型密度值,M(θ)为旋转矩阵,
本发明的有益效果是:将实测大气密度值引入轨道计算,构建了一种通过拟合、旋转实现大气密度模型值和实测值贴合的短时大气模型修正方法,通过算例应用和分析表明,针对某低轨道航天器,与传统的理论模型相比,该短时大气模型修正方法较好的提高了航天器轨道预报精度。
具体实施方式
下面结合实施例对本发明进一步说明,本发明包括但不仅限于下述实施例。
通常意义上讲,在太阳活动和地磁活动平静期,可以认为实测大气密度(考虑实测密度准确的情况下)比模型计算的大气密度更为贴近实际情况。主要是基于一段时间已知的实测大气密度,寻找已有大气密度模型同实测密度模型的关系,建立短时修正模型。这里,我们选取MSIS2000大气密度模型,找寻一种模型匹配原则以满足密度模型与实测密度的贴近,具体方法如下:
(1)分别对大气模型密度值和实测密度值进行线性拟合:令ρmod为MSIS2000大气模型密度值,ρmod_fit为模型密度线性拟合值,ρreal_fit为实测大气密度线性拟合值,a1和b1是MSIS2000大气模型线性拟合的常数项和一次项系数,其中a1的值域范围与密度值相当,b1直线斜率值域为[-∞,+∞],t为时间;a2和b2分别是实测大气密度线性拟合的常数项和一次项系数,意义同a1和b1。
(2)计算出实测大气密度拟合曲线和模型密度拟合曲线的夹角θ,然后逐点进行旋转,目的是将模型密度贴近实测大气密度。
其中,t′是旋转后的时间值,ρ′mod为MSIS2000大气模型密度旋转后的值。M(θ)为旋转矩阵,
这样,经过对密度模型值的拟合、旋转,就完成模型同实测密度的近似匹配。同时,这种匹配方法也适用于只有部分实测大气密度数据支撑的短时大气模型修正计算。
下面以具体实例说明,基于以上的模型修正方法,计算了某在轨600km高度航天器四个数据弧段的修正模型轨道预报(3天)。并同时分析了DTM2000、MSIS2000模型预报3天的精度(基准轨道为3天GPS数据定轨结果)。
由步骤(1)和(2)计算的实测密度拟合曲线与模型密度拟合曲线夹角θ见表1。
基于实测大气密度的短时大气模型修正方法计算的大气模型修正值用于航天器轨道预报结果见表2。
表1修正模型相关参数取值
表2给出不同大气密度模型预报3天RMS统计:
表2不同大气密度模型预报3天统计结果
从上表可以看出,利用修正的大气模型进行轨道预报,其预报精度大多数情况是最好的。该实测大气密度短时修正模型与传统的理论模型相比,对低轨航天器轨道预报精度有较大的提升。
Claims (1)
1.一种基于实测大气密度的短时大气模型修正方法,其特征在于包括下述步骤:
(1)选定大气密度模型,采集计算设定时长的实测大气密度,分别对大气模型密度值和实测密度值进行线性拟合,令得到实测大气密度拟合曲线和模型密度拟合曲线,其中,ρmod_fit为模型密度线性拟合值,ρreal_fit为实测大气密度线性拟合值,a1和b1是MSIS2000大气模型线性拟合的常数项和一次项系数,a2和b2是实测大气密度线性拟合的常数项和一次项系数,t为时间;
(2)计算实测大气密度拟合曲线和模型密度拟合曲线的夹角然后对模型密度拟合曲线逐点进行旋转,得到旋转后的时间值t′和旋转后的选定大气模型密度值ρ′mod,其中,ρmod为选定大气模型密度值,M(θ)为旋转矩阵,
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710299136.0A CN107132155A (zh) | 2017-05-02 | 2017-05-02 | 基于实测大气密度的短时大气模型修正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710299136.0A CN107132155A (zh) | 2017-05-02 | 2017-05-02 | 基于实测大气密度的短时大气模型修正方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107132155A true CN107132155A (zh) | 2017-09-05 |
Family
ID=59716478
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710299136.0A Pending CN107132155A (zh) | 2017-05-02 | 2017-05-02 | 基于实测大气密度的短时大气模型修正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107132155A (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107958102A (zh) * | 2017-11-07 | 2018-04-24 | 北京临近空间飞行器系统工程研究所 | 一种用于高超声速气动热预测的偏差大气参数确定方法 |
CN109033026A (zh) * | 2018-07-23 | 2018-12-18 | 中国人民解放军63920部队 | 一种大气密度探测数据的标校方法和设备 |
CN109323698A (zh) * | 2018-12-03 | 2019-02-12 | 西安四方星途测控技术有限公司 | 一种空间目标陨落多模型跟踪引导技术 |
CN112711736A (zh) * | 2021-01-28 | 2021-04-27 | 中国人民解放军63920部队 | 大气密度探测数据的标校方法、装置、存储介质和处理器 |
CN114180103A (zh) * | 2022-02-16 | 2022-03-15 | 中国人民解放军32035部队 | 一种基于空间环境变化的大气摄动参数动态修正方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103708045A (zh) * | 2014-01-16 | 2014-04-09 | 中国人民解放军国防科学技术大学 | 一种探月飞船跳跃式再入的在线参数辨识方法 |
CN103793618A (zh) * | 2014-02-28 | 2014-05-14 | 北京航天飞行控制中心 | 一种热层大气密度修正方法 |
CN204789249U (zh) * | 2015-07-15 | 2015-11-18 | 中国工程物理研究院材料研究所 | 一种高准确度β射线法在线测量大气颗粒物浓度的装置 |
CN106525651A (zh) * | 2016-10-24 | 2017-03-22 | 中国科学院国家空间科学中心 | 基于x射线掩日观测反演临近空间大气密度的方法 |
-
2017
- 2017-05-02 CN CN201710299136.0A patent/CN107132155A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103708045A (zh) * | 2014-01-16 | 2014-04-09 | 中国人民解放军国防科学技术大学 | 一种探月飞船跳跃式再入的在线参数辨识方法 |
CN103793618A (zh) * | 2014-02-28 | 2014-05-14 | 北京航天飞行控制中心 | 一种热层大气密度修正方法 |
CN204789249U (zh) * | 2015-07-15 | 2015-11-18 | 中国工程物理研究院材料研究所 | 一种高准确度β射线法在线测量大气颗粒物浓度的装置 |
CN106525651A (zh) * | 2016-10-24 | 2017-03-22 | 中国科学院国家空间科学中心 | 基于x射线掩日观测反演临近空间大气密度的方法 |
Non-Patent Citations (2)
Title |
---|
杜凯 等: "基于实测大气密度的航天器轨道计算技术", 《第七届中国卫星导航学术年会论文集》 * |
杜凯 等: "基于实测数据的大气密度模型修正方法", 《工程技术》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107958102A (zh) * | 2017-11-07 | 2018-04-24 | 北京临近空间飞行器系统工程研究所 | 一种用于高超声速气动热预测的偏差大气参数确定方法 |
CN107958102B (zh) * | 2017-11-07 | 2021-06-11 | 北京临近空间飞行器系统工程研究所 | 一种用于高超声速气动热预测的偏差大气参数确定方法 |
CN109033026A (zh) * | 2018-07-23 | 2018-12-18 | 中国人民解放军63920部队 | 一种大气密度探测数据的标校方法和设备 |
CN109323698A (zh) * | 2018-12-03 | 2019-02-12 | 西安四方星途测控技术有限公司 | 一种空间目标陨落多模型跟踪引导技术 |
CN109323698B (zh) * | 2018-12-03 | 2021-05-11 | 中科星图(西安)测控技术有限公司 | 一种空间目标陨落多模型跟踪引导方法 |
CN112711736A (zh) * | 2021-01-28 | 2021-04-27 | 中国人民解放军63920部队 | 大气密度探测数据的标校方法、装置、存储介质和处理器 |
CN114180103A (zh) * | 2022-02-16 | 2022-03-15 | 中国人民解放军32035部队 | 一种基于空间环境变化的大气摄动参数动态修正方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107132155A (zh) | 基于实测大气密度的短时大气模型修正方法 | |
CN103323888B (zh) | Gnss大气探测数据中对流层延迟误差的消除方法 | |
EP3364151B1 (en) | Precision calibration method for attitude measurement system | |
CN107870001B (zh) | 一种基于椭球拟合的磁力计校正方法 | |
CN107024674A (zh) | 一种基于递推最小二乘法的磁强计现场快速标定方法 | |
Niwa et al. | A three-dimensional icosahedral grid advection scheme preserving monotonicity and consistency with continuity for atmospheric tracer transport | |
CN113569201B (zh) | 地磁Ap指数预报方法、装置及电子设备 | |
CN114912301B (zh) | 一种低速风洞全机模型测力试验数据处理与修正系统 | |
Liléo et al. | Investigation on the use of NCEP/NCAR, MERRA and NCEP/CFSR reanalysis data in wind resource analysis | |
CN108154271A (zh) | 一种基于空间相关性和曲面拟合的地面气温质量控制方法 | |
Martin et al. | On the use of the synoptic vertical velocity in trajectory model: validation by geochemical tracers | |
CN106885892A (zh) | 对森林生态系统蒸散量拆分的方法及装置 | |
CN116297068A (zh) | 基于地表反射率优化的气溶胶光学厚度反演方法及系统 | |
CN109884666B (zh) | 一种基于数据同化技术的对流层延迟改正方法 | |
CN111765879A (zh) | 一种三维电子罗盘装置及实用校准方法 | |
Yang et al. | A brief introduction to BNU-HESM1. 0 and its earth surface temperature simulations | |
CN110110482A (zh) | 探空特性层温度偏差订正及效果检验算法 | |
CN111126466B (zh) | 一种多源pwv数据融合方法 | |
CN105158137B (zh) | 一种基于最小二乘支持向量机的水松纸透气度检测方法 | |
CN113991711B (zh) | 一种光伏电站储能系统容量配置方法 | |
Zhang et al. | Particle filter with unknown statistics to estimate liquid level in the silicon single-crystal growth | |
CN106483079A (zh) | 水汽中氢氧稳定同位素比率的测量数据的处理方法和装置 | |
Wang et al. | Testing and improving the performance of the Common Land Model: A case study for the Gobi landscape | |
CN110705924B (zh) | 一种基于风向扇区的测风塔测风数据处理方法及装置 | |
CN114880933A (zh) | 一种基于再分析资料的无探空站点地基微波辐射计大气温湿廓线反演方法及系统 |
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 | ||
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20170905 |