data=EEG.data;
dataavg=mean（data（17::）3）;
for?i=1:size（data3）
datanoise（i:）=data（17:i）-dataavg;
datanoise2=sum（datanoise.^22）/size（data2）;
pdata（i:）=sum（data（17:i）.^2）/size（data2）;
end
pnoise=sum（datanoise2）/（size（data3）-1）;
pdataavg=mean（pdata）
psoure=pdataavg-pnoise/size（data3）;
Q=sqrt（psoure/pnoise）;
SNR=log10（Q）;

data=MLRd;
dataavg=mean（