五點滑動平均法平滑處理
通過數(shù)據(jù)采集系統(tǒng)采樣得到的振動信號數(shù)據(jù)往往疊加有不規(guī)則的隨機干擾信號。由于隨機干擾信號的頻帶較寬，有時高頻成分所占比例還很大，是的采集到的數(shù)據(jù)繪成的振動曲線上呈現(xiàn)許多毛刺，為了消除干擾信號的影響，提高振動曲線的光滑度，常常需要對采樣數(shù)據(jù)進行平滑處理。利用最小二乘法原理對離散數(shù)據(jù)進行線性平滑的方法稱為直線滑動平均法。
MATLAB源程序如下:
- function y = mean5( x,n )
- %五點滑動平均分法平滑處理
- % 輸入
- % x 待平滑信號
- % n 平滑次數(shù)
- % 輸出
- % y 平滑后信號
- %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('原始信號')
- % subplot(312)
- % plot(t,y)
- % legend('加入噪聲信號')
- % z = mean5(y,10);
- % subplot(313)
- % plot(t,z)
- % legend('平滑處理后信號')
- L = length(x);%信號長度
- 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
所有資料51hei提供下載:
mean5.rar
(505 Bytes, 下載次數(shù): 42)
2018-12-21 15:30 上傳
點擊文件名下載附件
|
