Jump to content
Molecupedia is a parody encyclopedia written by Molecular Reality Corporation. It is not Wikipedia and is not affiliated with the Wikimedia Foundation. About

Sampling (signal processing)

Page protected
From Molecupedia, the encyclopedia of very small things
This is an old revision of this page. Molecupedia does not store old revisions, so the text below is the current version.
Measuring a continuous signal at a series of instants

Sampling is the measurement of a continuous signal at a series of instants, usually equally spaced. The number of measurements each second is the sampling rate, fs, and the time between them is 1/fs. A nanopore instrument samples the ionic current through a pore: at 200 kHz, the rate of the amplifier on the company's live molecular stream, a measurement is made every 5 µs, and a molecule that blocks the pore for 50 µs is seen in ten of them.[1] The question that sampling theory answers is when those ten numbers are enough to say what the current did in between.

The theorem

A signal that contains no frequencies above B is fixed completely by samples taken more than 2B times a second. Half the sampling rate, fs/2, is the Nyquist frequency. The result has several fathers. E. T. Whittaker studied the interpolation series in 1915; Harry Nyquist showed in 1928 that a telegraph channel of bandwidth B can carry no more than 2B independent pulses a second;[2] Vladimir Kotelnikov proved the theorem in 1933; and Claude Shannon stated it as a tool for communication engineers in 1949, in the form usually taught.[3] Shannon's statement includes the recipe for getting the signal back: place a curve of the form sin(πx)/(πx) at each sample, scale it by the sample, and add them up.

Aliasing

A frequency above the Nyquist frequency is not lost. It appears in the samples at a lower, false frequency, because a fast sine wave and a slow one can pass through the same row of points. A component at f shows up at |f − k fs|, with k the whole number that brings it into the band from zero to fs/2. At 200 kS/s a component at 150 kHz is recorded as one at 50 kHz, and nothing in the record tells the two apart. The Analog-to-digital converter article shows the case of 70 kHz sampled at 100 kS/s, which arrives at 30 kHz.

The cure is to remove what is too fast before the converter sees it, with an analog low-pass filter. Because a filter does not cut off sharply, the cutoff is set well below the Nyquist frequency. Papers on solid-state pores report filter cutoffs of 0.1 to 0.4 of the sampling rate, with a Bessel filter the usual choice because it does not distort the shape of a pulse.[4][5] A step in current, which is what an event is, contains frequencies at every height, so no sampling rate is high enough to record a pore's events without a filter. The filter decides what is kept, and the sampling rate has to be chosen to suit it.

Choosing a rate

Two panels showing the same filtered 50 microsecond dip in current, one sampled every 5 microseconds with ten points inside the dip and one every 10 microseconds with five
A blockade of 50 µs, filtered with a single pole at 70 kHz and sampled every 5 µs (top) and every 10 µs (bottom). Calculated; no noise is drawn.

Three things pull in different directions. Faster sampling resolves shorter events and lets a gentler filter pass them. It also produces more data: a 16-bit record at 100 kS/s is 200 kB a second, 720 MB an hour. And a wider filter admits more noise (see Signal-to-noise ratio). Commercial sequencers sample slowly, because a motor protein holds the strand: a 2018 chemistry advanced the strand at 450 bases a second and was sampled at 4 kHz, about nine samples to a base.[6] A solid-state pore passes DNA thousands of times faster and has no such luxury.

Sampling faster than the filter requires has one use: the extra samples can be averaged. Averaging blocks of 100 samples turns a 100 kS/s record into a 1 kS/s one with white noise lower by a factor of 10, which is how the web display of the MR1 Nanopore Reader is fed.[7] Oversampling in this way spreads the noise of the converter's rounding over a wider band, and each fourfold increase in rate buys about one bit of resolution after averaging.

At Molecular Reality

The company's instruments sample at 200 kHz on the live stream, at 100 kS/s in the MR1 Nanopore Reader and, in the first batch of amplifiers for the demonpore 64, at "up to 125,000 samples of current per second".[7] That last rate carries a signal only to 62.5 kHz. Elsewhere the figure 125 kHz is given as the bandwidth of the amplifier, which would call for more than 250,000 samples a second, a point worked through in Teabiscuits. The README of the MR1 repository labels the link between the amplifier board and its converter "SPI Bus (100 kHz)", the sample rate and not the clock of the bus, which for 16-bit samples must run at 1.6 MHz or more.[8]

The company's critique of its own plans for machine learning compares the rates of two instruments it proposes to combine: a camera at 30 frames a second sees the flow cell once every 33,000 µs, and the amplifier samples it every 5 µs, a ratio of about 6,700 to 1 (The Nanopore World Model).[7]

See also

References

  1. ^ "Keltar's Field Guide to Molecular Reality". Molecular Reality Corporation. August 2026.
  2. ^ Nyquist, H. (1928). "Certain topics in telegraph transmission theory". Transactions of the American Institute of Electrical Engineers 47 (2): 617-644. doi:10.1109/T-AIEE.1928.5055024.
  3. ^ Shannon, C. E. (1949). "Communication in the presence of noise". Proceedings of the IRE 37 (1): 10-21. doi:10.1109/JRPROC.1949.232969.
  4. ^ Kwok, H.; Briggs, K.; Tabard-Cossa, V. (2014). "Nanopore fabrication by controlled dielectric breakdown". PLOS ONE 9 (3): e92880. doi:10.1371/journal.pone.0092880.
  5. ^ Smeets, R. M. M.; Keyser, U. F.; Dekker, N. H.; Dekker, C. (2008). "Noise in solid-state nanopores". Proceedings of the National Academy of Sciences 105 (2): 417-421. doi:10.1073/pnas.0705349105.
  6. ^ Rang, F. J.; Kloosterman, W. P.; de Ridder, J. (2018). "From squiggle to basepair: computational approaches for improving nanopore sequencing read accuracy". Genome Biology 19: 90. doi:10.1186/s13059-018-1462-9.
  7. ^ a b c "Company records". Molecular Reality Corporation. October 2026.
  8. ^ "MR1 Nanopore Reader". GitHub. October 16, 2026.
Last edited 1 hour ago by Aperture tube