mopen.m

人体脑信号波形检测与分析

clear all;
[filename, pathname] = uigetfile( ...
       {'*.dat';'*.ch0';'*.ch1';'*.*.ch2';'*.*'}, ...
        '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;
            
            %lb=length(chb);
   %lm=length(chb);
   %chb=chb(1:lm);
   %cho=cho(1:lm);
   % Rxy=corrcoef(chb,cho);
  
   %rxy=num2str(Rxy(1,2));
        elseif (afx=='Ch1'|afx=='CH1'|afx=='ch1')
            cho=fread(fid,'int32');
            fclose(fid);
            plot(cho);
       elseif (afx=='Ch2'|afx=='CH2'|afx=='ch2')
            chb=fread(fid,'int32');
            fclose(fid);
            plot(chb);
    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));




    sf=1;
    lm=length(chb);
    i=0:1:lm-1;
    
    set(gcf,'Name','Show Singal:');
    plot(i/sf,cho,'b');
    hold on;
    plot(i/sf,chb,'r');
    hold off;
    ylabel('show Singal');
    xlabel('time (s)');
    title('show Singal');
    grid on;
    set(gca,'NextPlot','replacechildren');
    
end