Questions tagged [signal-processing]
Digital signal processing of sound, video or other single or multidimensional data
254 questions
2
votes
1
answer
253
views
Mathematica implementation of Mason's gain formula
Mason's gain formula (MGF) is a method for finding the transfer function of a linear signal-flow graph (SFG).
I'm trying to port the MATLAB code to Mathematica code, but the Mathematica code outputs ...
- 1
3
votes
2
answers
225
views
Low-pass filter for InterpolationFunction
Problem outline
Let's suppose we have an interpolating function (an InterpolatingFunction[] object) that depends on a single argument, interpreted as time (obtained ...
- 474
1
vote
1
answer
100
views
How to estimate and replace values in a signal? [closed]
I have this data:
j={{0.,0.46154},{0.02,0.96675},{0.04,BLINK},{0.06,BLINK},{0.08,BLINK},{0.1,BLINK},{0.12,BLINK},{0.14,0.95638},{0.16,BLINK},{0.18,BLINK},{0.2,BLINK},{0.22,0.97462},{0.24,BLINK},{0.26,...
- 119
2
votes
1
answer
185
views
How can I get the power spectrum and mean power of my signal?
I have measured this signal over time (where x is time in seconds and y is AC in D):
...
- 119
5
votes
2
answers
413
views
Frequency identification via Fourier when wave period is larger than dataset length
The following code is taken from the manual
...
- 2,072
0
votes
0
answers
157
views
Comparing analytical and numerical results of 2D Fourier transform
I would like to be able to transform back and forth with Fourier and InverseFourier, however, I can't quite match up the results....
- 2,011
1
vote
0
answers
220
views
How to Model a Parametric Fast Fourier Transform in Mathematica?
I'm using this code which evaluates the FFT of my original signal (which is a time series). Using the Manipulate function of Mathematica it is possible to vary the ...
1
vote
3
answers
183
views
1
vote
0
answers
115
views
Find knee by intersecting line method
I have an x-y data list. The data is here. It looks like below:
I want to mark the knee location. I mean I need to find the point of greatest slope change and get the x coordinate of this point. ...
- 319
2
votes
1
answer
169
views
Implementing Izhikevich's spiking neuron model
I am trying to implement Izhikevich's simple spiking neuron model which is given by
V'(t) = 0.04 V^2+5V+140-u+I
u'(t) = a(bv-u)
if V=30 mV then V=c, u=u+d
...
- 61
1
vote
1
answer
452
views
Periodogram and frequency
On the x axis of periodogram, one should have either a period or frequency. What frequency is used in Periodogram --- angular or ...
- 4,066
0
votes
1
answer
117
views
`Periodogram` returning only 1/0 error
I have a large dataset and would like to get a Periodogram (I saved the data as CloudObject):
...
- 4,066
0
votes
1
answer
81
views
How to find the difference equation from system response or primitive equation?
If I have a system response or the primitive equation , how can i find its difference equation by mathematica?
eg:
...
2
votes
1
answer
308
views
Calculation of CIC/FIR filter is too slow
I'm investigating the propeties of Cascoded Integrator Comb (CIC) filter. (for sigma-delta modulator application). I have 2 quesitons and below is what I have done:
a). Any alternative method to ...
- 157
2
votes
2
answers
457
views
Signal processing from Discrete data (Discrete Fourier Transform)
I am new one to Mathematica. So excuse me if my question is not absolutely accurate. I am trying to find all isolated frequencies in Fourier domain. And then to make inverse Fourier. And finally to ...
- 23
6
votes
1
answer
500
views
Problem with Spectrogram
Some time ago I wrote a program to analyse a short segment of music and plot as a spectrogram in Mathematica. It allows one to find the notes in a weird chord for example. The example below is from ...
1
vote
0
answers
62
views
Phase unwrapping of transfer function model for calculating group delay [duplicate]
Suppose that I have a transfer function model
...
- 61
3
votes
1
answer
91
views
What do we get from OutputResponse[sys {u1,u2}]?
I understand what OutputResponse is doing in the first two examples below.
...
- 7,174
4
votes
1
answer
170
views
What is this this application of KalmanEstimator doing?
I struggle to understand the Kalman filter and the documentation for KalmanEstimator. Please help me understand the first Application on that documentation page. The code in question is copied here:
<...
- 7,174
2
votes
2
answers
363
views
Why does this PhaseRange not work?
Any idea why this PhaseRange does not work?
BodePlot[-(1/(s + 1000)), PlotLayout -> "Phase",
PhaseRange -> {0, 2 \[Pi]}]
- 2,612
0
votes
1
answer
99
views
How to predict the delay of a higher order digital filter
My last question on this was a minimal example. I start with the following:
freq=0.1;
input=Table[Sin[freq t],{t,1,250}];
I use results above in two digital ...
- 7,174
3
votes
2
answers
206
views
How to predict the delay of a digital filter?
I define a simple digital filter.
transferFunction = TransferFunctionModel[{{{0.13793103448275862*z}}, {{-0.9 + z}}}, z,SamplingPeriod -> 1]
Then I put samples ...
- 7,174
1
vote
1
answer
471
views
How to excite a transfer function with continuous white noise
I'm trying to get the output response of linear system represented by transfer function, once the white noise is given as an input. Here's the code:-
For generating continuous noise, I'm getting help ...
- 13
0
votes
0
answers
116
views
Why can't I use my discrete lowpass filter?
I am trying demonstrate a discrete IIR lowpass filter using the feature of RecurrenceFilter described here.
This is my filter.
...
- 7,174
4
votes
1
answer
279
views
How to conciliate results between Fourier and FourierTransform?
I can't match the analytical results of a continuous-time Fourier transformation with the results I get from a discrete Fourier transform using Fourier[] in ...
- 123
2
votes
1
answer
205
views
Play two sounds at the same time in mathematica. For example the inventions of Bach NO.8
I use mathematica to play the Inventions No8 by Bach. I want to play the two parts (by right and left hands )at the same time. But I don't know how to do it. Can you help me? Thank you!
...
- 47
2
votes
0
answers
170
views
Finding recurrence relation from an inexact numerical sequence
Consider the recurrence relationship
d[0] = 0;
d[i_] := (1 + Sqrt[1 + 4 d[i - 1]^2])/2 /; i >= 1
which generates the following sequence (generated by using <...
- 592
13
votes
1
answer
920
views
3
votes
1
answer
438
views
Audio Digital Signal Processing
Resource and strategy request
I want to identify resources and strategies to implement a Mathematica based solution for audio digital signal processing (DSP) and digital to analog conversion (DAC).
...
- 14.7k
1
vote
1
answer
214
views
IIR Digital Filter: samples not matching the frequency response
So I am a small project in Mathematica involving a Chebyshev 4th order highpass filter, and to obtain its digital IIR filter with the bilinear transformation.
Then I implement the differences equation ...
1
vote
1
answer
220
views
Microfluctuations analysis and Power spectrum
I measured a parameter over time and obtained the values below:
0,627896 0,205004 0,259237 1,059125 0,832184 0,587992 0,565537 0,527323 0,460228 0,471958 0,26696 0,75367 ...
- 119
1
vote
1
answer
131
views
Why do I get this error message from LowpassFilter?
I have a list of numbers data of length 10^7.
I have values
rate = 100000.;
cutoff = 2000.;
I try to low pass filter with
...
- 2,093
1
vote
3
answers
477
views
0
votes
1
answer
225
views
Where to get the "SignalSystem" mathematica package?
I am reading this book: Filter Design for Signal Processing Using MATLAB and Mathematica.It can find on the following link:
https://www.wolfram.com/books/profile.cgi?id=3841
It's mathmatica examples ...
- 11
0
votes
1
answer
152
views
Building complex signals in Mathematica
How to build a particular signal in Mathematica such that we have a sine curve for half cycle, then have zeros for some time and this continues?
It is simple to have the zeros for the second half ...
1
vote
0
answers
245
views
How to compute Total Harmonic Distortion (THD) for the following signal?
The signal is as follows
Sin[x - j]
with $j$ being
...
- 55
4
votes
2
answers
436
views
Sampling and filtering data
I'm trying to teach myself about sampling and filtering — and I'm not sure if I'm misunderstanding the concepts or misapplying Wolfram Language. I started by creating a simple test signal with ...
- 43
3
votes
1
answer
172
views
How to code a Sum Block when working with Transfer Functions?
I am working using transfer functions models with Mathematica and i am missing a basic feature like the ability to use a Sum Block.
How could the model above be modeled?
Considering:
...
- 189
5
votes
1
answer
545
views
Why doesn't HighpassFilter fully remove the DC component? [duplicate]
The high pass filter is not completely removing the DC component. See the picture. Blue is the signal before filtering, green is after.
Here is some background on the subject matter. I'm studying ...
- 467
5
votes
0
answers
198
views
How to train each layer in a Neural Network so they optimize different loss functions in an adversarial network?
Example NetGraph to illustrate the idea (Input is an online signal with a value between -1 and 1, Noise is Gaussian Centered at 0 and standard deviation 0.1, EvilNet is constrained to output a value ...
- 536
2
votes
0
answers
142
views
Sound steganography with `InverseSpectrogram`
Here is my code of how to hide a message inside a sound:
...
- 32.8k
0
votes
0
answers
121
views
How can I deconvolve a spectrum to get the information about each of the deconvolved spectra?
I have a spectrum with two peaks. I want to deconvolve it by gaussian so that I get the parametres (e.g. peak position, area, FWHM, Height) of each deconvolved spectrum.
I have tried the following:
<...
- 552
24
votes
2
answers
1k
views
Can we use Mathematica to design an electronic active filter?
How can we use Mathematica to design an electronic implementation of an active filter?
- 5,434
1
vote
0
answers
140
views
How to convolve two audio signal? [closed]
Can we obtain reverb effect by convolve two audio signal?
Here are some example.
https://www.bilibili.com/video/BV1vE411h7W2
- 11
0
votes
1
answer
257
views
Are there some programs about inverse Fourier and Laplace transfroms?
To be more exact, I have a function F[w_,s_], where $w$ is the Fourier transform of $x$ and $s$ is
the Laplace transform of $t$.
Now I want to perform the double inverse transforms $s\to t$ and $w\to ...
- 173
1
vote
0
answers
102
views
SpectrogramArray timescale
I have a 600 second @ 8000 Hz WAV file which need to be analyzed in detail in terms of events (red parts in plot) defined in time and frequency that need to be rather exact, in particular the ...
- 125
1
vote
0
answers
103
views
Is there some programs about the fast inverse descrete time Fourier transform
At first, I give the definition of the inverse discrete-time Fourier transform
$$\phi(s)=\frac{1}{2\pi}\int_{-\pi}^{\pi}\exp(iks)f(k)dk$$
Here what I use is
...
- 173
1
vote
1
answer
564
views
DSP filter problem in the simulation of a Lock-In amplifier
I wanted To explain to students how a Lock-In amplifier works. The best way is to simulate the behaviour of one. However, in my program, I get unexpected behaviour, which I suspect has to do with ...
- 81
1
vote
1
answer
218
views
Noise filtering in audio
I was given audio with white Gaussian noise and i have to implement two types of filter: Rectangle and Triangle that filter out redundant frequencies and reduce the amount of noise
Here is the code ...
- 41
0
votes
1
answer
293
views
obtain the phase between two signals [duplicate]
I would like to calculate the phase between two signals like
f= Table[Sin[x-0.8],{x,-100,100,0.1}]
g= Table[Sin[x],{x,-100,100,0.1}]
For doing that I use the <...
- 601