mopen.asv

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

clear all;
[filename, pathname] = uigetfile( ...
       {'*.dat';'*.ch0';'*.ch1';'*.CH0';'CH1';'*.*'}, ...
        'Pick a file to open');                     %%%%%-----Get file 
fname=[pathname filename];
    
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));
        if (afx=='dat'|afx=='DAT')
            dat=fread(fid,'uint16');
            fclose(fid);
            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=='Ch1'|afx=='CH1'|afx=='ch1')
            cho=fread(fid,'int32');
            fclose(fid);
       elseif (afx=='Ch2'|afx=='CH2'|afx=='ch2')
            chb=fread(fid,'int32');
            fclose(fid2);            
    end
    
       %plot(cho,'b');subplot(2,2,1);
    
      %plot(chb,'r');subplot(2,2,2);
   lb=length(chb);
   lm=length(chb);
   chb=chb(1:lm);
   cho=cho(1:lm);
    Rxy=corrcoef(chb,cho);
  
    rxy=num2str(Rxy(1,2));


end