TFchirp: Time Frequency Decomposition Toolbox for Chirp Signals
Project description
Time Frequency Transform for Chirp Signals
Step 1: Quadratic chirp signal
import numpy as np
import scipy
import matplotlib.pyplot as plt
# Generate a quadratic chirp signal
dt = 0.0001
rate = int(1/dt)
ts = np.linspace(0, 5, int(1/dt))
data = scipy.signal.chirp(ts, 10, 5, 300, method='quadratic')
Step 2: S Transform Spectrogram
import TFchirp
# Compute S Transform Spectrogram
spectrogram = TFchirp.sTransform(data, sample_rate=rate)
plt.imshow(abs(spectrogram), origin='lower', aspect='auto')
plt.title('Original Spectrogram')
plt.show()
Step 3: Quick Inverse S Transform
# Quick Inverse S Transform
inverse_ts = TFchirp.inverse_S(spectrogram)
plt.plot(inverse_ts-data)
plt.title('Time Series Reconstruction Error')
plt.show()
Recovered spectrogram:
# Compute S Transform Spectrogram on the recovered time series
inverseSpectrogram = TFchirp.sTransform(inverse_ts, sample_rate=rate)
plt.imshow(abs(inverseSpectrogram), origin='lower', aspect='auto')
plt.title('Recovered Specctrogram')
plt.show()
Project details
Release history Release notifications | RSS feed
Download files
Download the file for your platform. If you're not sure which to choose, learn more about installing packages.
Source Distribution
TFchirp-0.0.1.tar.gz
(3.3 kB
view hashes)
Built Distribution
TFchirp-0.0.1-py3-none-any.whl
(14.9 kB
view hashes)