% COMPLEX PLANE-PROJECTIVE CHIRP
%
% planechirp given by project2 acting on cis(2pi ux) * cis(2pi vy)
%
% CALL FORMAT = planechirp(u,v,M,N,angle,f1,f2);
%
% u,v in cycles per picture for non-chirped wave
%
% u and v must be scalars, as Matlab cannot have three dimensional matrices 
%
% DEFAULTS: f2 = f1; equal focal lengths if not specified
%
% EXAMPLE:  dispmatb(planechirp(2,1,200,100,pi/4,1))
%   produces a 200 by 100 matrix of a chirped 1 by -2 Hz planewave
%
%
% Y
% ^
% |                              X size given by N
% |                              Y size given by M
% |
% |
% ---------------> X

function planechirp = dummy_name(u,v,M,N,angle,f1,f2,X0,Y0)

i=sqrt(-1);

if nargin<9
  Y0=0;
end%if
if nargin<8
  X0=0;
end%if
if nargin < 7
  f2 = f1;
end%if

[X,Y] = meshdom( (0:N-1)*u/N, (0:M-1)*v/M);
[X2, Y2] = project2(X-X0,Y-Y0,angle,f1,f2);

planechirp = exp(2*pi*i*(X2+Y2)); % 

%%%dispmatb(planechirp);   % display real part as byte image