首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 15 毫秒
1.
We construct a simple, regularized estimator for the dark energy equation of state by using the recently introduced linear response approximation. We show that even a simple regularization substantially improves the performance of the free-form fitting approach. The use of the linear response approximation allows an analytical construction of the maximum likelihood estimator, in a convenient and easy to use matrix form. We show that, in principle, such regularized free-form fitting can give us an unbiased estimate of the functional form of the equation of state of dark energy. We show the efficacy of this approach on simulated SuperNova Acceleration Probe class data, but it is easy to generalize this method to include other cosmological tests. We provide a possible explanation for the sweet spots seen in other reconstruction methods.  相似文献   

2.
We present a new general procedure for determining a given set of quantities. To this end, we define a certain statistic, which we call 'modified  χ2' (χ2M)  , because of its similarity to the standard  χ2  . The terms of this  χ2M  are made up of the fluctuations of an unbiased estimator of some statistical quantities and certain weights. Only the diagonal terms of the covariance matrix appear explicitly in our statistic, while the full covariance matrix (and not its inverse) is included implicitly in the calculation of the weights. Choosing these weights, we may obtain, through minimizing  χ2M  , the estimator that provides the minimum rms, either for those quantities or for the parameters on which these quantities depend. In this paper, we describe our method in the context of cosmic microwave background experiments, in order to obtain either the statistical properties of the maps or the cosmological parameters. The test here is constructed out of some estimator of the two-point correlation function at different angles. For the problem of one-parameter estimation, we show that our method has the same power as the maximum-likelihood method. We have also applied this method to Monte Carlo simulations of the COBE -DMR data, as well as to the actual 4-yr data, obtaining consistent results with previous analyses. We also provide a very good analytical approximation to the distribution function of our statistic, which could also be useful in other contexts.  相似文献   

3.
We present a new unbiased minimal variance (UMV) estimator for the purpose of reconstructing the large-scale structure of the Universe from noisy, sparse and incomplete data. Similar to the Wiener filter (WF), the UMV estimator is derived by requiring the linear minimal variance solution given the data and an assumed a priori model specifying the underlying field covariance matrix. However, unlike the WF, the minimization is carried out with the added constraint of an unbiased reconstructed mean field. The new estimator does not necessitate a noise model to estimate the underlying field; however, such a model is required for evaluating the errors at each point in space. The general application of the UMV estimator is to predict the values of the reconstructed field in unsampled regions of space (e.g. interpolation in the unobserved Zone of Avoidance), and to dynamically transform from one measured field to another (e.g. inversion of radial peculiar velocities to over-densities). Here, we provide two very simple applications of the method. The first is to recover a 1D signal from noisy, convolved data with gaps, for example CMB time-ordered data. The second application is a reconstruction of the density and 3D peculiar velocity fields from mock SEcat galaxy peculiar velocity catalogues.  相似文献   

4.
In this paper we analyse the relations between a previously described oblate Jaffe model for an ellipsoidal galaxy and the observed quantities for NGC 2974, and obtain the length and velocity scales for a relevant elliptical galaxy model. We then derive the finite total mass of the model from these scales, and finally find a good fit of an isotropic oblate Jaffe model by using the Gauss-Hermite fit parameters and the observed ellipticity of the galaxy NGC 2974. The model is also used to predict the total luminous mass of NGC 2974, assuming that the influence of dark matter in this galaxy on the image, ellipticity and Gauss-Hermite fit parameters of this galaxy is negligible within the central region, of radius 0.5R e.  相似文献   

5.
An unbiased method for improving the resolution of astronomical images is presented. The strategy at the core of this method is to establish a linear transformation between the recorded image and an improved image at some desirable resolution. In order to establish this transformation only the actual point spread function and a desired point spread function need be known. No image actually recorded is used in establishing the linear transformation between the recorded and improved image.
This method has a number of advantages over other methods currently in use. It is not iterative, which means it is not necessary to impose any criteria, objective or otherwise, to stop the iterations. The method does not require an artificial separation of the image into 'smooth' and 'point-like' components, and thus is unbiased with respect to the character of structures present in the image. The method produces a linear transformation between the recorded image and the deconvolved image, and therefore the propagation of pixel-by-pixel flux error estimates into the deconvolved image is trivial. It is explicitly constrained to preserve photometry and should be robust against random errors.  相似文献   

6.
Photoclinometry is the most common method used to obtain high-resolution topographic maps of planetary terrain. We derive the likelihood function of photoclinometric surface slope from (1) the probability distribution of the measured photon count of natural sunlight through a Charge-Coupled Device (CCD) including uncertainty due to camera shot noise, camera read noise, small-scale albedo fluctuation and atmospheric haze, and (2) a photometric model relating photocount to surface orientation. We then use classical estimation theory to determine the theoretically exact biases and errors inherent in photoclinometric surface slope and show when they may be approximated by asymptotic expressions for sufficiently high sample size. We show how small-scale albedo variability often dominates biases and errors, which may become an order of magnitude larger than surface slopes when surface reflectance has a weak dependence on surface tilt. We provide bounds on the minimum possible error of any unbiased photoclinometric surface slope estimate, and compute the sample sizes necessary to constrain errors within desired design thresholds.  相似文献   

7.
We consider the Sitnikov problem; from the equations of motion we derive the approximate Hamiltonian flow. Then, we introduce suitable action–angle variables in order to construct a high order normal form of the Hamiltonian. We introduce Birkhoff Cartesian coordinates near the elliptic orbit and we analyze the behavior of the remainder of the normal form. Finally, we derive a kind of local stability estimate in the vicinity of the periodic orbit for exponentially long times using the normal form up to 40th order in Cartesian coordinates.  相似文献   

8.
We estimate the distribution of intrinsic shapes of APM galaxy clusters from the distribution of their apparent shapes. We measure the projected cluster ellipticities using two alternative methods. The first method is based on moments of the discrete galaxy distribution while the second is based on moments of the smoothed galaxy distribution. We study the performance of both methods using Monte Carlo cluster simulations covering the range of APM cluster distances and including a random distribution of background galaxies. We find that the first method suffers from severe systematic biases, whereas the second is more reliable. After excluding clusters dominated by substructure and quantifying the systematic biases in our estimated shape parameters, we recover a corrected distribution of projected ellipticities. We use the non-parametric kernel method to estimate the smooth apparent ellipticity distribution, and numerically invert a set of integral equations to recover the corresponding distribution of intrinsic ellipticities under the assumption that the clusters are either oblate or prolate spheroids. The prolate spheroidal model fits the APM cluster data best.  相似文献   

9.
We study the limits of accuracy for weak lensing maps of dark matter using diffuse 21-cm radiation from the pre-reionization epoch using simulations. We improve on previous 'optimal' quadratic lensing estimators by using shear and convergence instead of deflection angles. This is a generalization of the deflection estimator, and is more optimal for non-Gaussian sources. The cross-power spectrum of shear and convergence is an unbiased estimator of lensing power spectrum which does not require knowledge of the source four-point function. We find that non-Gaussianity provides a limit to the accuracy of weak lensing reconstruction, even if instrumental noise is reduced to zero. The best reconstruction result is equivalent to Gaussian sources with effective independent cell of side length  2.0  h −1 Mpc  . Using a source full map from z = 10 to 20, this limiting sensitivity allows mapping of dark matter at a signal-to-noise ratio greater than 1 out to l ≲ 6000, which is better than any other proposed technique for large-area weak lensing mapping.  相似文献   

10.
We show with analytic models that the assumption of uncorrelated intrinsic ellipticities of target sources that is usually made in searches for weak gravitational lensing arising from large-scale mass inhomogeneities ('field lensing') is unwarranted. If the orientation of the galaxy image is determined either by the angular momentum or by the shape of the halo in which it forms, then the image should be aligned preferentially with the component of the tidal gravitational field perpendicular to the line of sight. Long-range correlations in the tidal field will thus lead to long-range ellipticity–ellipticity correlations that mimic the shear correlations arising from weak gravitational lensing. We calculate the ellipticity–ellipticity correlation expected if halo shapes determine the observed galaxy shape, and we discuss uncertainties (which are still considerable) in the predicted amplitude of this correlation. The ellipticity–ellipticity correlation induced by angular momenta should be smaller. We consider several methods for discriminating between the weak-lensing (extrinsic) and intrinsic correlations, including the use of redshift information. An ellipticity–tidal-field correlation also implies the existence of an alignment of images of galaxies near clusters. Although the intrinsic alignment may complicate the interpretation of field-lensing results, it is inherently interesting as it may shed light on galaxy formation as well as on structure formation.  相似文献   

11.
We develop a new method to estimate the redshift of galaxy clusters through resolved images of the Sunyaev–Zel'dovich effect (SZE). Our method is based on morphological observables which can be measured by actual and future SZE experiments. We test the method with a set of high-resolution hydrodynamical simulations of galaxy clusters at different redshifts. Our method combines the observables in a principal component analysis. After calibrating the method with an independent redshift estimation for some of the clusters, we show – using a Bayesian approach – how the method can give an estimate of the redshift of the galaxy clusters. Although the error bars given by the morphological redshift estimation are large, it should be useful for future SZE surveys where thousands of clusters are expected to be detected; a first preselection of the high-redshift candidates could be done using our proposed morphological redshift estimator. Although not considered in this work, our method should also be useful to give an estimate of the redshift of clusters in X-ray and optical surveys.  相似文献   

12.
We revisit the issue of the recent dynamical evolution of clusters of galaxies using a sample of Abell, Corwin & Olowin (ACO) clusters with   z < 0.14  , which has been selected such that it does not contain clusters with multiple velocity components nor strongly merging or interacting clusters, as revealed in X-rays. We use as proxies of the cluster dynamical state the projected cluster ellipticity, velocity dispersion and X-ray luminosity. We find indications for a recent dynamical evolution of this cluster population, which however strongly depends on the cluster richness. Poor clusters appear to be undergoing their primary phase of virialization, with their ellipticity increasing with redshift with a rate  dε/d z ≃ 2.5 ± 0.4  , while the richest clusters show an ellipticity evolution in the opposite direction (with  dε/d z ≃−1.2 ± 0.1  ), which could be due to secondary infall. When taking into account sampling effects due to the magnitude-limited nature of the ACO cluster catalogue we find no significant evolution of the cluster X-ray luminosity, while the velocity dispersion increases with decreasing redshift, independent of the cluster richness, at a rate  dσ v /d z ≃−1700 ± 400 km s−1  .  相似文献   

13.
Radio interferometers are used to construct high resolution images of the sky at radio frequencies and are the key instruments for accessing the statistical properties of the evolution of neutral hydrogen over cosmic time. Here we use simulated observations of the model sky to assess the efficacy of different estimators of the large-scale structure and power spectrum of the sky brightness distribution. We find that while the large-scale distribution can be reasonably estimated using the reconstructed image from interferometric data, estimates of the power spectrum of the intensity fluctuations calculated from the image are generally biased. This bias is found to be more pronounced for diffuse emission. The visibility based power spectrum estimator, however, gives an unbiased estimate of the true power spectrum. This work demonstrates that for an observation with diffuse emission the reconstructed image can be used to estimate the large-scale distribution of the intensity, while to estimate the power spectrum, visibility based methods should be preferred.With the upcoming experiments aimed at measuring the evolution of the power spectrum of the neutral hydrogen distribution, this is a very important result.  相似文献   

14.
We discuss Hill stability in the general three-body problem. The Hill curves in the general problem are the same as in the planar problem. We show that the bifurcation points correspond to the five equilibrium solutions, and derive the criterion for stability in the general case. Application of this criterion to 19 natural satellites of the Solar system leads to the result that, apart from Neptune 1, all the other 18 satellites are unstable in the sense of Hill. The dominant factor in producing this result is the finite eccentricity of the planetary orbits around the Sun.  相似文献   

15.
In this work we study the performance of linear multifilters for the estimation of the amplitudes of the thermal and kinematic Sunyaev–Zel'dovich (SZ) effects. We show that when both effects are present, estimation of these effects with standard matched multifilters is intrinsically biased. This bias is due to the fact that both signals have basically the same spatial profile. We find a new family of multifilters related to the matched multifilters that cancel this systematic bias, hence we call them unbiased matched multifilters. We test the unbiased matched multifilters and compare them with the standard matched multifilters using simulations that reproduce the future Planck mission observations. We find that in the case of the standard matched multifilters the systematic bias in the estimation of the kinematic Sunyaev–Zel'dovich effect can be very large, even greater than the statistical error bars. Unbiased matched multifilters cancel this kind of bias effectively. In concordance with other works in the literature, our results indicate that the sensitivity and resolution of Planck will not be enough to give reliable estimations of the kinematic Sunyaev–Zel'dovich effects of individual clusters. However, as the estimation with the unbiased matched multifilters is not intrinsically biased, it can be possible to use them to study statistically any peculiar cosmological bulk flows via the kinematic SZ effect.  相似文献   

16.
  相似文献   

17.
Rodin  A. V. 《Solar System Research》2003,37(2):101-111
The previously proposed (Rodin, 2002) method for calculating the microphysical properties of spatially inhomogeneous rarefied aerosol media with mixing using the lowest-order moments of the size distribution is generalized to particle coagulation. We show that when the problem is formulated in terms of moments, all of the solutions admitted by the stochastic coagulation equation lie within a narrow range whose boundaries can be determined by means of quadratic programming. We discuss the choice of an optimal solution within this range and compare the moment method with the results of our computations by the classical finite-difference method using a model of photochemical aerosols in Titan's atmosphere as an example. The moment method allows the efficiency of microphysical computations to be significantly increased by using precomputed low-dimension interpolation tables. It can be used to construct self-consistent models for the globular circulation of planetary atmospheres.  相似文献   

18.
The theory of broad-absorption-line (BAL) QSOs is worked out on the basis of the radiation hydrodynamics equation solution for the two-phase media of active galactic nuclei. We suppose that the BAL QSO physics depends on some 'hidden' AGN parameters, such as the mass and size of the compact stellar system. We therefore approach this problem in the more general framework of the 'interacting subsystems theory', which includes these parameters. We compare the results of the numerical model calculations with the observed spectra and show that the BAL QSOs (and the radio-quiet quasars as well) contain massive compact stellar kernels in their central regions. We show that the line-locking effect is determined by the radiation pressure, and is also favoured by the drag force of the hot gas acting on the line-absorbing clouds.
We derive some general conclusions about the physics of AGN. In particular, we show that the radio-quiet versus radio-loud dichotomy can be explained by using two types of hot gas outstreams in quasars.  相似文献   

19.
We present a new method of studying quadruple lenses in elliptical power-law potentials parametrized by ψ ( x , y )∝( x 2+ y 2 q 2) β /2 β (0 β <2). For this potential, the moments of the four image positions weighted by signed magnifications (magnification times parity) have very simple properties. In particular, we find that the zeroth moment – the sum of four signed magnifications satisfies ≃2/(2− β ); the relation is exact for β =0 (point-lens) and β =1 (isothermal potential), independent of the axial ratio. Similar relations can be derived when a shear is present along the major or minor axes. These relations, however, do not hold well for the closely related elliptical density distributions. For a singular isothermal elliptical density distribution without shear, the sum of signed magnifications for quadruple lenses is ≈2.8, again nearly independent of the ellipticity. For the same distribution with shear, the total signed magnification is around 2–3 for most cases, but can be significantly different for some combinations of the axial ratio and shear where six or eight images can appear.  相似文献   

20.
It is an objective fact that there exists error in the satellite dynamic model and it will be transferred to satellite orbit determination algorithm, forming a part of the connotative model error. Mixed with the systematic error and random error of the measurements, they form the unitive model error and badly restrict the precision of the orbit determination. We deduce in detail the equations of orbit improvement for a system with dynamic model error, construct the parametric model for the explicit part of the model and nonparametric model for the error that can not be explicitly described. We also construct the partially linear orbit determination model, estimate and fit the model error using a two-stage estimation and a kernel function estimation, and finally make the corresponding compensation in the orbit determination. Beginning from the data depth theory, a data depth weight kernel estimator for model error is proposed for the sake of promoting the steadiness of model error estimation. Simulation experiments of SBSS are performed. The results show clearly that the model error is one of the most important effects that will influence the precision of the orbit determination. The kernel function method can effectively estimate the model error, with the window width as a major restrict parameter. A data depth-weight-kernel estimation, however, can improve largely the robustness of the kernel function and therefore improve the precision of orbit determination.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号