specifyCoefficients

Задайте коэффициенты в модели PDE

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

Синтаксис

specifyCoefficients(model,Name,Value)

specifyCoefficients(model,Name,Value,RegionType,RegionID)

CA = specifyCoefficients(___)

Описание

Коэффициенты УЧП

solvepde решает УЧП формы

$m \frac{\partial^{2} u}{\partial t^{2}} + d \frac{\partial u}{\partial t} - \nabla \cdot (c \nabla u) + a u = f$

solvepdeeig решает задачи о собственных значениях УЧП формы

$\begin{array}{l} - \nabla \cdot (c \nabla u) + a u = λ d u \\ или \\ - \nabla \cdot (c \nabla u) + a u = λ^{2} m u \end{array}$

specifyCoefficients задает коэффициенты m, d, c, a и f в модели PDE.

пример

specifyCoefficients(model,Name,Value) задает заданные коэффициенты в каждом Name каждому сопоставил Value, и включает их в model. Необходимо задать все эти имена: mDCA, и f. Этот синтаксис применяет коэффициенты к целой геометрии.

Примечание

Включайте геометрию в model перед использованием specifyCoefficients.

пример

specifyCoefficients(model,Name,Value,RegionType,RegionID) коэффициенты присвоений для заданной области геометрии.

пример

CA = specifyCoefficients(___) возвращает указатель на содействующий объект присвоения в model.

Примеры

свернуть все

Задайте уравнение Пуассона

Скрипт Open Live Script

Задайте коэффициенты для уравнения Пуассона $- \nabla \cdot \nabla u = 1$ .

solvepde уравнения адресов формы

$m \frac{\partial^{2} u}{\partial t^{2}} + d \frac{\partial u}{\partial t} - \nabla \cdot (c \nabla u) + a u = f$ .

Поэтому коэффициенты для уравнения Пуассона $m = 0$ , $d = 0$ , $c = 1$ , $a = 0$ , $f = 1$ . Включайте эти коэффициенты в модель PDE L-образной мембраны.

model = createpde();
geometryFromEdges(model,@lshapeg);
specifyCoefficients(model,'m',0,...
                          'd',0,...
                          'c',1,...
                          'a',0,...
                          'f',1);

Задайте нуль граничные условия Дирихле, поймайте в сети модель и решите УЧП.

applyBoundaryCondition(model,'dirichlet', ...
                             'Edge',1:model.Geometry.NumEdges, ...
                             'u',0);
generateMesh(model,'Hmax',0.25);
results = solvepde(model);

Просмотрите решение.

pdeplot(model,'XYData',results.NodalSolution)

Figure contains an axes object. The axes object contains an object of type patch.

Содействующий указатель для непостоянных коэффициентов

Скрипт Open Live Script

Задайте коэффициенты для уравнения Пуассона в 3-D с непостоянными характеристиками выброса и получите содействующий объект.

Коэффициенты уравнения $m = 0$ , $d = 0$ , $c = 1$ , $a = 0$ . Для непостоянных характеристик выброса взять $f = y^{2} \tanh (z) / 1000$ .

f = @(location,state)location.y.^2.*tanh(location.z)/1000;

Установите коэффициенты в 3-D прямоугольной геометрии блока.

model = createpde();
importGeometry(model,'Block.stl');
CA = specifyCoefficients(model,'m',0,...
                               'd',0,...
                               'c',1,...
                               'a',0,...
                               'f',f)

CA = 
  CoefficientAssignment with properties:

    RegionType: 'cell'
      RegionID: 1
             m: 0
             d: 0
             c: 1
             a: 0
             f: @(location,state)location.y.^2.*tanh(location.z)/1000

Установите нуль условия Дирихле на поверхности 1, поймайте в сети геометрию и решите УЧП.

applyBoundaryCondition(model,'dirichlet','Face',1,'u',0);
generateMesh(model);
results = solvepde(model);

Просмотрите решение на поверхности.

pdeplot3D(model,'ColorMapData',results.NodalSolution)

Задайте коэффициенты В зависимости от субдомена

Скрипт Open Live Script

Создайте скалярную модель PDE с L-образной мембраной как геометрия. Постройте метки субдомена и геометрия.

model = createpde();
geometryFromEdges(model,@lshapeg);
pdegplot(model,'FaceLabels','on')
axis equal
ylim([-1.1,1.1])

Figure contains an axes object. The axes object contains 4 objects of type line, text.

Установите c коэффициент к 1 во всех областях, но f коэффициент к 1 в субдомене 1, 5 в субдомене 2, и-8 в субдомене 3. Установите все другие коэффициенты на 0.

specifyCoefficients(model,'m',0,'d',0,'c',1,'a',0,'f',1,'Face',1);
specifyCoefficients(model,'m',0,'d',0,'c',1,'a',0,'f',5,'Face',2);
specifyCoefficients(model,'m',0,'d',0,'c',1,'a',0,'f',-8,'Face',3);

Установите нуль граничные условия Дирихле на все ребра. Создайте mesh, решите УЧП и постройте результат.

applyBoundaryCondition(model,'dirichlet', ...
                             'Edge',1:model.Geometry.NumEdges, ...
                             'u',0);
generateMesh(model,'Hmax',0.25);
results = solvepde(model);
pdeplot(model,'XYData',results.NodalSolution)

Figure contains an axes object. The axes object contains an object of type patch.

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

свернуть все

`model` — Модель PDE
`PDEModel` объект

Модель PDE в виде PDEModel объект.

Пример: model = createpde

Аргументы name-value

Пример: specifyCoefficients(model,'m',0,'d',0,'c',1,'a',0,'f',@fcoeff)

Примечание

Необходимо задать все эти имена: mDCA, и f.

`m` — Коэффициент производной времени второго порядка
скаляр | вектор-столбец | указатель на функцию

Коэффициент производной времени второго порядка в виде скаляра, вектор-столбца или указателя на функцию. Для получения дополнительной информации на размерах, и для деталей формы указателя на функцию коэффициента, см. m, d, или Коэффициент для specifyCoefficients.

Задайте 0 если термин не является частью вашей проблемы.

Пример: specifyCoefficients('m',@mcoef,'d',0,'c',1,'a',0,'f',1,'Face',1:4)

Типы данных: double | function_handle
Поддержка комплексного числа: Да

`d` — Коэффициент производной времени первого порядка
скаляр | вектор-столбец | указатель на функцию

Коэффициент производной времени первого порядка в виде скаляра, вектор-столбца или указателя на функцию. Для получения дополнительной информации на размерах, и для деталей формы указателя на функцию коэффициента, см. m, d, или Коэффициент для specifyCoefficients.

Примечание

Если m коэффициент является ненулевым, d должен быть 0 или матрица, и не указатель на функцию. См. d Коэффициент, Когда m будет Ненулевым.

Задайте 0 если термин не является частью вашей проблемы.

Пример: specifyCoefficients('m',0,'d',@dcoef,'c',1,'a',0,'f',1,'Face',1:4)

Типы данных: double | function_handle
Поддержка комплексного числа: Да

`c` — Коэффициент производной пробела второго порядка
скаляр | вектор-столбец | указатель на функцию

Коэффициент производной пробела второго порядка в виде скаляра, вектор-столбца или указателя на функцию. Для получения дополнительной информации на размерах, и для деталей формы указателя на функцию коэффициента, см. c Коэффициент для specifyCoefficients.

Пример: specifyCoefficients('m',0,'d',0,'c',@ccoef,'a',0,'f',1,'Face',1:4)

Типы данных: double | function_handle
Поддержка комплексного числа: Да

`a` — Коэффициент множителя решения
скаляр | вектор-столбец | указатель на функцию

Коэффициент множителя решения в виде скаляра, вектор-столбца или указателя на функцию. Для получения дополнительной информации на размерах, и для деталей формы указателя на функцию коэффициента, см. m, d, или Коэффициент для specifyCoefficients.

Задайте 0 если термин не является частью вашей проблемы.

Пример: specifyCoefficients('m',0,'d',0,'c',1,'a',@acoef,'f',1,'Face',1:4)

Типы данных: double | function_handle
Поддержка комплексного числа: Да

`f` — Исходный коэффициент
скаляр | вектор-столбец | указатель на функцию

Исходный коэффициент в виде скаляра, вектор-столбца или указателя на функцию. Для получения дополнительной информации на размерах, и для деталей формы указателя на функцию коэффициента, см. f Коэффициент для specifyCoefficients.

Задайте 0 если термин не является частью вашей проблемы.

Пример: specifyCoefficients('m',0,'d',0,'c',1,'a',0,'f',@fcoeff,'Face',1:4)

Типы данных: double | function_handle
Поддержка комплексного числа: Да

`RegionType` — Геометрический тип области
`'Face'` для 2D модели | `'Cell'` для 3-D модели

Геометрический тип области в виде 'Face' или 'Cell'.

Пример: specifyCoefficients('m',0,'d',0,'c',1,'a',0,'f',10,'Cell',2)

Типы данных: char | string

`RegionID` — Геометрический ID области
вектор из положительных целых чисел

Геометрический ID области в виде вектора из положительных целых чисел. Найдите идентификаторы области при помощи pdegplot.

Пример: specifyCoefficients('m',0,'d',0,'c',1,'a',0,'f',10,'Cell',1:3)

Типы данных: double

Выходные аргументы

свернуть все

`CA` — Содействующее присвоение
`CoefficientAssignment` объект

Содействующее присвоение, возвращенное как объект CoefficientAssignment Properties.

Больше о

свернуть все

`d` Коэффициент, когда `m` является Ненулевым

d коэффициент принимает специальную матричную форму когда m является ненулевым. Необходимо задать d как матрица конкретного размера, и не как указатель на функцию.

d представляет коэффициент демпфирования в случае ненулевого m. Задавать d, выполните эти два шага:

Вызовите results = assembleFEMatrices(...) для проблемы с вашими исходными коэффициентами и использованием d = 0. Используйте 'none' по умолчанию метод для assembleFEMatrices.
Возьмите d коэффициент как матрица размера results.M. Обычно d любой пропорционально results.M, или линейная комбинация results.M и results.K.

Смотрите динамику ослабленного консольного луча.

Советы

Для уравнений собственного значения коэффициенты не могут зависеть от решения u или его градиент.
Можно преобразовать дифференциальное уравнение с частными производными в необходимую форму при помощи Symbolic Math Toolbox™. pdeCoefficients (Symbolic Math Toolbox) преобразует УЧП в необходимую форму и извлекает коэффициенты в структуру, которая может использоваться specifyCoefficients.
pdeCoefficients функция также может возвратить структуру символьных выражений, в этом случае необходимо использовать pdeCoefficientsToDouble (Symbolic Math Toolbox), чтобы преобразовать эти выражения, чтобы удвоить формат прежде, чем передать их specifyCoefficients.

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

findCoefficients | PDEModel | pdeCoefficients (Symbolic Math Toolbox) | pdeCoefficientsToDouble (Symbolic Math Toolbox)

Темы

Введенный в R2016a

Документация

specifyCoefficients

Синтаксис

Описание

Примеры

Задайте уравнение Пуассона

Содействующий указатель для непостоянных коэффициентов

Задайте коэффициенты В зависимости от субдомена

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

`model` — Модель PDE
`PDEModel` объект

Аргументы name-value

`m` — Коэффициент производной времени второго порядка
скаляр | вектор-столбец | указатель на функцию

`d` — Коэффициент производной времени первого порядка
скаляр | вектор-столбец | указатель на функцию

`c` — Коэффициент производной пробела второго порядка
скаляр | вектор-столбец | указатель на функцию

`a` — Коэффициент множителя решения
скаляр | вектор-столбец | указатель на функцию

`f` — Исходный коэффициент
скаляр | вектор-столбец | указатель на функцию

`RegionType` — Геометрический тип области
`'Face'` для 2D модели | `'Cell'` для 3-D модели

`RegionID` — Геометрический ID области
вектор из положительных целых чисел

Выходные аргументы

`CA` — Содействующее присвоение
`CoefficientAssignment` объект

Больше о

`d` Коэффициент, когда `m` является Ненулевым

Советы

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

Темы

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

Поддержка

Документация

specifyCoefficients

Синтаксис

Описание

Примеры

Задайте уравнение Пуассона

Содействующий указатель для непостоянных коэффициентов

Задайте коэффициенты В зависимости от субдомена

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

model — Модель PDE PDEModel объект

Аргументы name-value

m — Коэффициент производной времени второго порядка скаляр | вектор-столбец | указатель на функцию

d — Коэффициент производной времени первого порядка скаляр | вектор-столбец | указатель на функцию

c — Коэффициент производной пробела второго порядка скаляр | вектор-столбец | указатель на функцию

a — Коэффициент множителя решения скаляр | вектор-столбец | указатель на функцию

f — Исходный коэффициент скаляр | вектор-столбец | указатель на функцию

RegionType — Геометрический тип области 'Face' для 2D модели | 'Cell' для 3-D модели

RegionID — Геометрический ID области вектор из положительных целых чисел

Выходные аргументы

CA — Содействующее присвоение CoefficientAssignment объект

Больше о

d Коэффициент, когда m является Ненулевым

Советы

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

Темы

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

Поддержка

`model` — Модель PDE
`PDEModel` объект

`m` — Коэффициент производной времени второго порядка
скаляр | вектор-столбец | указатель на функцию

`d` — Коэффициент производной времени первого порядка
скаляр | вектор-столбец | указатель на функцию

`c` — Коэффициент производной пробела второго порядка
скаляр | вектор-столбец | указатель на функцию

`a` — Коэффициент множителя решения
скаляр | вектор-столбец | указатель на функцию

`f` — Исходный коэффициент
скаляр | вектор-столбец | указатель на функцию

`RegionType` — Геометрический тип области
`'Face'` для 2D модели | `'Cell'` для 3-D модели

`RegionID` — Геометрический ID области
вектор из положительных целых чисел

`CA` — Содействующее присвоение
`CoefficientAssignment` объект

`d` Коэффициент, когда `m` является Ненулевым