正常恒星的演化模型计算

本文提供了正常恒星演化模型的一种计算方法.它对常用的计算方法作了改进.计算程序用Algol-60语言编成,并在国产TQ-16型扩体的计算机上对质量为2.82M_⊙的恒星进行了演化模型计算.计算结果给出的各种演化特性证明这种方法和计算程序都是正确可用的.     

Hansen系数及其导数的直接计算方法

回顾总结了7种Hansen系数及其导数的直接计算方法,比较分析了这些方法的计算效率和计算稳定性.研究表明:Hansen系数的递推关系可以用来判别计算结果的稳定性.最后指出, Wnuk方法(双精度计算)和McClain方法(4精度计算)是稳定的,可以用来计算人造卫星轨道摄动.由于大多数人造卫星采用小偏心率轨道,需要计算无奇点摄动,推荐使用McClain方法1 (4精度计算).     

频率稳定度实时评估系统的设计与实现

为了实时计算原子钟的频率稳定度,设计了一种频率稳定度实时评估系统.利用现有理论基础(稳定度测量原理及其算法等)以及软件技术开发了一种实时计算并以图形显示原子钟稳定度的测量分析软件.该软件使用VC++的多媒体定时器实现数据的实时不间断采集及计算,同时采用双缓存绘图技术显示计算结果.测试结果表明该软件工作稳定,计算结果可靠、可信.     

1965年11月23日的日环食

越南最近所能看到的日食是1965年11月23日的日环食,这一次日环食经过越南境内的路线作者曾进行了计算.上次计算是利用纽康(S.Newcomb)的日食表,用它算得的日食路线只是近似的数值.现在利用更精确的数据进一步做精密而详细的计算.计算太阳月亮坐标是根据纽康的太阳表和布朗的月亮表.布朗为了使计算得的月亮位置能与观测到的位置更好地符合,在月亮表中加入了经验项     

一种适用于大偏心率轨道密集星历精密计算的快速处理方法

空间目标探测和编目是航天器安全的重要课题, 空间目标数量巨大, 造成轨道计算的工作量十分繁重. 分析方法虽然计算速度快, 却不能适应高精度观测资料的处理工作, 因而数值方法将成为目标轨道计算的重要方法. 空间目标编目工作普遍涉及密集星历的产生和计算问题, 针对这一问题的大偏心率轨道数值计算目前尚未形成兼具计算精度和计算效率的有效技术手段, 难以满足编目工作的处理要求. 在此背景下, 提出一种适用于大偏心率轨道密集星历精密计算的快速处理方法, 并通过数值实验对模型参数进行优选, 最后通过计算实例证实了该计算方案的优越性.     

一种初轨计算的稳健方法

光学测角资料的初轨计算在空间目标搜索发现中具有重要作用,当观测资料存在野值时,基于最小二乘的经典初轨计算方法不够稳健.采用最小一乘方法建立了一种初轨计算的稳健方法,方法将初轨计算问题转换为线性规划问题求解,并通过bootstrap方法给出估计精度.数值计算结果表明方法稳健有效,并具有较高的崩溃点.     

计算倾角函数的改进Gooding方法

将Gooding方法中的Vklm(I)改写为便于计算的形式,并利用标准递推过程计算Aklm(I).在规定k和l同奇偶的前提下改写了Gooding的计算程序,使程序缩短了一半,提高了程序计算效率和可读性,计算时间缩短了41％,计算精度和稳定性也略有提高.     

倾角函数的四精度计算试验

倾角函数是天体力学分析理论中一种常用的函数.当把摄动方程展开成时间和根数的形式时需要用到.历史上提出了很多经典的倾角函数递推算法,并在双精度平台下开发了Fortran程序.进行了1次四精度计算倾角函数的试验,结果表明:L平面递推方法的四精度计算精度可达10-22,计算速度比双精度Jacobi方法快6倍.     

地月系潮汐演化长期趋势的估计

根据角动量守恒原理,计算了地月系经潮汐演化到达平衡状态时的旋转周期和地月距离.并根据当前与平衡状态时地月系的总能量差,计算了到达平衡状态的时间.进而估计了地月距离变化和地球自转速率变化的长期趋势.     

Kerr黑洞的角动量

本文指出,计算Kerr黑洞角动量应考虑散射电子的逆向吸积效应.同时提出用顺、逆向粒子的最小圆轨道半径作为黑洞俘截粒子的有效截面重要标志,从而大大简化了问题的计算。对K.S.Thorne的结果,本文给出了几种修正.     

基于遗传算法的极短弧定轨

空间目标的巡天观测获取了海量的极短弧观测数据,而经典初轨计算方法对于极短弧几乎不能获得合理的结果.将初轨计算问题转换为两个三变量的分层优化问题,采用遗传算法,针对具体问题选择了优化变量以及相应的遗传操作,建立了一种极短弧初轨计算方法.基于实测资料的数值实验表明,方法可为后续工作提供有效的初值.     

磨镜机上镜面支承的一种近似计算

本文研究了中等大小镜面的九点支承系统,推导了公式,分别计算出对平面,凹、凸镜面,有中孔及无中孔时支承点的位置. 编制了计算程序,可用于PC—1500袖珍计算器.     

倾角函数展开及其在分析法轨道预报中的应用

田谐项摄动是分析法轨道预报中的重要部分,其中包含大量倾角函数及其偏导数的计算.由于具有精度更高、速度更快的优点,倾角函数一般通过递推方法计算.以文献中提出的改进Gooding方法为基础,将其给出的程序稍加改进,在计算2–50阶倾角函数时缩短了约24%的计算时间.考虑到分析法预报过程中轨道平倾角变化很小,以泰勒展开式计算倾角函数,可极大提高计算速度,较大程度地减小分析法预报耗时,且引力场阶次越高,减小幅度越大,取50阶时预报耗时缩短了48%.另一方面,以2阶展开式计算倾角函数时,与改进Gooding法相比,分析法预报星历偏差很小.对于500 km高度的低轨卫星,分别以改进Gooding法和2阶泰勒展开式计算倾角函数,预报3天,当地球引力场阶次不高于50时,二者预报星历偏差RMS (Root Mean Square)低于1 mm,且随着轨道高度的增加,预报星历偏差RMS逐渐减小.     

倾角函数递推公式的误差传递分析

在倾角函数递推关系这个题目上,一些研究者做过工作,如Challe,A(1969),R.H.Googing(1971),J.A.Campbell(1972),G.E.O.Giacaglia(1976)和童付(1979)。最好的结果是由童付获得的,递推关系仅由二个或三个相邻的函数来表达。然而,当倾角较小时,童付的公式存在严重的计算误差传递。本文导出了一个分析公式来估计此种计算误差传递,并在分析计算误差传递机制的基础上,指出可避免计算误差传递的途径,导出的分析公式见文中第(9)式。     

利用Jacobi多项式计算倾角函数

研究了一种利用Jacobi多项式计算倾角函数的方法,该方法表达式非常简单,不存在k和l是否同奇偶和计算非整数阶乘的问题,也不存在k＜0和k≥0的转换问题,递推公式可使用标准的Jacobi多项式的递推公式.而且,计算精度和适应阶数可与Gooding方法相当,计算时间比Gooding方法省9％.     

Hansen系数递推的效率

探讨了偏心率函数的递推方法的效率, 给出了一种成批递推方法, 其计算效率明显优于直接计算方法, 而且递推是正向的. 采用该成批递推方法递推时, 偏心率函数的量级从小到大变化, 可以保证递推的精度.     

EGUP与黑洞熵的修正值

目前,人们对黑洞Bekenstein-Hawking熵的量子修正值产生了极大的兴趣,尤其是黑洞熵对数修正项的系数.在广义不确定关系(GUP)的基础上,通过引入了推广的广义不确定关系(EGUP),运用面积定理计算了3类时空的黑洞熵的修正值,得到的黑洞熵的修正项的系数是正的.这种计算方法不仅对单视界时空适用,而且对有内视界的黑洞时空依然成立,并且在EGUP基础上计算出黑洞熵的修正值.相比GUP基础上得到的黑洞熵,EGUP可以应用于大尺度时空下,所以应用范围更广.此计算方法简洁明了,物理意义明确,可为黑洞熵对数修正值系数的确定提供参考.     

利用TALYS计算太阳耀斑伽玛射线

太阳耀斑伽玛射线能谱是加速粒子与太阳大气介质原子碰撞的结果,它是研究太阳耀斑中加速粒子和高能电子最为直接的手段.通过分析伽玛射线能谱,可以获得耀斑过程中加速粒子的成分、能谱、角分布及太阳大气元素丰度等重要信息.TALYS程序是一套模拟核反应的软件,对核反应过程中的所有信息均能完整地描述.利用TALYS计算得到了完整的太阳耀斑伽玛射线的核反应截面数据,开发了一套新的耀斑伽玛射线谱计算程序.详细介绍了耀斑伽玛射线计算的理论模型,并简单探讨了耀斑伽玛射线的特性,为未来的耀斑伽玛射线能谱分析奠定了理论基础.     

初轨计算精度的bootstrap估计

从非参数统计角度,在仅有观测资料而没有其它信息的情况下,分别给出了初轨计算的精度和置信区间的估计方法.该方法基于bootstrap方法,仅依赖于观测资料,不需要精密定轨结果作为参考,也无需假设观测资料误差呈正态分布,并且对各种初轨计算方法均适用.数值计算结果表明该方法应用简便,可作为初轨的评估和后续应用的重要参考.     

由两块拼接反射镜组成的光学系统的衍射能量分布

LAMOST的光学系统由反射施密特改正板MA 和球面主镜MB 组成 ,MA 和MB 均为对角径 1 .1m的六角形子镜的拼镜面 .本文计算了星在子午圈上时望远镜观测不同的天区零视场和有视场时的光学系统衍射能量分布 ,还计算了MA 不变MB 改为整块的 6m镜以及MA 和MB 均为对角径 1 .5m的六角形子镜的拼镜面的衍射能量分布 .对这些计算结果进行分析 ,论证了六角形子镜尺寸选择的合理性