draw5.asv

3D加速信号处理

%h0=figure('position',[50 50 900 650],...
%    'name','原信号归一化以及去重力处理演示');
%b3+b4
%偏移量补偿
Temp11=T1(1:5);
Temp11(6:(size(T1,1)+5))=T1(:);
Temp11(size(Temp11):(size(Temp11)+2))=T1((size(T1)-2):size(T1));%
Temp12=T2(1:5);
Temp12(6:(size(T2,1)+5))=T2(:);
Temp12(size(Temp12):(size(Temp12)+2))=T2((size(T2)-2):size(T2));%
Temp13=T3(1:5);
Temp13(6:(size(T3,1)+5))=T3(:);
Temp13(size(Temp13):(size(Temp13)+2))=T3((size(T3)-2):size(T3));%
T1=filter(b4,1,Temp11);
T1(1:7)=[];
NorMK1_DH_b3b4=T1;%
T2=filter(b4,1,Temp12);
T2(1:7)=[];
NorMK2_DH_b3b4=T2;%
T3=filter(b4,1,Temp13);
T3(1:7)=[];
NorMK3_DH_b3b4=T3;%

%b7+b8
%偏移量补偿
Temp11=T1(1:11);
Temp11(12:(size(T1,1)+11))=T1(:);
Temp11(size(Temp11):(size(Temp11)+5))=T1((size(T1)-5):size(T1));%
Temp12=T2(1:11);
Temp12(12:(size(T2,1)+11))=T2(:);
Temp12(size(Temp12):(size(Temp12)+5))=T2((size(T2)-5):size(T2));%
Temp13=T3(1:11);
Temp13(12:(size(T3,1)+11))=T3(:);
Temp13(size(Temp13):(size(Temp13)+5))=T3((size(T3)-5):size(T3));%
T1=filter(b8,1,Temp11);
T1(1:16)=[];
NorMK1_DH_b7b8=T1;%
T2=filter(b8,1,Temp12);
T2(1:16)=[];
NorMK2_DH_b7b8=T2;%
T3=filter(b8,1,Temp13);
T3(1:16)=[];
NorMK3_DH_b7b8=T3;%
subplot(411);
plot(NorMK2_DH);grid on;title('归一化去重力信号');
axis([500 1500 -1 1]);

subplot(412);
plot(NorMK2_DH_b5);grid on;title('N=5的均值滤波');
axis([500 1500 -1 1]);

subplot(413);
plot(NorMK2_DH_b6);grid on;title('N=5的中值滤波');
axis([500 1500 -1 1]);