unlock/analysis/andres/Code/getUnlockSpecgram.m File Reference

Functions
Calculating the spectrogram if	isfield (SpecInfo,'nEpochs')&&SpecInfo.doEpochs%For epochs
Calculating the spectrogram if if isfield(unlockInfo.epoch,'nOdds')%eegData(	eegData (:,:,:, 2)
	save (saveFile,'dataSpec','SpecInfo','-v7.3')
Variables
	function [specData, SpecInfo]
Getting the size of	spectrogram [firstSpec,~,~] = mtspecgramc(squeeze(eegData(:,:,1,1)),SpecInfo.movingWin,SpecInfo.params)
Initial vble	eegDataSpec = nan(size(firstSpec,1),size(firstSpec,2),nChs,2)
Getting the spectrogram tic for	iTrial
end toc Super structure specData	data = mtspecgramc(squeeze(eegData(:,:,iTrial,1)),SpecInfo.movingWin,SpecInfo.params) mtspecgramc(squeeze(eegData(:,:,iTrial,2)),SpecInfo.movingWin,SpecInfo.params) eegDataSpec
specData	tspec = tspec
specData	fspec = fspec
SpecInfo	specData = {'T','Freq.','chs','epoch','stimFreq'}
SpecInfo	size = size(specData.data)
SpecInfo	FreqRes = specData.fspec(3) - specData.fspec(2)
SpecInfo	hz2spec = length(specData.fspec)/SpecInfo.freqResolut
SpecInfo	sec2spec = length(specData.tspec)/(length(eegData)/SpecInfo.freqResolut)
Saving matrix and info structure if exist(fullfile(BCIinfo.dirs.saveFilename(1:end-10),'Spectrograms'))%else%mkdir(BCIinfo.dirs.saveFilename(1 SpecInfo	saveFile = saveFile

---

Function Documentation

Calculating the spectrogram if if isfield (unlockInfo.epoch,'nOdds') % eegData( eegData	(	:	,
		:	,
		:	,
		2
	)

Calculating the spectrogram if isfield	(	SpecInfo	,
		'nEpochs'
	)

save	(	saveFile	,
		'dataSpec'	,
		'SpecInfo'	,
		'-v7.3'
	)

---

Variable Documentation

decoded given by the HSD Taken only the channels with real info unlockData data = mtspecgramc(squeeze(eegData(:,:,iTrial,1)),SpecInfo.movingWin,SpecInfo.params) mtspecgramc(squeeze(eegData(:,:,iTrial,2)),SpecInfo.movingWin,SpecInfo.params) eegDataSpec

Initial vble eegDataSpec = nan(size(firstSpec,1),size(firstSpec,2),nChs,2)

SpecInfo FreqRes = specData.fspec(3) - specData.fspec(2)

SpecInfo fspec = fspec

function[specData, SpecInfo]

Initial value:

 getUnlockSpecgram(eegData,SpecInfo)
%
% Calculates the spectrogram using chronux's multitaper functions 
%
% INPUT
% eegData:      matrix. This is the EEG raw data. Has to be in the form [samples x channels]
%
% OUTPUT
% 
% Andres. v1.0
% Created 15 Dec 2013
% Last modified 24 Dec 2013

nChs = SpecInfo.nChs

SpecInfo hz2spec = length(specData.fspec)/SpecInfo.freqResolut

Getting the spectrogram tic for iTrial

Initial value:

 1:SpecInfo.nOdd       % number of epochs
        fprintf('Processing epoch %i\n',iTrial)

Saving matrix and info structure if exist (fullfile(BCIinfo.dirs.saveFilename(1:end-10),'Spectrograms')) % else % mkdir(BCIinfo.dirs.saveFilename(1 SpecInfo saveFile = saveFile

SpecInfo sec2spec = length(specData.tspec)/(length(eegData)/SpecInfo.freqResolut)

SpecInfo size = size(specData.data)

SpecInfo specData = {'T','Freq.','chs','epoch','stimFreq'}

Getting the size of spectrogram[firstSpec,~,~] = mtspecgramc(squeeze(eegData(:,:,1,1)),SpecInfo.movingWin,SpecInfo.params)

end Spec analysis transformation factors SpecInfo tspec = tspec