首页 > 其他分享 >EMG信号的低通巴特沃斯滤波器

EMG信号的低通巴特沃斯滤波器

时间:2022-10-10 16:34:55浏览次数:49  
标签:滤波器 EMG 沃斯 随机数 巴特 低通 信号


使用matlab中自带的randn函数产生一组随机数,作为EMG信号,然后EMG信号的采样率为2048hz。这里随机数产生的随机数种子采用的机遇系统时钟的随机数种子。系统输入有两个,一个是仿真时间,单位为s,一个是干扰值,输出有两个,一个为EMG信号,一格式带60hz正弦干扰的EMG信号。系统的设计步骤如下所示:

首先使用randn产生一组随机数,然后设计一个低通的巴特沃斯滤波器,其截止频率为150hz,将建立的随机数EMG信号输入到巴特沃斯滤波器。再将滤波得到的信号添加一个60hz的sin信号作为干扰。从而实现函数一个主要功能。

首先通过matlab自带的函数randn函数参数随机数,这里采用基于系统时钟的随机数种子,在matlab中使用的代码语句为:

RandStream.setDefaultStream(RandStream('mt19937ar','seed',sum(100*clock)));

这里使用了RandStream命令实现随机数种子的产生,其中sum(100*clock)为获取系统的时钟,这里如果需要每次产生固定的随机数,只需要将clock修改为一个固定的值就可以了。在完成这个语句之后,执行:

y1 = randn([N,1]);

从而完成EMG随机数,其仿真结果如下所示:

EMG信号的低通巴特沃斯滤波器_滤波器

图1 EMG随机信号

    为了设计需要,这里还需要对输入的信号进行FFT变换,对于FFT,主要使用MATLAB

的自带函数fft进行分析,matlab的代码如下所示:

fy1 = fft(y1,NFFT);
f = fs/2*linspace(0,1,NFFT/2+1);
subplot(322);
plot(f,abs(fy1(1:NFFT/2+1)));
title('Single-Sided Amplitude Spectrum of y(t)')
xlabel('Frequency (Hz)')
ylabel('|Y(f)|')
grid on;
axis([0,NFFT/2,1.2*min(abs(fy1)),1.2*max(abs(fy1))]);

其仿真如下所示:

EMG信号的低通巴特沃斯滤波器_巴特沃斯_02

图2 EMG随机信号的频谱图

从仿真结果如下所示,随机产生的EMG信号,其频谱图在频域上的各个频率点的都有分布。

[n,Wn]=buttord(Wp,Ws, Rp,As);可得到满足性能的模拟巴特沃斯滤波器的最小阶数n及截止频率ωc,其中ωp为通带的拐角频率,ωs为阻带的拐角频率,ωp和ωs的单位均为rad/s;Rs为通带区的最大波动系数,Rp为Rs阻带区的最小衰减系数,Rs和Rp的单位都为dB。

通过[b, a]=butter(n,Wn,'low');[b,a]=butter(n,ωc,′s′)可设计截止频率为的n阶低通模拟巴特沃斯滤波器。将得到的b和a作为滤波器的系数输入,其代码如下:

y2=filter(b,a,y1);

其仿真结果如下所示:

EMG信号的低通巴特沃斯滤波器_随机数_03

 

图3 通过巴特沃斯滤波之后的波形

其截止频率为150hz,通过FFT后,EMG信号的频谱图如下所示:

EMG信号的低通巴特沃斯滤波器_matlab_04

 

图4 通过巴特沃斯滤波之后的频谱

从上图可以看到,通过巴特沃斯滤波器之后,信号的频谱只存在于150hz的频域区域,大于150hz的频谱分量被滤除。

标签:滤波器,EMG,沃斯,随机数,巴特,低通,信号
From: https://blog.51cto.com/u_15815923/5744737

相关文章

  • 基于systemgenerator的根号计算
    一、系统设计仿真结果以及硬件资源估计(用于复制到你的那个txt文件中去即可。)顶层框图: 整个系统的结构如下所示: 进入内部模块则有: 其主要由三大部分组成: 第一部分:输入信......
  • 跟我学Python图像处理丨傅里叶变换之高通滤波和低通滤波
    摘要:本文讲解基于傅里叶变换的高通滤波和低通滤波。本文分享自华为云社区《[Python图像处理]二十三.傅里叶变换之高通滤波和低通滤波》,作者:eastmount。一.高通滤波傅......
  • 跟我学Python图像处理丨傅里叶变换之高通滤波和低通滤波
    摘要:本文讲解基于傅里叶变换的高通滤波和低通滤波。本文分享自华为云社区《​​[Python图像处理]二十三.傅里叶变换之高通滤波和低通滤波​​》,作者:eastmount。一.高通滤......
  • 低通滤波,高通滤波,中值滤波
    低通滤波和高通滤波低通滤波和高通滤波需用到傅里叶变换知识,可参考这里。图像在频域里面,频率低的地方说明它是比较平滑的,因为平滑的地方灰度值变化比较小,而频率高的地方......