Своеобразие русской архитектуры: Основной материал – дерево – быстрота постройки, но недолговечность и необходимость деления...
Таксономические единицы (категории) растений: Каждая система классификации состоит из определённых соподчиненных друг другу...
Топ:
Когда производится ограждение поезда, остановившегося на перегоне: Во всех случаях немедленно должно быть ограждено место препятствия для движения поездов на смежном пути двухпутного...
Основы обеспечения единства измерений: Обеспечение единства измерений - деятельность метрологических служб, направленная на достижение...
История развития методов оптимизации: теорема Куна-Таккера, метод Лагранжа, роль выпуклости в оптимизации...
Интересное:
Аура как энергетическое поле: многослойную ауру человека можно представить себе подобным...
Берегоукрепление оползневых склонов: На прибрежных склонах основной причиной развития оползневых процессов является подмыв водами рек естественных склонов...
Что нужно делать при лейкемии: Прежде всего, необходимо выяснить, не страдаете ли вы каким-либо душевным недугом...
Дисциплины:
2023-02-03 | 29 |
5.00
из
|
Заказать работу |
|
|
Разберем идею разностного метода решения .краевых задач на примере взаимодействия световых пучков (см. рис. 12.1), переобозначив в системе(12.1-1) интенсивность излучения вправо на Y, а интенсивность излучения влево на у (просто в целях удобства, чтобы не писать индекс).
рис 12.1
Суть метода заключается в покрытии расчетного интервала сеткой из N точек. Тем самым определяются (N-i) шагов (рис. 12.7). Затем надо заменить дифференциальные уравнения исходной краевой задачи аппроксимирующими их уравнениями в конечных разностях, выписав соответствующие разностные уравнения для каждого 1-го шага. В нашем случае достаточно просто заменить первые производные из (12.1-1) их разностными аналогами (такой метод называется еще методом Эйлера):
Рис. 12.7. Сетка, покрывающая расчетный интервал
Примечание
Существует множество способов аппроксимации дифференциальных уравнений разностными. От выбора конкретного варианта зависит не только простота, быстрота и удобство вычислений, но и сама возможность получения правильного ответа.
Получилась система (по числу шагов) 2-(м-1) разностных линейных алгебраических уравнений с 2-N неизвестными YI и yi. Для того чтобы она имела единственное решение, надо дополнить число уравнений до 2-м. Это можно сделать, записав в разностном виде оба граничных условия:
Y0=1O, yN=R-YN.
Сформированная полная система алгебраических уравнений называется разностной схемой, аппроксимирующей исходную краевую задачу. Обратите внимание, что правые части разностных уравнений системы на каждом шаге записаны для левой границы шага. Такие разностные схемы называют явными, т. к. все значения YI+I и yi+i находятся в левой части уравнений. Полученную явную разностную схему легко записать в матричной форме
|
A.Z=B,
где z - неизвестный вектор, получающийся объединением векторов Y и у. Решив систему (3), мы получим решение краевой задачи.
Примечание
На самом деле, все несколько сложнее, поскольку, вообще говоря, необходимо еще доказать, что, во-первых, разностная схема действительно аппроксимирует дифференциальные уравнения и, во-вторых, при м-*» разностное решение действительно сходится к дифференциальному.
22. Устойчивые и асимптотически устойчивые решения разностных уравнений
Рассмотрим дифференциальное уравнение
Пусть некоторое фиксированное решение y = φ(x) этого уравнения существует при всех x≥ x0 .
Решение y = φ(x) уравнения называется устойчивым по Ляпунову при x ≥ x0 , если для любого ε > 0 существует число δ > 0 (зависящее, вообще говоря, от ε ) такое, что:
— решение y = y(x) задачи Коши с начальным условием y(x0 ) , |y(x0) − φ(x0) | < δ , существует при всех x ≥ x0 ;
— для всех таких решений справедливо неравенство |y(x) − φ(x) | < ε , при всех x > x0 .
Геометрически это означает, что интегральные кривые y = y(x), близкие в момент x = x0 к интегральной кривой y = φ(x), остаются близкими к ней и на всем промежутке [x0, ∞) .
На рисунке красным изображено устойчивое решение задачи Коши y' = − y, y(1) = 1.
Видно, что все интегральные кривые, близкие к этому решению в начальный момент x = 1, остаются вблизи него и при x > 1 .
Решение y = φ(x) называется неустойчивым по Ляпунову при x ≥ x0 , если существует число ε > 0 такое, что для любого δ > 0
найдутся решения y = yδ(x) и значение x1 = x1(δ) > x0 такие, что хотя | yδ( x0) − φ( x0) | < δ , но |y( x1) − φ( x1) | ≥ ε .
На рисунке красным изображено неустойчивое решение y = 0 задачи Коши y' = sin2 y, y(0) = 0.
Видно, что интегральные кривые, близкие к y = 0 в начальный момент x0 = 0, удаляются от y = 0 с ростом x > 1 .
|
2f(n+2)-2f(n+1)+f(n)=0
2л^2-2л+1=0
D= - 4
л(1)=(2+2i)/4=(1+i)/2
л(2)=(2-2i)/4=(1-i)/2
y=(e)^0.5x*(C1*cos0.5x+C2*sin0.5x)
(0.5)^n*(C1*cos(n*pi)/4+C2*sin(n*pi)/4)
л=a+ib
л=sqr(a*a+b*b)>=1
1) Если л1 и л2 <1, то решение устойчиво
2) Если л1 и л2 >1, то решение неустойчиво
3) Если л1 и л2 =1, то решение устойчиво (но не ассимптотически)
4) Если хотя бы одно л>1, то неустойчиво
Размещено на Allbest.ru
|
|
История создания датчика движения: Первый прибор для обнаружения движения был изобретен немецким физиком Генрихом Герцем...
Архитектура электронного правительства: Единая архитектура – это методологический подход при создании системы управления государства, который строится...
Археология об основании Рима: Новые раскопки проясняют и такой острый дискуссионный вопрос, как дата самого возникновения Рима...
Эмиссия газов от очистных сооружений канализации: В последние годы внимание мирового сообщества сосредоточено на экологических проблемах...
© cyberpedia.su 2017-2024 - Не является автором материалов. Исключительное право сохранено за автором текста.
Если вы не хотите, чтобы данный материал был у нас на сайте, перейдите по ссылке: Нарушение авторских прав. Мы поможем в написании вашей работы!