提取VMF1中的Tm

clc;clear;
name=['y2018','.vmf1_g'];
fid = fopen(name);
a=[];a1=[];A1=[];
if fid>0
    i=0;
    while feof(fid) == 0
        i = i+1;
        tline{i} = fgetl(fid);
        a=tline(i);
        if i>22
            a1=cell2mat(a);
            if (a1(1)=='C'&&a1(2)=='H'&&a1(3)=='A'&&a1(4)=='N')||(a1(1)=='c'&&a1(2)=='h'&&a1(3)=='a'&&a1(4)=='n')
            A1=cat(1,A1,a1);
            end
            name1=['chan','.txt'];
            dlmwrite(name1, A1, 'delimiter','','precision',6)
        end
    end
end
fclose(fid);