Научная литература
booksshare.net -> Добавить материал -> Физика -> Рабинович М.И. -> "Введение в теорию колебаний и волн." -> 716

Введение в теорию колебаний и волн. - Рабинович М.И.

Рабинович М.И. Введение в теорию колебаний и волн. — НИЦ, 2000. — 564 c.
Скачать (прямая ссылка): vvedenievteoriuvoln2000.djvu
Предыдущая << 1 .. 710 711 712 713 714 715 < 716 > 717 718 719 720 721 722 .. 942 >> Следующая

(r)(z) и 0(2) не являются непрерывными функциями (на промежутке 2 е
[0,4995; 0,5015] мы рассматривали верхнее решение (c)). Из сравнения
рисунков видно, что профили решений y(z) и @(z) изменились.
6.1.3. Метод многократной стрельбы
Метод стрельбы, описанный в предыдущем пункте, иногда не позволяет
получить удовлетворительные результаты. Так, решение соответствующих
задач Коши (включая дифференциальные уравнения в вариациях) может
оказаться практически невозможным при наличии сильной чувствительности к
начальным условиям. В таких случаях часто оказывается удобным
использовать метод многократной стрельбы (см., например, [6.34]).
Опишем кратко идею этого метода на примере задачи (6.1.1),
(6.1.2) с граничными условиями ГУ2 вида (6.1.4). По аналогии с п.
6.1.1 выберем на промежутке ze[0, 1] сетку узловых точек 2о = 0, zi, ...,
z"=l, z/+i > Zi, которая, вообще говоря,
19 М. Холодниок и др.
290
Глава 6
может не быть эквидистантной. В практических задачах эта сетка узловых
точек выбирается гораздо менее плотной, чем при использовании метода
конечных разностей.
На каждом подынтервале [zf, zi+1], i = 0, ..., п-1 уравнения (6.1.1),
(6.1.2) интегрируются независимо. Для этого нам необходимо задать
начальные условия в точках z,-, i = 0,
..., п-1; эти начальные значения мы обозначим как тр. В точке го = 0 уже
заданы два условия (6.1.4а), и поэтому здесь, аналогично тому, как это
делалось в п. 6.1.2, мы выбираем два дополнительных условия вида
(6.1.24). При этом вектор т)о будет иметь только две составляющих rjoi
=х(0) и 1102 = = у{0). Остальные векторы тр, t=l, ..., п-1 будут иметь по
четыре составляющих
4i = (X (Zi), х' (zt), у (г,), у' (Z[)). (6.1.45)
Обозначим решение соответствующей задачи Коши на каждом подынтервале [z*,
z,+i] через x(z) = x(z; гг.тр) (и аналогично для x'(z), y{z), y'(z)).
Метод многократной стрельбы состоит в следующем: так подобрать векторы
тр, г = 0, ... ..., п-1, чтобы функции x(z), x'(z), y(z), у (z),
полученные "частями" на отдельных подынтервалах [zi, zi+\], оказались
непрерывными и чтобы при этом были удовлетворены граничные условия
(6.1.4Ь).
Для этого в точках zi, i= 1, ..., п- 1, должны выполняться условия
fa fa; zt_ь ч,.,), х'(zt; zf_b i|i-i), У fa; z,_b Пг-i), /e ,
<и-,))=ч" (6Л'46)
а в точке zn= 1 (в соответствии с условиями (6.1.4b)) должны быть
выполнены соотношения
x'(U гп_и Ч"-1) = 0, y'(U zn-1, 4n-i) = 0. (6.1.47)
Уравнения (6.1.46) и (6.1.47) представляют собой систему из 4(п-1) + 2
нелинейных уравнений относительно 4(п-1)+2 неизвестных - составляющих
векторов ц0,Чь •••.'Пп-ь Эти уравнения можно решать с помощью любого
подходящего метода (например, с помощью метода Ньютона). Для вычисления
матрицы Якоби можно опять использовать вариационные дифференциальные
уравнения (см. п. 6.1.2). Указанный алгоритм легко модифицируется на
случай ГУ1 и ГУЗ; при этом в случае ГУЗ мы выбираем две неизвестных,
например х(0), у(0), а затем из уравнений (6.1.5а) находим остальные
неизвестные, т. е. x'(0), z/'(0). Далее, с помощью условий (6.1.5Ь) мы
легко получаем соотношения, аналогичные (6.1.47).
6.2. Зависимость стационарных решений от параметра
291
6.2. ЗАВИСИМОСТЬ СТАЦИОНАРНЫХ РЕШЕНИИ ОТ ПАРАМЕТРА
Так же как и для "сосредоточенных" систем (описываемых обыкновенными
дифференциальными уравнениями, см. § 5.2),
Рис. 6.2. Диаграмма стационарных решений задачи 17, Re = 625.
существует целый ряд методов и подходов, позволяющих находить зависимость
стационарных решений от параметров в случае "распределенных" систем,
описываемых уравнениями
19*
292
Глава 6
с частными производными. Набор методов оказывается здесь даже более
широким, поскольку здесь есть еще и проблема аппроксимации задачи в
бесконечномерном пространстве посредством задачи в пространстве конечной
размерности.
Наиболее простым подходом является, конечно, последовательное
использование какого-либо метода предыдущего параграфа для каждого
значения параметра. Слова "наиболее простым" означают, что при этом мы не
должны строить специальный алгоритм продолжения по параметру. Недостатком
подобного подхода является, однако, необходимость осуществлять переходы
через точки бифуркации с помощью вмешательства человека. В результате
время, которое необходимо затратить в этом случае на диалог "человек-
машина", существенно возрастает. Иногда, тем не менее, предпочтение
отдается именно этому подходу, особенно в случаях, когда нам нужно
построить лишь небольшую совокупность решений (задан небольшой диапазон
значений параметра) или же в случаях, когда частое появление
бифуркационных значений параметра не ожидается.
Результаты последовательного использования метода Ньютона для системы
разностных уравнений (6.1.23), которая возникает в задаче 17,
представлены на диаграмме стационарных решений, показанной на рис. 6.2. В
данном случае уже само решение системы (6.1.23) для больших п являлось
достаточно серьезной проблемой.
В последующих пунктах мы рассмотрим методы продолжения по параметру
Предыдущая << 1 .. 710 711 712 713 714 715 < 716 > 717 718 719 720 721 722 .. 942 >> Следующая

Реклама

c1c0fc952cf0704ad12d6af2ad3bf47e03017fed

Есть, чем поделиться? Отправьте
материал
нам
Авторские права © 2009 BooksShare.
Все права защищены.
Rambler's Top100

c1c0fc952cf0704ad12d6af2ad3bf47e03017fed