ترغب بنشر مسار تعليمي؟ اضغط هنا

Efficient Calibration of Radio Interferometers Using Block LDU Decomposition

60   0   0.0 ( 0 )
 نشر من قبل Ahmad Mouri Sardarabadi
 تاريخ النشر 2018
  مجال البحث فيزياء
والبحث باللغة English




اسأل ChatGPT حول البحث

Having an accurate calibration method is crucial for any scientific research done by a radio telescope. The next generation radio telescopes such as the Square Kilometre Array (SKA) will have a large number of receivers which will produce exabytes of data per day. In this paper we propose new direction-dependent and independent calibration algorithms that, while requiring much less storage during calibration, converge very fast. The calibration problem can be formulated as a non-linear least square optimization problem. We show that combining a block-LDU decomposition with Gauss-Newton iterations produces systems of equations with convergent matrices. This allows significant reduction in complexity per iteration and very fast converging algorithms. We also discuss extensions to direction-dependent calibration. The proposed algorithms are evaluated using simulations.



قيم البحث

اقرأ أيضاً

The calibration of modern radio interferometers is a significant challenge, specifically at low frequencies. In this perspective, we propose a novel iterative calibration algorithm, which employs the popular sparse representation framework, in the re gime where the propagation conditions shift dissimilarly the directions of the sources. More precisely, our algorithm is designed to estimate the apparent directions of the calibration sources, their powers, the directional and undirectional complex gains of the array elements and their noise powers, with a reasonable computational complexity. Numerical simulations reveal that the proposed scheme is statistically efficient at low SNR and even with additional non-calibration sources at unknown directions.
In order to meet the theoretically achievable imaging performance, calibration of modern radio interferometers is a mandatory challenge, especially at low frequencies. In this perspective, we propose a novel parallel iterative multi-wavelength calibr ation algorithm. The proposed algorithm estimates the apparent directions of the calibration sources, the directional and undirectional complex gains of the array elements and their noise powers, with a reasonable computational complexity. Furthermore, the algorithm takes into account the specific variation of the aforementioned parameter values across wavelength. Realistic numerical simulations reveal that the proposed scheme outperforms the mono-wavelength calibration scheme and approaches the derived constrained Cramer-Rao bound even with the presence of non-calibration sources at unknown directions, in a computationally efficient manner.
Searching for dispersed radio pulses in interferometric data is of great scientific interest, but poses a formidable computational burden. Here we present two efficient, new antenna-coherent solutions: The Chirpolator and The Chimageator. We describe the equations governing both techniques and propose a number of novel optimizations. We compare the implementation costs of our techniques with classical methods using three criteria: the operations rates (1) before and (2) after the integrate-and-dump stage, and (3) the data rate directly after the integrate-and-dump stage. When compared with classical methods, our techniques excel in the regime of sparse arrays, where they both require substantially lower data rates, and The Chirpolator requires a much lower post-integrator operations rate. In general, our techniques require more pre-integrator operations than the classical ones. We argue that the data and operations rates required by our techniques are better matched to future supercomputer architectures, where the arithmetic capability is outstripping the bandwidth capability. Our techniques are, therefore, viable candidates for deploying on future interferometers such as the Square Kilometer Array.
The Tianlai Pathfinder is designed to demonstrate the feasibility of using a wide field of view radio interferometers to map the density of neutral hydrogen in the Universe after the Epoch of Reionizaton. This approach, called 21~cm intensity-mapping , promises an inexpensive means for surveying the large-scale structure of the cosmos. The Tianlai Pathfinder presently consists of an array of three, 15~m $times$ 40~m cylinder telescopes and an array of sixteen, 6~m diameter dish antennas located in a radio-quiet part of western China. The two types of arrays were chosen to determine the advantages and disadvantages of each approach. The primary goal of the Pathfinder is to make 3D maps by surveying neutral hydrogen over large areas of the sky %$20,000 {rm deg}^2$ in two different redshift ranges: first at $1.03 > z > 0.78$ ($700 - 800$~MHz) and later at $0.21 > z > 0.12$ ($1170 - 1270$~MHz). The most significant challenge to $21$~cm intensity-mapping is the removal of strong foreground radiation that dwarfs the cosmological signal. It requires exquisite knowledge of the instrumental response, i.e. calibration. In this paper, we provide an overview of the status of the Pathfinder and discuss the details of some of the analysis that we have carried out to measure the beam function of both arrays. We compare electromagnetic simulations of the arrays to measurements, discuss measurements of the gain and phase stability of the instrument, and provide a brief overview of the data processing pipeline.
We describe techniques concerning wavelength calibration and sky subtraction to maximise the scientific utility of data from tunable filter instruments. While we specifically address data from the Optical System for Imaging and low Resolution Integra ted Spectroscopy instrument (OSIRIS) on the 10.4~m Gran Telescopio Canarias telescope, our discussion is generalisable to data from other tunable filter instruments. A key aspect of our methodology is a coordinate transformation to polar coordinates, which simplifies matters when the tunable filter data is circularly symmetric around the optical centre. First, we present a method for rectifying inaccuracies in the wavelength calibration using OH sky emission rings. Using this technique, we improve the absolute wavelength calibration from an accuracy of 5 Angstroms to 1 Angstrom, equivalent to ~7% of our instrumental resolution, for 95% of our data. Then, we discuss a new way to estimate the background sky emission by median filtering in polar coordinates. This method suppresses contributions to the sky background from the outer envelopes of distant galaxies, maximising the fluxes of sources measured in the corresponding sky-subtracted images. We demonstrate for data tuned to a central wavelength of 7615~$rmAA$ that galaxy fluxes in the new sky-subtracted image are ~37% higher, versus a sky-subtracted image from existing methods for OSIRIS tunable filter data.
التعليقات
جاري جلب التعليقات جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
mircosoft-partner

هل ترغب بارسال اشعارات عن اخر التحديثات في شمرا-اكاديميا