Сгенерируйте 2D mesh и найдите ее область.

Создайте модель PDE.

model = createpde;

Включайте геометрию встроенной функции lshapeg. Постройте геометрию.

geometryFromEdges(model,@lshapeg);
pdegplot(model,'FaceLabels','on')

Сгенерируйте mesh и постройте ее.

mesh = generateMesh(model);
figure
pdemesh(model)

Вычислите область целой mesh.

ma = area(mesh)

ma = 3.0000

Сгенерируйте 2D mesh и найдите область каждого элемента.

Создайте модель PDE.

model = createpde;

Включайте геометрию встроенной функции lshapeg. Постройте геометрию.

geometryFromEdges(model,@lshapeg);
pdegplot(model,'FaceLabels','on')

Сгенерируйте mesh и постройте ее.

mesh = generateMesh(model);
figure
pdemesh(model)

Вычислите область целой mesh и область каждого отдельного элемента mesh. Отобразите области первых 5 элементов.

[ma,mi] = area(mesh);
mi(1:5)

ans = 1×5

    0.0047    0.0054    0.0053    0.0048    0.0061

Найдите объединенную область элементов сопоставленной с конкретной поверхностью 2D mesh.

Создайте модель PDE.

model = createpde;

Включайте геометрию встроенной функции lshapeg. Постройте геометрию.

geometryFromEdges(model,@lshapeg);
pdegplot(model,'FaceLabels','on')

Сгенерируйте mesh и постройте ее.

mesh = generateMesh(model);
figure
pdemesh(model)

Найдите элементы сопоставленными с поверхностью 1 и вычислите общую площадь этих элементов.

Ef1 = findElements(mesh,'region','Face',1);
maf1 = area(mesh,Ef1)

maf1 = 1.0000

Найдите, сколько из общей области mesh принадлежит этим элементам. Возвратите результат как процент.

maf1_percent = maf1/area(mesh)*100

maf1_percent = 33.3333