MATLAB编程处理肌电信号

clear all
clc
close all
x1=load（‘D:\JD\1.txt‘）;
x2=load（‘D:\JD\2.txt‘）;
x3=load（‘D:\JD\3.txt‘）;
x4=load（‘D:\JD\4.txt‘）;
ori=load（‘D:\JD\std.txt‘）;
ori=ori（:4）;
x1=x1（:4）;
x1=（x1-mean（ori））.*3.3/1024*1000;
x11=x1;
x2=x2（:4）;
x2=（x2-mean（ori））.*3.3/1024*1000;
x21=x2;
x3=x3（:4）;
x3=（x3-mean（ori））.*3.3/1024*1000;
x31=x3;
x4=x4（:4）;
x4=（x4-mean（ori））.*3.3/1024*1000;
x41=x4;
%=====%part3%===========%带通滤波器%========%
Wp = [30 200]/500; Ws = [10 250]/500;
Rp =0.1; Rs =30;
[nWn] = buttord（WpWsRpRs）;
[b1a1] = butter（nWn）;
%=====%part3%===========%带通滤波器%========%
%=====%part4%===========%50hz陷波器%=======%
Wp = [10 250]/500; Ws = [48 52]/500;
Rp =0.1; Rs = 30;
[nWn] = buttord（WpWsRpRs）;
[b2a2] = butter（nWn‘stop‘）;
%=====%part4%===========%50hz陷波器%=======