gridxyraydata.m

zemax与matlab调用程序

function RayDataOut = gridXYRayData(x, y, xname, yname, radius)
% gridXYRayData - Generate a 2D grid of data in a struct with the given field names.
%
% RayDataOut = gridXYRayData(x, y, xname, yname, radius);
%
% x and y are are numeric vectors, xname and yname are strings, radius is a numeric scalar.
% Creates a vector struct (RayDataOut) with two fields called xname and yname. The elements of the struct have a
% meshgrid of values taken from the x and y vector inputs. Only elements with a hypotenuse (sqrt(x^2+y^2)) of less
% than or equal to the given value are included in the struct. This function is useful for generating data
% compatible for input to zArrayTrace. The "radius" input is used to include only points within a circular
% pupil for example to generate data for a mode 0 call to zArrayTrace with rays spread evenly over the pupil
% in a rectangular grid ...
%
% x = -1:0.1:1; y = -1:0.1:1;
% RayData = gridXYRayData(x, y, 'z', 'l', 1.0);
%
% The "radius" input is only used if specified greater than 0.
% See also zArrayTrace, genRayDataMode0

% MZDDE - The ZEMAX DDE Toolbox for Matlab.
% Copyright (C) 2002-2004 Defencetek, CSIR
% Contact : [e-mail]
% 
% This file is part of MZDDE.
% 
%  MZDDE is free software; you can redistribute it and/or modify
%  it under the terms of the GNU General Public License as published by
%  the Free Software Foundation; either version 2 of the License, or
%  (at your option) any later version.
%
%  MZDDE is distributed in the hope that it will be useful,
%  but WITHOUT ANY WARRANTY; without even the implied warranty of
%  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
%  GNU General Public License for more details.
%
%  You should have received a copy of the GNU General Public License
%  along with MZDDE (COPYING.html); if not, write to the Free Software
%  Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
%


% Copyright 2003, Defencetek, CSIR
% $Revision: 1.0 $


[X,Y] = meshgrid(x,y);
if radius > 0.0
    r = sqrt(X.*X + Y.*Y); % Find the distance of each point from the origin
    list = find(r <= radius); % Get a list of indices for which the distance is less than radius
    % Create a struct with two fields having the given names 
    RayDataOut = struct(xname, num2cell(X(list)), yname, num2cell(Y(list)));
else
    RayDataOut = struct(xname, num2cell(reshape(X,1,numel(X))), yname, num2cell(reshape(Y,1,numel(Y))));
end