cardg.m

MATLAB科学计算与工程分析源代码源程序6

function [x,y]=cardg(bs,s) 
% 此函数用来编写心形线所围成的区域
% 编写日期：2007-7-1

nbs=4; 

if nargin==0  %如果没有输入参数
  x=nbs;   
  return 
end
dl=[  0      pi/2   pi       3*pi/2
      pi/2   pi     3*pi/2   2*pi;
      1      1      1        1
      0      0      0        0];

if nargin==1   %如果只有一个输入参数
  x=dl(:,bs);   
  return 
end 

x=zeros(size(s)); 
y=zeros(size(s)); 
[m,n]=size(bs); 
if m==1 & n==1,   
  bs=bs*ones(size(s)); % 扩展bs 
elseif m~=size(s,1) | n~=size(s,2),   
  error('bs must be scalar or of same size as s'); 
end 

nth=400; 
th=linspace(0,2*pi,nth); 
r=2*(1+cos(th)); 
xt=r.*cos(th); 
yt=r.*sin(th); 
th=pdearcl(th,[xt;yt],s,0,2*pi); 
r=2*(1+cos(th)); 
x(:)=r.*cos(th); 
y(:)=r.*sin(th);