
- •Основы вычислительной математики Учебно-методическое пособие
- •Оглавление
- •Предисловие
- •I. Введение
- •Контрольные вопросы
- •II. Системы линейных алгебраических уравнений
- •II.1. Основные понятия и определения
- •II.2. Прямые методы решения систем линейных уравнений
- •II.3. Метод Гаусса решения слау
- •II.4. Решение слау с помощью -разложения
- •II.5. Метод итераций решения слау
- •II.6. Метод Зейделя решения слау
- •II.7. Сходимость итерационных процессов для систем линейных уравнений
- •III. Методы решения нелинейных скалярных уравнений
- •III.1. Локализация корней
- •III.2. Метод простых итераций
- •III.3. Приведение нелинейного уравнения к виду , допускающему сходящиеся итерации
- •III.4. Метод Ньютона (Метод касательных)
- •III.5. Контрольные вопросы
- •III.6. Задачи для самостоятельного решения
- •IV. Приближенное решение систем нелинейных уравнений
- •IV.1. Метод Ньютона
- •IV.2. Метод итераций
- •IV.3. Контрольные вопросы
- •IV.4. Задачи для самостоятельного решения
- •V. Интерполяция и аппроксимация
- •V.1. Интерполяция
- •V.2. Интерполяционный многочлен в форме Лагранжа
- •V.3. Интерполяционный многочлен в форме Ньютона
- •V.4. Приближение методом наименьших квадратов
- •V.5. Контрольные вопросы
- •V.6. Задания для самостоятельного решения
- •VI. Равномерное приближение функций
- •V.1. Полиномы Чебышева
- •VI.2. Понятие о равномерном приближении функций
- •VI.3. Контрольные вопросы
- •VII. Приближенное дифференцирование
- •VII.1. Постановка вопроса
- •VII.2. Формулы приближенного дифференцирования, основанные на первой интерполяционной формуле Ньютона
- •VII.3. Формулы приближенного дифференцирования, основные на формуле Стирлинга
- •VII.4. Контрольные вопросы
- •VII.5. Задачи для самостоятельного решения
- •VIII. Численное интегрирование функций
- •VIII.1. Общие замечания
- •VIII.2.Квадратурные формулы Ньютона-Котеса
- •VIII.3. Формула трапеций и ее остаточный член
- •Общая формула трапеций (правило трапеций).
- •VIII.4. Формула Симпсона и ее остаточный член
- •Общая формула Симпсона (параболическая формула).
- •VIII.5. Формулы Ньютона-Котеса высших порядков
- •VIII.6. Контрольные вопросы
- •VIII.7. Задачи для самостоятельного решения
- •IX. Приближенное решение обыкновенных дифференциальных уравнений
- •IX.1. Интегрирование дифференциальных уравнений с помощью степенных рядов
- •IX.2. Метод последовательных приближений
- •IX.3. Метод Эйлера
- •IX.4. Метод Эйлера с итерациями
- •IX.5. Метод Рунге-Кутта
- •IX.6. Контрольные вопросы
- •IX.7. Задачи для самостоятельного решения
- •X. Преобразование фурье
- •X.1. Разложение периодических функций в ряд Фурье
- •X.2. Эффект Гиббса
- •X.3. Контрольные вопросы
- •IX.4. Задачи для самостоятельного решения
- •Список литературы
X.2. Эффект Гиббса
Эффект Гиббса возникает при решении задачи синтеза для разрывных функций. Изучим основные его проявления на примере прямоугольной импульсной функции
которая
имеет скачек, равный 1 в точке разрыва
.
Формальное разложение в ряд Фурье
находится подстановкой функции
в формулы (10.8), вычислением интегралов,
оказывающихся равными
и подстановкой выражений для коэффициентов в (10.9):
Использование конечного числа членов ряда приводит к тому, что частичные суммы ряда, являются функциями, в которых присутствуют периодические составляющие. Их период равен периоду последнего удержанного члена или первого оставленного.
К. Ланцош предложил способ уменьшения эффекта пульсаций за счет сглаживания усеченного ряда
(8.18)
интегрированием (усреднением) по периоду последнего оставленного члена.
При
выборе в качестве интервала сглаживания
периода последнего члена усеченного
ряда Фурье сглаженное значение
получим как среднее от
:
(10.19)
На практике берется как номер первого отбрасываемого члена ряда.
Подставляя (10.18) в (10.19), получим:
Применяя известные тригонометрические формулы, окончательно находим:
(10.20)
где
– так называемые сигма-факторы:
(10.21)
Таким
образом, сглаженный ряд Фурье есть
исходный ряд Фурье с коэффициентами,
умноженными на соответствующие
сигма-факторы. Отметим, что эффект
сглаживания достигается за счет того,
что при
сигма-фактор
-го
члена
оказывается равным нулю.
С
физической точки зрения формулу (10.20)
можно трактовать как наблюдение исходной
функции
через узкое просвечивающее прямоугольное
окно шириной
.
Наблюдаемая функция
есть средняя интенсивность света от
исходной функции
.
Она оказывается более яркой, когда
функция
больше, и менее яркой, когда функция
меньше.
X.3. Контрольные вопросы
Ортогональные функции.
Ряд Фурье для периодической Функции.
Различные формы записи для ряда Фурье.
Эффект Гиббса.
IX.4. Задачи для самостоятельного решения
Разложить заданную периодическую функцию в ряд Фурье. Построить на одном чертеже график исходной функции и ее аппроксимирующих полиномов .
1.
2.
3.
4.
5.
6.
7.
8.
9.
10.
11.
12.
13.
14.
15.
Список литературы
Пирумов У.Г. Численные методы: Учеб. Пособие для студ. втузов. – М.: Дрофа, 2003. – 224 с.
Вержбицкий В.М. Основы численных методов: Учебник для вузов. – М.: Высш. шк., 2002. – 840 с.
Поршнев С.В. Вычислительная математика. Курс лекций. – СПб.: БХВ-Петербург, 2004. – 320 с.
Бахвалов Н.С., Жидков Н.П., Кобельков Г.М. Численные методы. – М.: Лаборатория Базовых Знаний, 2002. – 632 с.
Демидович Б.П., Марон И.А., Шувалова Э.З. Численные методы анализа: Учеб. Пособие для студ. втузов. – М.: Изд-во «Наука», 1967. – 368 с.
Демидович Б.П., Марон И.А. Основы вычислительной математики: Учеб. Пособие для студ. втузов. – М.: Изд-во «Наука», 1970. – 664 с.