We're sorry but this page doesn't work properly without JavaScript enabled. Please enable it to continue.
Feedback

Monitoring the ionosphere altitude variation with a sound card

00:00

Formal Metadata

Title
Monitoring the ionosphere altitude variation with a sound card
Subtitle
software defined radio processing of DCF-77 signals
Title of Series
Number of Parts
611
Author
License
CC Attribution 2.0 Belgium:
You are free to use, adapt and copy, distribute and transmit the work or content in adapted or unchanged form for any legal purpose as long as the work is attributed to the author in the manner specified by the author or licensor.
Identifiers
Publisher
Release Date
Language
Production Year2017

Content Metadata

Subject Area
Genre
Abstract
[DCF-77] is a German very low frequency (VLF) emitter locked on the Cs clocks of PTBused for synchronizing radiofrequency-disciplined clocks found, for example,in low-cost weather stations. VLF signals propagate through the waveguidewhose boundary conditions are defined on the one hand by ground, and on theother hand by the ionosphere (charged) layer altitude. In addition to theamplitude modulation of the Cs-locked DCF-77 frequency standard, a [phasemodulation] is imprinted on the carrierfor time of flight measurement, allowing for "precise" time of flightmeasurements. The topic of the presentation is 1. reception of the DCF-77 VLF signal using a coil antenna, with enough power to feed a lock-in amplifier and extract phase and magnitude information. When the reference signal of the lock-in is referred to a local Cs primary standard, day/night ionosphere altitude variations are readily observed -- here using the carrier phase analysis, 2. extraction of the time of flight through CDMA processing of the phase output from the lock-in amplifier, recorded on a (low frequency) oscilloscope, emphasizing the time resolution improvement of the phase modulation with respect to the amplitude modulation gained from the increased signal bandwidth, 3. replace the lock-in amplifier with software-defined radio processing of signals recorded with a personal computer sound card or the analog to digital converter (RTL2832U) of a DVB-T receiver dongle, 4. refer the second channel of the stereo sound card/DVB-T to a low-cost GPS 1-PPS time reference for local oscillator drift compensation, and hence measure the ionosphere altitutde variation through time of flight variation by comparing the phase-encoded cross-correlation peak position with the 1-PPS position, yielding results consistent with the lab-grade instrumentation. Although the whole processing chain is trivial, various issues making itspractical implementation challenging will be emphasized. Most significantly,the processing chain is very similar to those applied to [GPSsignals], yet easier tograsp with the strong VLF signal. We aim at analyzing the time stability ofthe phase modulated signal -- allowing for sub-100 us timing resolution -- andcorrelated observed phase fluctuations with [ionospherealtitude]
Bounded variationPlastikkarteRadio-frequency identificationExecution unitConvex hullTrigonometryMassMathematicsSpecial unitary groupReflection (mathematics)Charge carrierRandomizationPresentation of a groupRaw image formatFrequencyPropagatorWaveFunction (mathematics)InformationStability theorySpectrum (functional analysis)QuicksortWellenwiderstand <Strömungsmechanik>TouchscreenPhase transitionSound cardMultiplication signSign (mathematics)MultilaterationMusical ensembleAdaptive behavior2 (number)Digital electronicsPosition operatorSound effectCalculus of variationsPhysical systemWorkstation <Musikinstrument>CodePhysicalismResonatorSource codeAtomic numberSlide ruleSpacetimeObject (grammar)Matter waveBitApproximationNoise (electronics)VelocityAdditionMeasurementOrder (biology)OscillationBounded variationPhysical quantitySurfaceAtomic clockRadical (chemistry)Software-defined radioContext awarenessMetrologieSoftware testingImage resolutionContent (media)FreewareVacuumOrder of magnitudeTime zoneGreatest common divisorMetreDaylight saving timeSphereXML
Amsterdam Ordnance DatumConvex hullGamma functionMaxima and minimaFunction (mathematics)FrequencyInverter (logic gate)Row (database)Raw image formatPhase transitionMultiplication signSynchronizationOrder of magnitude2 (number)Source codeBounded variationDrop (liquid)Different (Kate Ryan album)Pulse (signal processing)Connectivity (graph theory)Sign (mathematics)OscillationCharge carrierInformationCASE <Informatik>Sound cardTheoremSampling (statistics)SineQuicksortGoodness of fitTowerResultantSoftware-defined radioWeb pageProtein foldingPseudozufallszahlenThomas BayesCodeCross-correlationAtomic clockReal numberBitNoise (electronics)1 (number)SequenceDegree (graph theory)OctaveProcess (computing)Arithmetic meanMetreLogic gateLocal ringPhysikalisch-Technische BundesanstaltElectronic visual displayComputer animationProgram flowchart
Twin primeSample (statistics)FrequencyConvex hullHill differential equationExecution unitClique-widthValuation (algebra)Data modelMultiplication signPhase transitionSoftware protection dongle2 (number)Process (computing)Sampling (statistics)BitFilter <Stochastik>DecimalInformationElectronic visual displayOrder of magnitudePoint (geometry)FrequencySoftware-defined radioCoefficientArithmetic meanFunction (mathematics)Hacker (term)Row (database)TimestampBand matrixMultiplicationData conversionDirection (geometry)Cross-correlationFunctional (mathematics)Block (periodic table)QuantificationPlastikkarteWaveAverageSpring (hydrology)MeasurementLinear regressionSpecial unitary groupResultantSinc functionFourier transformSound cardQuicksortThermal fluctuationsDifferent (Kate Ryan album)MedianEstimatorPotenz <Mathematik>Stability theoryOscillationImage resolutionINTEGRALMathematicsBounded variationAlgorithmMetrologieStandard deviationException handlingMusical ensembleVirtualizationoutputOpen sourceSet (mathematics)AreaFinite impulse responseTowerFreewareGradientGreatest element
Image resolutionSampling (music)Data modelLibrary (computing)Maxima and minimaStability theoryCross-correlationMultiplication signWaveMeasurementMoment (mathematics)Content (media)Population densityChemical equationAtomic clockSound effectSpring (hydrology)Scripting languageOctaveSatelliteSound cardSource codeLaptopINTEGRALSoftware protection dongleSoftware-defined radioResultantWebsiteEmailPresentation of a groupPhase transitionPosition operatorPoint (geometry)Group actionBoundary value problemTunisData miningThermal fluctuationsConsistencyTerm (mathematics)TelecommunicationEndliche ModelltheorieSinc functionData transmissionFrequencyCurve3 (number)Flow separationSurfaceSlide ruleInformationFunctional (mathematics)Bounded variationQuicksortMicrocontrollerShape (magazine)BitPhysikalisch-Technische BundesanstaltAutocorrelationFormal languageoutputFilter <Stochastik>CodePseudozufallszahlenBinary fileField (computer science)Process (computing)Band matrixComputer animation
Scripting languageProcess (computing)ImplementationOctaveInterface (computing)PlastikkarteSoftware-defined radioPseudozufallszahlenRandom number generationSource codeMappingAutocorrelationLine (geometry)Point (geometry)Cross-correlationNumberOscillationGodINTEGRALLocal ringBitTheory of relativity2 (number)Phase transitionComputer animationLecture/Conference
Computer animation
Transcript: English(auto-generated)
Last talk. So Martin decided to put this in academia. And actually, this is more of a hobby than my actual work. So what we wanted to work with my colleague, Enrico Rubiola, is using a very high stability clock
source, which is the DCF77 low frequency, very low frequency emitter in Germany, for measuring ionosphere altitude variations. And that was actually decided during one of our barbecue parties in the lab when Enrico was discussing with me about his PhD on measuring Lorentz signals and being
able to separate the ionospheric wave from the surface ground propagating wave. So basically, what we wanted to say is, OK, he did this like 25, 30 years ago as a PhD. Can we do it nowadays with a sound card? So that was the original insight into this topic. And maybe some of you remember that a couple of years ago,
I discussed GPS reception. Now, GPS is a bit challenging because it's below terminal noise. You cannot see it on a spectrum analyzer. But the basic insight was that if you could collect the raw carrier signal from the GPS, well, GPS is basically a collection of atomic clocks in the space. And you could do some nice physics.
So a bit later now, Laurenti Starkey here from TESS has published in one of his reviews what physics you can do on GPS usage. And this is all based on GNSS DR. So basically, the insight in this work is if you have a very high stability clock, and because frequency is a physical quantity
that you can measure with the highest resolution, then maybe you can do something fun. And so in this context, what we wanted to look at is looking at DCF77 as a signal source that is locked or disciplined on the PTBs of the German National Metrology Center atomic clocks, cesium clocks.
And actually, that was all started during the EFTS, the European Frequency and Time Seminar, that is organized in Boussinesq every summer, when Andreas Bauer, actually the acknowledgement is here below, showed us one of these slides where he was presenting the code system that was implemented in DCF77.
So you might know that DCF77 is one of these radio frequency emitters that allows you to synchronize some of these low cost weather station clocks. And this is just done by amplitude modulation. I'll show you a bit later. But what Andreas showed during his presentation was that if you just look at the carrier with this amplitude modulation, you have a narrow band signal, which is basically just
the timing signal once every second. And you're not going to do good positioning with this because you have a very high stability frequency, but your timing capability is very poor because nothing looks more similar from one sign period than another sign period. So you have no timing capability. And what they did at PTB, the people running DCF77,
is that they added a sort of random phase noise over the carrier, which is amplitude modulated. And when you've been working for a year on GPS and reading all the documentation of GPS, that just jumps on you as a typical spectrum spreading using phase modulation and sort of random phase modulation.
So when we had a look at this with Enrico, we said we have to try to do this and do this as easy as possible. So what we want to do here is to have access to the raw signal from DCF77 and implement not the basic amplitude modulation detection that everyone has been doing maybe, but the more fancy phase modulation detection.
And why did we do this? Well, if you look at the literature, it's been very well known. And I think that people are forgetting this a little bit, actually. I have never been introduced to this before, this work, is that in the 60s, people were very much interested in very low frequency propagation. They were looking at ionosphere altitude because these waves are bouncing over ionosphere. And if you look into all the literature, you see that between nighttime and daytime,
ionosphere changes altitude because ionization due to cosmic rays from the sun will change the electronic density in the ionosphere. And basically, you have a 20 kilometer altitude variation in addition to reflection capability change due to electron density.
So what we want to do here in this work is we have this very high-stability clock source located in Manflingen, Germany. And we are located in Bezanson, 370 kilometers away. And we have this beautiful cesium-locked oscillator that generates a carrier wave with a 50 kilowatt emitter.
And we wonder whether we could do something more interesting than simply synchronizing a weather station clock by looking at the time of flight of this wave. And this requires, of course, accurate timing because carrier frequency is one story. But if you want to do some timing, you need to have some sort of time of flight measurement. And this is where the carrier phase
measurement will be involved. So just to get an order of magnitude, what are we looking for? We are considering DCF 77, 370 kilometer away from me, a ground wave propagating around 100 microseconds. And if you make this wave bounce over ionosphere, it will take another 100 microseconds additional
to go to the ionosphere and back to me. And if the ionosphere moves by about 20 kilometer, you could expect some variation of day to night, variation of about an additional 95 microseconds. So the question is, can we time the signal generated by DCF 77 received in Bezanson by better than,
or much better than 100 microseconds? That's a question with a few approximations. Will you consider that the wave is propagating at velocity in a vacuum, and you assume that the ground wave delay is negligible with respect to what we're looking at? So that's really the objective of this work.
So the first thing I would like to introduce, I know this is not a top session about antennas. But again, for someone who's been raised with VHF antenna, microwave antenna, for me, an antenna is a 50 ohm or 75 ohm piece of wire, which is usually resonant. And when you go back to this three kilometer wavelength signal, any antenna will be much, much shorter
than the wavelength. And if you go back in the literature, you know, or you might find, that such antennas are necessarily high quality factor antennas, which for someone raised in VHF and SHF is usually a poor thing. You don't like high resonance antenna because they are narrow band. And usually we want wide band signals to have a lot of information content, like Shannon told us.
And so on the one hand, we have these narrow band antennas. But actually narrow band is good for SDR because narrow band in SDR will allow you to select only one signal and get rid of all the very low frequency signals that might prevent us from receiving DCS-77. So we have this high quality factor antenna,
which is simply a coil which is made resonant. The other thing that is not usual to us, or at least to me, is that these antennas are very high impedance. If you calculate the impedance at the anti-resonance of the system, it's 175, 125 kilo ohm. So you need an adaptation circuit.
So just fill the effect transistor that has a very high impedance that will generate a low impedance output. And this you can feed directly to the sound card of your PC. Now if you do this in the lab, well, you get a lot of signals. Basically because your cathode ray screen is generating signal. Your switching supply is generating. That was the first discussion this morning.
So you need to take the antenna outdoor. So you go through the emergency exit outdoor that gets out of my lab, put the antenna outside. And if I put the antenna like three, five kilometer away from the lab, then I get my signal from DCS-77. So at least with respect to GPS, with DCS-77,
I can take a spectrum analyzer and just have a look at the signal, and I see it. So now I have a signal. I want to decode it. So I steal a lock-in amplifier from one of my colleagues. I take a synthesizer. And I simply generate 77.5 kilohertz, the DCS-77 carrier frequency. And surely enough, on the oscilloscope,
I see my amplitude modulated signal. These are the amplitude drops once every second. And I have the phase information, which comes out of the lock-in. And because there's no reason for this receiver to be locked to a cesium clock, the local oscillator of this synthesizer has no reason to being locked on the cesium.
Of course, you have a phase, which depends on the thing we want to decode, plus a time-varying signal, which is a frequency offset, as we've just seen. But because I'm in a time and frequency lab, I'm lucky enough to have access to a cesium clock. So now I take a cesium clock. I lock the synthesizer on a cesium clock. And surely enough, PTB has a cesium clock,
which is locked to the Bezant-Sou observatory clock. Of course, these two clocks are used to generate time. So these are synchronized, or at least they have a very slow drift. So basically, we can get amplitude modulation, phase modulation. And this is just my lock-in amplifier output. There's not much I can do out of this.
So I take this on an oscilloscope. I grab the oscilloscope signal. And by recording the raw signal from DCS-77, can I extract amplitude phase? And if I get the phase, can I do a cross-correlation so that I can recover my pseudorandom information? So here is the pseudorandom code
that has been added on top of DCS-77 carrier. So it's a pseudorandom. It's a 511-bit, so it's very similar to GPS. GPS is 10-bit. This one is 9-bit. Pseudorandom means that it doesn't repeat over this 511 sequence. And you have a plus or minus 13-degree phase variation,
as opposed to GPS, which is plus or minus 0, 180. This is just a very tiny phase variation, plus or minus 13-degree. So you generate this code. This is stolen from the DCS-77 web page from Wikipedia. You generate the code. You load the code.
You re-sample your signal to have as many 1's and 0's as you have samples in each bit. You cross-correlate. And surely enough, this is one of these Eureka days, where everything works nicely. Once every second, you've got one of these peaks coming out, cross-correlation peak. And you might say, no, this is just noise. Well, if you zoom into one of these peaks, you have a real cross-correlation signal
that is significant. So this is an example of getting the signal by cross-correlating my pseudorandom code over the phase information that I've recorded. And I can get a nice timing signal. And if you wonder why we're doing this,
this is the amplitude modulation. I zoomed into the amplitude modulation. This is one of my cross-correlation peaks. If you look at this over one minute, this is what you have once every second. You have a drop in amplitude. Once every second, you have a cross-correlation peak. And if you look at them, you see that obviously we're going to have a much better timing accuracy with such a signal than with this one. If you do this over one minute,
you see that you have a beautiful cross-correlation peak which is very narrow. And here you have this broad drop of the amplitude. So basically what has been published earlier is that you have at least a tenfold improvement in timing capability by using the phase cross-correlation as opposed to amplitude modulation. So far I've shown off with my fancy instruments.
I have an experiment which is a few thousand zero. And you need a cesium clock. Quite a few people here will not be able to perform it. So now let's go to software-defined radio. I want to do this now. I have my locking amplifier with my cesium locked source. And I want to compare the result with a sound card. And again, my sound card has a local oscillator
which has no reason of being locked on cesium. So I need to compare my sound card recording with some sort of time reference. And the obvious time reference is one PPS, one pulse per second coming from a GPS. I have one of these cheap U-blocks receivers. U-blocks is a Swiss company making very good GPS receivers.
This is for 80 euro, you can get a phase lock GPS receiver with your one PPS signal, which is measured to be plus or minus 30 nanoseconds. So you have a very good time reference. And so again, we do the same thing. We have our audio signal because I have a sound card that can record at 192 kilohertz. I am weaving Nyquist's theorem of two samples
for a 77.5 kilohertz carrier frequency. I do a frequency translating fear filter. And I just record magnitude and phase out of this signal. And if I do this, I will try to tune the frequency translating fear filter so that I
compensate for this 192 kilohertz offset from the wanted frequency. Now, as I just mentioned during the session panel, one of the things that always baffles me with this translating fear filter is, should I put plus 77 or minus 77? I always do it wrong and I have to try twice. What is the sign of the frequency offset that you have to put in the frequency filter?
And somehow, in this case, it worked. Whatever I put, plus or minus 77.5. And this is obvious once you realize that you've recorded a real signal. This is not an IQ output, a complex output. It's a real signal. A real signal is even. So you have the same component on the positive frequency and the negative frequency. So whether you put plus 77 or minus 77,
you always end up having one of the information components into the baseband. And so for this, once we can do whatever we want, plus or minus sign, and you end up having a nice signal. So you do this. And the first time you do this, you get once every second the amplitude drop. And this is your phase from the recording.
And you try playing with this delta f that you cannot reach here. So originally, I'm slightly off frequency with an offset of zero hertz. I tried to move it closer. So let's say one hertz. So removing this phase offset, we're getting better. Whoops, I went too far. And at the end, you can compensate for the phase offset due to the difference
between local oscillator of a sound card and DCF77. So 15 PPM is not too bad. It's much better than most DVB-T receivers. So now I can get my phase information. And again, same story. I have my one PPS, one pulse per second from GPS.
I have my amplitude signals. Somehow, I don't know why, but if I do have two sources, I get out of sync sometimes in the graphical display. I don't know why. And the green is the phase. So I'm going to show you I've been running this since October. So I don't want to tune manually the frequency offsets every time. So what I did here is a very ugly.
What I'm doing here is I'm recording using new radio. And then I process in Matlab. I have easier thought in Okta rather than working directly in C++. So I find it easier to record the data and then process them. So what I'm doing here is I'm recording the raw data. I'm having a local oscillator, 77.5 kilohertz. And I'm shifting by a discrete time my DCF77 signal.
And I low-pass filter, then I decimate. So basically here, I reproduce the translating filter. That's easy enough. And then what I do is I make a coarse frequency offset measurement by taking the magnitude of the Fourier transform. And I say, OK, my Fourier transform is going to tell me by how much.
It's a very rough estimate, but by how much my local oscillator is offset with respect to DCF77. And I create a local oscillator, which is exponential of J2 pi, the frequency offset times the discrete time. And I remove this offset. And then once I've removed this coarse frequency offset,
I make a linear fit of the phase, the resulting phase. And once again, I remove the average value of this linear fit of a phase, which is an accurate or precise frequency difference estimate. And if I do this, I have my offset-corrected phase information. If I cross-correlate the signal, I get my beautiful cross-correlation peaks
once every second. So I can do the same thing we did with a lock-in amplifier with a sound card. And all this will be referenced to the 1 PPS of GPS by recording on a stereo sound card on one channel. I record DCF77. On the second channel, I record GPS. And I can measure here the difference
between GPS and 1 PPS. At some point, you sometimes have a slight offset of the phase. Sometimes you lose some of these cross-correlation peaks. And your timing will go away. But if you run a median filter on this thing, you will find the median value. And it will give you a pretty good idea of what the time stamp is, so the time delay is.
And this is an example of a three-day measurement, where I have on top the lock-in amplifier, so the high-quality, high-grade lock-in amplifier. And on the bottom, this is the output of my sound card. So you see that the results are pretty consistent. You have here daytime, nighttime.
So indeed, I can measure some sort of fluctuation of the atmosphere during nighttime with respect to daytime. And what you cannot see here below the text is the timing resolution is about 9 microseconds. If you remember, I was claiming that the daytime to nighttime variation is about 100 microseconds. So we have 10 times better time resolution
by using this cross-correlation technique. So this is what I want to say here. This was done in September. This was my very first measurement in September. I had to leave for a three-week trip in end of September, October. And unfortunately for you, I will not be able to show you any better result than this,
because since October, it's been winter. In winter, ionosphere is not stable. And I have to wait for the spring to come back to have, again, these measurements. I was hoping that yesterday maybe the sun might be helpful. But no, not yet. We have to wait a few more weeks before the ionosphere shows, again, this high stability during daytime.
Now, one of the thing here is that not everyone has access to a 192 kilohertz sound card. So of course, we all see that we have these DVBT receivers. And these DVBT receivers, they also allow me to assemble at two mega samples per second. So with two mega samples per second, I can timestamp the GPS1 PPS with a much better accuracy,
because now I have 500 nanosecond resolution. However, I'm losing in quantification resolution, because I'm going from a 16-bit sound card to an 8-bit A to D converter. So will the 8-bit A to D converter from the RTL2832U be enough to get my air wave out of my ground wave?
And we've seen that the antenna is a high quality factor bandpass filter. So will I be able to get rid of a jammer here? So we did this. And actually, I'm putting this as a provocation for the awesome SDR people in the room, because I couldn't get the direct sampling function
to work with the new radio block. And so what I did is this very ugly hack. I'm sorry, but at least it works. Is I just put in the LibRTLSDR the fact that I'm using an E4000. E4000 is a tuner that is not using an intermediate frequency. So by setting your RTLSDR with such a setting,
you just tell it, record whatever sample you have, and don't put any frequency offset. And somehow this works. So if there is a better way to do it, I think this is a very powerful way of using this RTL2832U. And it would be nice to either make sure that they work, or at least find some way of getting people using them.
And so once you get this working, you have one I input, which is DCF77, the QR input, which is the GPS signal. And again, we do the same thing. I record for one minute, once every minute. I record for duration of one minute my sample GPS DCF77. I process my data. And I will save them by doing the same processing.
I will not get into the details. But it's exactly the same story, except that you need to do, well, by processing a 2 megahertz signal with a very narrow band, you are targeting a 700 kilohertz bandwidth. You realize why you need to do multiple filter with decimation. And again, we go back to this whole DSP thing. You can take all the courses on DSP that you want.
If you haven't realized that if you do directly 2 megahertz to 700 hertz, you have like 10,000 coefficients. And it's not going to run efficiently when you haven't understood what DSP means. And so by doing this with practical signals, you will realize why you have to cascade your filters with decimation steps in between. And so you do this again with your RTLSDR.
And at some point, you should have a chart displayed. But there's a lot of points on this next chart. So recording the data, you extract phase and magnitude. Again, I have my once every second amplitude information. My phase information, my GPS timestamp,
the phase which doesn't look very good, but at least cross-correlation peaks come once every second here because I have a bit of leftover phase drift. I lose my cross-correlation. And I get here my cross-correlation peaks. So again, I can do my timing experiments. And I will be able to get a precise time of time stamping
by using DVBT. So even some of you who do not have 192 kilohertz sampling audio sound cards, you can reproduce this with 10 euro worth of equipment with an RTLSDR dongle. And again, by doing this, we have a collection of timing. So this is a conclusion of this
has been running since November. This is 6th of November. This was done yesterday. So what you see here, here I changed my algorithm. So this is just a slight change in my algorithm. So what you see here is that phase. Of course, we are referring GPS to a clock which is used for defining UTC.
So even though GPS is not officially linked to UTC, to universal time, well, we have good reasons to think that GPS will be synchronized with UTC. So what you see here is over three months, you do not see any observable drift, which is well known from the metrology community. And every one of these dips here is another day. So day, night, day, night, each one of these like this.
And if I plot an island deviation for those of you who are more into the time and frequency, so this is basically the stability of my timing capability as a function of integration time. So this is going from one minute up to three months. You see that you are at the 10 to the minus 7 down to 10 to the minus 11, which
is about 1,000 times worse than what has been published by PTB. But I think I've spent much less than 1,000 times less than what they did in their experiment. And just to give you some hints, the red curve here is a tuning fork. This is basically the stability of your watch. I took a tuning fork from a watch,
and I measured the stability. So what this chart tells you is that a tuning fork is very stable in the short term, but on the long term, it's going to drift if only due to temperature variation when I wear or I do not wear my watch. So over time, the tuning fork is going to drift, whereas these two clocks that we're using here, the cesium-locked clock generated by DCF,
is going to always demonstrate better stability the longer you integrate. And at some point, you have this intersection course which tells you that below this, your quartz watch will be better than DCF77. And above some time, you should look at DCF77. This is basically what cesium atomic clocks are doing.
You have high stability from the quartz oscillator on the short term, and you lock your quartz oscillator on the atomic transition. So basically, what you can do with the sound card is get acquainted to what people are doing when they try to design atomic clocks. So having various time sources, and these time sources will give the best stability over various integration
durations. So basically, what we showed here is that we have the ability to receive DCF77 signal using active antenna. We have been able to extract amplitude and phase information using electron amplifier. We've been able to reproduce these experiments with SDR, either with an audio sound card or with a DVBT dongle.
And we have consistent results. I'm eagerly waiting for a spring to get, again, a beautiful, stable signal. I try to correlate fluctuations in the DCF77 signal with x-ray density measured by geosynchronous satellites from NOAA. And there is no obvious correlation. I believe that 77.5 kilohertz is a bit too high,
and we're probing layers of atmosphere that are too deep to detect what people are usually measuring with ionosuems, which are lower frequency devices, 10, 15 kilohertz, where you really see the effect of x-ray density variation on the ionic content of ionosphere. I do have a few questions, though,
if anyone wants to pursue this thing. We're supposed to be told that you have two waves propagating on DLF, a ground wave and an ionospheric wave. That was the initial insight that low-round sea allows you to distinguish these two waves. And if, as published by PTB, the ground wave is much stronger than the ionospheric wave, how come that my unique signal that I'm measuring
is varying from day to night? This is a question for me. And I wonder whether you have to use a waveguide model where you have varying boundary conditions or whether you have truly two waves propagating. At the moment, I did everything on a PC. Whether it's worth spending time putting all of this in a microcontroller is questionable, but it's kind of fun, and it's not too difficult to do.
So the next thing I would like to look at is basically, at the moment, I've just been looking at the position of a cross-correlation peak, but I have not investigated the shape of a cross-correlation peak. And what I would expect is that if an air wave is slightly shifted with respect to the ground wave, what you should see
is a broadening of a cross-correlation peak because you have two waves slightly time-shifted. You should remember that the phase bit is 1.5 millisecond apart. It's 700 hertz. These are 100 microseconds time delays, so you will not be able to see the two cross-correlation peak separate, but you should be able to see some sort of cross-correlation peak broadening
as a function of this air wave and ground wave separation. So this is still open to questions. So basically, what I wanted to show you is that you can get involved into these electrical phenomena and atmosphere involved in this topic while investigating software-defined radio. So I've put this, of course, the slides on my website
and on the FOSDEM website. There has been a very extensive literature on this topic, but I haven't seen anyone trying to tackle the issue of, can anyone reproduce these experiments? Usually these are very fancy experiments with rather fancy equipment. And I believe that this is a topic that
would be worth investigating by more people, especially because very low frequency communication has lost some of its appeal since everyone is looking at high bandwidth transmissions. And nevertheless, low frequency is a very amusing topic. So with this, I thank you for staying all the way to the end of the talk.
And if there is any question, I would be happy to answer. Thank you. If you have any questions, just on your way out, look around you, pick up all the trash. Glass goes in bins in the hallways, and trash goes here. Hopefully it's not overflowing.
Take it over. The design of your RDL, was it a commercial, or was it sales field in your dongle? The dongle is, actually, I heard 20 euro. I buy mine for 8 euros, actually. This is just one of these 8 euro dongles. And you just have to unsolder the two capacitance and put in a semi-connector instead.
Have you got some designs somewhere that you plan on sharing it somewhere? It's like the previous presenter. I believe that in academia, we are not careful enough to make something that other people can use. I would be very happy to put on the website all my Okta script.
It might be that you have to send me an email to tell you how it works, because they're definitely not plug and play. But I would be very happy to send them to you. Your academics need to get that out of you and just put the code out on GitHub. We'll do so. Yeah? I'm surprised that you can't receive DCF on a sound card,
because although you need to language carrier, there should be input filters that and even digital processing in your sound card doesn't allow send and send all into your code. I have a 192 kilohertz sound card.
I'm well below half the sampling frequency. Why wouldn't my sound card allow me to record this? Well, the fact is that it doesn't. Luckily enough, it doesn't. OK. Well, I tried it on the lab. It's fitted with Dell laptops.
I tried it on two Dell laptops. It works. And on Panasonic laptops, it works. Any other question? Yeah? So you cross-correlated your receive signals with the p-RNA, with the pseudorandom? Pseudorandom number.
So what about, I mean, and then you do a frequency correction to amount for the LO offset of that peak. So wouldn't you just use the autocorrelation of your receive signal, because, oh.
The first time I did this, before I was even sure that my pseudorandom number generator was correct, the first thing I did was I autocorrelated to see if I could get my peaks once every second. The question about using directly the autocorrelation, because then once I got my autocorrelation, I wanted the cross-correlation to work.
Like with GPS, you can only get the cross-correlation if you correct for a phase. And then I forgot returning back to the autocorrelation. I wonder whether what you can grab out of the autocorrelation could be. But at some point, I want to compare this autocorrelation with the GPS tag. And then I need to have my local copy
that I will phase align with my GPS. So at some point, by taking external, I think the autocorrelation will tell me whether I got the right signal and whether I get a peak once every second. But if I want to synchronize with an external source, which is the one PPS, then I need to run a local copy of my oscillator or my pseudorandom number generator and phase synchronize it as we do in GPS.
But the autocorrelation is beautiful. Actually, this is completely off topic. But I've been spending the last three weeks working on RDS because Bastian tells us that RDS is so easy. It took me three weeks to get it. And the first thing I got with this stupid RDS is indeed the autocorrelation. And once I got the autocorrelation,
I needed another three weeks to figure out how to get the bits out, how to get the correlation things. Autocorrelation is really something that should be the first thing tested whenever we want to check whether the signal integrity is correct. So autocorrelation is a very powerful tool indeed.
That's it. Thank you. Thank you. Thank you.