任务描述:我有一个站点观测降水的时间序列,需要判断降水观测是否连续出现10天以上的缺省值
写成了函数,直接调用即可
%% 本函数用于判断某数连续出现的次数超过某次,
function result = default_test(A,B,count)
&变量说明,A是时间序列,B是目标值,由于站点的缺省值默认为32700,因此B=32700,count是设定次数,我在使用时取了10
result = 0;
k= 0; %k是计数器
for i = 1:length(A)
if A(i) == B %等于目标数则开始计数
k = k+1;
if k>=count
result = 1; %说明连续出现次数超过设定次数
end
else
k = 0; %不是目标数则计数器归零且跳过
continue
end
end
以default_test(A,B,10)为例,如果返回结果是0,说明没有连续十天以上缺省;如果返回值是1,说明存在连续十天以上缺省
Hope this would bring you some inspiration!