Greetings,
I understand that you want to quantify the phase shift between the original and the reconstructed signals. The phase shift occurs as some information may be lost during the deconstruction and reconstruction process.
Please refer the code below to quantify the phase shift:
Fs = 1000; % Sampling frequency (Hz)
t = 0:1/Fs:1; % Time vector (1 second duration)
f = 10; % Frequency of the signal (Hz)
signal = sin(2*pi*f*t);
[c,l] = wavedec(signal,14,'db2');
lfnoise = wrcoef('a',c,l,'db2',10);
final = signal - lfnoise;
% Plot the original signal and the final filtered signal
figure;
subplot(2, 1, 1);
plot(t, signal);
title('Original Signal');
xlabel('Time (s)');
ylabel('Amplitude');
subplot(2, 1, 2);
plot(t, final);
title('Final Filtered Signal');
xlabel('Time (s)');
ylabel('Amplitude');
% Quantifying the Phase Shift
original_phase = angle(hilbert(signal));
filtered_phase = angle(hilbert(final));
phase_difference = filtered_phase - original_phase;
size(phase_difference)
Please find more information on the functions used in the code above using the following documentation links.
- hilbert - https://in.mathworks.com/help/signal/ref/hilbert.html
- angle - https://in.mathworks.com/help/matlab/ref/angle.html
Hope this helps.
