Commit 39c8d82c authored by Franco Buffa's avatar Franco Buffa
Browse files

source and data

parent af91204e

dev/unwr3.m

+3 −15
Original line number Original line Diff line number Diff line
@@ -33,19 +33,7 @@ function unwr3(fconf) 
  run('conf.m');

  run('conf.m');

 end 

 end 

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

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

% 

 load phase.dat

 load phase.dat

%   keyboard;

 psi=e2;

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

 o=zeros(mp,mp);

 o(blk)=1;

 for i=1:mp

  m0=mean(psi(:,i));

  o(:,i)=o(:,i)*m0;

 end

 psi(blk)=0;

 psi=psi+o;

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

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

 I=eye(mp);

 I=eye(mp);

%

%
@@ -55,7 +43,7 @@ function unwr3(fconf) 
 A=B=L'*L;

 A=B=L'*L;

 [U,S,V]=svd(B);

 [U,S,V]=svd(B);

 s=diag(S);

 s=diag(S);

 C=-(round(psi/(2.0*pi)*L')*L+L'*round(L*psi/(2.0*pi)));

 C=-(round(e2/(2.0*pi)*L')*L+L'*round(L*e2/(2.0*pi)));

 D=C*U;

 D=C*U;

 K=zeros(mp,mp);

 K=zeros(mp,mp);

 for i=1:mp

 for i=1:mp
@@ -64,7 +52,7 @@ function unwr3(fconf) 
 end

 end

%

%

 K=round(K*U');

 K=round(K*U');

 u0=psi+2.0*pi*K;

 u0=e2+2.0*pi*K;

 u0(blk)=0;

 u0(blk)=0;

% 

% 

 figure(1);

 figure(1);

manual/hop-user-guide.pdf

+714 B (1.27 MiB)

File changed.

No diff preview for this file type.

View original file
View changed file