Exponenta Banner
exponenta event banner

pdeInterpolant

Interpolant для узловых данных к выбранным местоположениям

расширьте все на странице

pdeInterpolant и [p,e,t] представление FEMesh данные не рекомендуются. Используйте interpolateSolution и evaluateGradient интерполировать решение для УЧП и его градиент к произвольным точкам, не переключаясь на [p,e,t] представление.

Описание

interpolant позволяет вам оценивать решение для УЧП в любой точке в геометрии.

Решатели Partial Differential Equation Toolbox™ возвращают значения решения в узлах, означая точки mesh. Чтобы оценить интерполированное решение в других точках в геометрии, создайте pdeInterpolant объект, и затем вызывает evaluate функция.

Создание

Синтаксис

F = pdeInterpolant(p,t,u)

Описание

пример

F = pdeInterpolant(p,t,u) возвращает interpolant F на основе точек данных p, элементы t, и значения данных в точках, u.

Используйте meshToPet получить p и t данные для интерполяции с помощью pdeInterpolant.

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

развернуть все

Местоположения точки данных, заданные как матрица с двумя или тремя строками. Каждый столбец p 2D или 3-D точка. Для получения дополнительной информации смотрите Данные о Mesh.

Для 2D проблем создайте p использование initmesh функция или экспорт в меню Mesh приложения PDE Modeler. Для 2D или 3-D геометрии с помощью PDEModel возразите, получите p использование meshToPet функция на model.Mesh. Например, [p,e,t] = initmesh(g) или [p,e,t] = meshToPet(model.Mesh).

Элементы триангуляции, указанные как матрица. Для получения дополнительной информации смотрите Данные о Mesh.

Для 2D проблем создайте t использование initmesh функция или экспорт в меню Mesh приложения PDE Modeler. Для 2D или 3-D геометрии с помощью PDEModel возразите, получите t использование meshToPet функция на model.Mesh. Например, [p,e,t] = initmesh(g) или [p,e,t] = meshToPet(model.Mesh).

Значения данных, чтобы интерполировать, заданный как вектор или матрица. Как правило, u решение проблемы УЧП, возвращенной assempde, parabolic, hyperbolic, или другой решатель. Например, u = assempde(b,p,e,t,c,a,f). Можно также экспортировать u в меню Solve приложения PDE Modeler.

Размерности матричного u зависьте от проблемы. Если np количество столбцов p, и N количество уравнений в системе УЧП, затем u имеет N*np 'Строки' . Первый np строки соответствуют уравнению 1, следующий np строки соответствуют уравнению 2 и т.д. Для параболических или гиперболических проблем, u имеет один столбец в течение каждого раза решения; в противном случае, u вектор-столбец.

Функции объекта

evaluateИнтерполируйте данные к выбранным местоположениям

Примеры

свернуть все

Скрипт Open Live Script

В этом примере показано, как создать pdeInterpolant от решения до скалярного УЧП.

Решите уравнение -Δu=1 на единичном диске с нулем условия Дирихле. 

g0 = [1;0;0;1]; % circle centered at (0,0) with radius 1
sf = 'C1';
g = decsg(g0,sf,sf'); % decomposed geometry matrix
problem = allzerobc(g); % zero Dirichlet conditions
[p,e,t] = initmesh(g);
c = 1;
a = 0;
f = 1;
u = assempde(problem,p,e,t,c,a,f);

Создайте interpolant для решения. 

F = pdeInterpolant(p,t,u);

Оцените interpolant в четырех углах квадрата. 

pOut = [0,1/2,1/2,0;
    0,0,1/2,1/2];
uOut = evaluate(F,pOut)
uOut = 4×1

    0.2485
    0.1854
    0.1230
    0.1852

Значения uOut(2) и uOut(4) почти равны, как они должны быть для симметричных точек в этой симметричной проблеме.

Смотрите также

|

Темы

Введенный в R2014b

Документация Partial Differential Equation Toolbox
Поддержка