Дифференциальные уравнения с частными производными

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

Два дифференциальных уравнения в частных производных, которые часто возникают при анализе физических явлений, — это уравнение Пуассона:

и его однородная форма — уравнение Лапласа.

Mathcad имеет две функции для решения этих уравнений в области с квадратной границей. Используйте функцию relax, если известны значения, принимаемые неизвестной функцией u(x, y), на всех четырех сторонах квадрата.

Если функция u(x, y) равна нулю на всех четырех сторонах квадрата, можно использовать функцию multigrid. Эта функция зачастую решает задачу быстрее, чем relax. Обратите внимание, что, если граничные условия одинаковы на всех четырех сторонах квадрата, можно достаточно просто преобразовать уравнение к эквивалентному виду, в котором граничное значение является нулем на всех четырех сторонах.

Функция relax возвращает квадратную матрицу, в которой:

· Расположение элемента в матрице соответствует расположению точки внутри квадрата, и

· Значение элемента аппроксимирует решение задачи в этой точке.

Эта функция использует метод релаксации для нахождения приближенного решения. Уравнение Пуассона на квадратной области представляется в виде:

a_j,ku_j+1,k + b_j,ku_j-1,k + c_j,ku_j,k+1 + d_j,ku_j,k-1 + e_j,ku_j,k = f_j,k

Аргументы функции relax следующие:

relax(a, b, n, d, e, f, u, rjac) relax function

Если граничное условие — ноль на всех четырех сторонах квадрата, используйте функцию multigrid вместо функции relax. Пример приведен на Рисунке 9.

multigrid(M, ncycle)

Рисунок 9: Использование функции multigrid для решения уравнения Пуассона в квадратной области.