power_data = eeg_power(ravgcsdm) extracts power values from avgcsdm matrix
0001 function power_data = eeg_power(ravgcsdm) 0002 %power_data = eeg_power(ravgcsdm) 0003 % 0004 %extracts power values from avgcsdm matrix 0005 % 0006 ch_pair_indices; 0007 if size(ravgcsdm,2) == 8385 0008 NChan = 129; 0009 else 0010 error('unknown number of channels') 0011 end; 0012 power_data = zeros(size(ravgcsdm,1),NChan); 0013 for ichan = 1:NChan 0014 power_data(:,ichan) = ravgcsdm(:,chpair(ichan,ichan)); 0015 end; 0016 0017 testreal = isreal(sum(power_data)); 0018 0019 if testreal ~= 1 0020 error('complex numbers in power array') 0021 end; 0022 0023