Commit 8511da61 authored by Franco Buffa's avatar Franco Buffa
Browse files

source and data

parent ac802c29

dev/errmap.m

+13 −11
Original line number Diff line number Diff line
@@ -57,6 +57,10 @@ function errmap(fconf) 
 RHO2(blk)=[];

 DELTA2(blk)=[];

 W(blk)=[];

%%%%%%%%%  

 a0=[ones(length(X),1) X Y RHO2 RHO2.*X RHO2.*Y];

 c0=a0\W;

%%%%%%%%%  

 if (~exist('aber','var'))

  aber=0;

 end
@@ -95,17 +99,15 @@ function errmap(fconf) 
 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;

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

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

 zf=c0(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;
@@ -121,14 +123,14 @@ function errmap(fconf) 
   jy=round((jy)*r1/R)+mp/2;

   jy=mp+1-jy;

   for k=1:length(jx)

    u1(jx(k),jy(k))=tol3;

    u1(jx(k),jy(k))=-tol3;

   end  

  case(2) 

   fname=[datapath,'panel.dat'];

   load(fname);

   p1=map0.*(map0==1);

   ip1=find(p1==0);

   u1(ip1)=tol3;

   u1(ip1)=-tol3;

   u1(blk)=0;

 end

dev/radrms.m

+10 −3
Original line number Diff line number Diff line 
function [x y u]=radrms(fconf)

function radrms(fconf)

%function [x y u]=radrms(fconf)

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

%

% Function:
@@ -54,14 +55,20 @@ function [x y u]=radrms(fconf) 
 y(end-d:end)=y(end-d);

 u(end-d:end)=u(end-d);

 

 figure(9);

 figure(1);

 plot(x,y,'*-');

 grid minor on;

 axis([4.5 32]);

 xlabel ('Antenna radius [m]');

 ylabel ('RMS deflession [mm]');



 figure(10);

 figure(2);

 plot(x,u,'*-');

 grid minor on;

 axis([4.5 32]); 

 xlabel ('Antenna radius [m]');

 ylabel ('Mean deflession [mm]');

 

 save ('radrms.dat','x', 'y', 'u'); 

end

dev/zern.m

+2 −2
Original line number Diff line number Diff line
@@ -131,14 +131,14 @@ function zern(fconf) 
   jy=round((jy)*r1/R)+mp/2;

   jy=mp+1-jy;

   for k=1:length(jx)

    u4(jx(k),jy(k))=tol4;

    u4(jx(k),jy(k))=-tol4;

   end  

 case(2) 

   fname=[datapath,'panel.dat'];

   load(fname);

   p1=map0.*(map0==1);

   ip1=find(p1==0);

   u4(ip1)=tol4;

   u4(ip1)=-tol4;

   u4(blk)=0;

 end