1 %# ssq=ssquare2(beta,f,g) returns sin(pi*beta*(f^2+g^2))
2 %# beta is equal to lambda*D/l^2, i.e. the Fresnelnumber
3 %# for the complete spectrum (from 0 to fsample)
5 %# only f and g from 0 to fsample/2 should be given, the other part is
6 %# calculated using the periodicity in frequency-space
8 %# complete frequency range is given
9 %# allows for astigmatism (5th argument)
11 function ssq=ssquare2(betah,f,g,varargin)
18 f_isall = varargin{1};
21 f_isall = varargin{1};
25 ssq=sin((pi*betah)*f.^2+(pi*betav)*g.^2);
31 ssq=cat(2,cat(1,ssq,flipdim(ssq(2:n,:,:),1)),cat(1,flipdim(ssq(:,2:m,:),2),flipdim(flipdim(ssq(2:n,2:m,:),1),2)));