Power_spectral_density_analysis (1)_alonedk2_matlab脑电功率_EEGlab_E

frex = 1:0.5:80; pain = zeros(length(frex),63); window = hann(EEG.srate*2); [pain,frex] = pwelch(EEG.data([1:31,33:end],:)',window,[],1:0.5:80,EEG.srate); delta_power = mean(mean(pain(dsearchn(frex,1):dsearchn(frex,4),:),2)); theta_power = mean(mean(pain(dsearchn(frex,4):dsearchn(frex,8),:),2)); alpha_power = mean(mean(pain(dsearchn(frex,8):dsearchn(frex,13),:),2)); beta_power = mean(mean(pain(dsearchn(frex,13):dsearchn(frex,30),:),2)); gamma_power = mean(mean(pain(dsearchn(frex,30):dsearchn(frex,80),:),2)); xlswrite('Power spectral density.xlsx',[delta_power theta_power alpha_power beta_power gamma_power])