Tuesday, January 15, 2019

KiwiSDR IQ data streams with >20.25 kHz bandwidth

Recently I got to implement a way of recording data from KiwiSDRs at more than 12 kHz (or 20.25 kHz) by combining streams of IQ data using a polyphase filter bank (PFB) synthesizer. Luckily the two most complicated blocks, the PFB synthesizer and delay blocks, are already available in GNURadio.

What was missing was a block reading KiwiSDR IQ data with GNSS timestamps and a block determining the delays of between the different IQ sources. As all of these come from the same KiwiSDR running the same ADC clock, the offsets are multiples of samples.

For testing I started with a CODAR signal centered around 4820kHz using the AB1BD KiwiSDR. IQ wav files with GNSS timestamps were recorded on 4808, 4820, 4832 kHz, each with 12kHz bandwidth. In order to avoid wrapping around one of the inputs, a PFB synthesizer with 4 channels and 3 inputs was used:

GNURadio flowgraph

The screenshot below shows the three input signals samples with 12kHz and the reconstructed signal @48kHz:
  • Gain vs. frequency for the input signals is not flat. This is a known feature of the KiwiSDR FPGA code (missing CIC filter compensation filters). As a consequence there are dips in the combined spectrum
  • Nevertheless the CODAR sweeps can be seen to be continuous in the combined waterfall display

Top: waterfall diagrams for three IQ streams @12kHz; bottom: combined IQ stream @48kHz


Another test is to look at the instantaneous frequency in the combined IQ data stream @48kHz which looks fine:
Instantaneous frequency vs. time for the combined IQ data @48kHz


Obviously this has more applications than CODAR data analysis, as it opens up the possibility to analyze signals with up to 36kHz bandwidths. When a PFB synthesizer with 6 channels is used, up to 5 KiwiSDR IQ streams can be combined which enlarges the bandwidth to 60 kHz. Instead of using prerecorded wav files, the KiwiSDR GNURadio block can be used to combine streams of IQ samples in realtime.

The code for this will be made available as part of gr-kiwisdr after some cleanup.

Sunday, January 13, 2019

Radar 4980 kHz US East Coast (2)

Here is a small update on the HF radar signals which can be received on KiwiSDRs located in the US east-coast area
  • The center frequency seems to be close to 4900 kHz and the bandwidth is 50-60kHz
  • The signals are active every few minutes for about 40 seconds
  • Up to four trains of pulses can be seen. These pulses repeat every 0.01 sec but the autocorrelation function peaks at 0.02 sec indicating that the pulses are phase modulated
  • The phase of all pulses is flipped by 180 degrees in each 0.02 sec slot
  • Each pulse train consists of at least 5 pulses
  • The fast pulse repetition rate suggests this is a radar for fast-moving targets, as opposed to, e.g., CODAR


abs(IQ) in 0.02 sec frames


abs(IQ) in 0.02 sec frames (zoom)


pulse trains w.r.t. GPS time


pulse trains w.r.t. GPS time (zoom)



arg(IQ) zoom in 0.01 sec frames, ×2 resampling

abs(IQ) zoom in 0.01 sec frames, ×2 resampling

arg(IQ) for two consecutive 0.01 sec frames, ×2 resampling

Tuesday, January 8, 2019

Radar 4980 kHz US East Coast

While monitoring WBHF signals on 4950 kHz, a radar-like signal was spotted around 4980 kHz. Note that its bandwidth exceeds the maximum bandwidth of the used KiwiSDR (20250kHz), so the full characteristics of this signal cannot be determined. It was active every few minutes but I was not patient enough to determine its schedule. Within the limited bandwidth the following can be seen:

1) The auto-correlation has peaks at multiples of 405/20250Hz = 0.02 sec = 1/50Hz:

autocorrelation

2) Plotting abs(IQ) in frames of 1/50Hz reveals a pair of pulses, each apparently consisting of three separate sub-pulses:

abs(IQ) in terms of 1/50Hz frames

3) TDoA analysis locates this is a real signal somewhere in upstate NY (and demonstrates that it is not local QRM):
KiwiSDR TDoA 4980 kHz
(This is an example of successfully using KiwiSDRs with different sampling rate in TDoA, cf.,  http://valentfx.com/vanilla/discussion/1424/possible-tdoa-problem-when-using-kiwis-in-3-ch-mode)

If anyone knows what kind of HF radar signal this is I would be grateful to know.