mncirc.m

来自「matlab仿真的文件包然后写出其具体功能(至少要20个字)。尽量不要让站长把时」· M 代码 · 共 25 行

function [h1,h2]=mncirc(M,N)
if nargin==1,  
   N=[1/3,1/4,1/6,1/12]*pi; 
   N=[-N(length(N):-1:1),N];
elseif nargin==0
   M=[1.3,1.4,1.5,2,3]; M=[M, 1./M];
   N=[1/3,1/4,1/6,1/12]*pi; 
   N=[-N(length(N):-1:1),N];
end
Mx=[]; My=[]; Nx=[]; Ny=[];
t=[0:pi/30:2*pi];
for i=1:length(M)
   r=M(i)/(1-M(i)^2); x0=r*M(i); 
   Mx=[Mx; r*cos(t)+x0]; 
   My=[My; r*sin(t)];
end
for i=1:length(N)
   r=sqrt(1+N(i)^2)/(2*N(i)); 
   x0=-1/2; y0=1/(2*N(i));
   Nx=[Nx; r*cos(t)+x0]; 
   Ny=[Ny; r*sin(t)+y0];
end
h1=line(Mx',My'); h2=line(Nx',Ny');
set([h1;h2],'Color',[1,0,0],'LineStyle',':');