How do I fix this code to get a proper moving average plot?
显示 更早的评论
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 个评论
DINUMOL VARGHESE
2020-10-19
Can you attach your moving average file here?
采纳的回答
Shubham Gupta
2019-9-25
0 个投票
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 个评论
Shubham Gupta
2019-9-25
Further, as I have mentioned that there are other modes for movmean, one of such is backward moving average which only takes previous elements. This can be implemeted much easier.
Try replacing your function with :
function avg = running_avg(signal, avg_window)
avg = zeros(1, length(signal));
for i = 1:length(signal)
if i<avg_window
windw = 1:i;
else
windw = i-avg_window+1:i;
end
avg(1, i) = mean(signal(windw));
end
end
Now, when you compare your plot with
movmean(signal,[avg_window,0])
You should get similar results. I hope it helps !
Ayugma Acharya
2019-9-25
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
更多回答（0 个）
类别
在 帮助中心 和 File Exchange 中查找有关 Whos 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
