Question

我使用下面的代码根据时间步数每五分钟生成一次文件名。但它没有正常工作。如果你打开precipFileNames，你会看到中途，代码每5分钟停止一次，而是5分1秒生成一个这样的文件名：

E:\MRMS\2004\PRECIPRATE.20040402.051959.tif

我该如何正确地做到这一点？

timeSteps = 417;
pathstr = 'E:\MRMS\2004';
startTime = 7.320395312500000e+05;
peakTime2 = 7.320400104166666e+05;
precipFileNames=cell(timeSteps,1);

for l = 1:timeSteps
       %precipFileNames{m} = strcat(fileparts(refFile), filesep, datestr(startTime, 'yyyy'), filesep,'PRECIPRATE.',datestr(startTime, 'yyyymmdd.hhMMss'), '.tif');
       precipFileNames{l} = strcat(pathstr(1:end-4), datestr(startTime, 'yyyy'), filesep, 'PRECIPRATE.',datestr(peakTime2, 'yyyymmdd.hhMMss'), '.tif');
       peakTime2 = addtodate(peakTime2, -5, 'minute');  %No. of times we go back in time from peak time
 end

Answer 1

使用浮点数在内部存储日期/时间。每次循环时，您都会将一个非常小的值（5分钟，0.0035）添加到一个相对较大的值（7e05 - ish），这会导致floating point arithmetic errors的累积。这些错误表现为与您的预期值略有不同。

因为您在循环中反复执行此添加（到peakTime2），所以一次迭代的浮点错误会被放大，因为下一次迭代取决于前一次迭代的结果。

而不是不断更新peakTime2我会改变 delta 值，并在每次循环时将其应用于原始日期时间对象。这样就不会产生错误累积，只需执行一次减法即可获得特定值。

for k = 1:timeSteps
    % Compute how many minutes to shift this iteration
    shift = -5 * (k - 1);

    % Apply the shift to the reference time
    thistime = addtodate(peakTime2, shift, 'minute');

    % Create the filename
    precipFileNames{k} = strcat(pathstr(1:end-4), ...
                                datestr(startTime, 'yyyy'), ...
                                filesep, ...
                                'PRECIPRATE.', ...
                                datestr(thistime, 'yyyymmdd.hhMMss'), ...
                                '.tif');
end

作为旁注，为了让人们阅读您的代码，我强烈建议不要将l用作变量，因为它看起来很像1。

如何正确生成带有5分钟时间戳的文件名？

1 个答案: