get(pluse1)
A=pluse1.Time;
B=pluse1.Data; 
t=zeros(10001,1);
t(1,1)=0;
for i=2:length(A)
    t(i,1)=0+(i*pi/200);
end
 x0=zeros(10001,1); 
  for k=1:length(A)   
     T(k)=cos(t(k,1))/sin(t(k,1));   
     x0(k)=B(k,1).*T(k);
  end
plot(A(1:400),B(51:450,1))
xlabel('t /s')
%polarplot(t(1:201,1),B(351:551,1))
% rlim ([-2,2])