Questions tagged [signal-processing]

AKA digital signal processing (DSP). A signal is an information-carrying changing attribute of an entity, but in the digital sense, 'signal' refers to either received or transmitted streams/blocks of data, commonly representing real-world quantities such as audio levels, luminosity, pressure etc over time or distance. 'Processing' is the act of altering, analyzing or characterizing the data to retrieve/modify information inherent in the signal in question.

Filter by
Sorted by
Tagged with
0 votes
0 answers
23 views

How to transform an array using its first spectral components

I am trying to transform an array based on its first spectral components, but I get an output feature which is so different from the original one, even if I keep the whole components. Let's consider ...
Kyv's user avatar
  • 667
0 votes
0 answers
12 views

CMSIS-DSP: Getting multiple definition of `sinTable_f32'

I am trying to compile part of the CMSIS-DSP and want to use the function arm_sin_cos_f32 in file ControllerFunctions/arm_sin_cos_f32.c Additionally the folder CommonTables is inclued. All others are ...
user19471767's user avatar
0 votes
0 answers
11 views

How to achieve a gradual change in frequencies for a sweep (chirp) signal on STM32F4?

I am generating a sweep signal (a sine wave whose frequency increases over time) on STM32F4 Discovery board by using a look-up table and changing the Prescaler value of TIM2. Unfortunately, the ...
Bro Bro's user avatar
  • 11
-2 votes
0 answers
53 views

Basic noise reduction using python [closed]

The following is a program that should take a some environment's ambient noise and create a profile that can be used to essentially reduce the amount of background noise in a recording. There is a ...
Noah Lee's user avatar
-1 votes
0 answers
31 views

How to a plot 2D spectrogram

I have a numpy array that is NxM.. where N is the time and M are sensors.. for example, sensor M1 records N1_1........................N1_36000 for 1 hour. sensor M100 records N100_1.......................
Reinhardt Rading's user avatar
0 votes
1 answer
31 views

Implementation differences between feedforward and backward comb filters

Following https://en.wikipedia.org/wiki/Comb_filter , I implemented the following feedforward comb filter in Python. def combfilter(x, delay=2, gain=-3): memory = [] y = np.zeros_like(x) ...
lussco's user avatar
  • 3
-2 votes
1 answer
11 views

FFmpeg - ducking music for voiceover (sidechain)

I am trying to duck the background music from a voiceover input file using FFmpeg (I tried doing this using pydub with no success). The below command is predominantly working, except that it is ...
kaushal's user avatar
  • 49
0 votes
0 answers
24 views

I 'am working on a problem of modulation and demodulation using CPFSK. but the Demodulated output stream is not correct even without adding any noise?

this is the matlab code I'm working on to get the demodulated signal without adding any noise or considering the channel % Adaptive CPFSK Modulation and Demodulation % Parameters bitrate = 1000; % ...
Fatima Haider's user avatar
0 votes
0 answers
32 views

Normalised cumulative energy difference

I need to recreate these two equations in Matlab. (https://i.stack.imgur.com/01vSV.png) The goal is to use them to detect peaks for a spike sorting algorithm using using a normalised cumulative energy ...
NICHOLAS LAMB's user avatar
0 votes
1 answer
135 views

Calculating time between high state

In python I want calculations of high / low states of a signal in chunks of a given size (samples_to_process). The two required calculations are number index between rising edges (...
Al Grant's user avatar
  • 1,967
1 vote
1 answer
85 views

Saturate 16-bit signed integer to 12-bits signed

I'm working with an SDR that has a 12-bit signed ADC/DAC that stores in 16-bit IQ samples. I want to ensure that after all the DSP is done the samples saturate at 12 bits instead of getting truncated ...
Joseph Freivald's user avatar
0 votes
1 answer
21 views

What is wrong with my FM synth implementation?

I'm trying to get an FM synth to work, the synth uses wavetable oscillators, I am adapting the formula for FM synthesis: To use wavetable oscillators as carrier and modulator, so I am doing something ...
TiagoLr's user avatar
  • 2,814
1 vote
1 answer
23 views

Librosa: concerning the spectrogram is the sound pressure level a SPL value or a weighted dB(A)?

I am new to Python and am trying to analyze an audio with librosa. Concerning the spectrogram is the sound pressure level a SPL value or a weighted dB(A)? And does anyone know where I can read about ...
Waffelmike's user avatar
0 votes
0 answers
43 views

Power Spectral Density and Total Power Calculation for Signal

I am currently working on analyzing a signal with the autocorrelation function given by: R(τ) = 100 * cos( 10000 * τ*pi) .* (Λ(2000τ))^2 I need assistance in plotting the power spectral density Gf​(ω) ...
Javad Ibrahimli's user avatar
0 votes
1 answer
24 views

Pole Zero plot in Scilab

While trying to generate the pole-zero plot for a discrete time all-pass system tha has a pole at 3/4 and a zero at 4/3, scilab code that I used didnt work. I used he code z=poly(0,'z') n = [z-4/3]; ...
lakshminarayanan raghavendran's user avatar
0 votes
1 answer
24 views

Is there any solution how to play a sound without declaring the duration of it?

I have a code which basically generates a wave(sine, triangle,...) but only for certain duration. Also I have midi implementation and when I press a button the sound starts to generate(and play). Is ...
Tom hala's user avatar
0 votes
0 answers
35 views

Is there any MATLAB filter or method to fix distortion in the color gradation of images?

I am working on an algorithm to reduce the size of images, but there is a problem, which is that the resulting image contains distortions in color gradations, as you can see below. Is there a filter ...
Arkan_Abd's user avatar
-1 votes
1 answer
82 views

How can I convert 16 bytes of data into 2 pieces of 8 bytes of data?

I am using the TMS320F28x card. I am sending a struct worth 16 bytes from the STM32 card. I want to receive this data on the TMS320F28x card. After receiving the data, I want to divide it into 2 parts ...
SerdarSoylemez's user avatar
0 votes
0 answers
37 views

Search for sound tags in the audio signal

I have an audio recording that was recorded through a telephone channel, the beginning and end of which contain an audio tag. I have a "clean" version of this audio signal, its spectrogram ...
Gawain's user avatar
  • 133
1 vote
0 answers
22 views

Using Python to fourier transform Interofrogram to spatial frequency

so essentially i have a set of laser fringes projected onto a camera. I captured this image and converted it to greyscale. I perform a 2D fourier transform and i get the following image: This is what ...
user2279603's user avatar
0 votes
0 answers
22 views

problem with GSAFT (backprojection) algorithm for nearfield image reconstruction on my simple planar array

I have a planar array containing 16 transmitters and 88 receivers located on the same plane with dimension 30cm-by-30cm. I have an object of size 6cm-by-6cm. Object and antenna array plane are in ...
mohammad rezza's user avatar
0 votes
0 answers
54 views

GNU Radio and Wildlife DSP

I have a RTL-SDR and want to count the numbers of beeps per minute for a wildlife tracker around 150Mhz. There is no information encoded/modulated in the beeps - it all comes down to beeps per minute -...
Al Grant's user avatar
  • 1,967
-1 votes
0 answers
21 views

RangeDopplerResponse flipped in Python

I'm trying to use the Matlab engine to create a RangeDopplerResponse graph as follows: rngdopresp = eng.phased.RangeDopplerResponse('DopplerFFTLength',2048, 'DopplerOutput', 'Speed', 'RangeMethod', '...
stressedutstudent's user avatar
0 votes
0 answers
28 views

How to run call back functions while playing audio in python?

Just a little background to what I am trying to achieve. I am using Audioowl python library for beat detection and music analyses. After analyzing a music file, the audioowl analyse_samples() function ...
Vibhore Jain's user avatar
2 votes
2 answers
115 views

How to implement multi-band-pass filter with scipy.signal.butter

Based on the band-pass filter here, I am trying to make a multi-band filter using the code bellow. However, the filtered signal is close to zero which affects the result when the spectrum is plotted. ...
Thoth's user avatar
  • 1,041
0 votes
1 answer
29 views

How to increase a value in a logarithmic range by 100th proportionally?

I'm working on a logarithmic slider/knob and I am not able to increase its value by a proportional step. Its easy in linear terms, given a range [a,b] to increase a value by a 100th step: increment = (...
TiagoLr's user avatar
  • 2,814
1 vote
0 answers
9 views

How to implement Bernaola Galvan segmentation algorithm with matlab to complete time sequence segmentation

enter image description here I want to use the bg algorithm to remove the part of the processing data in the figure that tends to 0 (that is, the stage when the tool and the workpiece are not in ...
Zhiqiang Liu's user avatar
0 votes
1 answer
34 views

Scipy sosfilt ValueError ndarray is not C-contiguous

sosfilt from the scipy=1.9.3 library is giving me a strange ValueError when inputting a numpy=1.23.4 array. MWE: import numpy as np from scipy.signal import sosfilt fs=48000 rng = np.random....
Mike's user avatar
  • 23
0 votes
1 answer
71 views

Plotting the PSD of EEG signals - python

I have this line of code to compute and plot the psd of the raw EEG signals 1st sample i tried, i took Sub-S01 from: https://openneuro.org/datasets/ds004367/versions/1.0.2 2nd sample i tired i took ...
NMA's user avatar
  • 1
-1 votes
1 answer
68 views

Fast fourier filtering treshold / forex tick data denoising

I have a big dataset of currency pairs tick data and i want to denoise it using fourier. The problem is that i want to automatically find a treshold for filtering, it will be a part of a larger system ...
Vinicento's user avatar
0 votes
0 answers
25 views

How to generate Full spectrum analysis in vibration using Python

Mathematical procedure for obtaining a full spectrum the spectrum obtained from the data from two probes I need some help. I aquired vibration signals using Adre 408 & Pulse 3560C. When using Adre ...
user20918008's user avatar
1 vote
0 answers
65 views

How to implement a symmetric autocorelation function?

I am trying to compute these values: where * is the complex conjugate,z[n]= exp( 1i* 2* pi* cumsum(wav_file_signal[n])/ wav_fs), t' = t + int(wav_fs/4), N is the signal length and t takes values in ...
Thoth's user avatar
  • 1,041
0 votes
0 answers
42 views

FFT Implementation issues (GPU Shader)

I am trying to implement an FFT in a shader to run on the GPU. I grabbed the algorithm from here: https://github.com/jbouny/fft-ocean ... nShader.js However, when implemented and compared with other ...
ndivuyo's user avatar
0 votes
0 answers
35 views

How to find time of beat (R-event) or RR interval from ECG signal using wfdb toolbox?

I am trying to find the time of beat from the physionet data https://physionet.org/content/prcp/1.0.0/. I have written the following code for one record 12726: To read data record = wfdb.rdrecord('...
bipvan's user avatar
  • 89
0 votes
0 answers
26 views

Savitzky-Golay Filter parameters for visualizing corresponding velocity profile

I have some Gaze position data(in pixels) that I have smoothened using a Savitzky-Golay filter with a first difference method. I intend to visualize the velocity of this eye data in order for further ...
I Paul Roy's user avatar
0 votes
1 answer
50 views

Persistent time shift after cross correlation

I have two data sets each containing data on the x and y positions of an object. The data is sampled with different sensors. One of the data sets is sampled at 125 Hz the other at 500Hz. I have ...
Simon Jeppesen's user avatar
0 votes
0 answers
58 views

Algorithm for calculating distance/similarity between sets(!) of colors

I would like to compare two images based on a color distribution. I've implemented an algorithm (based on kNN) that extracts 3 most dominant colors from an image, together with weights for each color. ...
mrpostman889's user avatar
0 votes
1 answer
34 views

Do these two DC filter algorithms achieve the same thing, and is one better?

I have two different DC filter algorithms and I'm wondering what their main differences are and if one is better than the other. Here's a DC filter algorithm from dsprelated.com: y = x - xm1 + 0.995 * ...
D Blues's user avatar
  • 33
0 votes
0 answers
9 views

Why is the phase negative in the APSK trigonometric notation, but positive in the constellation diagram?

Let's write a trigonometric record of QAM: If I=1 and Q=1, then the phase in the cosine will be equal to -pi/4, since the arctangent will be equal to pi/4. At the same time, this point is located in ...
Vladimir Korshunov's user avatar
0 votes
0 answers
15 views

Methodology to time series data pattern recognition

I did some researches on SO and found many related questions/answers about pattern recognition. I was hesitating whether to post this question because this could be duplicated... But since I'm a noob ...
HQW.ang's user avatar
  • 119
0 votes
0 answers
20 views

FFT Algorithm leaves discontinuities in input (audio processing)

I am trying to program an audio plugin which uses an FFT algorithm to convert audio input into frequency data, does some processing on that frequency data (I'll explain those in a second), and then ...
Kallyn's user avatar
  • 1
0 votes
1 answer
58 views

Numpy signal processing: Efficiently summarize subarrays of consecutive nonzero elements

so, I'm working with FFT data that when thoroughly cleaned should look like many many zeros and an occasional large number. What I have at the moment is many zeros and an occasional short subarray of ...
Alosapien's user avatar
1 vote
1 answer
42 views

Conjugate symmetric: 3D Fourier transform dimension

I have a real-valued input 3D array with the shape of (H,W,D)=[8,8,20], where H, W, and D represent height, width and depth(z dimension), respectively. When computing the DFT, what will be the ...
S.EB's user avatar
  • 2,094
0 votes
0 answers
34 views

What is equivalent of smooth(t,ECG_signal,0.1,'rloess') in python?

What is python equivalent of this code to remove baseline wandering. It is using the local polynomial regression smoother (LOESS). t = 1:length(ECG_signal); yy2 = smooth(t,ECG_signal,0.1,'rloess'); ...
Fateme Nazari's user avatar
0 votes
1 answer
8 views

Get Frequency Response of Decimated Filter with unit impluse function

I am working with CIC-decimator with sigma-delta implementation. Now I want to check it's frequency response, Let us say: I implementted them with python, I have no filter coefficients. So I simulate ...
马慧超's user avatar
  • 193
0 votes
0 answers
24 views

how to plot frequency domain of only a portion of wav and do variable smoothing?

premise: I need to plot the frequency domain of an audio sample, I know nothing about signal processing and very little of python I followed few examples and it is kinda working using fft and fftfreq. ...
davidear's user avatar
  • 129
0 votes
2 answers
112 views

How to Install the `thinkdsp` Package for Python

Below is a table summarizing the steps I have successfully completed so far: Description Windows 10 powershell command (1) run powershell PS C:\Users\WORKINGSTATION> (2) create new project ...
Marco Moldenhauer's user avatar
2 votes
1 answer
59 views

How to programmatically identify the first and second minima and the peak in a bell-shaped curve?

I'm working with a dataset that forms a bell-shaped (normal-like) distribution. I'm trying to find three specific points in this distribution: The first minimum point before the curve ascends. The ...
Joe's user avatar
  • 597
1 vote
0 answers
51 views

How to get surface topology from fringe pattern in python

I'm trying to recover the surface shape of a body from fringe pattern. There's a interferometer projecting circular fringes on the body. A camera takes a picture of the reflection of these fringes. I ...
cqRinO's user avatar
  • 23
0 votes
1 answer
55 views

Why does the list returned by librosa.feature.rms() have so few values?

I have audio data from a file that lasts a second : audiodata, rate = librosa.load(datafile, sr=None, mono=True) Having a sample rate of 44100 Hz, the resulting list has 44100 values. But when ...
Alice Mabille's user avatar

1
2 3 4 5
118