Extracting multiple envelopes from a signal
4 次查看(过去 30 天)
显示 更早的评论
I have a signal which when observed carefully shows a three individual sinusoidal signals (figures attached). I want these three individual signals to be extracted in terms of amplitude, frequency and phase. Please help!
2 个评论
John BG
2017-8-11
would you please be so kind to supply the signal?
in a .mat file
thanks in advance, awaiting answer
采纳的回答
Cedric
2017-10-17
Quick cheap trick gets 2/3. Was going to bed, no time for more but that could be a starting point.
function main
load( 'MultipleSignalEnvelope', 'ITDc0' ) ;
plot( ITDc0, 'b.' ) ;
grid on ;
hold on ;
t = 1 : numel( ITDc0 ) ;
A = 1.2 ;
P = fminsearch( @objective, zeros(3,1) ) ;
plot( A*sin((t+P)/360*2*pi).' ) ;
function mindif = objective( P )
mindif = sum( min( abs( ITDc0-A*sin((t+P)./360*2*pi )))) ;
end
end
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/178689/image.png)
%
9 个评论
Cedric
2017-10-20
编辑:Cedric
2017-10-20
You will have to try. Define P as a vector with 6 components instead of 3, the first three being amplitudes and the next three being "phases". In the objective function, replace A with P(1:3) and P with P(4:6). Also, update the initial conditions. Instead of zeros(3,1), try to start e.g.. with three unit amplitudes and three null phases: [ones(3,1);zeros(3,1)].
Note that you could use FMINCON and define constraints, e.g. upper and lower bounds for the amplitude and phases.
Also note that you'd better try to see if you could get help from someone with experience in digital signal processing, because this approach is just a quick trick. If you reword your question in a new thread and nobody answers, you may be able to hire someone for not that much on sites that market freelancers.
更多回答(1 个)
Image Analyst
2017-8-11
编辑:Image Analyst
2017-8-11
Look at the fft.
10 个评论
Image Analyst
2017-10-17
Well that's a change from what you said before. Too bad, that would have been so easy. I don't know how to solve it without spending way more time on it than I allow myself to spend on questions here, so . . . good luck.
另请参阅
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!