tan8.asv

3D加速信号处理

FileName='running22.txt';
MK=importdata(FileName);      % 将生成的MKMATLAB.txt文件导入工作空间，变量名为MK，实际上它不显示出来
%Y=0:300/(size(MK,1)-1):300;

MK1=MK(:,1);MK2=MK(:,2);MK3=MK(:,3);

%坐标选取
MK1min=min(MK1);MK2min=min(MK2);MK3min=min(MK3);
MK1max=max(MK1);MK2max=max(MK2);MK3max=max(MK3);
Xmin=0;
Xmax=(round(size(MK,1)/500)+1)*500;
Y1min=round(MK1min/50)*50;
Y2min=round(MK2min/50)*50;
Y3min=round(MK3min/50)*50;
Y1max=(round(MK1max/50)+1)*50;
Y2max=(round(MK2max/50)+1)*50;
Y3max=(round(MK3max/50)+1)*50;

%滤波器的设计
b1=[0.2 0.2 0.2 0.2 0.2];  %均值
b2=[0.1 0.1 0.6 0.1 0.1];  %中值
b3=[0.073 0.118 0.618 0.118 0.073];  %中值(黄金分割)
b4=[-0.1 -0.1 1.4 -0.1 -0.1];
%b4 锐化算法
%b5 量化
b=b2;
%n=1000;
%window=blackman(n);
%[h,w]=freqz(window,1);
[B,A]=butter(4,0.05);
N=10;
MK4=filter(B,A,MK1);%MK4(1:5)=mean(MK1(1:5));
MK4=uint8(MK4);
MK5=filter(window,1,MK2);%MK5(1:5)=mean(MK2(1:5));
MK5=uint8(MK5);
MK6=filter(window,1,MK3);%MK6(1:5)=MK3(1:5);
MK6=uint8(MK6);

h0=figure('toolbar','none',...
    'position',[50 50 900 650],...
    'name','3D加速信号的预处理');
subplot(811);
plot(MK1,'r');grid on;
axis([Xmin Xmax Y1min Y1max]);
subplot(812);
plot(MK4,'r');grid on;
%axis([Xmin Xmax Y1min Y1max]);

subplot(813);
plot(MK2,'b');grid on;
axis([Xmin Xmax Y1min Y1max]);
subplot(814);
plot(MK5,'b');grid on;
%axis([Xmin Xmax Y2min Y2max]);

subplot(815);
plot(MK(:,3),'m');grid on;
axis([Xmin Xmax Y3min Y3max]);
subplot(816);
plot(MK6,'m');grid on;
%axis([Xmin Xmax Y3min Y3max]);