mreadfile.m

来自「人体脑信号波形检测与分析」· M 代码 · 共 51 行

function [cho,chb,rxy,lm]=mreadfile(fname)

fid=fopen(fname);
if fid<=0
    chb=0;
    cho=0;
    rxy='0';
    lm=0;
    ERRORDLG('File Open Error'); 
else
    afx=fname(length(fname)-2:length(fname)-1);
    fclose(fid);
    if afx=='da'
            dat=fread(fid,'uint16');
            l=length(dat);
            cho=dat(1:4:l);
            chb=dat(2:4:l);
            cho=cho-mean(cho);
            chb=chb-mean(chb);
            chb=myzeros(chb);
            cho=myzeros(cho);
            lb=length(chb);
            lo=length(cho);
            omax=max(cho);
            bmax=max(chb);
            lm=min(lo,lb);
            chb=chb(1:lm)./bmax;
            cho=cho(1:lm)./omax;
        elseif afx=='Ch'
            fname0=fname(1:length(fname)-3);
            fname1=[fname0,'Ch1'];
            fname2=[fname0,'Ch2'];
            fid1=fopen(fname1);
            cho=fread(fid1,'int32');
            fclose(fid1);
            fid2=fopen(fname2);
            chb=fread(fid2,'int32');
            fclose(fid2);            
    end
    
       
%    lb=length(chb);
 %   lm=length(chb);
   % chb=chb(1:lm);
  %  cho=cho(1:lm);
  %  Rxy=corrcoef(chb,cho);
 %   rxy=num2str(Rxy(1,2));
 %   lm=length(chb);
end