- •Методические указания к лабораторным работам по дисциплине "Вычислительная математика "
- •Лабораторная работа №1
- •(Количество часов -2 час. )
- •Определение псевдообратной матрицы для произвольной матрицына основе метода ортогонализации Грамма-Шмидта (гшо)
- •Построение псевдообратной матрицы на основе метода ортогонализации Грамма-Шмидта.
- •1.1. Сокращение интервала неопределенности методом золотого сечения
- •1.2. Сокращение интервала неопределенности методом квадратичной аппроксимации
- •1.3. Минимизация многомерной функции при наличии линейных ограничений на основе метода Давидона-Флетчера-Пауэлла
- •Порядок выполнения работы
- •2. Учебный пример выполнения лабораторной работы
- •2.1. Минимизация многомерных функций
- •Содержание отчета
- •4.2. Быстрое преобразование Уолша
- •4.3. Быстрое преобразование Хаара
- •Порядок выполнения работы
- •Порядок выполнения работы
- •3. Метод Рунге-Кутта 4-го порядка
- •Порядок выполнения работы
- •Порядок выполнения работы
- •Построение графика функции одного аргумента
- •Просмотр участков двумерных графиков
- •Порядок выполнения работы
- •1.2. Абсолютные и относительные погрешности
- •1.5. Погрешности произведения и частного приближенных чисел
- •3. Содержание отчета по лабораторной работе
- •Контрольные вопросы
- •Лабораторная работа № 10 решение нелинейных уравнений
- •1. Теоретическая часть
- •1.1. Метод половинного деления
- •1.2. Метод хорд
- •1.3. Метод Ньютона (касательных)
- •1.4. Метод простой итерации
- •2. Порядок выполнения работы
- •3. Содержание отчета по лабораторной работе
- •4. Контрольные вопросы
- •5. Исходные данные
- •Интерполирование функций
- •1. Теоретическая часть
- •1.1. Формула Лагранжа
- •1.2. Интерполирование по схеме Эйткена
- •1.4. Формула Ньютона с разделенными разностями
- •Интерполяционный многочлен Ньютона, использующий разделенные разности, имеет вид:
- •1.6. Интерполяция сплайнами
- •2. Порядок выполнения работы
- •3. Содержание отчета по лабораторной работе
- •4. Контрольные вопросы
- •5. Исходные данные
4.2. Быстрое преобразование Уолша
Функции
Уолша обладают интересными свойствами,
привлекающими к ним все большее внимание.
Они принимают всего два значения {+1 или
-1} и потому удобны для вычислений на
ЭВМ. Существует различное упорядочение
функций Уолша .Рассмотрим одну из
возможных систем функций Уолша – систему
Уолша-Адамара. Элементарная матрица
Адамара, состоящая из одного элемента
(N=1),
имеет вид H1=[1].
Для N=2
элементами матрицы Адамара будут
элементарные матрицы (
H1):
![]()
.
Для
N
= 4 элементами матрицы Адамара будут
матрицы (
H2):
![]()
.
Для
N=2n
имеем
.
Дискретное преобразование Уолша, как и дискретное преобразование Фурье, представляется в матричной форме:
FY=HN*S/N; S= (HN)T*FY.
Здесь FY-коэффициенты спектрального разложения по функциям Уолша, S-дискретные временные отсчеты сигнала.
Быстрое преобразование Уолша (БПУ) можно получить из формулы для БПФ (10),исключив фазосдвигающие составляющие.
Матрица преобразований упрощается и имеет вид:
.
Коэффициенты спектрального разложения по функциям Уолша имеют прямую последовательность в отличие от БПФ, имеющей последовательность с двоично-инверсными номерами.
4.3. Быстрое преобразование Хаара
Преобразование
Хаара основывается на использовании
ортогональной матрицы Хаара. Приведем
пример матрицы Хаара четвертого порядка:
и матрицы
Хаара восьмого порядка:
![]()
-
1
1
1
1
1
1
1
1
1
1
1
1
-1
-1
-1
-1
21/2
21/2
-21/2
-21/2
0
0
0
0
0
0
0
0
21/2
21/2
-21/2
-21/2
2
-2
0
0
0
0
0
0
0
0
2
-2
0
0
0
0
0
0
0
0
2
-2
0
0
0
0
0
0
0
0
2
-2
Для записи алгоритма БПХ введем следующие обозначения матриц:
U02=[1
1]; U12=[1
-1]; I2=
;
,
где ,
∑ - знаки прямой суммы.
.
Таким образом, матрица Хаара представлена в виде произведения «n» слабо заполненных матриц.
Быстрое преобразование Хаара - самое быстрое из используемых ортогональных функций.
Дискретное преобразование Хаара, как и дискретное преобразование Фурье, представляется в матричной форме:
;
.
.
Здесь использовано правило для транспонирования произведения матриц: [M1*M2*M3]Т=[М3Т*М2Т*М1Т].
Коэффициенты спектрального разложения по функциям Хаара имеют прямую последовательность в отличие от БПФ, имеющей последовательность с двоично-инверсными номерами.
