i want to draw this DOP figure with these specific parametres
4 views (last 30 days)
Show older comments
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/1485212/image.png)
0 Comments
Accepted Answer
Star Strider
on 17 Sep 2023
This probably represents a multi-exponential decay, because a single exponential doesn’t faithfully reproduce it.
I would need the signal itself rather than an image of it, to estimate all its parameters.
Fs = 1E+9;
L = 2E-6;
t = linspace(0, L*Fs, L*Fs+1)/Fs;
k = log(0.6)/0.2E-6;
s = -4*(exp(k*t) .* sin(2*pi*2.5/5E-7*t));
figure
plot(t, s.*heaviside(t-0.5E-6))
grid
ylim([-1 1.2])
.
2 Comments
Star Strider
on 17 Sep 2023
O.K.
In that instance, you need to calculate the fft of the original signal to understand where the noise is with respect to the signal peak (that should be centred about
).
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/1485257/image.png)
A relatively easy way to do this is to either calculate the fft itself, or use the pspectrum function to estimate it.
Then use the bandpass function to filter it, eliminating the low-frequency baseline variation and the high-frequency noise. For bandpass and its friends, use the 'ImpulseResponse','iir' name-value pair to design an efficient elliptic filter to filter it. It may be necessary to threshold it, so initially eliminating the (relatively) constant section from 0 to 0.6 µs before filtering it. Another way to approach it would be to first use a lowpass filter to eliminate the high-frequency noise, and then a highpass filter to eliminate the baseline variation.
That is how I would approach it, at least initially, changing as necessary to get the most optimal result.
You can then experiment with the approach in Curve fitting to a sinusoidal function to estimate the parameters. It would be necessary to change the ‘fit’ function in that code to match the
function, so you get the appropriate information from it.
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/1485262/image.png)
More Answers (0)
See Also
Categories
Find more on Digital Filter Analysis in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!