meshcircle2.m

基于有限单元法的三角剖分圆的MATLAB程序

function meshcircle2
% meshcircle2  a practice for circle mesh
% synoposis:meshcircle2
% input:
% output:
%        cheng wf 2007 3
%
clear,clc,clf,hold on
theta=0:pi/20:2*pi;
x=0:10;y=0:10;r=2;
xr=3+2.*cos(theta);yr=3+2.*sin(theta);
n=0;xc=zeros(1,19);yc=zeros(1,19);
for i=1:19
    n=n+1;
       if i<=12
       xc(i)=3+r*cos((2*i*pi)/12);
       yc(i)=3+r*sin((2*i*pi)/12);
       end
      if i>12&i<19
       xc(i)=3+0.5*r*cos((4*i*pi)/12);
       yc(i)=3+0.5*r*sin((4*i*pi)/12);
       end
    if i==19
        xc(i)=3;
        yc(i)=3;
    end
end
X=xc';
Y=yc';
TRI=delaunay(X,Y)
[m,n]=size(TRI);
%[VX,VY]=voronoi(X,Y,TRI)
triplot(TRI,X,Y)
%plot(VX,VY)
axis ([0 10 0 10])
axis square
hold on;
%plot(X,Y,'or')
for i=1:19
text(X(i),Y(i),['( ',int2str(i), ')'])
end
for i=1:m
    x=(X(TRI(i,1))+X(TRI(i,2))+X(TRI(i,3)))/3;
    y=(Y(TRI(i,1))+Y(TRI(i,2))+Y(TRI(i,3)))/3;
    text(x,y,int2str(i))
end
shg