Create Moving Average filter WITHOUT filter() function

7 次查看（过去 30 天）
Terry Carney
Terry Carney 2021-5-6
评论： Terry Carney 2021-5-6
Hello,
I can't seem to plot a desired output using a Hanning Moving Average equation. My goal is to create a moving average filter with a 3-point moving average with a created signal x. However, I don't notice a difference between the unfiltered signal and the filtered signal. May someone give advice on how to enter a moving average equation and it's coefficients without the use of the filter() function? Or is the 3-point moving average so small that there isn't a difference?
Hanning Moving Average filter:
y[n] = 1/4(x[n] + 2x[n-1] + x[n-2]).
fs = 200; %sampling frequency
Ts = 1/200; %sampling time
t = 0:Ts:1;
x = sin(2*pi*2*t) + sin(2*pi*10*t) + sin(2*pi*90*t) %signal
plot(t,x); %unfiltered signal
title('Signal with 2, 10, and 90 Hz');
xlabel('time (s)');
%Difference equation of Hanning Moving Average filter:
n = 3
yn = (.25*x(n)) + (.5*x*(n-1)) + (.25*(n-2)) %Attempt at creating the moving average equation.
plot(t,yn)
Thank You
  5 个评论
显示 3更早的评论
Terry Carney
Terry Carney 2021-5-6
编辑：Terry Carney 2021-5-6
May you give me an example?
What I'm thinking of is:
y[1] = 1/4*(.25x + .25(-x))
y[2] = 1/4*(.5x + .5x)
y[3] = 1/4*(.75x + 1x +.25x)
Jonas
Jonas 2021-5-6
编辑：Jonas 2021-5-6
you could calculate the main part e.g. by
y=1/4*(x(3:end) + 2*x(2:end-1) + x(1:end-2))
but then you jave to think about the first two samples of y because the normal formula will have negative index. maybe you want to assume 0 before the actual x(n) starts. similar to that you have to think about two values beyond x because the formula for y contains up to n-2. this way y will be longer than your original x by two samples

请先登录，再进行评论。

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

采纳的回答

Mathieu NOE
Mathieu NOE 2021-5-6
hello
try this
fs = 200; %sampling frequency
Ts = 1/200; %sampling time
t = 0:Ts:1;
x = sin(2*pi*2*t) + sin(2*pi*10*t) + sin(2*pi*90*t) %signal
plot(t,x); %unfiltered signal
title('Signal with 2, 10, and 90 Hz');
xlabel('time (s)');
%Difference equation of Hanning Moving Average filter:
yn = zeros(size(t));
%%%%%%% main loop %%%%%%%%%%%%
yn(1) = 0.25*x(1); % first iteration
yn(2) = 0.25*x(2) + 0.5*x(1); % second iteration
ind = 3:length(t); % for n>= 3
yn(ind) = (0.25*x(ind)) + (0.5*x(ind-1)) + (0.25*x(ind-2));
%Attempt at creating the moving average equation.
plot(t,x,'b',t,yn,'r')
  2 个评论
Jonas
Jonas 2021-5-6
providing just the solution is boring :p providing help to find the solution would be the better way ;)
Terry Carney
Terry Carney 2021-5-6
@Jonas, I agree. But I was up for a couple nights trying to solve this. I do thank you for your guidance as well.

请先登录，再进行评论。

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Smoothing and Denoising 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by