makegriddat.m

来自「模式识别工具箱,希望对大家有用!」· M 代码 · 共 48 行

function [griddat,x,y] = makegriddat(minx,maxx,miny,maxy,nrstepx,nrstepy)
%MAKEGRIDDAT Make uniform 2D grid.
%
% [griddat,x,y] = makegriddat(minx,maxx,miny,maxy,nrstepx,nrstepy)
%
% Make a uniform 2D grid of objects and store it in a 2xN array. This
% array is not directly a PRTools dataset, but can be used for making
% a 2D plot of a user function:
%
%   grid = makegriddat(0,10,0,10);
%   out = f(grid);
%   plotg(grid,out);
%
% see also: gendatgrid, plotg

% Copyright: D. Tax, R.P.W. Duin, [e-mail]
% Faculty of Applied Physics, Delft University of Technology
% P.O. Box 5046, 2600 GA Delft, The Netherlands


global GRIDSIZE;

if nargin<6
  if isempty(GRIDSIZE)
    nrstepy = 30;
  else
    nrstepy = GRIDSIZE;
  end
end
if nargin<5
  if isempty(GRIDSIZE)
    nrstepx = 10;
  else
    nrstepx = GRIDSIZE;
  end
end

stepx = (maxx-minx)/(nrstepx-1);
stepy = (maxy-miny)/(nrstepy-1);

x = minx:stepx:maxx;
y = miny:stepy:maxy;

[gx,gy] = meshgrid(x,y);
griddat=[reshape(gx,nrstepx*nrstepy,1), reshape(gy,nrstepx*nrstepy,1)];

return