Що залишився після т кроків виключення стовпець (або стовпці) визначає невязки від істинного рішення (які повинні бути рівні нулю при т п) тільки теоретично. Обчислені невязки через помилки округлення не є істинними невязке, однак вони мають необхідний порядок величини. Якщо дійсні значення векторів невязки досить великі, то залишився стовпець (або стовпці) може служити для них гарною апроксимацією. [47]
Позначимо через G (x) замкнуту область допустимих значень. Тоді потрібно знайти кусочно-безперервну функцію 5 (х), що мінімізувала функціонал (14.1) при відповідних обмеженнях, пов'язаних з граничними станами по стійкості положення (проти спливання), міцності, деформативності і загальної поздовжньої стійкості. Обмеження використовуються у вигляді вектора нев'язок між розрахунковими параметрами напружено-деформованого стану і гранично допустимими. [48]
Значні труднощі пов'язані з оцінкою і забезпеченням необхідної точності рівнянь і особливо систем рівнянь. Точність рішення рівнянь оцінюють за величиною нев'язки - різниці правої і лівої частин рівнянь при підстановці обчислених значень коренів. Точність рішення систем рівнянь оцінюють по нормі вектора нев'язок рівнянь системи. [49]
Виконанню цієї процедури повинна передувати процедура unsymdet. У тілі процедури unsym асі solve використовується процедура unsymsol. Якщо матриця А погано обумовлена по відношенню до операції звернення, то передбачений вихід з процедури. Необхідно відзначити, що вектори невязок обчислюються з накопиченням скалярного твори або з подвійною точністю. [50]
Ці обурення можуть виникнути під час запису чисел А15 і bj в пам'ять машини або в силу невизначеностей, або властивих цим числам, або виникають через помилки експерименту. Останні невизначеності дуже великі для багатьох реальних задач в порівнянні з машинними помилками округлення. Розв'язана задача породжує подальші обурення через різного роду невідповідностей у машинної арифметики або в алгоритмі рішення машинної завдання. Вектор невязки г представляє загальний ефект цих невідповідностей. [51]
Сторінки: 1 2 3 4