相模變換

%V change before
Va=V.signals(1,1).values(1500:3000);
Vb=V.signals(1,2).values(1500:3000);
Vc=V.signals(1,3).values(1500:3000);
figure(1)
plot(V.time(1500:3000),Va,'r-')
hold on
plot(V.time(1500:3000),Vb,'b:')
hold on
plot(V.time(1500:3000),Vc,'g--')
hold off
legend('Va','Vb','Vc')
figure(2)
subplot(311)
plot(V.time(1500:3000),Va)
subplot(312)
plot(V.time(1500:3000),Vb)
subplot(313)
plot(V.time(1500:3000),Vc)
%---------------------------------------------
% % % % Va1=Va+16000;
% % % % Vb1=Vb+16000;
% % % % Vc1=Vc+16000;
% % % % Va2=Va1/32000;
% % % % Vb2=Vb1/32000;
% % % % Vc2=Vc1/32000;
% % % % Va3=roundn(Va2,-3);
% % % % Vb3=roundn(Vb2,-3);
% % % % Vc3=roundn(Vc2,-3);
Va1=Va;
Vb1=Vb;
Vc1=Vc;
Va2=Va1/16000;
Vb2=Vb1/16000;
Vc2=Vc1/16000;
Va3=roundn(Va2,-3);
Vb3=roundn(Vb2,-3);
Vc3=roundn(Vc2,-3);
Va4=Va3*1000;
Vb4=Vb3*1000;
Vc4=Vc3*1000;
Va5=roundn(Va4,0);
Vb5=roundn(Vb4,0);
Vc5=roundn(Vc4,0);
Va6=Va5*2;
Vb6=Vb5*2;
Vc6=Vc5*2;
figure(3)
plot(V.time(1500:3000),Va6,'r-')
hold on
plot(V.time(1500:3000),Vb6,'b:')
hold on
plot(V.time(1500:3000),Vc6,'g--')
hold off
legend('Va','Vb','Vc')
%---------------------------------------------
Xv1before=Va6-Vb6;
Xv10=Xv.signals.values(1500:3000);
figure(4)
plot(Xv.time(1500:3000),Xv10,'r--')
hold on
plot(Xv.time(1500:3000),Xv1before,'b--')
hold off
 legend('Xv1','Xv1before') %Va6,Vb6做爲有符號的AB測試輸入數據
%---------------------------------------------
Xv11=Xv10/6000;
Xv12=roundn(Xv11,-3);
Xv13=Xv12*1000;
Xv14=roundn(Xv13,0);
Xv15=Xv14*2;
figure(5)
plot(Xv.time(1500:3000),Xv15,'r--')
hold on
plot(Xv.time(1500:3000),Xv1before,'b--')
hold off
legend('Xv1after','Xv1before')
%---------------------------------------------
Xv110=Xv10-4000;
Xv120=Xv110/2000;
Xv130=roundn(Xv120,-3);
Xv140=Xv130*1000;
Xv150=roundn(Xv140,0);
Xv160=Xv150*2;
figure(6)
plot(Xv.time(1500:3000),Xv160,'r--')
hold on
plot(Xv.time(1500:3000),Xv1before,'b--')
hold off
legend('Xv1after','Xv1before')  %Xv160做爲有符號的相模變換數據後數據輸入