pre.m

来自「研究羊群效应的例程」· M 代码 · 共 35 行

temp=[];
[i,I]=xlsread('I');
I=I(3:end,1);
for m=[1:677]
[a,b]=xlsread(num2str(m));
b=b(3:end,1);
a=a(3:end,2:end);
for i=1:size(I,1)
    if i>size(b,1)
        b=[b;I(i)];      
        a=[a;a(end,:)];  %尾巴接上
    else
       x=I(i);
       y=b(i);
       if length(x{1,1})~=length(y{1,1})
           flag=1;
       else    
           z=x{1,1}-y{1,1};
           flag=length(find(z~=0));
       end
       if flag>0  %日期不匹配
           if i==1
              b=[I(1);b];
              a=[a(1,:);a]; %头补上
           else
              b=[b(1:i-1);I(i);b(i:end)];
              a=[a(1:i-1,:);a(i-1,:);a(i:end,:)];  %中间插入
           end
       end   
    end
end
temp=[temp a(:,4)];
end