Как добиться одного масштаба при разных длин волн? Голограмма Френеля

clear;close all;
I1=imread('D:\Study\disser\hologramm\hologram.bmp');
Ih1=double(I1)-mean2(double(I1));
%Ih1=I1(:,:,1);
%xout = zeros(size(Ih1)+96);
%xout(1:end-96,1:end-96) = Ih1;
[N1,N2]=size(Ih1);
K=min(N1,N2);
I=Ih1(1:K,1:K);
 
h=0.000532;
d=295;
L=5.6;
%L=pix*N;
secim=0;
 
%e=0.000732;
while h:0.000632
 
pix=L/K;
 
%pix=0.0152;
%disp(pix);
 
%Z1=zeros(K,0);
%Z2=zeros(0,K);
%Xp=[Z1,[Z2;Ih1;Z2],Z1];
%I=double(Xp);
%pix=abs(d)*h/L;
n=-K/2:K/2-1;
x=n*pix;y=x;
[xx,yy]=meshgrid(x,y);
 
k=2*pi/h;
%zn=K+1;
%e=exp(-1i*2*pi*(xx.^2/K+yy.^2/K));
 
w=exp(1i*pi/(h*d)*(xx.^2+yy.^2));
%Fresnel=w.*e;
f2=I.*w;
Uf=fft2(f2,K,K);
Uf=fftshift(Uf);
%figure,imagesc(abs(Uf)),colormap('gray');
nn=-K/2:K/2-1;
vksi=1/(K*pix);
x=nn*vksi;
y=x;
[xx1,yy1]=meshgrid(x,y);
ph=exp(1i*k*d)/(1i*h*d)*exp(1i*pi*h*d*(xx1.^2+yy1.^2));
U0=Uf.*ph;
 
%figure,imagesc(abs(U0)),colormap('gray');
csr=pix/vksi;
dmin=(K*pix.^2)/h;
 
%фильтр усреднения
e1=fspecial('average');
X2 = imfilter(U0,e1,'replicate');
h=h+0.000005;
 
figure,imagesc(abs(X2)),colormap('gray');
 
%sovm=(X2+secim);
%secim=sovm;
%ft=fftshift(fft(sovm));
%figure,plot(abs(ft));
%figure,imagesc(abs(sovm)),colormap('gray');
 
%disp(sovm(256));
 
end