How to align multiple signals?

21 次查看(过去 30 天)
Hi Matlab experts,
I'm trying to align these signals but have yet to be successful.
The signals are attached. They are periodic, noisy sin signals. I used the alignsignals function as well as xcorr to align the signals. Could somebody help me out with that?
Moreover, how can I compute the frequency of these signals?
Thanks in advance!

采纳的回答

Sulaymon Eshkabilov
Here is one of the viable solutions to compute resonant frequencies of each signal and find out if their are coherent:
S1 = load('s1.mat').s1;
S2 = load('s2.mat').s2;
S3 = load('s3.mat').s3;
S4 = load('s4.mat').s4;
S5 = load('s5.mat').s5;
S6 = load('s6.mat').s6;
fs = 1024;
S = [S1; S2; S3; S4; S5; S6];
N = length(S1);
freq = fs*(0:(N/2))/N; % Freq values
for ii=1:6
Y = fft(S(ii,:), N); % FFT of the signal s(t)
Yp = abs(Y/N); % Absolute value or magnitude of the signal spectrum
semilogy(freq,Yp(1:N/2+1), 'linewidth', 2)
hold all
end
title('FFT analysis of a signal')
xlabel('Frequency, [Hz]')
ylabel('Magnitude of signal = |S(f)|')
ylim([0, 0.025])
legend({'s1', 's2', 's3', 's4', 's5', 's6'})
grid on
You can check signal coherences using wcoherence() function:
wcoherence(S1,S2)
wcoherence(S2,S3)
wcoherence(S3,S4)
wcoherence(S4,S5)
wcoherence(S5,S6)
It looks like signal s5 and s6 are pure noise.
  2 个评论
Susan
Susan 2023-2-8
编辑:Susan 2023-2-8
@Sulaymon Eshkabilov Thank you so much for your response. I genuinely appreciate it.
Please tell me how I interpret the first figure's results regarding the relationship between the time domain signals.
I went through the wcoherence documentation but am still trying to figure out how to interpret the results. Any help on that would be greatly appreciated.
Finally, how should I do that if I want to align these signals in the time domain?
Sulaymon Eshkabilov
编辑:Sulaymon Eshkabilov 2023-2-8
Most welcome. What do mean to align the signals in the time domain? Put them together or what? to make their peak values correspond?
You can try alignsignals() fcn:
[S2a, S3a] = alignsignals(S2, S3, [], "truncate");
nexttile
plot(S2), hold on
plot(S3)
nexttile
plot(S2a), hold on
plot(S3a)
hold off
See - DOC.

请先登录,再进行评论。

更多回答(1 个)

Sulaymon Eshkabilov
One fundamental question for signal processing is a "must", that is "what is the data sampling time or sampling frequency?"
Pl., provide this info.
  1 个评论
Susan
Susan 2023-2-7
编辑:Susan 2023-2-7
@Sulaymon Eshkabilov Thanks for your response. The sampling rate of all signals equals 1024.

请先登录,再进行评论。

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by