
- •Содержание
- •6 Решение систем линейных алгебраических уравнений прямыми методами
- •6.1 Метод Гаусса с выбором главного элемента по столбцу
- •6.2 Метод Холецкого
- •6.3 Метод прогонки
- •Задание для самостоятельной работы
- •Вопросы
- •7 Решение систем линейных алгебраических уравнений итерационными методами
- •7.1 Метод Якоби
- •7.2 Метод Зейделя
- •7.3 Метод простой итерации.
- •Задание для самостоятельной работы
- •Вопросы
- •8 Приближение функций
- •8.1 Постановка задачи приближения функции по методу наименьших квадратов
- •8.2 Определение параметров эмпирической зависимости
- •8.3 Многочлены Бернштейна
- •Задание для самостоятельной работы
- •Вопросы
- •9 Интерполяция функций
- •9.1 Постановка задачи интерполяции функций
- •9.2 Оценка погрешности интерполяции
- •Задание для самостоятельной работы
- •Вопросы
- •10 Интерполяция сплайнами
- •10.1 Глобальная и кусочно-полиномиальная интерполяция
- •10.2 Интерполяция сплайнами
- •Задание для самостоятельной работы
- •Вопросы
- •11 Численное дифференцирование
- •11.1 Первая производная. Двухточечные методы
- •11.2 Вычисление первых производных по трёхточечным схемам
- •11.3 Вычисление производных второго порядка
- •11.4 Вычисление производных третьего порядка
- •Задание для самостоятельной работы
- •12 Численное интегрирование
- •Задание для самостоятельной работы
- •13 Решение задачи Коши одношаговыми методами
- •13.1 Постановка задачи Коши для дифференциального уравнения первого порядка
- •13.2 Численное решение задачи Коши методом Эйлера
- •13.3 Оценка погрешности метода Эйлера
- •13.4 Модификации метода Эйлера
- •Задание для самостоятельной работы
- •Варианты заданий
- •15 Преобразование Фурье
- •15.1 Аппроксимация функции по Фурье
- •2. Преобразование Фурье
- •15.2 Быстрое преобразование Фурье.
- •16 Решение дифференциальных уравнений в частных производных
- •16.1 Уравнение Лапласа (эллиптическое уравнение)
- •17 Решение дифференциальных уравнений в частных производных (продолжение)
- •17.1 Уравнение теплопроводности (параболическое уравнение)
- •Задание для самостоятельной работы
- •Синютин сергей алексеевич
- •Вычислительная математика (часть 2)
Задание для самостоятельной работы
Найти для x = 0.1m, m=0 – 10, t=0.02 решение уравнения ut = (1/a) uxx , 0 < x < 1, 0 < t <0.02, удовлетворяющее условиям u(0,t) =e at , u(1,t) = e a (t -1), 0 < t < 0.02; u(x,0) = e- ax , 0 < x < 1; a = 2 + 0.3 k, k = - 4.
а) по явной разностной схеме, взяв h = 0.1, t < ah2/2;
б) по неявной разностной схеме, взяв h = 0.1, t = 0.02.
Сравнить полученные решения.
Используя метод сеток, найти u(x, t) – решение задачи ut = uxx + (1 + ax)4 – 12 a 2 t(1 +ax)2, 0 < x < 1, 0 < t < 0.1; u(0,t) =t, u(1,t) = t (1 + a)4, 0 < t < 0.1; u(x,0) = 0, 0 < x < 1; a =
Найти решение задачи ut = uxx + (4x/(x 2 + 2t + 1)) ux , 0 < x < 1, 0 < t; u(0,t) = 1/(2t +1), u(1,t) = 1/(2(t + 1)), 0 < t; u(x,0) = 1/(x 2 + 1), 0 < x < 1 в точках (xm, T) xm = 0.1m (m = 1 – 10, T = 0.1) с точностью e = 0.005.
Найти решение задачи ut = uxx + (4x/(x 2 + t + 1)) ux + x/(x 2 + t + 1)2, 0 < x < 1, 0 < t; u(0,t) = 0, u(1,t) = 1/(t + 2), 0 < t; u(x,0) = x/(x 2 + 1), 0 < x < 1 в точках (xm, T) xm = 0.1m (m = 1 – 10, T = 0.4) с точностью e = 0.01.
Используя метод сеток, найти u(x, t) – решение задачи ut = uxx + (-a 2 t + 1)e-ax , 0 < x < 1, 0 < t < 0.01; u(0,t) =t, u(1,t) = t e-a , 0 < t < 0.01; u(x,0) = 0, 0 < x < 1; a= 0.5 k, k = 2 при t = 0.01 с точностью 0.001, если !utt ! < 0.02, !uxxxx ! < 6.
Найти решение задачи ut = (1/2) uxx - 4x 2/(x 2 + t + 1)3, 0 < x < 1, 0 < t; u(0,t) = 1/(t + 1), u(1,t) = 1/(t + 2), 0 < t; u(x,0) = 1/(x 2 + 1), 0 < x < 1 в точках (xm, T) xm = 0.1m (m = 1 – 10, T = 0.3) с точностью e = 0.005.
Найти для x = 0.1m, m=0 – 10, t=0.02 решение уравнения ut = (1/a) uxx , 0 < x < 1, 0 < t < 0.02, удовлетворяющее условиям u(0,t) =e at , u(1,t) = e a ( t -1), 0 < t < 0.02; u(x,0) = e- ax , 0 < x < 1; a = 2 + 0.3 k, k = 4:
а) по явной разностной схеме, взяв h = 0.1, t < ah2/2;
б) по неявной разностной схеме, взяв h = 0.1, t = 0.02.
Сравнить полученные решения.
Найти решение задачи ut = uxx + (4x/(x 2 + t + 1)) ux, 0 < x < 1, 0 < t; u(0,t) = 1/(2t + 1), u(1,t) = 1/(2(t + 1)), 0 < t; u(x,0) = 1/(x 2 + 0.1), 0 < x < 1 в точках (xm, T) xm = 0.1m (m = 1 – 10, T = 0.2) с точностью e = 0.005.
Найти решение задачи ut = (1/2) uxx + 4x 2/(x 2 + t + 1)3, 0 < x < 1, 0 < t; u(0,t) = 1/(t + 1), u(1,t) = 1/(t + 2), 0 < t; u(x,0) = 1/(x 2 + 1), 0 < x < 1 в точках (xm, T) xm = 0.1m (m = 1 – 10, T = 0.2) с точностью e = 0.005.
Используя метод сеток, найти u(x, t) – решение задачи ut = uxx + (-a 2 t + 1)e-ax , 0 < x < 1, 0 < t < 0.01; u(0,t) =t, u(1,t) = t e-a , 0 < t < 0.01; u(x,0) = 0, 0 < x < 1; a = 0.5 k, k = 4 при t = 0.01 с точностью 0.001, если !utt ! < 0.02, !uxxxx ! < 6.
Найти решение задачи ut = uxx + (4x/(x 2 + t + 1)) ux + x/(x 2 + t + 1)2, 0 < x < 1, 0 < t; u(0,t) = 0, u(1,t) = 1/(t + 2), 0 < t; u(x,0) = x/(x 2 + 1), 0 < x < 1 в точках (xm, T) xm = 0.1m (m = 1 – 10, T = 0.3) с точностью e = 0.01.
Используя метод сеток, найти u(x, t) – решение задачи ut = uxx + (1 + ax)4 – 12 a 2 t (1 +ax)2, 0 < x < 1, 0 < t < 0.1; u(0,t) =t, u(1,t) = t (1 + a)4, 0 < t < 0.1; u(x,0) = 0, 0 < x < 1; a = 1 + 0.4 k, k = 4 при t = 0.01 с точностью 0.001, если !utt ! < 0.02, !uxxxx ! < 6.
Найти для x = 0.1m, m=0 – 10, t=0.02 решение уравнения ut = (1/a) uxx , 0 < x < 1, 0 < t <0.02, удовлетворяющее условиям u(0,t) =e at , u(1,t) = e a (t -1), 0 < t < 0.02, u(x,0) = e- ax , 0 < x < 1; a = 2 + 0.3 k, k = - 2:
а)по явной разностной схеме, взяв h = 0.1, t < ah2/2;
б) по неявной разностной схеме, взяв h = 0.1, t = 0.02.
Сравнить полученные решения.
Найти решение задачи ut = (1/2) uxx + 4x 2/(x 2 + t + 1)3, 0 < x < 1, 0 < t; u(0,t) = 1/(t + 1), u(1,t) = 1/(t + 2), 0 < t; u(x,0) = 1/(x 2 + 1), 0 < x < 1 в точках (xm, T) xm = 0.1m (m = 1 – 10, T = 0.1) с точностью e = 0.005.
Используя метод сеток, найти u(x, t) – решение задачи ut = uxx + (-a 2 t + 1)e-ax , 0 < x < 1, 0 < t < 0.01; u(0,t) =t, u(1,t) = t e-a , 0 < t < 0.01; u(x,0) = 0, 0 < x < 1; a = 0.5 k, k = 6 при t = 0.01 с точностью 0.001, если !utt ! < 0.02, !uxxxx ! < 6.
Найти решение задачи ut = uxx + (4x/(x 2 + t + 1)) ux , 0 < x < 1, 0 < t; u(0,t) = 1/(2t + 1), u(1,t) = 1/(2(t + 1)), 0 < t; u(x,0) = 1/(x 2 + 0.1), 0 < x < 1 в точках (xm, T) xm = 0.1m (m = 1 – 10, T = 0.3) с точностью e = 0.005.
Используя метод сеток, найти u(x, t) – решение задачи ut = uxx + (1 + ax)4 – 12 a 2 t(1 +ax)2, 0 < x < 1, 0 < t < 0.1; u(0,t) =t, u(1,t) = t (1 + a)4, 0 < t < 0.1; u(x,0) = 0, 0 < x < 1; a = 1 + 0.4 k, k = 3 при t = 0.01 с точностью 0.001, если !utt ! < 0.02, !uxxxx ! < 6.
Найти решение задачи ut = uxx + (4x/(x 2 + t + 1)) ux + x/(x 2 + t + 1)2, 0 < x < 1, 0 < t; u(0,t) = 0, u(1,t) = 1/(t + 2), 0 < t; u(x,0) = x/(x 2 + 1), 0 < x < 1 в точках (xm, T) xm = 0.1m (m = 1 – 10, T = 0.2) с точностью e = 0.01.
Найти для x = 0.1m, m=0 – 10, t=0.02 решение уравнения ut = (1/a) uxx , 0 < x < 1, 0 < t <0.02, удовлетворяющее условиям u(0,t) =e at , u(1,t) = e a (t -1), 0 < t < 0.02; u(x,0) = e- ax , 0 < x < 1; a = 2 + 0.3 k, k = 3:
а) по явной разностной схеме, взяв h = 0.1, t < ah2/2;
б) по неявной разностной схеме, взяв h = 0.1, t = 0.02.
Сравнить полученные решения.
Используя метод сеток, найти u(x, t) – решение задачи ut = uxx + (-a 2 t + 1)e-ax , 0 < x < 1, 0 < t < 0.01; u(0,t) =t, u(1,t) = t e-a , 0 < t < 0.01; u(x,0) = 0, 0 < x < 1; a = 0.5 k, k = 1 при t = 0.01 с точностью 0.001, если ½utt ê < 0.02, ½uxxxx ê < 6.