
- •Основы вычислений
- •Типы данных
- •Cтроковые переменные
- •Работа с массивами
- •Принцип программирования в Mathcad
- •Логические операторы
- •Матричные операторы
- •Операторы выражения
- •Функции
- •Элементарные функции
- •Функция вывода текущего времени
- •Вычисление рядов и произведений
- •Вычисление предела
- •Простейшие матричные операции Транспонирование
- •Сложение и вычитание матриц
- •Векторы являются частным случаем, поэтому для них справедливы все те операции, что и для матриц. Но есть и специфические операторы.
- •Определитель квадратной матрицы
- •Ранг матрицы
- •Обращение квадратной матрицы
- •Сортировка элементов матрицы
- •Вывод размера матрицы
- •Нелинейные алгебраические уравнения
- •Численное решение уравнений
- •Системы линейных уравнений
- •Решение слау с использованием вычислительного блока Given/Find
- •Алгоритм исключения Гаусса
- •Произвольные системы линейных уравнений
- •Переопределённые системы
- •Недоопределённые системы
- •Вырожденные системы линейных уравнений
- •Регуляризация
- •Матричные разложения
- •Разложение Холецкого
- •Обыкновенные дифференциальные уравнения
- •Задачи Коши для оду
- •Дифференциальные уравнения n-го порядка
- •Схемы Рунге-Кутта
- •Система n дифференциальных уравнений
- •Решение систем оду в одной заданной точке
- •Жёсткие системы оду
- •Примеры динамических моделей
- •Модель "хищник-жертва"
- •Модель генератора автоколебаний
- •Модель Лоренца
- •Краевые задачи для оду
- •Алгоритм стрельбы
- •Двухточечные краевые задачи
- •Краевые задачи с условием во внутренней точке
Решение слау с использованием вычислительного блока Given/Find
Дана система линейных уравнений. Найти её решение.
3x
+ 2x
+ x
+ 11x
=
25
7x + 4x + 3x + 5x =37
x + 10x + 2x + 19x = 44
6x - 8x + 14x + 15x = 31
Алгоритм выполнения:
Ввести начальные значения корней
Ввести ключевое слово Given
Через логическое равно ввести все уравнения системы\
Записать Find(x)
Установить знак =
Если система содержит большое количество уравнений, удобнее применять матричную форму.
Алгоритм выполнения:
1.Ввести матрицу коэффициентов при неизвестных
2. Ввести вектор свободных членов
3. Ввести начальные значения корней через вектор
4. Ввести ключевое слово Given
5. Записать систему в матричной форму
6. Применить метод Find
7. В вычислительном блоке логическое равно.
Самый простой способ решения почти всякой несингулярной системы является алгоритм исключения Гаусса.
Алгоритм исключения Гаусса
Алгоритм последовательных исключений Гаусса основан на преобразовании матрицы A линейной системы Ax=b к треугольному виду (т.е. к форме, когда все элементы ниже главной диагонали матрицы являются нулевыми). Точнее, СЛАУ Ax=b заменяется эквивалентной системой с другой матрицей A* и другим вектором правых частей b*, но имеющей то же решение, что и исходная система.
Алгоритм состоит в следующем:
1. Проводится прямой ход исключения неизвестных путем подстановки одних уравнений в другие. Используются следующие формулы: aij = - aij / ajj aik = aik + aijaik bj = bj + aijbj j = 1,2, ... N-1 i=j+1, j+2, ... N k=j+1, j+2, ... j+N
2. При помощи обратного хода определяются все неизвестные х. Для этого хN сразу определяется из последнего уравнения, в которое не входят другие х (матрица системы теперь является трехдиагональной). Затем хN подставляется в предыдущее уравнение, из которого сразу определяется хN-1 и т.д.
Этот метод реализован во встроенной функции lsolve.
При использовании этого метода система уравнений должна быть записана в матричной форме.
Алгоритм выполнения:
Ввести матрицу с коэффициентами при неизвестных
Ввести столбец свободных членов
Записать функция lsolve
При численном решении после функции ставится знак =, а
п ри символьном
М
ожно
провести проверку:
Произвольные системы линейных уравнений
Классические задачи решения систем линейных уравнений предполагают равное количество уравнений и неизвестных . Рассмотрим теперь системы, в которых матрица коэффициентов при неизвестных является неквадратной либо плохо обусловлена.