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陷波器%=======