您现在的位置:首页 > >

信号处理除去加性噪音的一个方法

发布时间:



信号处理除去加性噪音的一个方法

原理:在信号采集时,经常有工业噪音的干扰,比如打桩机引起的震动等,这里噪音一般为高频分量,可以利用傅里叶变换,在巴特沃斯低通滤波器下进行滤波处理,并且需要设置相应的通带截止频率和阻带截止频率。

下图为对信号处理后及处理前的频域和时域分布
下面为matlab仿真代码

a = xlsread(‘C:UserslegionDocumentsTencent Files949800083FileRecvMobileFileincX_1.xlsx’,‘sheet2’);

y = a( 394:end,4);

fs =1000;

t = 1/fs;

l = 42484;

t = ( 0:l-1)
t;
Y = fft(y,l);
n = 0:l-1;
k = n
fs/l;

subplot(324);

plot(k,abs(Y)/5)

title(‘原频率信号’);

axis([0,5,0,20]);

wp = 2.1/fs;

ws = 6/fs;

[n,wn] = buttord(wp,ws,1,50);

[a,b] = butter(n,wn);

[h,f] = freqz(a,b,‘whole’,fs);

f = (0:length(f)-1)’*fs/length(f);

subplot(325);

plot(f(1:length(f)/2),abs(h(1:length(f)/2)));

title(‘滤波器’);

axis([0,10,0,10]);

sf = filter(a,b,y);

subplot(321);

plot(t,sf);

title(‘滤波后信号’);

subplot(322);

plot(t,y);

title(‘原信号’);

SF = fft(sf);

subplot(323);

plot((1:length(SF)/2)
fs/length(SF),2abs(SF(1:length(SF)/2))/length(SF));

title(‘滤波后频率信号’);

axis([0,1,0,1]);


热文推荐
猜你喜欢
友情链接: