(Не рекомендуемый) Треугольные данные о геометрии
pdetrg не рекомендуется. Использование area вместо этого.
[ возвращает области отдельных треугольников mesh как векторный ar,g1x,g1y,g2x,g2y,g3x,g3y] = pdetrg(p,t)ar и компоненты градиента треугольной основы функционируют как векторы g1x,g1y,g2x,g2y,g3x,g3y.
Найдите область каждого элемента 2D mesh. Кроме того, возвратите половину отрицательного котангенса каждого угла и компоненты градиента треугольных основных функций для каждого элемента.
Сгенерируйте и постройте mesh для модульной геометрии квадрата.
[p,e,t] = initmesh(@squareg); pdemesh(p,e,t)
Вычислите область каждого отдельного элемента mesh и половины отрицательного котангенса каждого угла. Отобразите первые 5 элементов для каждого результата.
[ar,a1,a2,a3] = pdetrg(p,t); ar(1:5)
ans = 1×5
0.0126 0.0148 0.0144 0.0156 0.0118
a1(1:5)
ans = 1×5
-0.2819 -0.3905 -0.5332 -0.1812 -0.5237
a2(1:5)
ans = 1×5
-0.5124 -0.2842 -0.1613 -0.4616 -0.3267
a3(1:5)
ans = 1×5
-0.1329 -0.2061 -0.2362 -0.2588 -0.0928
Найдите область самого маленького и самый большой элемент mesh.
min(ar)
ans = 0.0061
max(ar)
ans = 0.0216
Используйте синтаксис с семью выходными аргументами, чтобы вычислить компоненты градиента треугольных основных функций для каждого элемента.
[ar,g1x,g1y,g2x,g2y,g3x,g3y] = pdetrg(p,t);