Удалите контуры между субдоменами
[ удаляет контуры между субдоменами, заданными в dl1,bt1] = csgdel(dl,bt,bl)bl. При удалении контуров в bl делает анализируемую матрицу геометрии противоречивой, затем csgdel удаляет дополнительные сегменты границы (сегменты ребра между субдоменами), чтобы сохранить непротиворечивость.
Удаление контуров обычно изменяет идентификаторы ребра остающихся контуров.
csgdel не удаляет граничные сегменты (внешние границы).
Удалите ребра в 2D геометрии, созданной в приложении PDE Modeler и экспортируемой в MATLAB® рабочая область.
Создайте геометрию в приложении PDE Modeler путем ввода следующих команд в командное окно MATLAB:
pdecirc(0,0,1,'C1') pdecirc(0,0,0.5,'C2') pderect([-0.2 0.2 0.2 0.9],'R1') pderect([0 1 0 1],'SQ1')
Уменьшайте геометрию до первого квадранта путем пересечения его с квадратом. Для этого введите (C1+C2+R1)*SQ1 в поле Set formula.
Из приложения PDE Modeler экспортируйте матрицу описания геометрии, установите формулу и матрицу пространства имен к рабочему пространству MATLAB путем выбора Export Geometry Description, Set Formula, Labels из меню Draw.
В командном окне MATLAB используйте decsg функция, чтобы разложить экспортируемую геометрию на минимальные области. Это создает AnalyticGeometry объект dl. Постройте dl.
[dl,bt] = decsg(gd,sf,ns); pdegplot(dl,'EdgeLabels','on','FaceLabels','on') xlim([-0.1 1.1]) ylim([-0.1 1.1])
Удалите ребра 1, 2, и 13 использований csgdel функция. Задайте ребра, чтобы удалить как вектор из идентификаторов ребра. Постройте получившуюся геометрию.
[dl1,bt1] = csgdel(dl,bt,[1 2 13]); pdegplot(dl1,'EdgeLabels','on','FaceLabels','on') xlim([-0.1 1.1]) ylim([-0.1 1.1])
Теперь удалите все контуры между субдоменами и постройте получившуюся геометрию.
[dl1,bt1] = csgdel(dl,bt); pdegplot(dl1,'EdgeLabels','on','FaceLabels','on') xlim([-0.1 1.1]) ylim([-0.1 1.1])
