domaincreator.m

本代码为基于matlab的求解偏微分方程的利用径向基函数的无网格方法的源码

figure(2)
cla
axis([-1 1 -1 1])
hold on
title('please mark boundary points by clicking the left mouse button');
xlabel('please press the right mouse button for the last boundary point ');
vertex=[];
u=1;
while u==1
    [x,y,u]=ginput(1);
    vertex=[vertex;x y];
    plot(vertex(:,1),vertex(:,2)  )
    plot(vertex(:,1),vertex(:,2) ,'+')
end
plot([vertex(:,1); vertex(1,1)],...
       [vertex(:,2) ;vertex(1,2)] )
    plot([vertex(:,1); vertex(1,1)],...
       [vertex(:,2) ;vertex(1,2)] ,'+')
intpoints=insidepolygon(RBFstepsize,vertex);
bndpoints=onpolygon(RBFstepsize,vertex);
test=intbnddistances(intpoints,bndpoints);
[i1,i2]=find(test<0.001);
intpoints(unique(i1),:)=[];
size(intpoints)
plot(intpoints(:,1),intpoints(:,2),'+')
plot(bndpoints(:,1),bndpoints(:,2),'o');
hold off