How do I fix this code to get a proper moving average plot?

Question

Ayugma Acharya 2019-9-25

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/481972-how-do-i-fix-this-code-to-get-a-proper-moving-average-plot

评论： DINUMOL VARGHESE 2020-10-19

This is a homework question. I was given a set of wind speed data, and was asked to calculate moving averages as well as other information using MATLAB. Specifically, for the moving average, I am not allowed to use movmean, conv, filter, smooth or any other built in functions. So I made my own function, but its not giving the proper output. I know my plot using my function is wrong since it does not give me the same plot when I check with movmean.

This is the plot I get using movmean:

This is the plot I get using my own function:

I have attached both the data set and the .m file for the code. Any help would be very much appreciated. Thank you!

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

DINUMOL VARGHESE 2020-10-19

Can you attach your moving average file here?

请先登录，再进行评论。

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

Answer 1

Shubham Gupta 2019-9-25

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/481972-how-do-i-fix-this-code-to-get-a-proper-moving-average-plot#answer_393386

Your basic function 'running_avg' seeems to be working perfectly. However, If you want answers exactly as you would get using movmean then you should understand how movmean works.

movmean does center averaging ( there are other modes for averaging too but default is center ).

Center averaging is basically to find average by keeping the current element at the center and then take previous and future elements to make the array of size specified by window. Now,

When the window is even, average is centered about the current and previous elements.

When the window is odd, average is centered about the current element.

Note : When the window needs to be shorten due to array's indexing at the endpoints, the average is taken over only the elements that can fit in the window. For e.g. if we 3 previous elements but there are only 2 previous elements then, movmean will take only 2 previous elements and shortens the window to be averaged by 1.

So, you will have to change the code for your function 'running_avg' accordingly to get output 'exactly' as movmean.

4 个评论
显示 2更早的评论隐藏 2更早的评论

Shubham Gupta 2019-9-25

You mean at the start ? That's because you have set "avg" to 0 for "i" smaller than window size. So now, the moving mean will be zero for indeces smaller than window size. Instead you can shorten the window size and increase it with each increment of "i" till it becomes equal to window size.

Ayugma Acharya 2019-9-25