五點(diǎn)均值法濾波濾波方法(平滑處理)MATLAB程序

ID:451569 · 發(fā)表于 2018-12-21 15:30

五點(diǎn)滑動(dòng)平均法平滑處理
通過數(shù)據(jù)采集系統(tǒng)采樣得到的振動(dòng)信號(hào)數(shù)據(jù)往往疊加有不規(guī)則的隨機(jī)干擾信號(hào)。由于隨機(jī)干擾信號(hào)的頻帶較寬，有時(shí)高頻成分所占比例還很大，是的采集到的數(shù)據(jù)繪成的振動(dòng)曲線上呈現(xiàn)許多毛刺，為了消除干擾信號(hào)的影響，提高振動(dòng)曲線的光滑度，常常需要對(duì)采樣數(shù)據(jù)進(jìn)行平滑處理。利用最小二乘法原理對(duì)離散數(shù)據(jù)進(jìn)行線性平滑的方法稱為直線滑動(dòng)平均法。

MATLAB源程序如下:

function y = mean5( x,n )
%五點(diǎn)滑動(dòng)平均分法平滑處理
% 輸入
% x 待平滑信號(hào)
% n 平滑次數(shù)
% 輸出
% y 平滑后信號(hào)
%exanple:
% N = 1000;
% t=0:0.1:(N-1)*0.1;
% x = 5*sin(t);
% a = rand(1,1000)*0.5;
% y = x+a;
% subplot(311)
% plot(t,x);
% legend('原始信號(hào)')
% subplot(312)
% plot(t,y)
% legend('加入噪聲信號(hào)')
% z = mean5(y,10);
% subplot(313)
% plot(t,z)
% legend('平滑處理后信號(hào)')
L = length(x);%信號(hào)長度
a = x;
for k = 1:n
b(1) = (3*a(1)+2*a(2)+a(3)-a(4))/5;
b(2) = (4*a(1)+3*a(2)+2*a(3)+a(4))/10;
for j = 3:L-2
b(j) = mean(a(j-2:j+2));
end
b(L-1) = (a(L-3)+2*a(L-2)+3*a(L-1)+4*a(L))/10;
b(L) = (-a(L-3)+a(L-2)+2*a(L-1)+3*a(L))/5;
a = b;
end
y = a';
end

復(fù)制代碼

所有資料51hei提供下載:

mean5.rar (505 Bytes, 下載次數(shù): 42)

ID:393578 · 發(fā)表于 2019-6-26 08:31

挺好的

帳號(hào)		自動(dòng)登錄	找回密碼
密碼			立即注冊(cè)