data=load('中圈120.wtf');
t=data(:,1);
x=data(:,2);
y=data(:,3);
z=data(:,4);
% figure(1);%画原始图形
% plot(t,x);
% xlabel('t/s');
% ylabel('幅值');
subplot(2,2,1);plot(t,x);
xlabel('t/s');ylabel('x');
subplot(2,2,2);plot(t,y);
xlabel('t/s');ylabel('y');
fs=5000;%采样率
n=500; %滤波器阶数
Wn=0.0001;%数字截止频率
b=fir1(n,Wn);%滤波系数
X=filter(b,1,x);%用Matlab自带的filter函数进行滤波
% figure(2);%画滤波后的波形
% plot(t,X);
% xlabel('t/s');
% ylabel('幅值');
subplot(2,2,3);plot(t,X);
xlabel('t/s');ylabel('X');
Y=filter(b,1,y);
subplot(2,2,4);plot(t,Y);
xlabel('t/s');ylabel('Y');
评论0