tan11.m

3D加速信号处理

FileName='running22.txt';
MK=importdata(FileName); 
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;

%选取归一化的零点并作归一化运算(选初始状态作为基准)
AverageT1=mean(MK1(1:100));
Tq1=find(abs(MK1-AverageT1)>=5);
Tq4=Tq1(10);
Average1=mean(MK1(1:(Tq4-50)));
Nor1=(max(MK1)-Average1)/0.9;
NorMK1=(MK1(:)-Average1)/Nor1;

%为消除使用filter函数制成的推移,对原向量前后添项
Temp1=NorMK1(1:200);
Temp1(201:(size(NorMK1,1)+200))=NorMK1(:);
Temp1(size(Temp1):(size(Temp1)+100))=NorMK1((size(NorMK1)-100):size(NorMK1));

window=ones(1,200)/200;
T=filter(window,1,Temp1);
T(1:100)=[];%删除增加项
T(size(T)-199:size(T))=[];
T1=NorMK1-T;

%bb=resample(MK1,1,4);
%MK1=bb;
%cc=resample(NorMK1,1,4);
%NorMK1=cc;
%输出对比图形
h0=figure('position',[50 50 900 650],...
    'name','原信号归一化以及去重力处理演示');
subplot(221);
plot(MK1);grid on;title('原信号');
axis([Xmin Xmax Y1min Y1max]);
subplot(222);
plot(NorMK1);grid on;title('初步归一化');
axis([Xmin Xmax -1 1]);
subplot(223);
plot(T);grid on;title('提取重力加速度');
%axis([Xmin Xmax -1 1]);
subplot(224);
plot(T1);grid on;title('去重力后的归一化信号');
%axis([Xmin Xmax -1 1]);
figure;
plot(T,'-');hold on;plot(NorMK1);hold off;grid on;
axis([1000 1500 -0.5 0.4]);