Skip to main content

Utilities for digital signal processing (DSP) fundamentals.

Project description

DSPFTW

Starting points and helper functions for learning digital signal processing.

Setup

If you haven't already, install Python 3.5 or greater, and the dspftw package.

python3 -m pip install dspftw --user

Intro

Decomposition

Superposition: The foundation of DSP

Sinusoids

Prefer to use the complex exponential form instead of real, as it's a more natural fit for fourier analysis and synthesis.

The following functions represent a complex sinusoid. They are equivalent, and take the same parameters:

  • A = Amplitude (y axis is amplitude).
  • f = Frequency (cycles per second) in Hertz. It is multiplied by 2π to get radians per second.
  • t = Time in seconds. These functions work in the time domain (x axis is time).
  • ϕ = Phase offset at t=0, in radians.
z(t) = A*exp(j*(2π*f*t+ϕ))

z(t) = A*cos(2π*f*t+ϕ)+j*A*sin(2π*f*t+ϕ)

We can define this in Python with the following.

import numpy as np

# Here we use the name "complex_sinusoid" instead of just "z".
def complex_sinusoid(A, f, t, phi): return A * np.exp(1j*(2*np.pi*f*t+phi))

In fact, this is defined in the dspftw package, so let's import that.

import dspftw

We can create our own sinusoid by defining everything except t.

def my_sinusoid(t): return dspftw.complex_sinusoid(A=5, f=5, t=t, phi=0)

We can get the signal at a bunch of times thanks to numpy arrays. We use numpy.linspace to generate the evenly spaced times.

import numpy as np
times = np.linspace(0, 1, num=25)  # 25 evenly spaced values between 0 and 1
my_signal = my_sinusoid(times)  # returns an array of complex values representing the signal

Now plot it out with dspftw.plot_complex(), which uses matplotlib.

import matplotlib.pyplot as plt
dspftw.plot_complex(my_signal)
plt.show()

Complex Exponential Signals

Roots of Unity

Roots of Unity Animation

Wolfram Mathworld

Delta Function

Conjugate

Complex Conjugate

Convolution

Convolution

Correlation

numpy.correlate

Kronecker Product

numpy.kron

Sample Signals

SigIDWiki sample signals

Loading Signals

Complex 8 bit

import numpy as np
signal = np.fromfile('filename', dtype='b')  # load the whole file
signal = np.fromfile('filename', dtype='b', count=1024)  # only load the first 1024 bytes of the file
signal = np.fromfile('filename', dtype='b', offset=1024)  # skip the first 1024 bytes of the file
signal = np.fromfile('filename', dtype='b', offset=1024, count=1024)  # skip 1024, then load 1024

This just loads the values as an array of real numbers, but we want it as complex. We have to interpret every other value as the imaginary component.

signal = signal[0::2] + signal[1::2]*1j

Complex 32 bit float, little-endian (x86)

signal = np.fromfile('filename', dtype='<f')

count and offset work as above, but note that offset is in bytes, so you must multiple by 4 since there are 4 bytes in 32 bits.

Complex 32 bit float, big-endian

signal = np.fromfile('filename', dtype='>f')

count and offset work as above, but note that offset is in bytes, so you must multiple by 4 since there are 4 bytes in 32 bits.

Project details


Download files

Download the file for your platform. If you're not sure which to choose, learn more about installing packages.

Source Distribution

dspftw-2021.132.1012.tar.gz (15.3 kB view details)

Uploaded Source

File details

Details for the file dspftw-2021.132.1012.tar.gz.

File metadata

  • Download URL: dspftw-2021.132.1012.tar.gz
  • Upload date:
  • Size: 15.3 kB
  • Tags: Source
  • Uploaded using Trusted Publishing? No
  • Uploaded via: twine/3.4.1 importlib_metadata/4.0.1 pkginfo/1.7.0 requests/2.25.1 requests-toolbelt/0.9.1 tqdm/4.60.0 CPython/3.8.10

File hashes

Hashes for dspftw-2021.132.1012.tar.gz
Algorithm Hash digest
SHA256 fc27b818521fe6931a019adf4c18d5dfabf34f316739852f5a9b0a6fa24b63b5
MD5 e21a8d78f4caa92370e980ca176d0fec
BLAKE2b-256 99c94a2d708e9733e6bab7dba80bb8bddad9569d4ca70b6f6a9e580830d72cd0

See more details on using hashes here.

Supported by

AWS Cloud computing and Security Sponsor Datadog Monitoring Fastly CDN Google Download Analytics Pingdom Monitoring Sentry Error logging StatusPage Status page