Commit 1e7d1e2a authored by Franco Buffa's avatar Franco Buffa
Browse files

code

parent 70464b90

dev/errmap.m

+20 −17
Original line number Original line Diff line number Diff line
@@ -34,7 +34,7 @@ function errmap(fconf) 
 

 

 load unwr.dat

 load unwr.dat

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% u0=sfactor*u0;

 u0=sfactor*u0;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

 dx=D/(mp-1)/k0;

 dx=D/(mp-1)/k0;
@@ -77,21 +77,6 @@ function errmap(fconf) 
   otherwise

   otherwise

    error("invalid aber value");

    error("invalid aber value");

  endswitch

  endswitch

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

 if(aber<3)

   printf('********************************\n');

   xf=c(5)*2*F^3*lambda/pi;

   yf=c(6)*2*F^3*lambda/pi;

   zf=c(4)*2*F^2*lambda/pi;

%  phi=atan2((c(3)*lambda/(2*pi)-yf/F),(c(2)*lambda/(2*pi)-xf/F))*180/pi  

%  theta=((c(3)*lambda/(2*pi)-yf/F)/sind(phi))*180/pi

   printf('Feed displacements:\n');

   printf('x=%f (mm)\n',1E3*xf);

   printf('y=%f (mm)\n',1E3*yf);

   printf('z=%f (mm)\n',1E3*zf);

   printf('********************************\n');

 end  

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

 else

 else

  a=[ones(length(X),1) X Y];

  a=[ones(length(X),1) X Y];

  c=a\W;

  c=a\W;
@@ -105,7 +90,25 @@ function errmap(fconf) 
%i3=find(abs(u1)>tol3);

%i3=find(abs(u1)>tol3);

%u1(i3)=0;

%u1(i3)=0;





 u1=sfactor*u1;

 u2=reshape(u1,mp*mp,1);

 u2(blk)=[];

 stdg=std(u2);

 printf('********************************\n');

 printf('global RMS = %f [mm]\n',stdg);

 if((aber>0) && (aber<3))

   xf=c(5)*2*F^3*lambda/pi;

   yf=c(6)*2*F^3*lambda/pi;

   zf=c(4)*2*F^2*lambda/pi;

%  phi=atan2((c(3)*lambda/(2*pi)-yf/F),(c(2)*lambda/(2*pi)-xf/F))*180/pi  

%  theta=((c(3)*lambda/(2*pi)-yf/F)/sind(phi))*180/pi

   printf('Feed displacements:\n');

   printf('x=%f (mm)\n',1E3*xf);

   printf('y=%f (mm)\n',1E3*yf);

   printf('z=%f (mm)\n',1E3*zf);

 end  

 printf('********************************\n');

 

 %   u1=sfactor*u1;





 save errmap.dat u1

 save errmap.dat u1