Matlab  |  Mathcad  |  Maple  |  Mathematica  |  Statistica  |  Другие пакеты Поиск по сайту
Internet-класс  |  Примеры  |  Методики  |  Банк задач  |  Консультации & Форум  |  Download  |  Ссылки  |  Конкурсы
Научно-практический журнал "Exponenta Pro. Математика в приложениях". Вышел 1/2004 номер журнала
 
Решение нелинейных и дифференциальных уравнений. Решение систем уравнений
Сугробов, Болотский

Архив работы RAR-архив работы (69 кб, mcd)

1. Решение нелинейных уравнений.

В данном разделе мы узнаем, каким образом в системе Mathcad решаются простейшие уравнения вида F(x)=0. Решить уравнение аналитически - значит найти все его корни, т.е. такие числа, при подстановке которых в исходное уравнение получим верное равенство. Графически - значит найти точки пересечения графика функции с осью OX. Для решения уравнений такого типа существует специальная функция:


Root(Выражение, Имя переменной)>


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


Задача решения уравнения в системе Mathcad разбивается на два основных этапа:
1) отделение корней и нахождение начального значения,
2) непосредственное решение уравнения с помощью функции root.


Рассмотрим на примере 1 выполнение этих этапов.


Пример 1



1-й этап:

Задание функции F(x)
Второй способ задания функции через разбиение ее на две части

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


Построив график функции F(x), мы убедились, что точки пересечения графика с осью OX приблизительно равны 0 и 1.2
Аналогично, преобразовав уравнение F(x) к виду f(x)=g(x) и построив два графика, найдем на плоскости точки их пересечения.


Абсциссы этих точек также приблизительно равны 0 и 1.2
Найдем далее непосредственное решение уравнения с помощью функции root, т.е. реализуем второй этап.


2-й этап:

Начальное приближение для первого корня
Нахождение первого корня
Вывод значения

Аналогично для второго корня :

Начальное приближение для второго корня
Нахождение второго корня
Вывод значения

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

Проверка

Анологично решению примера1, решите самостоятельно пример2.

Пример2:



2.Решение систем уравнений.


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


Способ 1(матричный):


Пусть нам дана система уравнений:




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


Для того ,чтобы задать матрицу щелкаем мышью ВИД, затем ПАНЕЛИ ИНСТРУМЕНТОВ, далее МАТЕМАТИКА. Появляется панель МАТЕМАТИКА с различными шаблонами. Щелкаем на шаблоне матрицы и появляется панель МАТРИЦЫ с математическими операциями над ними. При помощи мыши выбираем нужную операцию.Рассмотрим примеры:


Пример 3 : Умножить две матрицы.

Задаем матрицы A и B нужного размера
Находим на панели шаблон скалярного произведения

Пример 4 : Найти обратную матрицу.

Задаем матрицу A
Матрице В присваиваем обратное значение матрицы А (находим на панели шаблон инверсии )
Находим матрицу B

Пример для самостоятельного решения.

Пример 5 : Умножить две матрицы и найти обратную для первой матрицы


Используя примеры 3-4, решим систему уравнений матричным способом.

Пример 6 :
Присвоение переменной А матрицы,состоящей из коэффициентов при переменных в системе
Присвоение переменной В матрицы, состоящей из столбца свободных членов системы
Матричное решение X , выраженное в матричном виде
Корни системы
Проверка

Способ 2 ( использование функции lsolve(A,B) ):


В системе Mathcad введена встроенная функция lsolve (A,B), которая возвращает вектор X для системы линейных уравнений при заданной матрице коэффициентов А и векторе свободных членов В. На следующем примере видна реализация вышеприведенных методов в системе Mathcad для системы из трех уравнений.

Матрица коэффициентов системы линейных уравнений
Вектор свободных членов
Решение системы
Решение с применением функции lsolve
Результат решения
Проверка

Способ 3 ( использование двух функций Given и Find):


Заключается в использовании двух функций системы: Given и Find. Преимущество данного метода заключается в том, что система уравнений вводится в Mathcad без использования матриц, т. е. в "натуральном" виде. Предварительно необходимо только указать начальные значения неизвестных. Это могут быть любые числа, входящие в область определения. (Часто за них принимают столбец свободных членов).

Рассмотрим пример:

Задаем начальные значения неизвестных
Знак равенства между выражениями в уравнениях вызывается с клавиатуры клавишами "CTRL" + "=", или с математической панели БУЛЕВО

Решите самостоятельно пример 7 тремя способами.

Пример 7:




3. Решение дифференциальных уравнений.


Система Mathcad имеет специальную встроенную функцию для решения дифференциальных уравнений. Ее вид :

Odesolve ( x , b [ , steps ] )

Для решения задачи Коши необходимы так называемые начальные условия и указание конца интервала. Эти данные вместе с самим уравнением записываются в блок функции Given и лишь затем применяется сама функция odesolve. Функция имеет ряд особенностей. Если указано число шагов step, то решение выполняется с фиксированным шагом, иначе - адаптивным методом. Ниже приведен пример решения дифференциального уравнения с помощью функции odesolve.



Задано дифференциальное уравнение
Заданы начальные условия
Задано решение дифференциального уравнения
Вычисление производной от b(a)
График решения заданного дифференциального уравнения b(a) и производной от функции решения - c(a)

Пример для самостоятельного решения.

Пример 8 :

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




Карта сайта | На первую страницу | Поиск |О проекте |Сотрудничество |
Exponenta Pro | Matlab.ru

Наши баннеры


Copyright © 2000-2003. Компания SoftLine. Все права защищены.

Дата последнего обновления информации на сайте: 11.05.04
Сайт начал работу 1.09.00

Программное обеспечение Microsoft, Macromedia, VERITAS, Novell, Borland, Symantec, Oracle и др.