雨滴谱处理_ISXD30雨滴谱_雨滴谱_

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);