Курсовая_ИНФОРМАТИКА (МАШ, КТО, ПТМК)-2015

12

Рис. 8. Вычисление матрицы D в Excel

Операции линейной алгебры в EXEL выполняются в следующей последовательности:

1.Выделяется диапазон ячеек (матрица), в котором будет размещен результат. Также в Excel для обозначения матрицы иногда используется термин – «массив».

2.Применяется нужная процедура из категории «математические»:

«МУМНОЖ» - возвращает произведение матриц (матрицы хранятся в массивах). Результатом является массив с таким же числом строк, как массив №1, и с таким же числом столбцов, как массив №2. Обращаем внимание, что массивы №1 и №2 должны задаваться именно в том порядке (слева направо), как они записаны в формуле. Нарушение этого правила приведет не только к неправильному вычислению коэффициентов результирующей матрицы, но и к получению результирующей матрицы другого размера (кроме случая перемножения квадратных матриц).

«ТРАНСП» - транспортирование матрицы (возвращает вертикальный диапазон ячеек в виде горизонтального и наоборот).

«МОБР» - возвращает обратную матрицу для матрицы, хранящейся в массиве.

«МОПРЕД» - возвращает определитель матрицы (матрица хранится в массиве).

3.После появления в массиве результатов первого значения (соответствующего индексам 1,1) для появления остальных

15

Еще один пример приведен на рис. 12

Рис. 12. Решение системы линейных алгебраических уравнений с проверками выполнения отдельных шагов

Здесь реализованы несколько проверок: проверка существования решения (вычисление определителя (детерминанта матрицы А); проверка правильности получения обратной матрицы (перемножение исходной матрицы на обратную должно давать единичную матрицу); проверка правильности вычисления искомого вектора (перемножение исходной матрицы на искомый вектор должно давать вектор правых частей). Предпочтительно выполнение задания с матрицей D именно в таком (более полном) объеме со всеми перечисленными проверками.

Реализованная последовательность решения системы уравнений относится к прямым методам, таким как метод Гаусса. Но системы линейных алгебраических уравнений можно решать и

17

Решение системы уравнений в Mathcad приведено на рис. 10б. В Mathcad также можно при менять прямые (метод Гаусса) и итерационные методы решения систем алгебраических уравнений.

Вкачестве дополнительного задания (задание №6) может быть реализована подобная процедура для заданного варианта системы уравнений с матрицей D (с применением как Excel, так и Mathcad).

3.4.Приближение таблично заданной функции

Винженерной практике довольно часто встречается задача получения зависимости (формулы) между двумя величинами в том случае, когда имеется таблица значений (например, в результате

экспериментальных исследований получены значения величины Yi при некоторых конкретных значениях аргумента Xi, где i=1,2,…m). Но для последующей деятельности требуется формула, связывающая эти две величины, по которой можно было бы с некоторой точностью определить величину Y при любом значении величины X. Получение такой зависимости называется аппроксимацией таблично (точечно) заданной функции.

Частным случаем аппроксимации является интерполяция полиномом (многочленом) n-й степени. Из курса математики известно, что имея m точек, можно получить неизвестные коэффициенты полинома n=m-1 степени. Интерполяцией называется получение такой полиномиальной зависимости Y ( X ) , в соответствии

скоторой для всех пар точек (называемых узлами интерполяции)

имеется совпадение табличных с рассчитанными по зависимости Y ( X ) . Интерполяционная зависимость может быть построена единая для всего интервала определения таблично заданной функции. В этом случае степень полинома на единицу ниже количества несовпадающих узлов интерполяции. Можно построить кусочную интерполяцию, разделив все количество узлов на группы. В пределах каждого такого участка интерполяции получают полином степенью на единицу ниже количества узлов интерполяции в этой группе. Если в группу включать 2 узла, то получается полином первой степени (прямая линия). Такая интерполяция называется кусочно-линейной.

Кроме построения интерполяционных зависимостей, можно строить аппроксимирующие зависимости на основе различных функций взаимосвязи двух рассматриваемых величин. Обычно в этом

18

случае количество неизвестных параметров выбранной функции намного ниже количества пар значений в исходной таблице Y ( X ) . Очевидно, что в этом случае возможно бесконечно большое количество сочетаний конкретных значений параметров аппроксимирующей функции. Необходимо выбирать такие значения параметров, которые в максимальной степени обеспечивали бы близость аппроксимирующей функции к исходным табличным значениям. Обычно применяют критерии квадратичного приближения. А именно, подбирают такие значения параметров аппроксимирующей функции, чтобы сумма квадратов отклонений (разницы вычисленных по ней значений и значений из таблицы) была

аппроксимирующих зависимостей на основе этого принципа (метода наименьших квадратов) реализовано в программах Excel и Mathcad. Точность аппроксимирующих зависимостей разного вида может весьма существенно различаться. В ходе аппроксимации с помощью Mathcad могут быть использованы следующие функции:

«expfit» для определения коэффициентов A, B и C экспоненциальной функции вида C ;

«lnfit» для определения коэффициентов A и B логарифмической функции вида f x A ln x B ;

«pwfit» для определения коэффициентов A, B и C степенной функции вида f x A xB C .

3.4.1.Условие задания № 4

Дана таблично заданная функция - пары точек (xi,yi) (прил. 3), для которых необходимо выполнить следующее.

С помощью программы Mathcad провести кусочно-линейную интерполяцию и найти значения y для следующих значений x:

1,3, 2,6, 4,4, 5,9, 7,1, 8,75. Построить график.

С помощью программы Mathcad провести полиномиальную интерполяцию и найти значения y для следующих значений x: 1,3, 2,6, 4,4, 5,9, 7,1, 8,75. Построить график. Записать уравнение полинома (коэффициенты полинома указать с точностью 10 знаков после запятой).

Провести 2 вида аппроксимации согласно варианту (прил.3). Оба графика построить на одной координатной плоскости. В

19

обоих случаях определить сумму квадратов отклонений для узловых точек. Данное задание выполнить как в Excel, так и в

Mathcad.

3.4.2. Пример выполнения задания № 4

Рассмотрим выполнение данного задания для следующих точек.

В качестве видов аппроксимации будем использовать:

полином 3-й степени;

степенную функцию.

Решим все поставленные задачи с использованием Mathcad.

1. Проведем кусочно-линейную интерполяцию с помощью функции «linterp» для заданных точек (xi,yi) и определим значение функции для указанных значений аргумента

(рис. 14).

Рис. 14. Кусочно-линейная интерполяция в Mathcad

2.Проведем полиномиальную интерполяцию с помощью функции «regress». Учитывая, что задано 10 точек, то

20

полином должен быть 9-й степени. В результате получили полином, представленный на рис. 15. Построим его график, на котором отметим исходные точки (рис. 15). По графику можно убедиться, что полином проходит через все заданные точки.

Рис. 15. Полиномиальная интерполяция в Mathcad