基于网络的GNSS共视时间比对数据传输系统的设计与实现

为了实现远程时间同步中的共视数据自动、连续传输,提出了基于TCP/IP协议来实现C/S模式的数据传输方案。阐述了数据发送端与接收端的传输原理和实现方案。发送端通过串口采集数据并通过网络将其传输到接收端,接收端接收发送端的数据并对其检验并保存。设计了发送端与接收端自动重连接功能,实现了数据传输的连续性。测试表明:采用基于TCP/IP协议和C/S模式的数据传输系统能够实现共视数据自动、连续、可靠的传输。     

嵌入式TCP/IP协议在接收机杜瓦性能参数监测系统中的应用

介绍了射电天文望远镜接收机杜瓦温度与真空度监测系统的设计,具体讨论了监测系统的硬件结构与嵌入式TCP/IP协议在该系统中的实现.针对目前大量使用8位微控制器的嵌入式系统,硬件主控系统芯片采用ATMEL公司生产的8位单片机ATmega16,网络控制芯片采用美国MICROCHIP公司生产的带SPI接口的独立以太网微控制器ENC28J60,设计实现了低成本、高效率、高稳定性的性能参数监测系统.对该系统进行测试分析、数据比对,测试结果表明该系统能够高效与PC机进行通信,PC机所接收到的温度值与标准数据的误差在±0.3 K范围内,接收到的真空值绝对误差在±10 ubar范围内,能够达到系统设计要求.     

相干消色散脉冲星观测系统的研究

云南天文台射电天文研究团组利用从美国伯克利大学CASPER天文信号与电子学研究中心购买的现场可编程门阵列(Field-Programmable Gate Array,FPGA)接收机平台ROACH2,实现了512 MHz输入带宽,512 MHz带宽分为128通道输出(每个通道4 M带宽),8比特采样和双极化输入(1 024 MHz)的基带数据采集终端。海量数据传输方式通过SFP+万兆网口实现,利用编写的脚本文件调用DSPSR程序包实现数据的解码、相干消色散、偏振计算和折叠等处理。数据处理结果以PSRFITS格式存储。构建以ROACH2为基带数据采集终端和DSPSR为数据处理核心的脉冲星观测系统,相比于以VLBI观测终端为基础构建的观测系统,在观测模式、数据处理方法、运算效率和观测数据的通用性等方面具有更好的优越性。     

新疆天文台天文观测数据传输日志系统设计与实现

新疆天文台望远镜产生的观测数据通过数据专线传输到新疆天文台本部数据中心进行长期存储。采用rsync进行数据传输,传输过程中没有详细的日志记录信息,无法对传输情况进行统计分析。通过调研新疆天文台数据传输的需求,设计并开发了新疆天文台观测数据传输日志系统。系统以新疆天文台南山站、奇台站的数据存储服务器和新疆天文台本部长期存储为基础,设计了日志系统的整体架构,实现了数据传输过程的日志记录。系统基于shell多线程技术及Qt框架,开发了后台服务程序及日志系统软件界面,实现了日志分析、统计备份与数据传输的可视化页面。     

四通道可拓展射电数字接收机系统

针对目前射电观测设备对数字接收机高采样率、宽带宽、多通道幅相一致,以及高速率直接采样和时域数据存储的需求,通过调研分析多个射电观测设备的数字处理系统的技术架构和功能,提出了基于ZYNQ SOC和ADS54J60等核心器件实现的4通道数字接收机系统方案。单通道最高可实现1 GSPS采样率,具有灵活可拓展的优势,可以通过添加板卡实现采样通道数的增加,能够满足未来大规模可拓展射电干涉阵列的需求。系统由高速数据采集卡、光通信接收卡、服务器3部分设备组成,基于SerDes的高速串行接口技术目前实现16位量化精度、300 MSPS的直接采样功能,具备60 dB满量程信噪比、40 Gbps的SFP+数据传输带宽、1.5 GB/s的PCIe通信带宽,可对4.5～150 MHz范围内中频模拟信号进行采集,目前该系统已经完成软硬件设计和测试。在测试中对频率为10 MHz、幅度为125 mV的正弦波信号进行采样,得到4通道间幅度差小于1 mV,相位延时小于3.3 ns。系统集成可编程SFP+通信接口和多机同步机制,可适应多元阵列同步采集的需要,可同时对原始射电信号的数据进行存储,为射电研究提供更详细的时域数...     

基于以太网的射电天文制冷接收机的监视系统

制冷接收机是射电望远镜的核心设备,它是否正常工作直接决定望远镜观测的效果。制冷接收机的制冷温度和杜瓦真空度是反映接收机是否正常工作的最重要、最直接的指标之一。因此实现制冷接收机制冷温度和杜瓦真空度的远程实时监控,及时了解接收机是否正常工作对保证射电望远镜正常运行、提高观测效率有重要意义。为乌鲁木齐天文站25m射电望远镜1.3cm制冷接收机研制的一套基于单片机和以太网的数据采集和数据传输的远程监视系统,实现了制冷接收机制冷温度和杜瓦真空度的远程实时监控。该系统采用了美国ATMEL公司生产的8位单片机AVR ATmega16、Microchip Technology公司生产的ENC28J60以及MAXIM公司生产的MAX7219,实现对射电天文制冷接收机制冷温度和杜瓦真空度的数据采集,并利用以太网传输数据实现了远程实时监控接收机的制冷状态。阐述了以太网数据采集及远程监控电路的设计原理及其实现方法。该系统首次在国内大型射电望远镜上实现了对制冷接收机工作状态的远程实时监控,对于保证乌鲁木齐天文站25m射电望远镜1.3cm波段的观测效果有重要作用。     

1米太阳望远镜液压监控系统

1m太阳望远镜的方位驱动控制系统采用液压轴承,液压系统正常工作时望远镜才能实现平稳旋转。可通过油压、油温监测液压系统的工作状态。针对原有液压系统不能提供其工作状态实时监测和自动化程度较低的缺陷,设计了液压监控系统并将其集成到望远镜控制系统中,可实现油压油温数据采集、油压异常报警、油泵自动启停等功能,保证液压轴承的安全运转。介绍了该监控系统的结构和工作原理,阐述了油压、油温数据采集系统的软硬件实现并给出了相应的实验数据和运行结果。理论分析和实验结果表明系统工作正常,数据采集速度快、分辨率高,精度可满足要求,能长期稳定运行。     

H-alpha太阳全日面像数据管理系统的建立

云南天文台的H-alpha太阳色球望远镜每天的数据量约有2G bytes，大量数据带来了使用上的不便。因此，实现数据的查询和共享成为了一个相当重要的问题。为此，试图通过建立H-alpha数据的管理系统，实现对工作流程的改进，完成对数据的检索、在线浏览和远程下载。在WindwsXP操作系统和MysQL数据库管理系统的平台上，利用JAVA语言和JSP技术开发了一个基于浏览器／服务器模式的原型系统，可以完成对数据的压缩存储，实现数据的检索、实时浏览、删除(需要相应权限)和下载等功能。     

基于Mark5B＋GPU脉冲星观测系统

云南天文台40m射电望远镜进行的脉冲星观测数据量巨大,必须实现数据的实时处理,否则将会产生海量的数据积压.为实现这一目标,采用图形处理器架构,对Mark5B数据进行解码、消色散、折叠等处理.实验结果表明,对以1s8MB的实时采样,可以在0.51s内处理完成,从而实现了实时处理的要求.首先介绍这一观测系统各部分的图形处理器实现,然后相对于传统中央处理器构架,对各部分的运算速度进行了详细的对比.针对时间开销最大的消色散部分,分析了单次傅里叶变换的数据量大小对执行效率的影响.从系统最终的输出轮廓和柱状图上可以看到实时处理的结果符合要求.最后对存在的问题和未来的工作进行了讨论.     

2.4m望远镜摄谱仪狭缝监视系统

为了实现对云南天文台2.4m望远镜摄谱仪狭缝的远程监视,我们以CPLD为核心设计了CCD控制器,用PC-104嵌入式计算机采集图像,利用pcAnywhere软件实现远程数据传输,以较低的成本设计制作了一套狭缝监视系统。该系统目前已经通过测试,达到了狭缝监视的基本要求,并为下一步实现自动导星和Robotic观测打下了基础。本文详细介绍了该系统的软硬件设计和测试结果。     

Impact of the Quadrupole Moment of the Sun on the Dynamics of the Earth-Moon System

Range of values of the Sun's mass quadrupole moment of coefficient J₂ arising both from experimental and theoretical determinations enlarge across literature on two orders of magnitude, from around 10^-7 until to 10^-5. The accurate knowledge of the Moon's physical librations, for which the Lunar Laser Ranging data reach an outstanding precision level, prove to be appropriate to reduce the interval of J₂ values by giving an upper bound of J₂. A solar quadrupole moment as high as 1.1 10^-5 given either from the upper bounds of the error bars of the observations, or from the Roche's theory, is not compatible with the knowledge of the lunar librations accurately modeled and observed with the LLR experiment. The suitable values of J₂ have to be smaller than 3.0 10^-6. As a consequence, this upper bound of 3.0 10^-6 is accepted to study the impact of the Sun's quadrupole moment of mass on the dynamics of the Earth-Moon system. Such as effect (with J₂ = 5.5±1.3 × 10^-6) has been already tested in 1983 by Campbell & Moffat using analytical approximate equations, and thus for the orbits of Mercury, Venus, the Earth and Icarus. The approximate equations are no longer sufficient compared with present observational data and exact equations are required. As if to compute the effect on the lunar librations, we have used our BJV relativistic model of solar system integration including the spin-orbit coupled motion of the Moon. The model is solved by numerical integration. The BJV model stems from general relativity by using the DSX formalism for purposes of celestial mechanics when it is about to deal with a system of n extended, weakly self-gravitating, rotating and deformable bodies in mutual interactions. The resulting effects on the orbital elements of the Earth have been computed and plotted over 160 and 1600 years. The impact of the quadrupole moment of the Sun on the Earth's orbital motion is mainly characterized by variations of , , and . As a consequence, the Sun's quadrupole moment of mass could play a sensible role over long time periods of integration of solar system models. This revised version was published online in July 2006 with corrections to the Cover Date.     

On the role of the motion of the photospheric footpoints of the magnetic field lines

By means of a simple relation between the velocity v of the fluid particle and the velocity v_f of the photospheric footpoint of the magnetic field line v_z and B_z being respectively the components of v and the magnetic field B normal to the photospheric surface, it is shown formally that through the phtospheric surface the transport of all the quantities attributed to the magnetic field, such as the magnetic flux, the magnetic energy and the helicity, is independent of v_z, and v_f is the only kinematical quantity on which the transport depends. In addition, in the neighborhood of the neutral line the velocity v_l of the moving curve of constant B_z is found to be equal approximately to the component of v or v_f in the direction of v_l. Since v_l can be measured or extimated, so can the components of v and v_f near the neutral line.     

在后牛顿引力理论(PPN形式)中引力常数变化对地球自转产生的效应

介绍和论述了在后牛顿引力理论(PPN形式)中在优越参考系和非优越参考系中经过参数化后引力常数变化对地球自转产生的效应,其中特别重点介绍了年周期变化的效应。此外也将理论结果同观测结果相对比。     

On the theory of the oblateness of the Sun

In this paper we first emphasize why it is important to know the successive zonal harmonics of the Sun's figure with high accuracy: mainly fundamental astrometry, helioseismology, planetary motions and relativistic effects. Then we briefly comment why the Sun appears oblate, going back to primitive definitions in order to underline some discrepancies in theories and to emphasize again the relevant hypotheses. We propose a new theoretical approach entirely based on an expansion in terms of Legendre's functions, including the differential rotation of the Sun at the surface. This permits linking the two first spherical harmonic coefficients (J ₂ and J ₄) with the geometric parameters that can be measured on the Sun (equatorial and polar radii). We emphasize the difficulties in inferring gravitational oblateness from visual measurements of the geometric oblateness, and more generally a dynamical flattening. Results are given for different observed rotational laws. It is shown that the surface oblateness is surely upper bounded by 11 milliarcsecond. As a consequence of the observed surface and sub-surface differential rotation laws, we deduce a measure of the two first gravitational harmonics, the quadrupole and the octopole moment of the Sun: J ₂=−(6.13±2.52)×10⁻⁷ if all observed data are taken into account, and respectively, J ₂=−(6.84±3.75)×10⁻⁷ if only sunspot data are considered, and J ₂=−(3.49±1.86)×10⁻⁷ in the case of helioseismic data alone. The value deduced from all available data for the octopole is: J ₄=(2.8±2.1)×10⁻¹². These values are compared to some others found in the literature. Supplementary material to this paper is available in electronic form at http://dx.doi.org/10.1023/A:1005238718479     

Third Integral and the Dynamics of the Galaxy in the Neighborhood of the Sun

Observational data on the dynamics of stars in the neighborhood of the sun indicate the existence of a third integral besides the integrals of the angular momentum and energy. The Poincaré integral is proposed as a third integral. The consequences of this assumption are derived and compared with available astrophysical data.     

Accounting for the viscosity of the fluid core in the problem of the physical libration of the moon

A two-component theoretical model of the physical libration of the Moon in longitude is constructed with account taken of the viscosity of the core. In the new version, a hydrodynamic problem of motion of a fluid filling a solid rotating shell is solved. It is found that surfaces of equal angular velocity are spherical, and a velocity field of the fluid core of the Moon is described by elementary functions. A distribution of the internal pressure in the core is found. An angular momentum exchange between the fluid core and solid mantle is described by a third-order differential equation with a right-hand side. The roots of a characteristic equation are studied and the stability of rotation is proved. A libration angle as a function of time is found using the derived solution of the differential equation. Limiting cases of infinitely large and infinitely small viscosity are considered and an effect of lag of a libration phase from a phase of action of an external moment of forces is ascertained. This makes it possible to estimate the viscosity and sizes of the lunar fluid core from data of observations.     

射电日像仪灵敏度测量

灵敏度是射电望远镜的一个重要性能指标,它反映了望远镜监测弱信号的能力。基于明安图射电频谱日像仪(Mingantu Ultrawide Spectral Radioheliograph,MUSER)的调试观测,给出了日像仪灵敏度的测量方法,对天线系统以及整个阵列的灵敏度进行测量分析,得到了日像仪系统整体的灵敏度性能参数。测量同时给出了天线系统的效率以及接收机系统的增益,这将为下一步日像仪展开常规的科学观测提供参考。     

Investigation of the kinematics and structure of the Galaxy in the vicinity of the North galactic pole

In the framework of the programme of studying the meridional section of the Galaxy (MEGA) the absolute proper motions of more than 11000 stars with respect to 3000 galaxies and their stellar magnitudes in the B, V Johnson system are determined in two sky regions near the North Galactic Pole (NGP) by means of Tautenburg Schmidt plates. The limiting and completing apparent stellar magnitudes are B = 20.4 and 18.3 mag, the overall and the investigated sky areas are 16.4 and 14.6 square degrees, respectively. Distances have been determined using the stellar magnitudes, colours, proper motions and reduced proper motions. Stellar kinematics, eccentricities of Galactic orbits, spatial distribution and changes of these characteristics with Z-distance from the Galactic plane are obtained up to 15 kpc. Four subsystems distinguished in the NGP direction, respectively with semithicknesses of 0.25, 0.38, 0.67, 1.48 kpc and density ellipsoid axial ratios of 0.09, 0.20, 0.28, 0.49 show mean velocities in the Galactic rotation direction relative to the LSR of 5.6 ± 0.6, − 11.0 ± 0.6, − 62.5 ± 1.2, − 181.6 ± 4.4 km/s, and ages of 0.1, 0.4, 0.9, 1 of the Galaxy age.     

Rotation of the elastic Earth: the role of the angular-velocity-dependence of the elasticity-caused perturbation

We calculate the so-called convective term, which shows up in the expression for the angular velocity of the elastic Earth, within the Andoyer formalism. The term emerges due to the fact that the elasticity-caused perturbation depends not only on the instantaneous orientation of the Earth but also on its instantaneous angular velocity. We demonstrate that this term makes a considerable contribution into the overall angular velocity. At the same time the convective term turns out to be automatically included into the correction to the nutation series due to the elasticity, if the series is defined by the perturbation of the figure axis (and not of the rotational axis) in accordance with the current IAU resolution. Hence it is not necessary to take the effect of the convective term into consideration in the perturbation of the elastic Earth as far as the nutation is related to the motion of the figure axis.