function [begin_time_num,end_time_num]=get_nihe_num(q,z,year,month,time)
%,year,month為各期觀測的年、月數組，time為總的觀測次數
%,q為計算起始年、月,z為計算終止年、月，該函數的作用是獲得對應的起始和算終止年、月觀測期數
j=1;
for i=1:time
    if year(i)==q(1)
       temp(j)=i;
       j=j+1;
   end
end
for i=1:length(temp)
    if month(temp(i))==q(2)
        begin_time_num=temp(i);
    end
end
j=1;
for i=1:time
    if year(i)==z(1)
       temp1(j)=i;
       j=j+1;
   end
end
for i=1:length(temp1)
    if month(temp1(i))==z(2)
        end_time_num=temp1(i);
    end
end