Hello,
I am calculating the wavelet transform coherence betweem two signals. I am trying to extract the average coherence between 0.01-1Hz. I have tried using the f output of the wcoherence() function however, I'm confused by this output becasure it is a matrix of complex numbers that are both positive and negative. In the plot, I can clearly see that there are values within this frequency fange however, the current approach I'm using returns NaN or the values increase linearly with the increasing frequency band which does not seem correct.
How do I extract the average coherence within this frequency range?
Here is the code I am currently using:
[wtc, f, coi] = wcoherence(x, y,fs);
for k = 0.01:0.05:1-fbandStep
upperFrequencyBound = k+fbandStep;
freqIndices = find(f_trial >= lower & f_trial <= upper);
avgCoherence(k) = mean(wtc_trial(freqIndices));
