Мне нужно реализовать средний фильтр для набора данных, но у меня нет доступа к набору инструментов обработки сигналов. Есть ли способ сделать это без использования цикла for? Вот код, который у меня работает:
x=0:.1:10*pi;
noise=0.5*(rand(1,length(x))-0.5);
y=sin(x)+noise; %generate noisy signal
a=10; %specify moving window size
my=zeros(1,length(y)-a);
for n=a/2+1:length(y)-a/2
my(n-a/2)=mean(y(n-a/2:n+a/2)); %calculate mean for each window
end
mx=x(a/2+1:end-a/2); %truncate x array to match
plot(x,y)
hold on
plot(mx,my,'r')
РЕДАКТИРОВАТЬ:
После реализации решения merv встроенный метод фильтрации отстает от исходного сигнала. Есть ли способ обойти это?