- •Конспект лекцій Частина і з дисципліни “Числові методи і моделювання на еом”
- •Лекция 1 числові методи алгебри. Особливості алгоритмування обчислювальних задач. Елементи теорії похибок обчислень та аналізу помилок округлення. Порядок виконання операцій
- •1.1. Про наближені обчислення
- •1.2. Лінійні заміни змінних
- •1.3. Системи лінійних алгебраїчних рівнянь
- •2.1. Апроксимація функції по Фур'є.
- •2.1.1. Перетворення Фур'є
- •2.2. Зворотна матриця
- •3.1. Метод ділення відрізка навпіл для розв'язання рівнянь
- •3.2. Метод хорд для рішення рівнянь
- •3.3. Метод дотичних для розв'язання рівнянь
- •3.4. Методика рішення алгебраїчного рівняння
- •Метод простих ітерацій
- •Метод Зейделя
- •Метод ітерацій для рішення рівнянь
- •4.4. Метод ітерацій для рішення систем нелінійних алгебраїчних і
- •Лекция 5 звернення матриць. Подвійність у лінійному програмуванні. Одночасне рішення пари двоїстих задач лінійного програмування.
- •Лекція 6
- •6.1. Чисельне диференціювання функції однієї змінної.
- •6.2. Чисельне інтегрування функції однієї змінної.
- •6. 3. Постановка задачі про чисельне рішення звичайного диференціального рівняння.
- •6.5. Метод Рунге-Кутта чисельного рішення звичайного диференціального рівняння.
- •6.6. Підхід до чисельного рішення системи звичайних диференціальних
- •Лекция 7 методи розв’язку диференціальних рівнянь та їх систем. Розв'язання систем лінійних алгебричних рівнянь із допомогою жорданових виключень
- •Лекция 8 чисельне диференціювання та інтегрування. Основна задача лінійного програмування. Дослідження її окремих випадків. Модифікований варіант жордановых винятків
- •8.1. Постановка основної задачі лінійного програмування (озлп)
- •8.2. Екстремальні задачі, що зводяться до озлп заміною змінних
- •8.3. Лінійна заміна змінних і її використання в дослідженні основної
- •8.4. Модифікований варіант жордановых виключень як спосіб організації лінійної заміни змінних
- •Лекция 9 диференціювання інтерполяційних формул. Мова « n-мірних» точок. Геометрія задачі лінійного програмування. Опорне рішення й оптимальне рішення. Загальні установки симплекса-методу
- •9.1.Мова n-мірних точок.
- •9.2. Геометрія задачі лінійного програмування.
- •Опорне рішення й оптимальне рішення. Загальні установки симплекс-методу
- •Підготовка озлп до рішення симплекс-методом.
- •Список рекомендованої літератури
2.1.1. Перетворення Фур'є
Так називається дія, за допомогою якого по заданій в інтервалі функції будується система чисел (2.1.2). За традицією, саме ці числа також позначають словами «перетворення Фур'є» даної функції. При цьому числа називаються косинус- перетворенням Фур'є функції , а числа називаються синусом-перетворенням Фур'є функції . У перетворення Фур'є функції є множина властивостей і застосувань; зокрема,
відомий ряд прикладних питань, у яких ту або іншу інформацію про функцію вдається одержати тільки через її перетворення Фур'є, що у таких ситуаціях вдається одержати деякими непрямими засобами.
Розглянемо наступний окремий випадок. Функція розглядається на інтервалі
(-p ,p ) і притім тільки в його окремих точках
при деякому заздалегідь заданому й фіксованому числі . Значення функції в цих крапках уважаються відомими; позначимо
.
У рівності
=
покладемо . Одержимо
(2.2.1) .
Проаналізуємо співвідношення (2.2.1). Якщо довільне ціле ненегативне число розділити з остачею на число , то вийде співвідношення , де для цілих є лише наступні можливості:
.
З урахуванням періодичності косинуса й синуса у вираженні (2.2.1) можна привести подібні члени, у результаті чого вийде:
(2.2.2) ,
де
,
Відзначимо, що тепер всі суми в (2.2.2) - кінцеві. Відомий наступний факт про тригонометричні суми:
для всіх чисел мають місце рівності
Якщо обидві частини співвідношення (2.2.2) помножити на й потім просумировать по , те, з обліком тільки що сказаного, легко одержати, що
(2.2.3а) ;
а якщо обидві частини (7.2.2) помножити на й, з обліком того ж твердження про підсумовування косинусів і синусів, одержимо співвідношення
(2.2.3б)
причому в (2.2.3а) і (2.2.3б) . Числа , називаються дискретним перетворенням Фур'є функції . Якщо в рівності (2.2.2) замінити на довільний , то воно з точного стане наближеним. Його праву частину в цьому випадку називають тригонометричною інтерполяцією функції .
2.2. Зворотна матриця
Нехай А и В - дві матриці наступного виду:
і ;
оборотний увага на те, що в першій з них стовпців стільки ж, скільки в другий з них - рядків. У цьому випадку можна побудувати нову матрицю З
,
називану добутком А на В (пишуть ІЗ=АВ), за правилом
,
де i=1,...,m і j=1,...,r. Серед матриць прийнято виділяти й позначати буквою E
наступну матрицю, називану одиничної, -
;
пояснимо: одинична матриця завжди квадратна, по її діагоналі коштує 1, а всі інші елементи - нулі, тобто якщо Е=(еij), i,j=1,...,n, те
.
Говорять, що квадратна матриця А и квадратна матриця В того ж розміру є зворотними по відношенню друг до друга, якщо виконано хоч одне з матричних рівностей - або АВ=Е або ВА=Е. Можна перевірити, що якщо одне із цих матричних рівностей виконується, то виконується й інше. Про всякий випадок пояснимо, що слова «матрична рівність» означають рівність заелементне, тобто якщо має місце рівність матриць U=V, те це означає, що в них рівні кількості рядків і рівні кількості стовпців і при цьому uij=vij для всіх i і j. Зокрема, система лінійних алгебраїчних рівнянь із п.3 Лекції 1 може бути записана у вигляді AX=B, де
,
а А и В мають той же змив, що й у п.3 Лекції 1.
Неважко помітити, що якщо для квадратної матриці А потрібно знайти зворотну матрицю В, те треба вирішити n систем лінійних алгебраїчних рівнянь із однієї й тої ж матрицею: запишемо матричну рівність АВ=Е поелементно (тут усього вийде n2 рівностей), причому будемо послідовно дорівнювати стовпці - перший стовпець лівої частини до першого стовпця правої частини, другий стовпець лівої частини до другого стовпця правої частини й т.д. –
Із цього спостереження треба, що знайти зворотну матрицю можна вирішуючи одночасно n систем лінійних алгебраїчних рівнянь із однієї й тією же матрицею А.
Лекция 3
АЛГОРИТМИ ТА ПРОГРАМИ РОЗРАХУНКІВ ЗА ФОРМУЛАМИ ТА ТАБЛИЦЯМИ. ВИКОРИСТАННЯ ЛАНЦЮГОВИХ ДРОБІВ. СХЕМА ГОРНЕРА. МЕТОДИ ДІЛЕННЯ ВІДРІЗКА НАВПІЛ, ХОРД І ДОТИЧНИХ ДЛЯ УТОЧНЕННЯ КОРЕНЯ РІВНЯННЯ. ВІДШУКАННЯ РЕЧОВИННИХ КОРНІЙ АЛГЕБРАЇЧНИХ РІВНЯНЬ.