draw1.m

3D加速信号处理

figure;
subplot(411);
plot(NorMK2_DH,'b');grid on;
axis([500 1500 -1 1]);

subplot(412);
plot(NorMK2_DH_b2,'b');grid on;
axis([500 1500 -1 1]);

Temp22=NorMK2_DH(1:51);
Temp22(52:(size(NorMK2_DH,1)+51))=NorMK2_DH(:);
Temp22(size(Temp22):(size(Temp22)+25))=NorMK2_DH((size(NorMK2_DH)-25):size(NorMK2_DH));%
b10=0.01*ones(1,51);
b10(26)=0.5;
TT1=filter(b10,1,Temp22);
TT1(1:76)=[];

%bbb=NorMK2_DH_b5(:)-NorMK2_DH(:);
subplot(413);
plot(NorMK2_DH_b6,'b');grid on;
axis([500 1500 -1 1]);
subplot(414);
plot(TT1,'b');grid on;
axis([500 1500 -1 1]);