共查询到20条相似文献，搜索用时 62 毫秒
1.
2.
Stéphane Colombi rew Jaffe Dmitri Novikov Christophe Pichon 《Monthly notices of the Royal Astronomical Society》2009,393(2):511-526
A method to rapidly estimate the Fourier power spectrum of a point distribution is presented. This method relies on a Taylor expansion of the trigonometric functions. It yields the Fourier modes from a number of fast Fourier transforms (FFTs), which is controlled by the order N of the expansion and by the dimension D of the system. In three dimensions, for the practical value N = 3 , the number of FFTs required is 20.
We apply the method to the measurement of the power spectrum of a periodic point distribution that is a local Poisson realization of an underlying stationary field. We derive an explicit analytic expression for the spectrum, which allows us to quantify – and correct for – the biases induced by discreteness and by the truncation of the Taylor expansion, and to bound the unknown effects of aliasing of the power spectrum. We show that these aliasing effects decrease rapidly with the order N . For N = 3 , they are expected to be, respectively, smaller than ∼10−4 and 0.02 at half the Nyquist frequency and at the Nyquist frequency of the grid used to perform the FFTs. The only remaining significant source of errors is reduced to the unavoidable cosmic/sample variance due to the finite size of the sample.
The analytical calculations are successfully checked against a cosmological N -body experiment. We also consider the initial conditions of this simulation, which correspond to a perturbed grid. This allows us to test a case where the local Poisson assumption is incorrect. Even in that extreme situation, the third-order Fourier–Taylor estimator behaves well, with aliasing effects restrained to at most the per cent level at half the Nyquist frequency.
We also show how to reach arbitrarily large dynamic range in Fourier space (i.e. high wavenumber), while keeping statistical errors in control, by appropriately 'folding' the particle distribution. 相似文献
We apply the method to the measurement of the power spectrum of a periodic point distribution that is a local Poisson realization of an underlying stationary field. We derive an explicit analytic expression for the spectrum, which allows us to quantify – and correct for – the biases induced by discreteness and by the truncation of the Taylor expansion, and to bound the unknown effects of aliasing of the power spectrum. We show that these aliasing effects decrease rapidly with the order N . For N = 3 , they are expected to be, respectively, smaller than ∼10
The analytical calculations are successfully checked against a cosmological N -body experiment. We also consider the initial conditions of this simulation, which correspond to a perturbed grid. This allows us to test a case where the local Poisson assumption is incorrect. Even in that extreme situation, the third-order Fourier–Taylor estimator behaves well, with aliasing effects restrained to at most the per cent level at half the Nyquist frequency.
We also show how to reach arbitrarily large dynamic range in Fourier space (i.e. high wavenumber), while keeping statistical errors in control, by appropriately 'folding' the particle distribution. 相似文献
3.
针对包含饱和样本数据的频数幂律分布拟合,提出一个新的幂律分布指数的极大似然估计方法的修正公式.对比研究显示,修正公式适用于剔除异常饱和值的幂律频数拟合.如果不剔除饱和值,幂律指数的估计只能使用修正前的公式,其误差随幂律指数变化,指数较小时误差较大.由此建议,对于包含饱和样本的频数分布拟合,首先剔除异常的饱和值,然后对剩余不含饱和值的子集使用修正公式进行参数估计. 相似文献
4.
Chris Koen 《Monthly notices of the Royal Astronomical Society》2005,361(3):887-896
The theory of low-order linear stochastic differential equations is reviewed. Solutions to these equations give the continuous time analogues of discrete time autoregressive time-series. Explicit forms for the power spectra and covariance functions of first- and second-order forms are given. A conceptually simple method is described for fitting continuous time autoregressive models to data. Formulae giving the standard errors of the parameter estimates are derived. Simulated data are used to verify the performance of the methods. Irregularly spaced observations of the two hydrogen-deficient stars FQ Aqr and NO Ser are analysed. In the case of FQ Aqr the best-fitting model is of second order, and describes a quasi-periodicity of about 20 d with an e-folding time of 3.7 d. The NO Ser data are best fitted by a first-order model with an e-folding time of 7.2 d. 相似文献
5.
6.
多任务学习方法在机器学习、计算机视觉、人工智能领域已得到广泛关注,利用任务间的相关性,将多个任务同时学习的效果优于每个任务单独学习的情况.采用多任务Lasso回归法(Multi-task Lasso Regression)用于恒星光谱物理参量的估计,不仅可以获取不同物理参量间的共同的特征信息,而且也可以很好地保留不同物理参量的特有的补充信息.使用恒星大气模拟模型合成光谱库ELODIE中的光谱数据和美国大型巡天项目Sloan发布的SDSS实测光谱数据进行实验,模型估算精度优于相关文献中的方法,特别是对重力加速度(lg g)和化学丰度([Fe/H])的估计.实验中通过改变光谱的分辨率,施加不同信噪比(SNR)的噪声,来说明模型的稳定性强.结果表明,模型精度受光谱分辨率和噪声的影响,但噪声对其影响更大,可见,多任务Lasso回归法不仅操作简便,稳定性强,而且也提高了模型的整体预测精度. 相似文献
7.
8.
9.
Chris Koen 《Monthly notices of the Royal Astronomical Society》2001,325(4):1335-1340
The observation of flux sources near the limit of detection requires a careful evaluation of possible biases in magnitude determination. Both the traditional logarithmic magnitudes and the recently proposed inverse hyperbolic sine (asinh) magnitudes are considered. Formulae are derived for three different biasing mechanisms: the statistical spread of the observed flux values arising from e.g. measurement error; the dependence of these errors on the true flux; and the dependence of the observing probability on the true flux. As an example of the results, it is noted that biases at large signal-to-noise ratios R , at which the two types of magnitude are similar, are of the order of −( p +1)/ R 2 , where the exponent p parametrizes a power-law dependence of the probability of observation on the true flux. 相似文献
10.
Chris Koen 《Monthly notices of the Royal Astronomical Society》1999,309(3):769-802
The estimation of the frequency, amplitude and phase of a sinusoid from observations contaminated by correlated noise is considered. It is assumed that the observations are regularly spaced, but may suffer missing values or long time stretches with no data. The typical astronomical source of such data is high-speed photoelectric photometry of pulsating stars. The study of the observational noise properties of nearly 200 real data sets is reported: noise can almost always be characterized as a random walk with superposed white noise. A scheme for obtaining weighted non-linear least-squares estimates of the parameters of interest, as well as standard errors of these estimates, is described. Simulation results are presented for both complete and incomplete data. It is shown that, in finite data sets, results are sensitive to the initial phase of the sinusoid. 相似文献
11.
The auroras on Jupiter and Saturn can be studied with a high sensitivity and resolution by the Hubble Space Telescope ( HST ) ultraviolet (UV) and far-ultraviolet Space Telescope Imaging Spectrograph (STIS) and Advanced Camera for Surveys (ACS) instruments. We present results of automatic detection and segmentation of Jupiter's auroral emissions as observed by the HST ACS instrument with the VOronoi Image SEgmentation (VOISE). VOISE is a dynamic algorithm for partitioning the underlying pixel grid of an image into regions according to a prescribed homogeneity criterion. The algorithm consists of an iterative procedure that dynamically constructs a tessellation of the image plane based on a Voronoi diagram, until the intensity of the underlying image within each region is classified as homogeneous. The computed tessellations allow the extraction of quantitative information about the auroral features, such as mean intensity, latitudinal and longitudinal extents and length-scales. These outputs thus represent a more automated and objective method of characterizing auroral emissions than manual inspection. 相似文献
12.
Coryn A. L. Bailer-Jones Mike Irwin & Ted von Hippel 《Monthly notices of the Royal Astronomical Society》1998,298(2):361-377
We investigate the application of neural networks to the automation of MK spectral classification. The data set for this project consists of a set of over 5000 optical (3800–5200 Å) spectra obtained from objective prism plates from the Michigan Spectral Survey. These spectra, along with their two-dimensional MK classifications listed in the Michigan Henry Draper Catalogue, were used to develop supervised neural network classifiers. We show that neural networks can give accurate spectral type classifications (σ68 = 0.82 subtypes, σrms = 1.09 subtypes) across the full range of spectral types present in the data set (B2–M7). We show also that the networks yield correct luminosity classes for over 95 per cent of both dwarfs and giants with a high degree of confidence. Stellar spectra generally contain a large amount of redundant information. We investigate the application of principal components analysis (PCA) to the optimal compression of spectra. We show that PCA can compress the spectra by a factor of over 30 while retaining essentially all of the useful information in the data set. Furthermore, it is shown that this compression optimally removes noise and can be used to identify unusual spectra. This paper is a continuation of the work carried out by von Hippel et al. (Paper I). 相似文献
13.
N. Solomos 《Astronomische Nachrichten》2008,329(3):252-254
This contribution aims to introduce the idea that a well‐evolved HTN of the far future, with the anticipated addition of very large apertures, could also be made to incorporate the ability to carry out photonic astronomy observations, particularly Optical VLBI in a revived Hanbury‐Brown and Twiss Intensity Interferometry (HBTII) configuration. Such an HTN could exploit its inherent rapid reconfigurational ability to become a multi‐aperture distributed photon‐counting network able to study higher‐order spatiotemporal photon correlations and provide a unique tool for direct diagnostics of astrophysical emission processes. We very briefly review various considerations associated with the switching of the HTN to a special mode in which single‐photon detection events are continuously captured for a posteriori intercorrelation. In this context, photon arrival times should be determined to the highest time resolution possible and extremely demanding absolute time keeping and absolute time distribution schemes should be devised and implemented in the HTN nodes involved. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim) 相似文献
14.
望远镜调度是望远镜运行中的关键组成部分, 用于辅助科研人员进行合理的观测计划安排, 提高望远镜的运行效率, 获取高质量的观测数据. 然而, 由于不同观测项目的科学需求不同, 望远镜的调度过程十分复杂. 针对短周期多目标的观测项目, 考虑望远镜换源时转动时长、观测高度角等因素进行建模, 采用贪心算法对中国科学院新疆天文台南山26m望远镜脉冲星到达时间观测列表进行调度. 通过模拟表明, 使用算法完成的观测列表可以有效地减少观测过程中的平均转动时长, 提升观测数据的质量, 提高望远镜时间利用率, 减少科研人员对观测列表编排的负担. 相似文献
15.
T. D. Kitching A. Amara F. B. Abdalla B. Joachimi A. Refregier 《Monthly notices of the Royal Astronomical Society》2009,399(4):2107-2128
In the absence of any compelling physical model, cosmological systematics are often misrepresented as statistical effects and the approach of marginalizing over extra nuisance systematic parameters is used to gauge the effect of the systematic. In this article, we argue that such an approach is risky at best since the key choice of function can have a large effect on the resultant cosmological errors.
As an alternative we present a functional form-filling technique in which an unknown, residual, systematic is treated as such. Since the underlying function is unknown, we evaluate the effect of every functional form allowed by the information available (either a hard boundary or some data). Using a simple toy model, we introduce the formalism of functional form filling. We show that parameter errors can be dramatically affected by the choice of function in the case of marginalizing over a systematic, but that in contrast the functional form-filling approach is independent of the choice of basis set.
We then apply the technique to cosmic shear shape measurement systematics and show that a shear calibration bias of | m ( z )| ≲ 10−3 (1 + z )0.7 is required for a future all-sky photometric survey to yield unbiased cosmological parameter constraints to per cent accuracy.
A module associated with the work in this paper is available through the open source icosmo code available at http://www.icosmo.org . 相似文献
As an alternative we present a functional form-filling technique in which an unknown, residual, systematic is treated as such. Since the underlying function is unknown, we evaluate the effect of every functional form allowed by the information available (either a hard boundary or some data). Using a simple toy model, we introduce the formalism of functional form filling. We show that parameter errors can be dramatically affected by the choice of function in the case of marginalizing over a systematic, but that in contrast the functional form-filling approach is independent of the choice of basis set.
We then apply the technique to cosmic shear shape measurement systematics and show that a shear calibration bias of | m ( z )| ≲ 10
A module associated with the work in this paper is available through the open source icosmo code available at http://www.icosmo.org . 相似文献
16.
Chris Koen 《Monthly notices of the Royal Astronomical Society》2000,316(3):613-620
One of the tools used to identify the pulsation modes of stars is a comparison of the amplitudes and phases as observed photometrically at different wavelengths. Proper application of the method requires that the errors on the measured quantities, and the correlations between them, be known (or at least estimated). It is assumed that contemporaneous measurements of the light intensity of a pulsating star are obtained in several wavebands. It is also assumed that the measurements are regularly spaced in time, although there may be missing observations. The amplitude and phase of the pulsation are estimated separately for each of the wavebands, and amplitude ratios and phase differences are calculated. A general scheme for estimating the covariance matrix of the amplitude ratios and phase differences is described. The first step is to fit a time series to the residuals after pre-whitening the observations by the best-fitting sinusoid. The residuals are then cross-correlated to study the interdependence between the errors in the different wavebands. Once the multivariate time-series structure can be modelled, the covariance matrix can be found by bootstrapping. An illustrative application is described in detail. 相似文献
17.
A. Pagliaro V. Antonuccio-Delogu U. Becciani M. Gambera 《Monthly notices of the Royal Astronomical Society》1999,310(3):835-841
We present and discuss a method to identify substructures in combined angular-redshift samples of galaxies within clusters. The method relies on the use of the discrete wavelet transform (hereafter DWT) and has already been applied to the analysis of the Coma cluster. The main new ingredient of our method with respect to previous studies lies in the fact that we make use of a 3D data set rather than a 2D one. We test the method on mock cluster catalogues with spatially localized substructures and on a N -body simulation. Our main conclusion is that our method is able to identify the existing substructures provided that: (a) the subclumps are detached in part or all of the phase space, (b) one has a statistically significant number of redshifts, increasing as the distance decreases due to redshift distortions; (c) one knows a priori the scale on which substructures are to be expected. We have found that to allow an accurate recovery we must have both a significant number of galaxies (≈200 for clusters at z ≥0.4 or about 800 at z ≤0.4) and a limiting magnitude for completeness m B =16.
The only true limitation to our method seems to be the necessity of knowing a priori the scale on which the substructure is to be found. This is an intrinsic drawback of the method and no improvement in numerical codes based on this technique could make up for it. 相似文献
The only true limitation to our method seems to be the necessity of knowing a priori the scale on which the substructure is to be found. This is an intrinsic drawback of the method and no improvement in numerical codes based on this technique could make up for it. 相似文献
18.
Chris Koen 《Monthly notices of the Royal Astronomical Society》2006,365(2):489-495
It is assumed that O − C ('observed minus calculated') values of periodic variable stars are determined by three processes, namely measurement errors, random cycle-to-cycle jitter in the period, and possibly long-term changes in the mean period. By modelling the latter as a random walk, the covariances of all O − C values can be calculated. The covariances can then be used to estimate unknown model parameters, and to choose between alternative models. Pseudo-residuals which could be used in model fit assessment are also defined. The theory is illustrated by four applications to spotted stars in eclipsing binaries. 相似文献
19.
A statistical model is formulated that enables one to analyse jointly the times between maxima and minima in the light curves of monoperiodic pulsating stars. It is shown that the combination of both sets of data into one leads to analyses that are more sensitive. Illustrative applications to the American Association of Variable Star Observers data for a number of long-period variables demonstrate the usefulness of the approach. 相似文献
20.
T. Toutain Y. Elsworth W. J. Chaplin 《Monthly notices of the Royal Astronomical Society》2006,371(4):1731-1736
It is well known that in the power spectrum solar p modes have asymmetric profiles, which depart from a Lorentzian shape. We present a framework to explain the contribution of correlated background noise, from the acoustic source, to this asymmetry. An important prediction is that observed peak asymmetry may differ depending on the way the p-mode observations are made, and on how the data are prepared. Furthermore, if valid, the proposed framework may provide the basis for separating the contribution of the correlated noise from that of the source location and properties. 相似文献