GA_Zernike.m

Go to the documentation of this file.

zernike = 13;
n_fraction = 15;
n_exponent = 8;
bias = (2^(n_exponent-1)-1);
depth = 1+n_fraction+n_exponent;
population = 1;
 
 
 
n0 = round(rand(zernike,depth,population));
 
w = zeros(zernike,population);
s = zeros(zernike,population);
exponent = zeros(zernike,population);
fraction = zeros(zernike,population);
 
s(:,:) = (-1.*ones(zernike,population)).^reshape(n0(:,1,:),[zernike,population]);
 
 
for n=2:1+n_exponent
    exponent(:,:) = double(exponent(:,:) + 2^(n-2)*reshape(n0(:,n,:),[size(n0,1),size(n0,3)]));
end
for n=2+n_exponent:depth
    fraction(:,:) = double(fraction(:,:) + 2^-(n-10)*reshape(n0(:,n,:),[size(n0,1),size(n0,3)]));
end
 
for z=1:1:zernike    
    w=s.*(1+fraction).*2.^(exponent-bias)
end
 
img = Zernike(w);
imshow(Phase_Image(img));