3点GPS定位MATLAB仿真

自己写的3点确定目标点的MATLAB仿真文件。 包含MATLAB画圆的程序

function h=circle（rx0y0CNb） 

% CIRCLE adds circles to the current plot 

% 

% CIRCLE（rx0y0） adds a circle of radius r centered at point x0y0. 

% If r is a vector of length L and x0y0 scalars L circles with radii r 

% are added at point x0y0. 

% If r is a scalar and x0y0 vectors of length M M circles are with the same 

% radius r are added at the points x0y0. 

% If r x0y0 are vector of the same length L=M M circles are added. （At each 

% point one circle）. 

% if r is a vector of length L and x0y0 vectors of length M~=L L*M circles are 

% added at each point x0y0 L circles of radius r. 

% 

% CIRCLE（rx0y0C） 

% adds circles of color C. C may be a string （‘r‘‘b‘...） or the RGB value. 

% If no color is specified it makes automatic use of the colors specified by 

% the axes ColorOrder property. For several circles C may be a vector. 

% 

% CIRCLE（rx0y0CNb） Nb specifies the number of points used to draw the 

% circle. The default value is 300. Nb may be used for each circle individually. 

% 

% h=CIRCLE（...） returns the handles to the circles. 

% 

% Try out the following （nice） examples: 

% 

%% Example 1 

% 

% clf; 

% x=zeros（1200）; 

% y=cos（linspace（01200）*4*pi）; 

% rad=linspace（10200）; 

% cmap=hot（50）; 

% circle（radxy[flipud（cmap）;cmap]）; 

% 

%% Example 2 

% 

% clf; 

% the=linspace（0pi200）; 

% r=cos（5*the）; 

% circle（0.1r.*sin（the）r.*cos（the）hsv（40））; 

% 

% 

%% Example 3 

% 

% clf 

% [xy]=meshdom（1:101:10）; 

% circle（[0.50.30.1]xy[‘r‘;‘y‘]）; 

% 

%% Example 4 

% 

% clf 

% circle（1:1000[]3:12）; 

% 

%% Example 5 

% 

% clf; 

% circle（（1:10）[002020][020200]）; 

% rewritten by Din-sue Fon. Dept. of Bio-Industrial Mechatronics Engineering 

% National Taiwan University March 102001 

% dsfong@ccms.ntu.edu.tw 

% written by Peter Blattner Institute of Microtechnology University of 

% Neuchatel Switzerland blattner@imt.unine.ch 

  

% Check the number of input arguments 

switch nargin 

case 0 

   r=[];x0=[];y0=[];C=[];Nb=[]; 

case 1 

   x0=[];y0=[];C=[];Nb=[]; 

case 2 

   y0=zeros（1length（x0））;C=[];Nb=[]; 

case 3 

   C=[];Nb=[]; 

case 4 

   Nb=[]; 

end 

  

if length（x0）~=length（y0） 

   if length（y0）==1 

      y0=ones（1length（x0））*y0; 

   elseif length（x0）==1 

      x0=ones（1length（y0））*x0; 

   else 

      error（‘The lengths of x0 and y0 must be identical‘）; 

   end; 

end; 

  

% set up the default values 

  

if isempty（r）r=1;end; 

if isempty（x0）x0=0;end; 

if isempty（y0）y0=0;end; 

if isempty（Nb）Nb=300;end; 

if isempty（C）C=get（gca‘colororder‘）;end; 

if isstr（C）C=C（:）;end; 

  

% work on the variable sizes 

  

x0=x0（:）; 

y0=y0（:）; 

r=r（:）; 

Nb=Nb（:）; 

  

  

% how many rings are plottet 

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       3862  2010-03-20 09:50  circle.m

     文件       1231  2010-04-28 14:58  threePointLocate.m

----------- ---------  ---------- -----  ----

                 5093                    2