
- •Математичні моделі в сапр: методичні вказівки до лабораторних робіт
- •Передмова до методичних вказівок
- •1.Лабораторна робота EvalExpr на тему "Знаходження похибок обчислень"
- •Теоретичні питання.
- •Приклад виконання роботи EvalExpr
- •Варіанти індивідуальних завдань до роботи EvalExpr
- •2.Лабораторна робота DetKram на тему “Розв’язання слр за формулами Крамера”
- •Теоретичні питання
- •Приклад виконання роботи DetKram
- •Варіанти індивідуальних завдань до роботи DetKram
- •Допоміжні матеріали до роботи DetKram Інструкція користувачу до додатку InteractiveDet
- •Приклад змісту вхідного файла Det.Inp:
- •Вихідні файли програми.
- •Натискання клавіші “Введення з файлу”.
- •Натискання клавіші “Зменшити у рядку”.
- •Натискання клавіші “Зменшити у стовпці”.
- •Натискання клавіші “Зменшити розміри”.
- •Завершення роботи додатку.
- •3.Лабораторна робота InvMatr на тему “Розв’язання слр за методом оберненої матриці”
- •Теоретичні питання
- •Приклад виконання роботи InvMatr
- •Варіанти індивідуальних завдань до роботи InvMatr
- •Допоміжні матеріали до роботи InvMatr Програма AdjMatr
- •Інструкція користувачу до додатку InteractiveAdjMatr
- •4.Лабораторна робота Gauss на тему „Наближене розв`язання слр методами виключення невідомих”
- •Теоретичні питання
- •Приклад виконання роботи Gauss
- •2В. Обчислення компонент вектора нев’язки.
- •3A. Знаходження наближеної оберненої матриці за методом Гауса.
- •3B. Знаходження компонент наближеного розв’язку .
- •3C. Знаходження компонент вектору нев’язки
- •5. Знаходження міри обумовленості головної матриці системи.
- •Варіанти індивідуальних завдань до роботи Gauss
- •Додаткові матеріали до роботи Gauss Програма NormVect знаходження різних норм вектора
- •Програма NormMatr знаходження різних норм матриці
- •Процедура g3v33l4
- •Модуль GaussRow
- •Iнструкцiя користувачу до програми InteractiveGauss
- •5.Лабораторна робота LinSysIt на тему “Ітераційні методи розв’язання слр”
- •Теоретичні питання
- •Приклад виконання роботи LinSysIt
- •Варіанти індивідуальних завдань до роботи LinSysIt
- •Допоміжні матеріали до роботи LinSysIt Програма NormForm
- •Програма SysIter
- •Програма Zeidel
- •6.Лабораторна робота TransEq на тему "Розв'язання трансцендентних рівнянь з одним невідомим"
- •Теоретичні питання
- •Приклад виконання роботи TransEq
- •Розрахункова таблиця знаходження кореня
- •Розрахункова таблиця знаходження кореня
- •Розрахункова таблиця знаходження кореня
- •Уточнення відокремленого кореня за методом простих ітерацій.
- •Розрахункова таблиця знаходження кореня рівняння за методом простих ітерацій
- •Розрахункова таблиця знаходження кореня рівняння за методом простих ітерацій
- •Варіанти індивідуальних завдань до роботи TransEq
- •Додаткові матеріали до роботи TransEq Інструкція користувачу для роботи з додатком TranscEq
- •Робота з додатком TranscEq
- •7.Лабораторна робота NonLinSys на тему „Розв’язання систем нелінійних рівнянь”
- •Теоретичні питання
- •Приклад виконання роботи NonLinSys
- •Розрахункова таблиця знаходження роз'вязку
- •Варіанти індивідуальних завдань до роботи NonLinSys
- •Допоміжні матеріали до роботи NonLinSys Інструкція користувачу для роботи з програмою Newton2
- •Текст програми Newton2
- •Текст модуля Funct
- •Текст модуля NewtonS
- •8.Лабораторна робота SelfEq на тему ”Знаходження характеристичного рівняння”
- •Теоретичні питання
- •Приклад виконання роботи SelfEq
- •Варіанти індивідуальних завдань до роботи SelfEq
- •9.Лабораторна робота SelfVect на тему “Знаходження власних векторів”
- •Теоретичні питання.
- •Приклад виконання роботи SelfVect
- •Варіанти індивідуальних завдань до роботи SelfVect
- •Література
Допоміжні матеріали до роботи NonLinSys Інструкція користувачу для роботи з програмою Newton2
Тут наведено текст
програми Newton2.pas
для уточнення
відокремленого розв’язку
системи двох нелінійних рівнянь
Вихідний файл Newton2.inp. Вихідний файл повинен бути у поточної теці (в якої розміщено тексти програми Newton2.pas і модуля Funct.pas). Структура вхідного файлу.
----------------------- початок ----------------------------
0.2 : x0 - координата x початкового наближення
0.9 : y0 - координата y початкового наближення
4 : n - число ітерацій
9 : All - число всіх позицій для друкування числа
6 : dec - число цифр після коми для друкування числа
------------------------- кінець ----------------------------
Робота з заготівкою модуля Funct.
Перед початком роботи потрібно підготувати записи функцій
,
,
і
,
як рядків програми на мові Turbo Pascal у вигляді функція:= вираз функцїї (див. відповідні рядкі у тексті модуля Funct);
Скопіювати повністю теку Newton2 (під тим же ім’ям) на власний магнітний носій (комп’ютер).
Завантажити Turbo Pascal, знайти в теці Newton2.pas і завантажити текст головної програми Newton2.pas у текстовий редактор Turbo Pascal.
Відкрити для змін вихідний текстовий файл Newton2.inp і встановити значення вхідних параметрів (кожен у окремому рядку з прогалиною після параметру, а далі до кінця рядка будь-яка послідовність символів). Наприклад:
------------------------ початок -----------------------
0.8 х0- координата х початкового наближення,
1.5 у0- координата у початкового наближення,
5 n - число ітерацій,
9 all- число всіх позицій для друкування числа,
8 dec- число цифр після коми для друкування числа.
------------------------------ кінець -----------------------------
Відкрити для зміни текстовий файл Funct.pas і замінити відповідні тіла функцій F(x,y), Fx’(x,y), Fy’(x,y) i G(x,y), Gx’(x,y), Gy’(x,y) на підготовлені раніше вирази. Скомпілювати замінений модуль.
Скомпілювати і виконати програму Newton2.pas.
Результати роботи програми будуть знаходитись у текстовому файлі Newton2.out.
Результуючий
файл Newton2.out.
У першому
рядку наведені початкові дані (див. опис
вхідного файлу). Далі йдуть рядки
ітераційних кроків з номером кроку
,
попереднього ітераційного наближення
(у дужках)
плюс вектор ітераційного кроку
(у дужках). Підраховане нове ітераційне
наближення
міститься у наступному рядку (як попереднє
ітераційне наближення для наступного
кроку). Наведемо приклад вихідного
файлу.
---------------------- початок -----------------------------
S<Newton2>
x0= 0.800 y0= 1.500 n=5 all=9 dec=6
1:{ 0.800000 1.500000} + {-0.012964 -0.087106}
2:{ 0.787036 1.412894} + { 0.002920 -0.005103}
3:{ 0.789956 1.407791} + { 0.000025 -0.000014}
4:{ 0.789981 1.407777} + { 0.000000 -0.000000}
5:{ 0.789981 1.407777} + { 0.000000 -0.000000}
F<Newton2>
--------------------- кінець ----------------------------