I'm trying to do simple averaging in Simulink and need some help with it.

Question

Guillermo 2011-11-9

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/20700-i-m-trying-to-do-simple-averaging-in-simulink-and-need-some-help-with-it

回答： Martin Murtagh 2021-3-23

Assume AI'm acquiring a noisy signal at 1 mS interval and I want to get 1 point every 50 mS. I would like to collect the data for the 50 mS and spit out an average for logging. I'm not sure how to do it.

Any simple solutions would be appreciated.

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Alexander Bottema 2011-11-9

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/20700-i-m-trying-to-do-simple-averaging-in-simulink-and-need-some-help-with-it#answer_27352

在 MATLAB Online 中打开

function y = my_average(u) %# codegen
persistent i buf bufSum;
if isempty(buf)
    buf = zeros(1,50);
    bufSum = 0;
    i = 1;
end
bufSum = bufSum - buf(i);
buf(i) = u;
bufSum = bufSum + u;
i = i + 1;
if i > numel(buf)
    i = 1;
end
y = bufSum / numel(buf);

2 个评论
显示无隐藏无

Guillermo 2011-11-9

Thank you, I will try to make this a sub-system and I hope it works in simulink...

Guillermo 2011-11-9

It seems to have worked. It generates a warning at every iteration but I can live with that.

Thank you,

请先登录，再进行评论。

Answer 2

Fangjun Jiang 2011-11-9

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/20700-i-m-trying-to-do-simple-averaging-in-simulink-and-need-some-help-with-it#answer_27328

Use the "Weighted Moving Average" block from Simulink>Discrete, click help for details.

2 个评论
显示无隐藏无

Guillermo 2011-11-9

I don't have a Weighted Moving Average" block. I tried to create one with the following:

****************************************************

function Vector_out = fcn(Raw_In, Vector_in)

%#codegen

Vector_in = circshift(Vector_in, [1,1]);

Vector_in(1) = Raw_In;

Vector_out(:) = Vector_in(:);

******************************************************

Vector_in is initialized by a data store block

and then Vector_out would be summed with the Sum of elements block then divided by the dimension. In this way I would have a moving average which is more desirable but a block average would do too.

If this sparks ideas let me know. It didn't work for me because I could not use the data store block to define the vector_in or Out I don't recall which.

Fangjun Jiang 2011-11-9

Hah, it is obsolete! http://www.mathworks.com/help/toolbox/simulink/ref_obsolete_blocks/weightedmovingaverageobsolete.html

Use the new Discrete FIR Filter, http://www.mathworks.com/help/toolbox/simulink/slref/discretefirfilter.html

请先登录，再进行评论。

Answer 3

Walter Roberson 2011-11-9

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/20700-i-m-trying-to-do-simple-averaging-in-simulink-and-need-some-help-with-it#answer_27329

Feed the signal in to an integration block that is reset to 0 every 50 samples; take the 50th output and divide it by 50 to get the average for that period.

There might be a cleaner way to only trigger the output collection every 50 outputs and tie that in to the integration block being cleared right afterwards.

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

Guillermo 2011-11-9

I would like to try this, but I'm not yet very familiar with simulink. How do it set up the integration block?

The initial condition is set to 0, how about the saturation limits? and the Ignore limit and reset when linearizing check box?

Thanks again,

请先登录，再进行评论。

Answer 4

Guillermo 2011-11-9

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/20700-i-m-trying-to-do-simple-averaging-in-simulink-and-need-some-help-with-it#answer_27440

在 MATLAB Online 中打开

I got my code to work too, thanks to the sample by Alexander. The two codes give the same result and I'm not sure either which is more efficient. I tried passing the number of elements n to acquire so that I can pass a number of points to average but got an error...

function y = fcn(u)
%#codegen
persistent i n buf;
if isempty (buf)
    i = 0;
    n = 50;
    buf = NaN(50,1);
end
buf = circshift(buf, [1,1]);
buf(1) = u;
n = numel(~isnan(buf));
y = sum(buf)/n;

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

Answer 5

Martin Murtagh 2021-3-23

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/20700-i-m-trying-to-do-simple-averaging-in-simulink-and-need-some-help-with-it#answer_655757

I have replaced the moving average block in my models with the subsystem shown. By setting the delay length to the moving average value divided by the time-step and the gain value to one over the moving average value, this subsystem will calcuate the moving average in the same way as the moving average function with a sliding window.