daetools:: findDecoupledBlocksИщите разъединенные блоки в системах уравнений
Блокноты MuPAD® будут демонтированы в будущем релизе. Используйте live скрипты MATLAB® вместо этого.
Live скрипты MATLAB поддерживают большую часть функциональности MuPAD, хотя существуют некоторые различия. Для получения дополнительной информации смотрите, Преобразовывают Notebook MuPAD в Live скрипты MATLAB.
[eqsBlocks,varsBlocks] := daetools::findDecoupledBlocks(eqs,vars)
[eqsBlocks,varsBlocks] := daetools::findDecoupledBlocks(eqs,vars) идентифицирует подмножества (блоки) уравнений, которые могут использоваться, чтобы задать подмножества переменных. Количество переменных vars должно совпасть с количеством уравнений eqs.
i th блок является системой уравнений, определяющей переменные в vars[varsBlocks[i]]. Переменные в vars[varsBlocks[1],...,varsBlocks[i-1]] определяются рекурсивно предыдущими блоками уравнений. После того, как вы решаете первый блок уравнений для первого блока переменных, второй блок уравнений, данных eqs[eqsBlocks[2]], задает разъединенное подмножество уравнений, содержащих только подмножество переменных, данных вторым блоком переменных vars[varsBlock[2]], плюс переменные из первого блока (эти переменные известны в это время). Таким образом, если нетривиальная блочная декомпозиция возможна, можно разделить процесс решения для большой системы уравнений, включающей много переменных в несколько шагов, где каждый шаг включает меньшую подсистему.
Количеством блоков является nops(eqsBlocks). Это совпадает с nops(varsBlocks). Если nops(eqsBlocks) = nops(varsBlocks) = 1, то нетривиальная блочная декомпозиция уравнений не возможна.
Реализованный алгоритм требует, чтобы для каждой переменной в vars было по крайней мере одно соответствующее уравнение в eqs, включающем эту переменную. То же уравнение не может также быть соответствующим к другой переменной. Если система не удовлетворяет это условие, то daetools::findDecoupledBlocks выдает ошибку. В частности, daetools::findDecoupledBlocks требует того nops(eqs) = nops(vars).
Вычислите блок нижнее треугольное разложение символьной системы дифференциальных алгебраических уравнений (ДАУ).
Создайте следующую систему четырех дифференциальных алгебраических уравнений. Здесь, выражения x1(t), x2(t), x3(t) и x4(t) представляют переменные состояния системы. Система также содержит символьные параметры c1, c2, c3, c4 и выражения f(t,x,y) и g(t,x,y).
eqs := [c1*diff(x1(t),t) + c2*diff(x3(t),t) = c3*f(t,x1(t),x3(t)),
c2*diff(x1(t),t) + c1*diff(x3(t),t) = c4*g(t,x3(t),x4(t)),
x1(t) = g(t,x1(t),x3(t)),
x2(t) = f(t,x3(t),x4(t))]:
vars:= [x1(t), x2(t), x3(t), x4(t)]:Используйте daetools::findDecoupledBlocks, чтобы найти блочную структуру системы.
[eqsBlocks, varsBlocks] := daetools::findDecoupledBlocks(eqs, vars)
Первый блок содержит два уравнения в двух переменных.
eqs[eqsBlocks[1]]
vars[varsBlocks[1]]
После того, как вы решите этот блок для переменных состояния x1(t), x3(t), можно решить следующий блок уравнений. Этот блок состоит из одного уравнения.
eqs[eqsBlocks[2]]
Этот блок включает одну переменную.
vars[varsBlocks[2]]
После того, как вы решите уравнение от блока 2 для переменной состояния x4(t), остающийся блок уравнений, eqs[eqsBlocks[3]] задает остающуюся переменную vars[varsBlocks[3]].
eqs[eqsBlocks[3]]; vars[varsBlocks[3]]
Найдите перестановки, которые преобразовывают систему в блок нижняя треугольная форма.
eqsPerm := [op(eqsBlocks[i]) $ i = 1..nops(eqsBlocks)]; varsPerm := [op(varsBlocks[i]) $ i = 1..nops(varsBlocks)]
Преобразуйте систему в блок нижняя треугольная система уравнений.
eqs := eqs[eqsPerm]; vars := vars[varsPerm]
Найдите матрицу падения получившейся системы. Матрица падения показывает, что система переставленных уравнений имеет три диагональных блока размера 2-by-2, 1-by-1 и 1-by-1.
daetools::incidenceMatrix(eqs, vars)
|
Список или вектор уравнений или выражений в переменных состояния |
|
Список или вектор идентификаторов или выражений, таких как |
Вложенный список целых чисел, представляющих перестановки, требуемые преобразовывать исходную систему eqs,vars в блок нижняя треугольная форма.