relat_shaf.m 612 B

123456789101112131415161718192021222324
  1. function PAC0=relat_shaf(PAC)
  2. % 将原始数据PAC值与shaffle数据PAC做相对值。PAC为结构数据
  3. [m,n]=size(PAC.pacmat);
  4. relat=PAC.pacmat-PAC.shf_data_mean;
  5. relatperc=(PAC.pacmat-PAC.shf_data_mean)./PAC.shf_data_mean;
  6. relat_mi=(PAC.pacmat-PAC.shf_data_mean)./PAC.shf_data_std;
  7. for i=1:n
  8. for j=1:m
  9. if relat(j,i)<0
  10. relat(j,i)=0;
  11. end
  12. if relatperc(j,i)<0
  13. relatperc(j,i)=0;
  14. end
  15. if relat_mi(j,i)<0
  16. relat_mi(j,i)=0;
  17. end
  18. end
  19. end
  20. PAC.pac_relat=relat;
  21. PAC.pac_relatperc=relatperc;
  22. PAC.relat_mi=relat_mi;
  23. PAC0=PAC;
  24. end