Студопедия

Главная страница Случайная страница

Разделы сайта

АвтомобилиАстрономияБиологияГеографияДом и садДругие языкиДругоеИнформатикаИсторияКультураЛитератураЛогикаМатематикаМедицинаМеталлургияМеханикаОбразованиеОхрана трудаПедагогикаПолитикаПравоПсихологияРелигияРиторикаСоциологияСпортСтроительствоТехнологияТуризмФизикаФилософияФинансыХимияЧерчениеЭкологияЭкономикаЭлектроника






  • Как продвинуть сайт на первые места?
    Вы создали или только планируете создать свой сайт, но не знаете, как продвигать? Продвижение сайта – это не просто процесс, а целый комплекс мероприятий, направленных на увеличение его посещаемости и повышение его позиций в поисковых системах.
    Ускорение продвижения
    Если вам трудно попасть на первые места в поиске самостоятельно, попробуйте технологию Буст, она ускоряет продвижение в десятки раз, а первые результаты появляются уже в течение первых 7 дней. Если ни один запрос у вас не продвинется в Топ10 за месяц, то в SeoHammer за бустер вернут деньги.
    Начать продвижение сайта
  • Краевые задачи для обыкновенных дифференциальных уравнений






    Лабораторная работа 2

    Численные методы решения

    обыкновенных дифференциальных уравнений (4 часа)

     

    Краевые задачи для обыкновенных дифференциальных уравнений

    Примерами краевых задач для обыкновенных дифференциальных уравнений являются следующие:

    1) Найти функцию u(x), которая удовлетворяет на отрезке [a, b] уравнению

    u" (x)=-f(x), (20)

    а на концах отрезка условиям

    u(a)=u(b)=0. (21)

    Задача (20), (21) имеет следующее физическое состояние. Между точками x= a и x= b натянута упругая струна, находящаяся под действием внешней изгибающей нагрузки. f(x) – величина нагрузки, а u(x) – прогиб струны в безразмерных единицах.

    2) Двухточечная краевая задача для линейного дифференциального уравнения второго порядка:

    u" (x)+ p(x) u¢ (x)+ q(x) u(x)= f(x), (22)

    α 0 u(a)+ α 1 u'(a)= A,

    β 0 u(b)+ β 1 u'(b)= B. (23)

    Рассмотрим частный случай задачи (22), (23):

    u" (x)+ p(x) u(x)= f(x), x [0, X], (24)

    u(0)= a, u(X)= b. (25)

    Введем сетку xi= ih, i= 1, 2, …, N; h= X/N. Обозначим через ui приближенное значение u(x) в узлах сетки. Рассмотрим уравнение (24) во внутренних узлах xi, i= 1, 2, …, N-1 и заменим производную второго порядка разностной формулой

    u" (xi)= . (26)

    Тогда из (24), (25) получим для определения ui систему линейных уравнений

    (27)

    Система (27) при p(x)≥ 0 имеет решение.

    Система (27) представляет собой систему линейных уравнений с трехдиагональной матрицей

    (28)

    и для ее решения применим метод прогонки, который фактически является методом исключения неизвестных Гаусса.

    1. Прямой ход прогонки. Запишем первое уравнение (28) в виде

    u1= α 1u2 + β 1, α 1= , β 1= . (29)

    Подставив (29) во второе уравнение (28) и упростив выражения, увидим, что можно для ui получить формулы

    ui= α iui+1 + β i, ,

    , i= 2, …, N-2 (30)

    Из последнего уравнения (28), учитывая (30) при i= N – 2, получим

    uN-1= β N-1, α N-1= 0, β N-1= . (31)

    2. Обратный ход прогонки. После вычисления прогоночных коэффициентов можно найти значения решения задачи. Формула (31) дает значение uN-1:

    uN-1= β N-1. (32)

    Остальные значения вычисляем в обратном порядке:

    ui= α iui+1 + β i, , i= N-2, N-3, …, 2, 1. (33)

    Рассмотрим метод пристрелки на примере решения линейной краевой задачи (24), (25). Если известны частное решение u(x) неоднородного уравнения u" (x)- p(x) u(x)= f(x) и два линейно независимых решения u1(x), u2(x) однородного уравнения u" (x) - p(x) u(x)= 0, то общее решение неоднородного дифференциального уравнения u" (x)- p(x) u(x)= f(x) можно записать в виде

    u(x) + C1u1(x) + C2u2(x).

    Постоянные C1, C2 можно определить из краевых условий (25).

    В методе пристрелки используется следующий способ.

    Сначала находят частное решение u(x) неоднородного уравнения u" (x)- p(x)u(x)= f(x), удовлетворяющее условию u(0)= a, и частное решение u1(x) однородного уравнения u" (x) - p(x) u(x)= 0 удовлетворяющее условию u(0)=0.

    Затем общее решение u(x) неоднородного уравнения u" (x)- p(x) u(x)= f(x), удовлетворяющее условию u(0)= a, записывают в виде u(x) + Cu1(x). Остается найти постоянную С из условия u(x) + Cu1(x)= b.

    Приведем сеточный аналог метода пристрелки. Пусть краевая задача приведена к системе уравнений (27)

    (34)

    Найдем частные решения неоднородной (35) и однородной (36) систем уравнений:

    (35)

    (36)

    Выбирая произвольные значения для , (при этом должно быть ≠ 0), находим из (35) и (36) формулы для вычисления частных решений:

    (37)

    (38)

    Найдем С из условия и запишем решение:

    (39)

    (40)

    Алгоритм метода пристрелки заключается в том, чтобы выбрать шаг h и выполнить последовательно вычисления по формулам (37) - (40).






    © 2023 :: MyLektsii.ru :: Мои Лекции
    Все материалы представленные на сайте исключительно с целью ознакомления читателями и не преследуют коммерческих целей или нарушение авторских прав.
    Копирование текстов разрешено только с указанием индексируемой ссылки на источник.