<div dir="ltr">Hi EEGLablist,<div><br></div><div>The spectopo function calculates power spectral density values for a given signal. However I would like to try get the same exact values of the spectopo function but using MATLAB's pwelch() function instead since this function should give the power spectral density as well. However I am getting completely different values and I am not sure if the way I am calculating it is correct:</div><div><br></div><div><div>X  = EEG.data(channel_no,:,epoch_no);</div><div>         length_X  = EEG.pnts;                    % signal length</div><div>         nfft = 2^nextpow2(length_X);             % Next power of 2 from length of x</div><div>         </div><div>         no_of_overlapped_samples = 0;</div><div>         window = 512;</div><div>         </div><div>         % get psd using pwelch()</div><div>         [spectra,freqs] = pwelch(X, window, no_of_overlapped_samples, nfft, EEG.srate);</div><div><br></div><div>         spectra = 10*log(spectra);</div><div><br></div><div>Above I am calculating it using the pwelch() function, with spectopo it is as below:</div><div><br></div><div>[spectra,freqs] = spectopo(EEG.data(channel_no,:,epoch_no), 0, EEG.srate);<br></div><div><br></div><div><br></div><div>Is there something I am missing or am doing wrong?</div><div><br></div><div>Thank you!</div><div><br></div>-- <br><div class="gmail_signature"><b><i>James Bonello</i></b></div>
</div></div>