阵列光栅振动信号均值降噪
阵列光栅振动信号均值降噪
相关说明:
1、光栅点25个,采集的信号为1000HZ, 时间 1分钟;
2、每个光栅点一分钟的采样60000次(可计算出),每相邻的10个采样点取均值,作为一个采样点,因此采样次数可以为6000次;代码如下:
clear all;
clc;
a=fopen(‘C:\Users\02\Desktop\振动数据处理\0110\1-1-25-20200110153300.bin’,‘rb’); %一分钟的信号数据
data1=fread(a,‘float32’); % 读去二进制文件
fclose(a);
N=25; %光栅点的个数
data2=reshape(data1,N,length(data1)/N); %分为25个光栅点每一分钟采集的点数 每一个点采集60000次
data2=data2’; %转置 (25个光栅的点原始信号)
p=data2(:,1); %取第一个光栅的信号
p_1=p’;
for i=1:6000
p2(i)=sum(p_1(10i-9:10i));
end %每十个采样点进行求和
p2=p2’;
p2=p2/10; % 取均值
subplot(211)
plot(p2(:,1),‘r’);
subplot(212);
plot(data2(:,1),‘b’); %画图对比原始采集的信号,与降噪的信号
截图如下:matlab代码运行的情况: