clear % 清空工作區(qū)變量

clc % 清空命令窗口

filepath = 'D:\Crop_ET'; % 文件路徑

filelist = dir(fullfile(filepath, '*.txt')); % 獲取路徑下所有txt文件

for i = 1:length(filelist)

? ? % 獲取文件完整路徑名

? ? filename = fullfile(filelist(i).folder, filelist(i).name);

? ??

? ? % 讀取數(shù)據(jù)到矩陣E，并將-9999替換為NaN

? ? E = load(filename);

? ? E(E==-9999) = NaN;

? ??

? ? % 計(jì)算E的平均值

? ? Result_E = nanmean(nanmean(E));

? ??

? ? % 顯示當(dāng)前文件的平均值

? ? disp([filelist(i).name, ': ', num2str(Result_E)]);

end