1152 |
fluktuacia |
1 |
\chap Results obtained in the trial version
|
1142 |
kaklik |
2 |
|
1152 |
fluktuacia |
3 |
The trial version construction was tested for proper handling of sampling rates in the range of 5 MSPS to 15 MSPS, but it should work even above this limit. The system works on i7 8 cores computer with Ubuntu 12.04 LTS operating system. Data recording of input signal is impossible above the sampling rates of around 7 MSPS due to bottleneck at HDD speed limits, but it should be resolved by the use of SSD disk drive. However, such design has not been tested in our setup.
|
1142 |
kaklik |
4 |
|
|
|
5 |
\sec Measured parameters
|
|
|
6 |
|
|
|
7 |
Two prototypes of ADC modules were assembled and tested. The first prototype, labeled ADC1, has LTC2190 ADC chip populated with LT6600-5 front-end operational amplifier. It also has a 1kOhm resistors populated on inputs which give it an ability of an internal attenuation of the input signal. The value of this attenuation $A$ is calculated by the following formula \ref[ADC1-gain]
|
|
|
8 |
|
|
|
9 |
\label[ADC1-gain]
|
|
|
10 |
$$
|
|
|
11 |
A = {806 \cdot R_1 \over R_1 + R_2} \eqmark
|
|
|
12 |
$$
|
|
|
13 |
|
|
|
14 |
Where the letters stand for the following:
|
|
|
15 |
\begitems
|
|
|
16 |
* $A$ - Gain of an input amplifier.
|
|
|
17 |
* $R_1$ - Output impedance of signal source (usually 50 Ohm).
|
|
|
18 |
* $R_2$ - Value of serial resistors at operational amplifier inputs.
|
|
|
19 |
\enditems
|
|
|
20 |
|
|
|
21 |
We have $R_2 = 1000 \Omega$ and $R_1 = 50 \Omega$ which imply that $A = 0.815$. That value of A was further confirmed by the measurement.
|
|
|
22 |
In our measurement setup we have H1012 Ethernet transformer connected to inputs of ADC. We have used this transformer for signal symetrization from BNC connector at Agilent 33220A signal generator. Circuit diagram of the used transformer circuit is shown in picture \ref[balun-circuit] and circuit realization in photograph \ref[SMA2SATA-nest].
|
|
|
23 |
|
|
|
24 |
\midinsert
|
|
|
25 |
\clabel[balun-circuit]{Balun transformer circuit}
|
|
|
26 |
\picw=10cm \cinspic ./img/SMA2SATA.pdf
|
|
|
27 |
\caption/f Simplified balun transformer circuit diagram.
|
|
|
28 |
\endinsert
|
|
|
29 |
|
|
|
30 |
The signal generator Agilent 33220A which we used, does not have optimal parameters for this type of dynamic range measurement. Signal distortion and spurious levels are only -70 dBc according to Agilent datasheet \cite[33220A-generator]. We have managed to measure an ADC saturation voltage of 705.7 mV (generator output) with this setup, mostly due to an impedance mismatch and uncalibrated measurement setup, with 1V ADC range selected by sense pin. This is a relatively large error, but the main result of our measurement, seen as a FFT plot shown in image \ref[ADC1-FFT], confirms $>$80 dB dynamic range at ADC module input.
|
|
|
31 |
|
|
|
32 |
\midinsert
|
|
|
33 |
\clabel[ADC1-FFT]{ADC1 sine test FFT}
|
|
|
34 |
\picw=15cm \cinspic ./img/screenshots/ADC1_CH2_FFT.png
|
|
|
35 |
\caption/f Sine signal sampled by ADC1 module with LTC2190 and LT6600-5 devices.
|
|
|
36 |
\endinsert
|
|
|
37 |
|
|
|
38 |
|
|
|
39 |
Similar test was performed at ADC2 module. For ADC2 we have to use formula with a different constant \ref[ADC1-gain]. The ADC2 module has LT6600-2.5 amplifiers populated on it with a gain equal to $A = 2.457$ and uses the same $R_2$ resistors. We measured saturation voltage of 380 mV (generator output) at channel 1 on this ADC. It is well within the parameter tolerances of the used setup.
|
|
|
40 |
|
|
|
41 |
\label[ADC2-gain]
|
|
|
42 |
$$
|
|
|
43 |
A = {1580 \cdot R_1 \over R_1 + R_2} \eqmark
|
|
|
44 |
$$
|
|
|
45 |
|
|
|
46 |
Where the letters stand for the following:
|
|
|
47 |
\begitems
|
|
|
48 |
* $A$ - Gain of an input amplifier.
|
|
|
49 |
* $R_1$ - Output impedance of signal source (usually 50 Ohm).
|
|
|
50 |
* $R_2$ - Value of serial resistors at operational amplifier inputs.
|
|
|
51 |
\enditems
|
|
|
52 |
|
|
|
53 |
\midinsert
|
|
|
54 |
\clabel[ADC2-FFT]{ADC2 sine test FFT}
|
|
|
55 |
\picw=15cm \cinspic ./img/screenshots/ADC2_CH1_FFT.png
|
|
|
56 |
\caption/f Sine signal sampled by ADC2 module with LTC2271 and LT6600-2.5 devices.
|
|
|
57 |
\endinsert
|
|
|
58 |
|
|
|
59 |
Computed FFT spectra for measured signal are shown in the images \ref[ADC2-FFT] and \ref[ADC1-FFT]. Both images confirm that ADCdual01A modules have input dynamical range of at least 80 dB.
|
|
|
60 |
|
|
|
61 |
\midinsert
|
|
|
62 |
\clabel[SMA2SATA-nest]{Used balun transformer}
|
|
|
63 |
\picw=15cm \cinspic ./img/SMA2SATA_nest1.JPG
|
|
|
64 |
\caption/f Balun transformer circuit used for ADC parameters measurement. It is constructed from H1012 transformer salvaged from an old Ethernet card.
|
|
|
65 |
\endinsert
|
|
|
66 |
|
|
|
67 |
|
|
|
68 |
\sec Example of usage
|
|
|
69 |
|
1147 |
kaklik |
70 |
At current state the constructed radioastronomy digitization unit paired with SDRX01B receiver module could be used in several experiments. We describe overall ideas of these experiments and show preliminary results in cases where we obtain the data.
|
1142 |
kaklik |
71 |
|
1147 |
kaklik |
72 |
\secc Simple polarimeter station
|
|
|
73 |
|
|
|
74 |
If we use two antennas with different linear polarization (Crossed Yagi antennas for example), we should determine polarization state of received signal. Such kind of measurement is useful if we need an additional information about reflection to distinguish between targets. This configuration needs more complicated antenna configuration and we had no experience with this type of observation, so we have not implemented this experiment.
|
|
|
75 |
|
1142 |
kaklik |
76 |
\secc Basic interferometric station
|
|
|
77 |
|
|
|
78 |
Interferometry station was chosen to serve as the most basic experimental setup. We connected the new data acquisition system to two SDRX01B receivers. Block schematics of the setup used is shown in the image \ref[block-schematic]. Two ground-plane antennae were used and mounted outside the balcony at CTU building at location 50$^\circ$ 4' 36.102'' N, 14$^\circ$ 25' 4.170'' E.
|
|
|
79 |
Antennae were equipped with LNA01A amplifiers. All coaxial cables had the same length of 5 meters. Antennae were isolated by common mode ferrite bead mounted on cable to minimise the signal coupling between antennas. Evaluation system consisted of SDGPSDO local oscillator subsystem used to tune the local oscillator frequency.
|
|
|
80 |
|
|
|
81 |
\midinsert
|
|
|
82 |
\clabel[block-schematic]{Receiver block schematic}
|
|
|
83 |
\picw=\pdfpagewidth \setbox0=\hbox{\inspic ./img/Basic_interferometer.png }
|
|
|
84 |
\par\nobreak \vskip\wd0 \vskip-\ht0
|
|
|
85 |
\centerline {\kern\ht0 \pdfsave\pdfrotate{90}\rlap{\box0}\pdfrestore}
|
|
|
86 |
\caption/f Complete receiver block schematic of dual antenna interferometric station.
|
|
|
87 |
\endinsert
|
|
|
88 |
|
1145 |
kaklik |
89 |
Despite of the schematic diagram proposed at beginning of system description \ref[expected-block-schematic].
|
1147 |
kaklik |
90 |
We have used two separate oscillators -- one oscillator drives ENC signal to ADCs still through FPGA based divider and the other one drives it to SDRX01B mixer.
|
1142 |
kaklik |
91 |
The reason for this modification was an attempt to simplify the frequency tuning during the experiment. A single oscillator may be used only with a proper setting of FPGA divider and this divider may be modified only by recompilation of FPGA code and loading/flashing a new FPGA schema. Due to fact that the FPGA is connected to PCI express and kernel drivers with hardware must be reinitialized, reboot of PC is required every time a FPGA scheme is changed. Instead of this complicated procedure, we set the FPGA divider to a constant division factor of 30 and used another district oscillator for ADCdual01 sampling modules and for SDRX01B receiver.
|
|
|
92 |
We have used ACOUNT02A MLAB instrument for frequency checking of correct setup on both local oscillators.
|
|
|
93 |
|
|
|
94 |
\midinsert
|
|
|
95 |
\clabel[phase-difference]{Phase difference}
|
|
|
96 |
\picw=10cm \cinspic ./img/screenshots/phase_difference.png
|
|
|
97 |
\caption/f Demonstration of phase difference between antennae.
|
|
|
98 |
\endinsert
|
|
|
99 |
|
1147 |
kaklik |
100 |
For the simplest demonstration of phase difference between antennae, we have analysed part of the signal by complex conjugate multiplication between channels. Results of this analysis can be seen in the following picture \ref[phase-difference]. Points of the selected part of the signal create a clear vector, which illustrates the presence of the constant phase difference determined by RF source direction.
|
1142 |
kaklik |
101 |
|
|
|
102 |
\secc Simple passive Doppler radar
|
|
|
103 |
|
1147 |
kaklik |
104 |
If we use an existing transmitter with known carrier frequency and proper antenna, we can detect flying object as signals surrounding the transmitter carrier frequency. We planned this experiment with the same station configuration as was described in section \ref[expected-block-schematic]. The ISS \glos{ISS}{International Space Station} as object and GRAVES radar transmitter were selected as adequate testing objects (We know ISS reflections from previous experiments). This experiment could be realised by previously described interferometer station, but unfortunately we missed the suitable orbit pass due to technical lacks with station configuration.
|
1142 |
kaklik |
105 |
|
1147 |
kaklik |
106 |
\secc Meteor detection station
|
1142 |
kaklik |
107 |
|
1147 |
kaklik |
108 |
The same observational station configuration should be used for meteor detection system \cite[mlab-rmds]. We used the GRAVES radar as suitable signal source and monitored its carrier frequency. GRAVES radar is located in France therefore we could not see its direct carrier signal, but meteors reflect it signal and as consequence we could easily detect meteor presence as reflection appearance. One meteor detected by this method is shown in picture \ref[meteor-reflection].
|
1142 |
kaklik |
109 |
|
1147 |
kaklik |
110 |
\midinsert
|
|
|
111 |
\clabel[meteor-reflection]{Meteor reflection}
|
|
|
112 |
\picw=10cm \cinspic ./img/screenshots/observed_meteor.png
|
|
|
113 |
\caption/f Meteor reflection (the red spot in centre of image) received by an evaluation design.
|
|
|
114 |
\endinsert
|
1142 |
kaklik |
115 |
|
1147 |
kaklik |
116 |
|
1142 |
kaklik |
117 |
\chap Proposition of the final system
|
|
|
118 |
|
|
|
119 |
The construction of the final system, that is supposed to be employed for real radioastronomy observations will be described in this chapter. It is mainly a theoretical analysis of the data handling systems. Realization of the described ideas might be possible as a part of our future development after we fully evaluate and test the current trial design.
|
|
|
120 |
|
|
|
121 |
The system requires proper handling of huge amounts of data and either huge and fast storage capacity is needed to store the captured signal data, or enormous computational power is required for online data processing and filtering. Several hardware approaches currently exist and are in use for data processing problem handling. Either powerful multi gigahertz CPUs, GPUs, FPGAs, or specially constructed ASICs are used for this task.
|
|
|
122 |
|
|
|
123 |
\sec Custom design of FPGA board
|
|
|
124 |
|
|
|
125 |
In the beginning of the project, a custom design of FPGA interface board had been considered. This FPGA board should include PCI express interface and should sell at lower price than the trial design. It should be compatible with MLAB internal standards which are further backward compatible with the existing or improved design of ADC modules. For a connection of FPGA board to another adapter board with PCIe we expect a use of a PCIe host interface.
|
|
|
126 |
Thunderbolt technology standard was expected to be used in this PC to PCIe module communication which further communicates with MLAB compatible FPGA module. Thunderbolt chips are currently available on the market for reasonable prices \cite[thunderbolt-chips]. However, a problem lies in the accessibility to their specifications, as they are only available for licensed users and Intel has a mass market oriented licensing policy, that makes this technology inaccessible for low quantity production. As a consequence, an external PCI Express cabling and expansion slots should be considered as a better solution, if we need to preserve standard PC as a main computational platform.
|
|
|
127 |
|
|
|
128 |
However, these PCI express external systems and cables are still very expensive. Take Opal Kelly XEM6110 \cite[fpga-pcie] as an example, with its price tag reaching 995 USD at time of writing the thesis. Therefore, a better solution probably needs to be found.
|
|
|
129 |
|
|
|
130 |
An interface problem will by probably resolved by other than Intel ix86 architecture. Many ARM computers have risen on market due to an increased demand of embedded technologies, which require high computation capacity, low power consumption and small size -- especially smartphones. Many of those ARM based systems have interesting parameters of signal processing. These facts make Intel's ix86 architecture unattractive for future projects.
|
|
|
131 |
|
|
|
132 |
\sec Parralella board computer
|
|
|
133 |
|
|
|
134 |
Parallella is a new product created by Adapteva, Inc. \cite[parallella-board]. It represents a small supercomputer, that has been in development for almost two years with only testing series of boards produced until now (first single-board computers with 16-core Epiphany chip were shipped in December 2013) \cite[parallella-board]. The board has nearly ideal parameters for signal processing (as it provides around 50 GFLOPS of computational power). It is is equipped with Epiphany coprocessor which has 16 High Performance RISC CPU Cores, Zynq-7020 FPGA with Dual ARM® Cortex™-A9 MPCore™ and operating frequency of 866 MHz, 1GB RAM, 85K Logic Cells, 10/100/1000 Ethernet and OpenCL support \cite[parallella16-board]. In addition to this, the board consumes only 3 Watts of power if both Zynq and Epiphany cores are running simultaneously.
|
|
|
135 |
|
|
|
136 |
The main disadvantage of Parralella board is its unknown lead time and an absence of SATA interface or other interface suitable for data storage connection. Fast data storage interface would be useful and would allow bulk processing of captured data. Following that, the results of data processing may be sent over the Ethernet interface to data storage server.
|
|
|
137 |
|
|
|
138 |
\midinsert
|
|
|
139 |
\clabel[img-parallella-board]{Parallella board overview}
|
|
|
140 |
\picw=15cm \cinspic ./img/ParallellaTopView31.png
|
|
|
141 |
\caption/f Top view on Parallella-16 board \cite[parallella16-board].
|
|
|
142 |
\endinsert
|
|
|
143 |
|
|
|
144 |
If Parallella board will be used as a radioastronomy data interface, there would be a demand for new ADC interface module. The interface module will use four PEC connectors mounted on the bottom of the Parallella board. This daughter module should have MLAB compatible design and should preferably be constructed in the form of separable modules for every Parallella's PEC connector.
|
|
|
145 |
|
|
|
146 |
\sec GPU based computational system
|
|
|
147 |
|
|
|
148 |
A new GPU development board NVIDIA K1, shown in the following picture \ref[img-NVIDIA-K1], has recently been released. These boards are intended to be used in fields including computer vision, robotics, medicine, security or automotive industry. They have good parameters for signal processing for a relatively low price of 192 USD. Unfortunately, they are currently only in pre-order release stage (in April 2014).
|
|
|
149 |
|
|
|
150 |
\midinsert
|
|
|
151 |
\clabel[img-NVIDIA-K1]{NVIDIA Jetson TK1 Development Kit}
|
|
|
152 |
\picw=15cm \cinspic ./img/Jetson_TK1_575px.jpg
|
|
|
153 |
\caption/f The NVIDIA Jetson TK1 Development Kit \cite[nvidia-k1].
|
|
|
154 |
\endinsert
|
|
|
155 |
|
|
|
156 |
NVIDIA board differs from other boards in its category by a presence of PCI Experess connector. If we decide to use this development board in our radio astronomy digitalisation system, the PCI express should be used for FPGA connection. A new FPGA board with PCI Express direct PCB connector
|
|
|
157 |
|
|
|
158 |
% doplnit popis pripojeni FPGA desky s HDMI Kabelem.
|
|
|
159 |
|
|
|
160 |
\sec Other ARM based computation systems
|
|
|
161 |
|
|
|
162 |
Other embedded ARM based computers, for example ODROID-XU, lack a suitable high speed interface. Their highest speed interface is USB 3.0 which has currently unsettled development support and needs commercial software tools for evaluation and testing.
|
|
|
163 |
|
|
|
164 |
|
|
|
165 |
From the summary analysis mentioned above, the Parrallella board seems to be a best candidate for computational board in radioastronomy data acquisition system, as it is optimised for high data flow processing. On one hand, Parrallella does not have much memory to cache the processing data but on the other hand it has wide bandwidth data channels instead. Other boards might provide much more computational power -- 300 GFLOPS in case of NVIDIA K1, but they are optimised for heavy computational tasks on limited amount of data which represents a typical problem in computer graphics. However, in our application we do not need such extreme computation power at data acquisition system level.
|
|
|
166 |
As a result we should presumably wait until Parallella becomes widely available. Following that, a new ADCdual interface board should be designed and prepared to be used in new scalable radio astronomy data acquisition system. In the meantime, before suitable computing hardware become accessible, the required applications and algorithms should be optimised using the proposed trial version with FPGA development board on standard PC host computer (having a PCI Express interface to development board).
|
|
|
167 |
|