48 views (last 30 days)

My current code is this:

{for mdo=1:12

index_month = find(month==mdo);

month_mean = mean(avhrr.sst(index_month));

seasonal_cycle(mdo)=month_mean;

s_cycle_month(index_month) = month_mean;

end

hold on

plot(t, month_mean, 'r')}

When I run this, it produces a blank graph. I think the loop is correctly producing the seasonal cycle, how can I get it to plot the graph?

Cedric Wannaz
on 10 Jan 2018

Edited: Cedric Wannaz
on 10 Jan 2018

Cedric Wannaz
on 11 Jan 2018

Check every variable. Maybe one is an empty array. Try to compute intermediary terms manually, e.g.

index_month = find(month==1)

If it outputs nothing it means that month, which you suppose is a vector with some elements equal to 1, has in fact no element equal to 1. If the output looks fine, then evaluate:

avhrr.sst(index_month)

If the output is not a vector of floating point elements, then maybe avhrr.sst is not what you think it is.

Check every expression for a given month index at first, so you are sure that the inputs are what you think they are, and that the outputs are valid. Once this works, embed it in the loop and try to debug the loop.

Also, try plotting seasonal_cycle first with e.g.

plot( seasonal_cycle ) ;

If you tried plotting s_cycle_month against t, I guess that it means that t is not 1:12. This is a little more complex to debug, because s_cycle_month is more complex than seasonal_cycle (you build the former by re-injecting the mean at positions that correspond to the current month apparently).

Sign in to comment.

Chad Greene
on 6 Mar 2019

Sign in to comment.

Sign in to answer this question.

Opportunities for recent engineering grads.

Apply Today
## 1 Comment

## Direct link to this comment

https://nl.mathworks.com/matlabcentral/answers/376376-writing-a-loop-to-calculate-a-seasonal-cycle-and-then-plot-the-seasonal-cycle#comment_659186

⋮## Direct link to this comment

https://nl.mathworks.com/matlabcentral/answers/376376-writing-a-loop-to-calculate-a-seasonal-cycle-and-then-plot-the-seasonal-cycle#comment_659186

Sign in to comment.