clear all;
clc;
[NUM,TXT,RAW]=xlsread('4号降水现象仪原始数据.xlsx');
P=size(NUM);% 取矩阵行数
fp=fopen('20191004.txt','wt+');
for i=1:P(1)
Empty=zeros(22,20);%重新清空矩阵
c=num2str(NUM(i,2));
obsTime=strcat( c(9:10),':', c(11:12));
fprintf(fp,'%c',obsTime);
fprintf(fp,'\n');
%没有雨粒子情况,写入空矩阵
if(NUM(i,4)==0)
for m=1:1:22
for n=1:1:20
if n==20
fprintf(fp,'%5d\n',0);
else
fprintf(fp,'%5d',0);
end
end
end
else %有粒子情况
Count=NUM(i,4);%取多少种类型粒子
%针对每种类型粒子,搜索在22*20矩阵的位置,并将粒子数更新Empty空矩阵对应位置
for k=6:2:(5+2*Count)
temp2=RAW{i,k};
for u=1:22 %22个粒径,纵坐标
for v=u:22:(u+418)
temp=strcat('','ANU',num2str(v,'%04d'),'');
tf=strcmp(temp,temp2);
if tf==1
Empty(u,(round((v-u)/22))+1)=RAW{i,k+1};
end
end
end
end
%将更新完的Empty矩阵按照要求写入txt文件
for m=1:1:22
for n=1:1:20
if n==20
fprintf(fp,'%5d\n',Empty(m,n));
else
fprintf(fp,'%5d',Empty(m,n));
end
end
end
end
end
fprintf(fp,'%s',' CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC');
fclose(fp);