MATLAB.Exponenta
MATLAB и Simulink на русском
Технологии разработки и отладки
		сложных технических систем
 

Математика\Partial Differential Equations Toolbox

Список функций Partial Differential Equations Toolbox: Полезные функции

  В оглавление \ К следующему разделу \ К предыдущему разделу

PDEJMPS
Ошибка в процессе адаптивного переопределения конечноэлементной сетки

Синтаксис.

errf=pdejmps(p,t,c,a,f,u,alfa,beta,m)

Описание:

Вычисляет ошибку решения PDE для адаптивного переопределения конечноэлементной сетки. Столбцам матрицы errf соответствуют треугольники, строкам - уравнения в системе PDE.

Входные параметры.

  • p - массив узлов конечноэлементной сетки (столбцам соответствуют узлы):
    - первая строка - горизонтальные координаты узлов,
    - вторая строка - вертикальные координаты узлов;
  • t - матрица треугольных конечных элементов (столбцам соответствуют треугольники):
    t(1:3,ie) - глобальные номера узлов треугольника с номером ie,
    t(4,ie) - номер зоны, которой принадлежит треугольник с номером ie.
  • с - массив, описывающий распределение коэффициента c в расчётной области (см. уравнение (1) в assempde);
  • a - массив, описывающий распределение коэффициента a в расчётной области (см. уравнение (1) в assempde);
  • f - массив, описывающий распределение правой части PDE f в расчётной области (см. уравнение (1) в assempde).
  • u – узловое распределение решения PDE (строкам соответствуют узлы, столбцам - переменные PDE в системе).

errf(:,K) вычисляется по формуле:

errf(:,K)=alfa*L2K(H^m*(f - AU) + beta*sqrt(0.5*sum((L(J)^m*JMP(J))^2)),

где L2K, H - L2-норма и линейный размер треугольника с номером K; AU - оператор левой части PDE, применённый к u; L(J) - длина J-й стороны K-го треугольника; JMP(J) - нормальная производная решения PDE на J-й стороне K-го треугольника.

  • alfa, beta, m - параметры представленной формулы.

   В оглавление \ К следующему разделу \ К предыдущему разделу

 

Поиск по сайту:

Система Orphus

Яндекс.Метрика