读取指定文件夹下小时数据并计算面平均日降水数据

clc

clear

data = xlsread('D:\TY\大论文\黄河数据\气象数据\华县权重.xls');%泰森多边形站点和权重

sign = data(:,1);

power = data(:,2);

len = size(sign,1);

FilePath = 'D:\TY\大论文\黄河数据\气象数据\小时数据';

DATE1 = datenum(1956,1,1):datenum(2022,12,31);%修改年份

DATE1 = DATE1';

DATE1 = datevec(DATE1);

date(:,1:3) = DATE1(:,1:3);

len1 = size(date,1);

for k=1:len

files = dir(fullfile([FilePath,'\',num2str(sign(k,1)) ,'.xlsx']));

x = xlsread(fullfile(FilePath,files.name));

sign1 = find(x(:,5)>=999);

x(sign1,5)=0;

first = find(x(:,1)==1956&x(:,2)==1&x(:,3)==1);

new_data = x(first:end,:);

for j = 1:len1

final(j,k) = sum(new_data(24*j-23:24*j,5))*power(k,1);

end

end

finaldata(:,4) = sum(final,2);

finaldata(:,1:3) = date(:,1:3);

xlswrite(['D:\TY\大论文\黄河数据\气象数据\降水\1956-2022面平均日降水','\','华县','.xlsx'], finaldata)