How to fill in a matrix and compute its binary frequency?

Question

Hydro 2014-11-25

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/164278-how-to-fill-in-a-matrix-and-compute-its-binary-frequency

编辑： dpb 2014-11-26

Hello All,

I have thirty year of preciptation data in a vector form. I need to find the frequency of the dry and wet days and plot the frequency against days (1:365).

I am stuck as the following code place in the precipitation data however, i am not able to change it to a binary condition (Wet=1, dry=0). Can somebody help me out. I still have to count all the wet and dry days and make the frequency graph against the days (1:365).

Thanks

R=zeros(365,30); % Space for the data
P_data=P(:,4); % Preciptation data
n=length(P_data);
C=1;
for i=1:30;
    for j=1:365;
    R(j,i)=P_data(C);
    C=C+1;
    end
end
if R(j,i)>0; 
    R(j,i)=1; % imposing the wet conditon if the the value in the ith row and jth colum is larger then zero
else
    R(j,i)=0;
end

3 个评论
显示 1更早的评论隐藏 1更早的评论

Image Analyst 2014-11-25

Yes, but you still didn't attach P, or give code to read it in.

Hydro 2014-11-25

Data.mat

the data

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

dpb 2014-11-25

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/164278-how-to-fill-in-a-matrix-and-compute-its-binary-frequency#answer_160218

编辑：dpb 2014-11-26

在 MATLAB Online 中打开

Leap years not accounted for?

If not, as the above would indicate, then

R=reshape(P_Data,365,[])>0;

2 个评论
显示无隐藏无

Hydro 2014-11-25

yes, the leap year is ignored for this analysis. any thought about the code would be appreciated. I still don't see what i wanted to see (Probability of precipitation occurrence against duration in days (1:365).

dpb 2014-11-25

在 MATLAB Online 中打开

Not quite sure what you mean, precisely by Probability of precipitation occurrence against duration in days (1:365) but R above will be a logical array of T for rainy days arranged by year in columns. So, P(R|day) =

Prainbyday=sum(R,2)/size(R,2);

Clearly P(noRain) is the complement. Is that not what you're looking for?

请先登录，再进行评论。