matlab生成网格编号每个单元的节点号

根据边界方程自动划分网格生成网格编号和每个单元的节点号。可设置网格的长宽，基于matlab。

clear;
k=0;
xstp=1000;
ystp=1000; %没步的大小  
for i=0:xstp:2000
  ymax=1000; % ymax=-0.12*i+150;
  ymin=0;    % ymin=0.12*i;
    
    for j=0:ystp:1000
        if j>=ymin && j<=ymax
            k=k+1;
            coord（k:）=[ij];
            plot（ij ‘.‘）;
            text（i j num2str（k） ‘fontsize‘ 8）;
            hold on;
        end
        if abs（j-ymin）            k=k+1;
            coord（k:）=[iymin];
            plot（iymin ‘.‘）;
            text（iyminnum2str（k） ‘fontsize‘ 8）;
             hold on;
        end
         if abs（j-ymax）ymax
            k=k+1;
            coord（k:）=[iymax];
            plot（iymax ‘.‘）;
            text（iymaxnum2str（k） ‘fontsize‘ 8）;
             hold on;
        end
    end
end
length（coord）
plot（[02000][00]）;
plot（[02000][10001000]）;
axis equal;          
pp=0;
%{
    *

      1st column — The unique values of x
    *

      2nd column — The number of instances of each value
    *

      3rd column — The percentage of each value
%}
kk=tabulate（coord（:1））

lnods=setelment（kkystp）
disp（‘element number is‘）
bb=size（lnods）;
ofile1=fopen（‘Meshout.txt‘‘w‘）;
fprintf（ofile1‘Results\n‘）;
for i=1:bb（1） 
fprintf（ofile1‘%d %d %d %d %d\n‘ilnods（i:））;
end
fprintf（ofile1‘coords\n‘）;
dd=size（coord）;
for i=1:dd（1）
    fprintf（ofile1‘%d %9.2f %9.2f\n‘icoord（i:））;
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        1441  2020-09-20 08:36  meshcoord.m
     文件         167  2020-09-20 08:36  Meshout.txt