Главная / Информатика / Реферат на тему "Возможности математического пакета MathCAd. Приближенные решения уравнений и их систем."

Реферат на тему "Возможности математического пакета MathCAd. Приближенные решения уравнений и их систем."

ФГБОУ ВПО «МОРДОВСКИЙ ГОСУДАРСТВЕННЫЙ ПЕДАГОГИЧЕСКИЙ

ИНСТИТУТ ИМЕНИ М. Е. ЕВСЕВЬЕВА»

Факультет физико-математический


Реферат

на тему:

Возможности математического пакета MathCAd.

Приближенные решения уравнений и их систем.



Выполнила

студентка 2 курса группы МДМ-214____________А. И. Ежова

Направление подготовки 050100 «Педагогическое образование».

Профили подготовки «Математика» и «Информатика».

Реферат проверила ______________________Т. В. Кормилицына











Саранск 2016

Введение

уравнение мathcad векторный погрешность

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

Программирование занимает особое место в Mathcad. Для начального обучения изучать его совершенно не обязательно. Возможности Mathcad позволяют решить подавляющее количество задач без всякого программирования, да к тому же, как правило, несколькими способами. Однако есть класс задач, при решении которых без программирования не обойтись. Это задачи, в которых часть документа из нескольких или многих операторов надо выполнить многократно. В таких случаях документ должен состоять из отдельных подпрограмм, объединенных в единую «головную» программу. И в этом случае программирования можно было бы избежать, если бы создатели Mathcad предусмотрели оператор перехода к метке, который есть (был) во многих языках программирования, в частности в FORTRAN. позволяет вводить программы любой сложности. Образцом достаточно сложной программы является программа решения плоской упругопластической задачи методом конечных элементов.

Программирование в Mathcad чрезвычайно наглядно и понятно, так как программа представляет собой последовательность формул. Основные операторы программирования расположены на панели Programming (Программирование), вызываемой щелчком на кнопке Programming Toolbar (Панель программирования) математической панели. Далее программами будем называть не все документы Mathcad, а только те из них, которые написаны с использованием панели программирования.

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

Многие задачи, решаемые с помощью математических пакетов, сводятся к решению уравнений - алгебраических, степенных, тригонометрических, к поиску значений неизвестных, превращающих эти уравнения в тождества строго или приближенно. Успех в решении подобных задач зависит не только от мощности соответствующих инструментов, встроенных в Mathcad, но и от знания пользователем их особенностей, нюансов, сильных и слабых сторон.


1. Решение уравнений с одной переменной


Уравнения можно решать различными методами, такими как: метод половинного деления, метод простой итерации, метод хорд, метод касательных. В этой главе рассмотрим подробнее вышеперечисленные методы.


1.1 Постановка задачи


Наиболее общий вид нелинейного уравнения:


hello_html_m237e0a59.gifhello_html_m237e0a59.gif(1.1)


где функция hello_html_m2d2b6b5c.gifhello_html_m2d2b6b5c.gif определена и непрерывна на конечном или бесконечном интервале hello_html_21909f0b.gifhello_html_21909f0b.gif.

Определение 1.1. Всякое число hello_html_m529b14b2.gifhello_html_m529b14b2.gif обращающее функцию F(x) в нуль, называется корнем уравнения (1.1).

Определение 1.2. Число hello_html_27c6b16d.gifhello_html_27c6b16d.gif, называется корнем hello_html_m3756eb60.gifhello_html_m3756eb60.gif- ой кратности, если при hello_html_m197e6bc1.gifhello_html_m197e6bc1.gif вместе с функцией hello_html_m6f706b4d.gifhello_html_m6f706b4d.gif равны нулю ее производные до hello_html_m16219f94.gifhello_html_m16219f94.gif -го порядка включительно:


hello_html_m23b47f9b.gifhello_html_m23b47f9b.gif(1.2)


Определение 1.3. Однократный корень называется простым.

Определение 1.4. Уравнения hello_html_m67dbd355.gifhello_html_m67dbd355.gif и hello_html_mb89a757.gifhello_html_mb89a757.gifназываются равносильными (эквивалентными), если множества решений данных уравнений совпадают.

Нелинейные уравнения с одной переменной подразделяются на алгебраические и трансцендентные.

Определение 1.5. Уравнение (1.1) называется алгебраическим, если функция hello_html_m6f706b4d.gifhello_html_m6f706b4d.gif является алгебраической.

Путем алгебраических преобразований из всякого алгебраического уравнения можно получить уравнение в канонической форме:


hello_html_3da3a825.gifhello_html_3da3a825.gif(1.3)


где hello_html_7a1eb03d.gifhello_html_7a1eb03d.gif - действительные коэффициенты уравнения; hello_html_mcccec95.gifhello_html_mcccec95.gif - неизвестное.

Из алгебры известно, что всякое алгебраическое уравнение имеет, по крайней мере, один вещественный или два комплексно сопряженных корня.

Определение 1.6. Уравнение (1.1) называется трансцендентным, если функция hello_html_m6f706b4d.gifhello_html_m6f706b4d.gif не является алгебраической.

Определение 1.7. Решить уравнение (1.1) означает следующее:

установить имеет ли уравнение корни;

определить число корней уравнения;

найти значения корней уравнения с заданной точностью.


1.2 Отделение корней


Определение 1.8. Отделение корней - процедура нахождения отрезков, на которых уравнение (1.1) имеет только одно решение.

В большинстве случаев отделение корней можно провести графически. Для этого достаточно построить график функции hello_html_m6f706b4d.gifhello_html_m6f706b4d.gif и определить отрезки, на которых функция hello_html_m6f706b4d.gifhello_html_m6f706b4d.gif имеет только одну точку пересечения с осью абсцисс.

В сомнительных случаях графическое отделение корней необходимо подкреплять вычислениями. При этом можно использовать следующие очевидные положения: если непрерывная функция принимает на концах отрезка hello_html_73423030.gifhello_html_73423030.gifзначения разных знаков (т. е. hello_html_m2aaad9dd.gifhello_html_m2aaad9dd.gif), то уравнение (1.1) имеет на этом отрезке по меньшей мере один корень;

если функция hello_html_m2b6abc8f.gifhello_html_m2b6abc8f.gifк тому же и строго монотонна, то корень на отрезке единственный.


1.3 Метод половинного деления


Пусть уравнение (1.1) имеет на отрезке hello_html_m2bb96378.gifhello_html_m2bb96378.gif единственный корень, причем функция hello_html_m2d2b6b5c.gifhello_html_m2d2b6b5c.gif на данном отрезке непрерывна (рисунок 1.1).

Разделим отрезок hello_html_73423030.gifhello_html_73423030.gifпополам точкой hello_html_m329cd18e.gifhello_html_m329cd18e.gif. Если hello_html_77ab6cdc.gifhello_html_77ab6cdc.gif,то возможны два случая:

- функция hello_html_m6f706b4d.gifhello_html_m6f706b4d.gif меняет знак на отрезке hello_html_67fa9c1d.gifhello_html_67fa9c1d.gif;

функция hello_html_m6e63d0be.gifhello_html_m6e63d0be.gif меняет знак на отрезке hello_html_m74b6505c.gifhello_html_m74b6505c.gif.

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


hello_html_f2dd07e.gif

Рисунок 1.1 - К объяснению метода половинного деления


Пример 1.1

Решение в пакете Mathcad методом половинного деления уравнения


hello_html_m4160c39.gif


1. Задание функции:


hello_html_2ce09af3.gif


2. Построение графика функции (Рисунок 1.2).


hello_html_m5c70e14.gif

Рисунок 1.2 - График функции hello_html_m3f0f706d.gifhello_html_m3f0f706d.gif


3. Задание функции, реализующей метод половинного деления (Рисунок 1.3). Здесь аргументы функции: hello_html_2b160d01.gifhello_html_2b160d01.gif - имя функции, hello_html_1503d851.gifhello_html_1503d851.gif- левая и правая координаты концов отрезка; hello_html_7cbebf79.gifhello_html_7cbebf79.gif - точность вычисления корня.

4. Вычисление значения корня уравнения:

hello_html_5716a076.gif


5. Проверка найденного значения корня:


hello_html_m48b87bf3.gif


Для рассмотрения процесса нахождения корня уравнения в динамке необходимо сохранить значение корня на каждом шаге вычислительной процедуры и построить зависимость значения корня от номера шага. Функция, возвращающая значение корня на каждом шаге метода половинного деления, представлена на рисунке 1.4. Аргументы функции: hello_html_m5a42b9ec.gifhello_html_m5a42b9ec.gif - имя функции,hello_html_1503d851.gifhello_html_1503d851.gif- левая и правая координаты концов отрезка, hello_html_7cbebf79.gifhello_html_7cbebf79.gif -точность вычисления корня.


hello_html_m37b4f0d8.gif

Рисунок 1.3 - Функция, реализующая метод половинного деления


hello_html_21cf60c7.gif

Рисунок 1.4 - Функция, реализующая метод половинного деления и возвращающая значение корня уравнения на каждом шаге процесса вычислений


После создания функции необходимо дополнить описанный выше документ следующей последовательностью команд.

1. Вычисление матрицы, первый столбец которой содержит номер итерации, второй - значение корня:


hello_html_5beedae3.gif


2. Визуализация зависимости значения корня от номера шага вычислительной процедуры (рисунок 1.5).


hello_html_m7881930.gif

Рисунок 1.5 - Зависимость значения корня от номера шага вычислительной процедуры

1.4 Метод простой итерации


Заменим уравнение (1.1) равносильным уравнением


hello_html_m5ff94c2a.gifhello_html_m5ff94c2a.gif. (1.4)


Пусть hello_html_27c6b16d.gifhello_html_27c6b16d.gif - корень уравнения (1.4), а hello_html_m371fbcfd.gifhello_html_m371fbcfd.gif, полученное каким-либо способом нулевое приближение к корню hello_html_27c6b16d.gifhello_html_27c6b16d.gif. Подставляя hello_html_m371fbcfd.gifhello_html_m371fbcfd.gifв правую часть уравнения (1.4),получим некоторое число hello_html_6415214e.gifhello_html_6415214e.gif. Повторим данную процедуру с hello_html_m607d1b1a.gifhello_html_m607d1b1a.gifи получим hello_html_2ba76097.gifhello_html_2ba76097.gif. Повторяя описанную процедуру, получим последовательность называемую итерационной последовательностью.


hello_html_m23edc68a.gifhello_html_m23edc68a.gif(1.5)


Итерационная последовательность, вообще говоря, может быть как сходящейся, так и расходящейся, что определяется видом функции hello_html_m5a42b9ec.gifhello_html_m5a42b9ec.gif.

Теорема 1.1. Если функция hello_html_m5a42b9ec.gifhello_html_m5a42b9ec.gif непрерывна, а последовательность (1.5) сходится, то предел последовательности (1.5) является корнем уравнения (1.4).

Действительно, пусть hello_html_m17c99510.gifhello_html_m17c99510.gif. Перейдем к пределу в равенстве


hello_html_m63e9b9a1.gifhello_html_m63e9b9a1.gif

hello_html_30ceec6e.gifhello_html_30ceec6e.gif. (1.6)


Условие сходимости итерационного процесса определяется теоремой о достаточном условии сходимости итерационного процесса.

Геометрическая интерпретация данного алгоритма представлена на рисунке 1

hello_html_f2dd07e.gif

Рисунок 1.6 - К объяснению метода простой итерации


Теорема 1.2. Пусть уравнение hello_html_3a5a1685.gifhello_html_3a5a1685.gif имеет единственный корень на отрезке hello_html_5a5a7608.gifhello_html_5a5a7608.gif и выполнены условия:

. hello_html_m5a42b9ec.gifhello_html_m5a42b9ec.gif определена и дифференцируема на hello_html_m34884505.gifhello_html_m34884505.gif.


2. hello_html_m4fd7ede6.gifhello_html_m4fd7ede6.gif для всех hello_html_m15b82bce.gifhello_html_m15b82bce.gif.


3. Существует такое вещественное hello_html_73aa8e73.gifhello_html_73aa8e73.gif, что hello_html_3c9ba70c.gifhello_html_3c9ba70c.gif для всех hello_html_m15d364b3.gifhello_html_m15d364b3.gif.

Тогда итерационная последовательность hello_html_m6d76b9b5.gifhello_html_m6d76b9b5.gif hello_html_m71dffbb7.gifhello_html_m71dffbb7.gif сходится при любом начальном приближении hello_html_6ee3d7de.gifhello_html_6ee3d7de.gif

Доказательство. Построим итерационную последовательность вида (1.5) с любым начальным значением hello_html_m576163a4.gifhello_html_m576163a4.gif. В силу условия 2 теоремы 1.2 все члены последовательности находятся в отрезке hello_html_m34884505.gifhello_html_m34884505.gif.

Рассмотрим два последовательных приближения hello_html_m6d76b9b5.gifhello_html_m6d76b9b5.gif и hello_html_m23d51cbe.gifhello_html_m23d51cbe.gif.

По теореме Лагранжа о конечных приращениях имеем

hello_html_281e9417.gif


Переходя к модулям и принимая во внимание условие 3 теоремы 1.2, получим


hello_html_10aa94c7.gif


При hello_html_576cce03.gifhello_html_576cce03.gif имеем hello_html_m4ad4f727.gifhello_html_m4ad4f727.gif


hello_html_m6ab6fd17.gifhello_html_m6ab6fd17.gif(1.7)

hello_html_162f89d.gif


Рассмотрим ряд


hello_html_mcfb994d.gifhello_html_mcfb994d.gif(1.8)


Составим частичные суммы этого ряда


hello_html_5803bf61.gif


Заметим, что hello_html_m30d6c51.gifhello_html_m30d6c51.gifя частичная сумма ряда (1.8) совпадает с n - ым членом итерационной последовательности (1.5), т. е.


hello_html_m8dc5de2.gifhello_html_m8dc5de2.gif(1.9)


Сравним ряд (1.8) с рядом


hello_html_5432da84.gifhello_html_5432da84.gif(1.10)


Заметим, что в силу соотношения (1.7) абсолютные величины членов ряда (1.8) не превосходят соответствующих членов ряда (1.10). Но ряд (1.10) сходится как бесконечно убывающая геометрическая прогрессия (hello_html_mb71ef05.gifhello_html_mb71ef05.gifпо условию ). Следовательно, и ряд (1.8) сходится, т. е. его частичная сумма (1.9) имеет предел. Пусть hello_html_m17c99510.gifhello_html_m17c99510.gif. В силу непрерывности функции hello_html_62564e5f.gifhello_html_62564e5f.gifполучаем (по формуле 1.6):


hello_html_m536403fc.gif


т. е. hello_html_27c6b16d.gifhello_html_27c6b16d.gif - корень уравнения hello_html_m5ff94c2a.gifhello_html_m5ff94c2a.gif.

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


1.5 Оценка погрешности метода простой итерации


Пусть hello_html_m26239daa.gifhello_html_m26239daa.gif приближение к истинному значению корня уравнения hello_html_m5ff94c2a.gifhello_html_m5ff94c2a.gif. Абсолютная ошибка приближения hello_html_m26239daa.gifhello_html_m26239daa.gif, оценивается модулем


hello_html_m4f04b5da.gifhello_html_m4f04b5da.gif.


Принимая во внимание (1.8) и (1.9), имеем


hello_html_m2320419b.gifhello_html_m2320419b.gif(1.11)


Сравним (1.11) с остатком ряда (1.9):

hello_html_m37c223cf.gifhello_html_m37c223cf.gif(1.12)


Учитывая оценку (1.7), получаем


hello_html_mbd5dd4.gifhello_html_mbd5dd4.gif|hello_html_m519331f.gifhello_html_m519331f.gif|.


Таким образом, для оценки погрешности hello_html_m1d375f05.gifhello_html_m1d375f05.gif - го приближения получается формула


hello_html_8702241.gifhello_html_8702241.gif(1.13)


На практике удобнее использовать модификацию формулы (1.13).

Примем за нулевое приближение hello_html_1e9ede4a.gifhello_html_1e9ede4a.gif(вместоhello_html_m371fbcfd.gifhello_html_m371fbcfd.gif). Следующим приближением будет hello_html_m26239daa.gifhello_html_m26239daa.gif (вместо hello_html_m1291cca6.gifhello_html_m1291cca6.gif). Так как hello_html_7a151c67.gifhello_html_7a151c67.gif то


hello_html_m7177119a.gifhello_html_m7177119a.gif(1.14)


При заданной точности ответа hello_html_7cbebf79.gifhello_html_7cbebf79.gif итерационный процесс прекращается, если


hello_html_m10a7951b.gifhello_html_m10a7951b.gif.


1.6 Преобразование уравнения к итерационному виду


Уравнение hello_html_m6f706b4d.gifhello_html_m6f706b4d.gif=0 преобразуется к виду, пригодному для итерационного процесса, следующим преобразованием

hello_html_m1b304f7f.gif


где hello_html_m2e34ba6b.gifhello_html_m2e34ba6b.gifотличная от нуля константа.

В этом случае


hello_html_m3f6dbf75.gifhello_html_m3f6dbf75.gif(1.15)


Функция hello_html_23e924e7.gifhello_html_23e924e7.gifдолжна удовлетворять условиям теоремы 1.2. Дифференцируя (1.15), получим


hello_html_7796f2f3.gifhello_html_7796f2f3.gif(1.16)


Для выполнения условия 3 теоремы 1.2 достаточно подобрать hello_html_6e78dbd4.gifhello_html_6e78dbd4.gif так, чтобы для всех hello_html_m15b82bce.gifhello_html_m15b82bce.gif


hello_html_m1668df13.gifhello_html_m1668df13.gif(1.17)


1.7 Решение уравнений методом простой итерации в пакете Mathcad


Продемонстрируем использование метода простой итерации на примере нахождения корня уравнения


hello_html_m3151506b.gif


1. Задание функции, стоящей в правой части (1.1):


hello_html_m1c17a98f.gif

2. Задание функции в соответствии с (1.15):


hello_html_64bd1426.gif


3. Задание функции в соответствии с (1.16):


hello_html_m2512036.gifhello_html_m2512036.gif.


4. Построение графиков функций hello_html_m3f03ddb3.gifhello_html_m3f03ddb3.gif (рисунок 1.7).


hello_html_m5327c5e.gif

Рисунок 1.7 - Графики функций hello_html_3e4d59ac.gifhello_html_3e4d59ac.gif и hello_html_3fc3c0c0.gifhello_html_3fc3c0c0.gif


5. Задание функции, реализующей вычислительную схему метода простой итерации на каждом шаге итерационного процесса (рисунок 1.8).


hello_html_m52f37deb.gif

Рисунок 1.8 - Функция, реализующая вычислительную схему метода простой итерации


6. Задание функции, стоящей в правой части (1.15):


hello_html_m26947b5.gif


7. Задание начального приближения: hello_html_5686b0b0.gifhello_html_5686b0b0.gif

8. Вычисление значений корня уравнения на каждом шаге итерационного процесса:


hello_html_m6810283b.gif


9. Визуализация итерационного процесса (рисунок 1.9):


hello_html_4d861a18.gif


hello_html_191cae6c.gif

Рисунок 1.9 - Зависимость значения корня уравнения от номера шага итерационного процесса

10. Вывод точного значения корня: hello_html_7f34e505.gifhello_html_7f34e505.gif.

11. Вывод значения функции: hello_html_m516d5f1d.gifhello_html_m516d5f1d.gif.


1.8 Метод хорд


Пример 1.2

Решить уравнение hello_html_11992007.gifhello_html_11992007.gifметодом хорд с точностью hello_html_m26ab6fd9.gifhello_html_m26ab6fd9.gif

1. Отделяем корни. Этот этап решения осуществляется с помощью аналитического или графического метода. После того как корень, подлежащий уточнению, отделен, за начальное приближение может быть выбрана любая точка hello_html_6bff567f.gifhello_html_6bff567f.gif(начало отрезка, его середина и т. д.).

Воспользуемся графическим методом. Построим график функции и найдем точки пересечения его с осью Ох (рисунок 1.10).


hello_html_341d28d7.gif

Рисунок 1.10 - Отделение корней графически


Получили два интервала: hello_html_5f2b9cae.gifhello_html_5f2b9cae.gif Интервал, в котором мы будем уточнять корень - hello_html_3633fb40.gifhello_html_3633fb40.gif

2. Уточняем корни. Находим первую производную функции

hello_html_m3d87ac91.gifhello_html_m3d87ac91.gif:

hello_html_19ec818a.gif


3. Определяем знаки hello_html_m715d03e8.gifhello_html_m715d03e8.gif на отрезке hello_html_m3c05bde.gifhello_html_m3c05bde.gif:


hello_html_m703a82a4.gif


hello_html_2057f26b.gif

Рисунок 1.11 - Проверка критерия достижения заданной точности


Значит, на данном отрезке действительно существует корень нашего уравнения.

4. Строим последовательность значений с использованием рекуррентной формулы метода хорд и проанализируем результаты вычисленных значений последовательности hello_html_64ea96da.gifhello_html_64ea96da.gif(рисунок 1.11). Для этого рассмотрим значения функции hello_html_36a578.gifhello_html_36a578.gif - эта величина является критерием достижения заданной точности hello_html_m48f22f90.gifhello_html_m48f22f90.gif. Начиная с hello_html_28138281.gifhello_html_28138281.gif, значения hello_html_m34627d28.gifhello_html_m34627d28.gif удовлетворяют критерию достижения заданной точности hello_html_m38e8a93b.gifhello_html_m38e8a93b.gif значит,hello_html_2b1c6381.gifhello_html_2b1c6381.gif является решением нашего уравнения.

5. Создаем функцию, реализующую вычисления корня уравнения hello_html_11992007.gifhello_html_11992007.gif на отрезке hello_html_m3c05bde.gifhello_html_m3c05bde.gif с точностью hello_html_m48f22f90.gifhello_html_m48f22f90.gif методом хорд (рисунок 1.12). Решением будет являться число 1,927, получившееся на третьем шаге решения.


hello_html_3e114629.gif

Рисунок 1.12 - Функция, возвращающая значения корня уравнения методом хорд. Аргументы функцииhello_html_6a9c79ec.gifhello_html_6a9c79ec.gif - концы отрезка; hello_html_110402d3.gifhello_html_110402d3.gif - погрешность вычисленийhello_html_m490cf4c6.gifhello_html_m490cf4c6.gifhello_html_m2042d992.gifhello_html_m2042d992.gif - функция первой производной


6. Проверяем решение (рисунок 1.13)

7.

hello_html_14f04c0d.gif

Рисунок 1.13 - Проверка решения уравнения встроенными функциями Mathcad

Ответ: корень уравнения по методу хорд равен 1,927 с точностью 0,001, найденный на третьем шаге.


1.9 Метод касательных


Пример 1.3

Уточнить методом касательных корень уравнения


hello_html_m3d87ac91.gif


на отрезке hello_html_m3c05bde.gifhello_html_m3c05bde.gif с точностью hello_html_m48f22f90.gifhello_html_m48f22f90.gif.

1. Отделяем корни уравнения.

2. Определяем неподвижную точку.

Для этого определим знаки функции и второй производной на отделенном интервале hello_html_3633fb40.gifhello_html_3633fb40.gif Для этого составим функцию, проверяющую условие неподвижности точки (рисунок 1.14).


hello_html_m7d47db0a.gif

Рисунок 1.14 - Определение неподвижной точки


Тогда подвижной точкой будет точка hello_html_m5fe833a5.gifhello_html_m5fe833a5.gif.

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

hello_html_6c0e944f.gif

Рисунок 1.15 - Построение итерационной последовательности по методу касательных


Анализируя полученные значения для достижения критерия заданной точности, можно сказать, что решением уравнения будет значение


hello_html_m5ddc7ff6.gifhello_html_m5ddc7ff6.gifпри hello_html_fae9224.gifhello_html_fae9224.gif, т. к. hello_html_m2c52c217.gifhello_html_m2c52c217.gif.


4. Создаем функцию, реализующую метод касательных (аналогично методу хорд).

5. Проверяем полученные результаты.

Отметим, что в пакете Mathcad имеется еще несколько функций, позволяющих решать уравнения, например, функция solve, вызываемая с панели Symbolic (рисунок 1.16).


hello_html_m245ee49.gif

Рисунок 1.16 - Панель Symbolic


Пример использования команды solve представлен на рисунке 1.17.


hello_html_m2ea9cb0b.gif

Рисунок 1.17 - Решение уравнения с помощью команды solve


2. Методы решения систем нелинейных уравнений


Системы нелинейных уравнений можно решать различными методами, такими как: метод простых итераций, метод Ньютона, методами спуска. Рассмотрим вышеперечисленные методы.


2.1 Векторная запись нелинейных систем. Метод простых итераций


Пусть требуется решить систему уравнений:


hello_html_m26dabe78.gifhello_html_m26dabe78.gif(2.1)


где hello_html_m3d3e44e5.gifhello_html_m3d3e44e5.gif - заданные, вообще говоря, нелинейные вещественнознаные функции hello_html_m1d375f05.gifhello_html_m1d375f05.gif вещественных переменных hello_html_m1f7c26a6.gifhello_html_m1f7c26a6.gif.

Введя обозначения:


hello_html_1d8a8c88.gif

hello_html_684267ce.gifhello_html_684267ce.gif=hello_html_1564c3b7.gifhello_html_1564c3b7.gif

hello_html_m59f780c.gif


систему (2.1) можно заменить одним уравнением:

hello_html_m67dbd355.gifhello_html_m67dbd355.gif, (2.2)


относительно векторной функции hello_html_46335e99.gifhello_html_46335e99.gif векторного аргумента hello_html_mcccec95.gifhello_html_mcccec95.gif.

Таким образом, исходную задачу можно рассматривать как задачу о нулях нелинейного отображения hello_html_6af2cf88.gifhello_html_6af2cf88.gif. В такой постановке данная задача является прямым обобщением задачи о нахождении решения нелинейного уравнения для случая пространств большей размерности. Это означает, что можно строить методы ее решения как на основе обсужденных в предыдущей лекции подходов, так и осуществлять формальный перенос выведенных для скалярного случая расчетных формул. Однако не все результаты и не все методы оказывается возможным перенести формально (например, метод половинного деления). В любом случае следует позаботиться о правомерности тех или иных операций над векторными переменными и векторными функциями, а также о сходимости получаемых таким способом итерационных процессов. Отметим, что переход от hello_html_7ea9b176.gifhello_html_7ea9b176.gif к hello_html_m632e1941.gifhello_html_m632e1941.gifвносит в задачу нахождения нулей нелинейного отображения свою специфику, учет которой привел к появлению новых методов и различных модификаций уже имеющихся методов. В частности, большая вариативность методов решения нелинейных систем связана с разнообразием способов, которыми можно решать линейные алгебраические задачи, возникающие при пошаговой линеаризации данной нелинейной вектор - функции hello_html_m6f706b4d.gifhello_html_m6f706b4d.gif.

Начнем изучение методов решения нелинейных систем с метода простых итераций.

Пусть система (2.1) преобразована к следующей эквивалентной нелинейной системе:


hello_html_m2f2c32dc.gifhello_html_m2f2c32dc.gif(2.3)

или в компактной записи:


hello_html_m4c2632ac.gifhello_html_m4c2632ac.gif=hello_html_2700ee2b.gifhello_html_2700ee2b.gif. (2.4)


Для этой задачи о неподвижной точке нелинейного отображения hello_html_3f33e0f.gifhello_html_3f33e0f.gifзапишем формальное рекуррентное равенство:


hello_html_1601d2ec.gifhello_html_1601d2ec.gif(2.5)


где hello_html_m65bc0dfc.gifhello_html_m65bc0dfc.gifкоторое определяет метод простых итераций для задачи (2.3).

Если начать процесс построения последовательности hello_html_5bb3fc32.gifhello_html_5bb3fc32.gif с некоторого вектора hello_html_m79e1b4d9.gifhello_html_m79e1b4d9.gif и продолжить вычислительный процесс по формуле (2.5), то при определенных условиях данная последовательность со скоростью геометрической прогрессии будет приближаться к вектору

hello_html_1e090e6b.gifhello_html_1e090e6b.gif- неподвижной точке отображения hello_html_m125e5a75.gifhello_html_m125e5a75.gif.

Справедлива следующая теорема, которую мы приводим без доказательства

Теорема 2.1. Пусть функция hello_html_m125e5a75.gifhello_html_m125e5a75.gifи замкнутое множество hello_html_m112beb4f.gifhello_html_m112beb4f.gifтаковы, что:

1. hello_html_5b8f3841.gifhello_html_5b8f3841.gif для любого hello_html_ca2c977.gifhello_html_ca2c977.gif

2. hello_html_2774e9da.gifhello_html_2774e9da.gif для любого hello_html_m26956b14.gifhello_html_m26956b14.gif

Тогда hello_html_7983d3d7.gifhello_html_7983d3d7.gifимеет в М единственную неподвижную точку hello_html_1e090e6b.gifhello_html_1e090e6b.gif последовтельность hello_html_m4b52322.gifhello_html_m4b52322.gif определяемая (2.5), сходится при любом hello_html_61f97829.gifhello_html_61f97829.gif справедливы оценки

hello_html_m7d57648d.gifhello_html_m7d57648d.gif


для любого hello_html_172ae239.gifhello_html_172ae239.gif

Отметим низкую практическую ценность данной теоремы из - за неконструктивности ее условий. В случаях, когда выбрано хорошее начальное приближение hello_html_57649f4b.gifhello_html_57649f4b.gif решению hello_html_1e090e6b.gifhello_html_1e090e6b.gif, больший практический интерес представляет следующая теорема.

Теорема 2.2. Пусть hello_html_m125e5a75.gifhello_html_m125e5a75.gif дифференцируема в замкнутом шаре hello_html_m2c607ed4.gifhello_html_m2c607ed4.gifпричём hello_html_m72cce20.gifhello_html_m72cce20.gif Тогда если центр hello_html_57649f4b.gifhello_html_57649f4b.gif и радус r шара S таковы, чтоhello_html_5b3e27d4.gifhello_html_5b3e27d4.gif то справедливо заключение теоремы 2.1 hello_html_m1f507e60.gifhello_html_m1f507e60.gif.

Запишем метод последовательных приближений (2.5) в развернутом виде:


hello_html_23034d22.gifhello_html_23034d22.gif(2.6)


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


hello_html_m4039413a.gifhello_html_m4039413a.gif(2.7)

Заметим, что, как и для линейных систем, отдельные уравнения в (2.7) неравноправны, т. е. перемена местами уравнений системы (2.3) может изменить в некоторых пределах число итераций и вообще ситуацию со сходимостью последовательности итераций. Для того чтобы применить метод простых итераций (2.6) или его зейделеву модификацию (2.7) к исходной системе (2.1), необходимо сначала тем или иным способом привести эту систему к виду (2.3). Это можно сделать, например, умножив (2.2) на неособенную hello_html_1111c7b0.gifhello_html_1111c7b0.gif матрицу А и прибавив к обеим частям уравнения hello_html_5d6427c7.gifhello_html_5d6427c7.gif вектор неизвестных. Полученная система


hello_html_5ab98512.gifhello_html_5ab98512.gif, (2.8)


эквивалентна исходной и имеет вид, аналогичный уравнению в методе итераций в одномерном случае. Проблема состоит лишь в правильном подборе матричного параметра.


2.2 Метод Ньютона решения систем нелинейных уравнений


Для решения системы (2.3) будем пользоваться методом последовательных приближений. Предположим, известно hello_html_m42bfaea8.gifhello_html_m42bfaea8.gif е приближение


hello_html_m21d39abf.gifhello_html_m21d39abf.gif,


одного из изолированных корней hello_html_m2ba86a88.gifhello_html_m2ba86a88.gif векторного уравнения (2.2). Тогда точный корень уравнения (2.2) можно представить в виде:


hello_html_m49417c02.gifhello_html_m49417c02.gif, (2.9)

где hello_html_m53d90212.gifhello_html_m53d90212.gif - поправка (погрешность корня). Подставляя выражение (2.9) в (2.2), имеем:


hello_html_m530a5fd3.gifhello_html_m530a5fd3.gif. (2.10)


Предполагая, что функция hello_html_m6f706b4d.gifhello_html_m6f706b4d.gif непрерывно дифференцируема в некоторой выпуклой области, содержащей hello_html_49db667e.gifhello_html_49db667e.gif и hello_html_154ecfb0.gifhello_html_154ecfb0.gif, разложим левую часть уравнения (2.10) по степеням малого вектора hello_html_643feacc.gifhello_html_643feacc.gif, ограничиваясь линейными членами:


hello_html_m3eef92ec.gifhello_html_m3eef92ec.gif,(2.11)


или, в развернутом виде,


hello_html_715620ff.gifhello_html_715620ff.gif(2.12)


Из формул (2.11) и (2.12) видно, что под производной hello_html_m12d47be9.gifhello_html_m12d47be9.gifследует понимать матрицу Якоби системы функций hello_html_m3d3e44e5.gifhello_html_m3d3e44e5.gif относительно переменных hello_html_m1f7c26a6.gifhello_html_m1f7c26a6.gif,т. е.


hello_html_m6b6b5b47.gifhello_html_m6b6b5b47.gif,

или в краткой записи:


hello_html_m342785c.gif


поэтому формула (2.12) может быть записана в следующем виде:


hello_html_m7b9c9646.gifhello_html_m7b9c9646.gif. (2.13)

Если hello_html_1196cb2d.gifhello_html_1196cb2d.gif то

hello_html_7ed598d8.gifhello_html_7ed598d8.gif(2.14)


Отсюда видно, что метод Ньютона решения системы (2.1) состоит в построении итерационной последовательности:


hello_html_d8636d3.gifhello_html_d8636d3.gif(2.15)

где hello_html_2fab670.gifhello_html_2fab670.gif


Если все поправки становятся достаточно малыми, счет прекращается. Иначе новые значения hello_html_3324f04.gifhello_html_3324f04.gif. Используются как приближенные значения корней, и процесс повторяется до тех пор, пока не будет найдено решение или не станет ясно, что получить его не удастся.

Пример 2.1

Найти методом Ньютона приближенное положительное решение системы уравнений:


hello_html_1bef7d85.gif

исходя из начального приближения hello_html_2f8d2864.gifhello_html_2f8d2864.gif.

Полагая:


hello_html_daddbce.gif


имеем:


hello_html_45757673.gif


Отсюда:


hello_html_m2edd3495.gif


Составим матрицу Якоби:


hello_html_43eccd38.gif


Имеем:


hello_html_m7d051bee.gif

причем


hello_html_m17cf4e19.gif


Следовательно, матрица hello_html_m17d8c765.gifhello_html_m17d8c765.gif- неособенная. Вычисляем обратную ей матрицу:


hello_html_47343314.gif


По формуле (2.15) получаем первое приближение:


hello_html_7bafd71c.gif

hello_html_m26dcba7f.gif


Аналогично находятся дальнейшие приближения. Результаты вычислений приведены в таблице 2.1.


Таблица 2.1 - Последовательные приближения корней

i

x

y

z

0

0,5

0,5

0,5

1

0,875

0,5

0,375

2

0,78981

0,49662

0,36993

3

0,78521

0,49662

0,36992

Останавливаясь на приближении, будем иметь:


х = 0,7852; у = 0,4966; z = 0,3699.


Пример 2.2

Решить систему двух нелинейных уравнений


hello_html_m7851d7eb.gif


методом Ньютона.

Решение


hello_html_7fef4e48.gif

Рисунок 2.1 - Задание координатной сетки


1. Зададим координатную сетку и вычислим значения координат hello_html_45a7f150.gifhello_html_45a7f150.gifи hello_html_1371ea34.gifhello_html_1371ea34.gif в узлах сетки (рисунок 2.1).

hello_html_m1ce19099.gif

Рисунок 2.2 - График функции и карта линий уровня


2. Построим график функции и карты линий уровня (рисунок 2.2) (на которых наглядно видно, что данная система имеет решение, и причем единственное) с использованием панели Graph (рисунок 2.3).


hello_html_729e2d8a.gif

Рисунок 2.3 - Панель Graph


3. Точки пересечения линий одинакового уровня дают решение данной системы уравнений.

4. Зададим начальное приближение переменных: hello_html_161f304a.gifhello_html_161f304a.gif


hello_html_6a351826.gif

Рисунок 2.4 - Вектор-функция, задающая систему уравнений



Рисунок 2.5 - Функция, возвращающая решение системы методом Ньютона


6. Зададим функцию (рисунок 2.5), реализующую метод Ньютона (функция hello_html_46335e99.gifhello_html_46335e99.gif возвращает таблицу, содержащую значения координат hello_html_m16f13324.gifhello_html_m16f13324.gif на каждом шаге итерации и соответствующие значения координат вектор-функции).


hello_html_m71a5634e.gif


Запустив программу, получим итерационную последовательность, которая показывает, как находятся приближения. Здесь две первые строки - это значения hello_html_45a7f150.gifhello_html_45a7f150.gifи hello_html_m73bb104c.gifhello_html_m73bb104c.gifсоответственно, а последние две строки - значения данных функций при найденных значениях hello_html_45a7f150.gifhello_html_45a7f150.gifиhello_html_m5afd22cf.gifhello_html_m5afd22cf.gif. В ноль функции обращаются на седьмом шаге. Значит, решением будет являться пара чисел hello_html_m7c859ca.gifhello_html_m7c859ca.gif и hello_html_m5f250f0.gifhello_html_m5f250f0.gif.

7. Проверяем решение системы нелинейных уравнений с помощью блока Given...Minerr (рисунок 2.6).


hello_html_m39b7f412.gif

Рисунок 2.6 - Проверка численного решения при помощи встроенных функций пакета Mathcad


Ответ: hello_html_m7c859ca.gifhello_html_m7c859ca.gif hello_html_m5f250f0.gifhello_html_m5f250f0.gif.


2.3 Решение нелинейных систем методами спуска


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

hello_html_m678f296a.gifhello_html_m678f296a.gif(2.16)


Из функций hello_html_3594a7e4.gifhello_html_3594a7e4.gif,hello_html_53a79cdc.gifhello_html_53a79cdc.gif системы (2.16) образуем новую функцию:


hello_html_m13083364.gifhello_html_m13083364.gif(2.17)


Так как эта функция неотрицательная, то найдется точка (вообще говоря, не единственная) hello_html_5d092ad8.gifhello_html_5d092ad8.gif такая что:


hello_html_m39e35086.gif



Следовательно, если тем или иным способом удается получить точку

hello_html_m55e99ce8.gifhello_html_m55e99ce8.gif, минимизирующую функцию hello_html_1682af79.gifhello_html_1682af79.gif и если при этом окажется, что hello_html_5afbd75b.gifhello_html_5afbd75b.gif то точка hello_html_m57a813d3.gifhello_html_m57a813d3.gif - истинное решение системы (2.16), поскольку

hello_html_15bcadbf.gif


Последовательность точек hello_html_365f8453.gifhello_html_365f8453.gif - приближений к точке hello_html_m55e99ce8.gifhello_html_m55e99ce8.gif минимума функции hello_html_344e47ef.gifhello_html_344e47ef.gif - обычно получают по рекуррентной формуле:


hello_html_m2537ec92.gifhello_html_m2537ec92.gif(2.18)


где hello_html_m40a6543d.gifhello_html_m40a6543d.gif- вектор, определяющий направление минимизации, а hello_html_3b924a96.gifhello_html_3b924a96.gif- скалярная величина, характеризующая величину шага минимизации (шаговый множитель). Учитывая геометрический смысл задачи минимизации функций двух переменных hello_html_344e47ef.gifhello_html_344e47ef.gif - "спуск на дно" поверхности hello_html_11ef559.gifhello_html_11ef559.gif (рис. 2.4), итерационный метод (2.18) можно назвать методом спуска, если вектор hello_html_m207e140b.gifhello_html_m207e140b.gif при каждом hello_html_m3173feb0.gifhello_html_m3173feb0.gif является направлением спуска (т. е. существует такое hello_html_2f385099.gifhello_html_2f385099.gif, что hello_html_m4d5b1373.gifhello_html_m4d5b1373.gif) и если множитель hello_html_60ece9b.gifhello_html_60ece9b.gif подбирается так, чтобы выполнялось условие релаксации hello_html_ma3916e7.gifhello_html_ma3916e7.gif<hello_html_m6138b0bf.gifhello_html_m6138b0bf.gif означающее переход на каждой итерации в точку с меньшим значением минимизируемой функции.

Таким образом, при построении численного метода вида (2.18) минимизации функцииhello_html_m78799f03.gifhello_html_m78799f03.gif следует ответить на два главных вопроса: как выбирать направление спуска hello_html_210cde8d.gifhello_html_210cde8d.gif и как регулировать длину шага в выбранном направлении с помощью скалярного параметра - шагового множителя hello_html_3b924a96.gifhello_html_3b924a96.gif? Приведем простые соображения по этому поводу.

При выборе направления спуска естественным является выбор такого направления, в котором минимизируемая функция убывает наиболее быстро.


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


hello_html_mea5e80a.gif


- антиградиент функции hello_html_344e47ef.gifhello_html_344e47ef.gif Таким образом, из семейства методов (2.18) выделяем градиентный метод:


hello_html_22cd4c7.gifhello_html_22cd4c7.gif(2.19)


Оптимальный шаг в направлении антиградиента - это такой шаг, при котором значение hello_html_ma3916e7.gifhello_html_ma3916e7.gif - hello_html_6ba2d527.gifhello_html_6ba2d527.gifнаименьшее среди всех других значений hello_html_344e47ef.gifhello_html_344e47ef.gif в этом фиксированном направлении, т. е. когда точка hello_html_m5b83366f.gifhello_html_m5b83366f.gif является точкой условного минимума. Следовательно, можно рассчитывать на наиболее быструю сходимость метода (2.19), если полагать в нем


hello_html_m5fd1eb01.gifhello_html_m5fd1eb01.gif. (2.20)


Такой выбор шагового множителя, называемый исчерпывающим спуском, вместе с формулой (2.19) определяет метод наискорейшего спуска.

Геометрическая интерпретация этого метода хорошо видна из рисунка 2.4, 2.5. Характерны девяностоградусные изломы траектории наискорейшего спуска, что объясняется исчерпываемостью спуска и свойством градиента (а значит, и антиградиента) быть перпендикулярным к линии уровня в соответствующей точке.

Наиболее типичной является ситуация, когда найти точно (аналитическими методами) оптимальное значение hello_html_3b924a96.gifhello_html_3b924a96.gif не удается. Следовательно, приходится делать ставку на применение каких-либо численных методов одномерной минимизации и находить hello_html_3b924a96.gifhello_html_3b924a96.gif в (2.18) лишь приближенно.

Несмотря на то, что задача нахождения минимума функции одной переменной hello_html_m35119810.gifhello_html_m35119810.gif намного проще, чем решаемая задача, применение тех или иных численных методов нахождения значений hello_html_m74391159.gifhello_html_m74391159.gif с той или иной точностью требует вычисления нескольких значений минимизируемой функции. Так как это нужно делать на каждом итерационном шаге, то при большом числе шагов реализация метода наискорейшего спуска в чистом виде является достаточно высокозатратной. Существуют эффективные схемы приближенного вычисления квазиоптимальных hello_html_3b924a96.gifhello_html_3b924a96.gif, в которых учитывается специфика минимизируемых функций (типа сумм квадратов функций).

Зачастую успешной является такая стратегия градиентного метода, при которой шаговый множитель hello_html_3b924a96.gifhello_html_3b924a96.gif в (2.18) берется либо сразу достаточно малым постоянным, либо предусматривается его уменьшение, например, делением пополам для удовлетворения условию релаксации на очередном шаге. Хотя каждый отдельный шаг градиентного метода при этом, вообще говоря, далек от оптимального, такой процесс по числу вычислений функции может оказаться более эффективным, чем метод наискорейшего спуска.

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

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

Примечание. Порядком сходимости последовательности hello_html_m68c84bf.gifhello_html_m68c84bf.gif к hello_html_m6d075155.gifhello_html_m6d075155.gif называют такое число hello_html_m76a1eb1e.gifhello_html_m76a1eb1e.gif, что


hello_html_21dd8dee.gif


где hello_html_1d2261d5.gifhello_html_1d2261d5.gif при всех hello_html_m58bf9c73.gifhello_html_m58bf9c73.gif.

Разработан ряд методов решения экстремальных задач, которые соединяют в себе низкую требовательность к выбору начальной точки и высокую скорость сходимости. К таким методам, называемым квазиньютоновскими, можно отнести, например, метод переменной метрики (Дэвидона-Флетчера-Пауэлча), симметричный и положительно определенный методы секущих (на основе формулы пересчета Бройдена).

При наличии негладких функций в решаемой задаче следует отказаться от использования производных или их аппроксимаций и прибегнуть к так называемым методам прямого поиска Циклического покоординатного спуска, Хука и Дживса, Розенброка и т. п.).

Замечание 1. Для разных семейств численных методов минимизации могут быть рекомендованы свои критерии останова итерационного процесса. Например, учитывая, что в точке минимума дифференцируемой функции должно выполняться необходимое условие экстремума, на конец счета градиентным методом можно выходить, когда достаточно малой становится норма градиента. Если принять во внимание, что минимизация применяется к решению нелинейной системы, то целесообразнее отслеживать близость к нулю минимизируемой неотрицательной функции, т. е. судить о точности получаемого приближения по квадрату его евклидовой метрики.

Замечание 2. Для решения hello_html_m1d375f05.gifhello_html_m1d375f05.gif - мерной системы (2.1) следует свести задачу к решению экстремальной задачи:


hello_html_m455fe1c4.gif


Рассмотрим далее примеры реализации некоторых алгоритмов поиска экстремумов функций, зависящих от нескольких переменных в пакете Mathcad.





































Заключение


В результате выполнения данной курсовой работы были получены навыки по использованию системы компьютерной математики Mathcad на примерах решения нелинейных уравнений и их систем.

В ходе выполнения работы были рассмотрены следующие программные реализации методов решений уравнений и их систем:

Метод половинного деления;

Метод простой итерации;

Метод хорд;

Метод касательных;

Метод спуска.


Список используемых источников


1. Гурский, Д.А. Вычисления в MATCHCAD 12 / Д.А.Гурский, Е.С. Турбина. - СПб.: Питер, 2006. - 544с.

2. Поршнев, С.В. Численные методы на базе MATCHCAD /С.В. Поршнев, И.В. Беленкова. - СПб.: БХВ-Питербург, 2005. - 464с.

. Макаров, Е.Г. Инженерные расчёты в MATCHCAD 14 / Е.Г Макаров. -СПб.: Питер, 2007.- 592с.

. Очков, В. Mathcad 14 для студентов, инженеров и конструкторов / В. Очков. - BHV.: - Спб, 2007. - 368с.

. Шушкевич, Г. Компьютерные технологии в математике. Система Mathcad 14. Часть1 / Г. Шушкевич, С. Шушкевич. - Издательство Гревцова. 2010. - 288с.

. Максфилд, Б. Mathcad в инженерных расчётах/Б. Максфилд.- Корона-век, 2012. - 368с.

. Охорзин, В.А. Прикладная математика в системе Mathcad/ В.А.Охорзин.- Лань, 2009. - 352с.

. Копчёнова, Н.В. Вычислительная математика в примерах и задачах/Н.В.Копчёнова, И.А.Марон. - М.: - Наука, 1972. - 368с.

. Дьяконов, В. Mathcad 2000. Учебный курс / В. Дьяконов. - СПб.: Питер, 2001. - 592с.

. Березин, И.С. Методы вычислений / И.С. Березин. Н.П. Жидков.-М.: - Наука,1966. - 632с.

Реферат на тему "Возможности математического пакета MathCAd. Приближенные решения уравнений и их систем."
  • Информатика
Описание:

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

Программирование занимает особое место в Mathcad. Для начального обучения изучать его совершенно не обязательно. Возможности Mathcad позволяют решить подавляющее количество задач без всякого программирования, да к тому же, как правило, несколькими способами. Однако есть класс задач, при решении которых без программирования не обойтись. Это задачи, в которых часть документа из нескольких или многих операторов надо выполнить многократно. В таких случаях документ должен состоять из отдельных подпрограмм, объединенных в единую «головную» программу. И в этом случае программирования можно было бы избежать, если бы создатели Mathcad предусмотрели оператор перехода к метке, который есть (был) во многих языках программирования, в частности в FORTRAN. позволяет вводить программы любой сложности. Образцом достаточно сложной программы является программа решения плоской упругопластической задачи методом конечных элементов.

Автор Ежова Александра Ивановна
Дата добавления 12.05.2016
Раздел Информатика
Подраздел Другое
Просмотров 212
Номер материала MA-066886
Скачать свидетельство о публикации

Оставьте свой комментарий:

Введите символы, которые изображены на картинке:

Получить новый код
* Обязательные для заполнения.


Комментарии:

↓ Показать еще коментарии ↓