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

Foodband

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

Mathcad предлагает новый способ для решения обыкновенных дифференциальных уравнений, разрешенных относительно старшей производной. Для этих целей служит уже известный нам блок given совместно с функцией odesolve. Дифференциальное уравнение совместно с начальными или граничными условиями записывается в блоке given. Производные можно обозначать как штрихами (Ctrl+F7), так и с помощью знака производной . Пример использования функции для решения задачи Коши приведен ниже.

У искомой функции явно указан аргумент, знак производной стоит перед скобкой.

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

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

Решение уравнений в частных производных

Одним из методов решения дифференциальных уравнений в частных производных является метод сеток. Идея метода заключается в следующем. Для простоты, ограничимся случаем только функции двух переменных, и будем полагать, что решение уравнения ищется на квадратной области единичного размера. Разобьем область сеткой. Шаг сетки по оси x и по оси y, вообще говоря, может быть разный. По определению частная производная равна

Если рассматривать функцию только в узлах сетки, то частную производную можно записать в форме

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

Такое выражение называется левой конечной разностью. Можно получить центральную конечную разность, найдя среднее этих выражений.

Теперь получим выражения для вторых производных.

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

1. Уравнения гиперболического типа

В качестве примера рассмотрим решение волнового уравнения (уравнения гиперболического типа).

Уравнение будем решать методом сеток. Запишем уравнение в конечных разностях

Полученное уравнение позволяет выразить значение функции u в момент времени через значения функции в предыдущие моменты времени.

Такая разностная схема называется явной, так как искомая величина получается в явном виде. Она устойчива, если.

Зададим начальные условия: смещение струны U в начальный и последующий моменты времени описывается синусоидальной функцией.

(Совпадение смещений при j=0 и j=1 соответствует нулевой начальной скорости).

Зададим граничные условия: на концах струны смещение равно 0 в любой момент времени

Будем полагать коэффициент

Записываем уравнение в конечных разностях, разрешенное относительно

Представляем результат на графике

2. Уравнения параболического типа

Еще один пример использования конечных разностей — уравнение диффузии.

Это уравнение параболического типа. Явная разностная схема для этого уравнения имеет вид

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

и диапазон изменения пространственной и временной координат:

Задаем начальные и граничные условия

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

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

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

Рассмотрим применение неявной схемы на примере уравнения теплопроводности

Запишем неявную разностную схему для этого уравнения

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

или в матричной форме

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

Задаем значения параметров

и начальное распределение температуры в области

Формируем матрицы уравнения (7)

Находим решение системы

3. Решение уравнений Лапласа и Пуассона

Для решения уравнений Пуассона и Лапласа (частный случай, когда ) — уравнений эллиптического типа — предназначена функция relax (a, b, c, d, e, f, u, rjac), реализующая метод релаксации. Фактически, эту функцию можно использовать для решения эллиптического уравнения общего вида

которое может быть сведено к уравнению в конечных разностях

В частности, для уравнения Пуассона коэффициенты .

Идея метода релаксации заключается в следующем. Если нет источников (уравнение Лапласа), то значение функции в данном узле на текущем шаге определяется как среднее значение функции в ближайших узлах на предыдущем шаге k

При наличии источников разностная схема имеет вид

Метод релаксации сходится достаточно медленно, так как фактически он использует разностную схему (3) с максимально возможным для двумерного случая шагом .

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

Функция relax возвращает квадратную матрицу, в которой:

расположение элемента в матрице соответствует его положению внутри квадратной области,

это значение приближает решение в этой точке.

Эта функция использует метод релаксации для приближения к решению.

Вы должны использовать функцию relax, если Вы знаете значения искомой функции u (x, y) на всех четырех сторонах квадратной области.

a, b, c, d, e — квадратные матрицы одного и того же размера, содержащие коэффициенты дифференциального уравнения.

f — квадратная матрица, содержащая значения правой части уравнения в каждой точке внутри квадрата

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

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

Задаем правую часть уравнения Пуассона — два точечных источника

Задаем значения параметров функции relax

Задаем граничные условия и начальное приближение — нули во всех внутренних точках области

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

Если граничные условия равны нулю на всех четырех сторонах квадрата, можно использовать функцию multigrid.

Алгоритм метода достаточно громоздкий, поэтому рассматривать его мы не будем.

28. Тема 7. Решение дифференциальных уравнений и систем в MathCad. Краткие теоретические сведения

Для решения дифференциальных уравнений с начальными условиями система Mathcad имеет ряд встроенных функций:

Rkfixed – функция для решения ОДУ и систем ОДУ методом Рунге–Кутта четвертого порядка с постоянным шагом;

Rkadapt – функция решения ОДУ и систем ОДУ методом Рунге–Кутта с переменным шагом;

Odesolve – функция, решающая ОДУ блочным методом.

Ниже приведено описание стандартной функции Rkfixed с указанием параметров функции.

Rkfixed(y, x1, x2, p, D)

Y – вектор начальных условий из K элементов (k – количество уравнений в системе);

X1 и X2 – левая и правая границы интервала, на котором ищется решение ОДУ или системы ОДУ;

P – число точек внутри интервала (x1, x2), в которых ищется решение;

D – вектор, состоящий из K-Элементов, который содержит первую производную искомой функции или первые производные искомых функций, если речь идет о решении системы.

Результатом работы функции является матрица из p+1 строк, первый столбец которой содержит точки, в которых получено решение, а остальные столбцы – сами решения.

На рисунке 2.7.1 приведены конкретные примеры решения различных дифференциальных уравнений и систем ОДУ в MathCAD.

При решении дифференциального уравнения первого порядка нужно создать вектор начальных условий из одного элемента Y1, который затем используется при формировании вектора-функции правой части дифференциального уравнения. При обращении к функции Rkfixed Указывается имя вектора Y, границы интервала, на котором ищется решение уравнения, например, (0 ; 2), количество точек, в которых ищется решение – 100, вектор-функция, описывающая правую часть дифференциального уравнения – D. В результате получается матрица Z, в первом столбце которой содержатся значения аргумента искомой функции, во втором – значения самой результирующей функции. При построении графика функции первый столбец полученной матрицы указывается как аргумент, второй столбец – как функция.

При решении системы дифференциальных уравнений нужно создать вектор начальных условий из двух элементов, например, вектор V, который затем используется при формировании вектора-функции правой части дифференциального уравнения. При обращении к функции Rkfixed Указывается имя вектора V, и границы интервала, на котором ищется решение уравнения, например, (0 ; 5), количество точек, в которых ищется решение – 100, вектор-функция, описывающая правую часть дифференциального уравнения – D. В результате получается матрица S, в первом столбце которой содержатся значения аргумента искомых функций, во втором и третьем столбцах – значения самих функций при соответствующем значении аргумента. При построении графика можно воспользоваться первым столбцом полученной матрицы как аргументом, а вторым и третьим столбцами – как функциями.

На рисунке 2.7.2 приведен пример решения дифференциального уравнения второго порядка с использованием функции Rkfixed. Необходимо решить дифференциальное уравнение второго порядка с заданными начальными условиями вида:

Рисунок 2.7.2 – Пример решения дифференциальных уравнений второго порядка с помощью Rkfixed

Для решения уравнения с помощью функции Rkfixed нужно выполнить замену переменных и привести дифференциальное уравнение второго порядка к двум дифференциальным уравнениям первого порядка. Вид этих уравнений приведен ниже.

Документ формируется точно так же, как и при решении системы ОДУ.

На рисунке 2.7.2 показана возможность вычисления вектора второй производной найденной функции – вектора А, построены графики исходной функции, функций первой и второй производных.

Реферат на тему: "Возможности математического пакета MathCAd. Вычисления, производные, интегралы, уравнения и их системы. Обыкновенные дифференциальные уравнения и их системы. Символьное решение."

дифференциальные уравнения и их системы. Символьное решение.

Выполнила студентка группы МДИ-117

Содержание

Введение

Математические и научно-технические расчеты являются важной сферой применения персональных компьютеров. Часто они выполняются с помощью программ, написанных на языке высокого уровня, например Бейсике или Паскале. Сегодня эту работу нередко выполняет обычный пользователь ПК. Для этого он вынужден изучать языки программирования и многочисленные, подчас весьма тонкие капризные численные методы математических расчетов. Нередко при этом из под руки способного физика, химика или инженера выходят далёкие от совершенства программы.

Это не вполне нормальное положение может изменить к лучшему применение интегрированных программных систем автоматизации математических расчетов (Eureka, MathCAD, MatLab и др.). Рассматрим возможности и эволюция одной из таких систем — MathCAD.

Foodband

MATHCAD — универсальный математический пакет, предназначенный для выполнения инженерных и научных расчетов. Основное преимущество пакета — естественный математический язык, на котором формируются решаемые задачи. Объединение текстового редактора с возможностью использования общепринятого математического языка позволяет пользователю получить готовый итоговый документ. Пакет обладает широкими графическими возможностями, расширяемыми от версии к версии. Практическое применение пакета существенно повышает эффективность интеллектуального труда.

От других продуктов аналогичного назначения MATHCAD отличается ориентация на создание высококачественных документов (докладов, отчетов, статей) в режиме WYSIWYG (What You See Is What You Get). Это означает, что, внося изменения, пользователь немедленно видит их результаты и в любой момент может распечатать документ во всем блеске. Работа с пакетом за экраном компьютера практически совпадает с работой на бумаге с одной лишь разницей — она более эффективна. Преимущества MATHCAD состоит в том, что он не только позволяет провести необходимые расчеты, но и оформить свою работу с помощью графиков, рисунков, таблиц и математических формул. А эта часть работы является наиболее рутинной и мало творческой, к тому же она и время емкая и малоприятная.

MathCAD является интегрированной системой программирования, ориентированной на проведение математических и инженерно-технических расчетов.

Система MathCAD содержит текстовый редактор, вычислитель и графический процессор.

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

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

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

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

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

Элементарные вычисления

Константами называются поименованные объекты, хранящие некоторые значения, которые не могут быть изменены.

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

Переменные являются поименованными объектами, которым присвоено некоторое значение, которое может изменяться по ходу выполнения программы. Тип переменной определяется ее значением; переменные могут быть числовыми, строковыми, символьными и т. д.

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

— латинские буквы любого регистра;

— арабские цифры от 0 до 9;

— символ подчеркивания (_), символ процент (%) и символ (.);

— буквы греческого алфавита (набираются с использованием клавиши Ctrl или применяется палитра греческих букв).

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

Операторы, обозначающие основные арифметические действия, вводятся с панели Calculator (Калькулятор, Арифметика).

Вычислительные операторы вставляются в документы при помощи панели инструментов Calculus (Матанализ).

Вычислительные операторы сгруппированы на панели Evaluation (Вычисления):

— Численный вывод (Evaluate Numerically) =

— Символьный (аналитический) вывод (Evaluate Symbolically) →

— Глобальное присваивание (Global Definition) ≡.

hello_html_749c67fd.png

Рассмотрим самый простой пример вычисления:

С помощью Calculator записываем дробные числа и арифметический знак сложения.

hello_html_47a52398.png

Более сложный пример вычислений: так же решается с помощью Calculator .

hello_html_1695a062.png

Переменные должны быть предварительно определены пользователем, т. е. им необходимо хотя бы однажды присвоить значение. В качестве оператора присваивания используется знак :=, тогда как знак = отведен для вывода значения константы или переменной. Попытка использовать неопределенную переменную ведет к выводу сообщения об ошибке. В MathCAD различают: локальные и глобальные переменные.

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

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

hello_html_mc2b30ce.png

Производные

Для аналитического вычисления производной нужно выбирать кнопку d/dx на панели Calсulus. На рабочем листе в окошко после оператора производной нужно вписать вычисляемое выражение. Теперь ввести знак стрелки с панели, либо набрать на клавиатуре сочетание Ctrl+”.”. Нажать F9. Значение производной функции будет выдано в виде математического выражения.

hello_html_m52a2510a.png

Решение задачи нахождения производной в определенной точке осуществляйте по следующей схеме. Сначала некоторой новой функции нужно присвоить значение производной от заданной функции. Затем подставить значение известной точки в эту функцию. Правильным будет и другой вариант. Можно задать известное значение точки, а затем вычислить производную от нужной функции.

hello_html_m4c0e2785.png

Вычисление производных высших порядков выполняйте с помощью кнопки dn/dxn, расположенной также в панели Calculus. Важно помнить, что показатель порядка n должен быть обязательно натуральным числом.

hello_html_43132185.png

Интегралы

Оператор интегрирования в Mathcad предназначен для численного вычисления определенного интеграла функции по некоторому интервалу. Для этого нужно выполнить следующие действия:

Наберате знак &. Появится знак интеграла с пустыми полями для подинтегрального выражения, пределов интегрирования и переменной интегрирования.

Далее задаются верхний и нижний пределы интегрирования.

Щёлкните на поле между знаком интеграла и d и набрать выражение, которое нужно интегрировать.

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

hello_html_mf5ed9ed.png

Переменные пределы интегрирования

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

hello_html_m1939f96d.png

Криволинейные и двойные интегралы

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

hello_html_m4f4d4d88.png

В Mathcad можно также вычислять двойные или кратные интегралы. Чтобы ввести знак двойного интеграла, нужно набирать & дважды. Ввести подынтегральное выражение, пределы и переменные интегрирования для каждого интеграла.

hello_html_m2dd9f61c.png

Численное решение нелинейного уравнения

Алгоритм приближенного решения уравнения f(x)=0 состоит из двух этапов:

1. нахождения промежутка, содержащего корень уравнения (или начальных приближений для корня);

2. получения приближенного решения с заданной точностью с помощью функции root

hello_html_2a365bd2.png

Нахождение корней полинома

Для нахождения корней выражения, имеющего вид v 0 +v 1 x+… v n-1 x n-1 +v n x n , лучше использовать функцию polyroots, нежели root. В отличие от функции root, функция polyroots не требует начального приближения и возвращает сразу все корни, как вещественные, так и комплексные.

Функция Polyroots(v) — возвращает корни полинома степени n. Коэффициенты полинома находятся в векторе v длины n + 1. Возвращает вектор длины n, состоящий из корней полинома.

hello_html_m2f437b56.png

Решение систем уравнений с помощью функции Lsolve

Системы линейных уравнений удобно решать с помощью функции lsolve. Функция lsolve(А, b) — возвращает вектор решения x такой, что Ах = b.

hello_html_m159e2e7.png

Решение системы уравнений методом Гаусса

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

В MathCAD прямой и обратный ходы метода Гаусса выполняет функция rref(A).

hello_html_m2d52eed0.png

Решение систем уравнений с помощью функций Find

Для решения системы уравнений с помощью функции Find необходимо выполнить следующее:

1. Задать начальное приближение для всех неизвестных, входящих в систему уравнений. MathCAD решает систему с помощью итерационных методов;

2. Напечатать ключевое слово Given. Оно указывает MathCAD, что далее следует система уравнений;

3. Введите уравнения и неравенства в любом порядке. Используйте [Ctrl]= для печати символа =. Между левыми и правыми частями неравенств может стоять любой из символов <, >, ≥ и ≤;

4. Введите любое выражение, которое включает функцию Find, например: х:= Find(х, у). Ключевое слово Given, уравнения и неравенства, которые следуют за ним, и какое — либо выражение, содержащее функцию Find, называют блоком решения уравнений.

hello_html_624d658e.png

Символьное решение уравнений

Имеются некоторые задачи, для которых возможности MathCAD позволяют находить решения в символьном (аналитическом) виде.

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

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

• если нужно найти все комплексные корни полинома со степенью меньше или равной 4, символьное решение даст их точные значения в одном векторе или в аналитическом или цифровом виде.

Команда Символы → Переменные → Вычислить позволяет решить уравнение относительно некоторой переменной и выразить его корни через остальные параметры уравнения.

Чтобы решить уравнение символьно, необходимо:

1. Напечатать выражение (для ввода знака равенства используйте комбинацию клавиш Ctrl + =);

2. Выделить переменную, относительно которой нужно решить уравнение, щелкнув на ней мышью;

3. Выбрать пункт меню Символы → Переменные → Вычислить.

Нет необходимости приравнивать выражение нулю. Если MathCAD не находит знака равенства, он предполагает, что требуется приравнять выражение нулю.

Чтобы решить систему уравнений в символьном виде, необходимо выполнить следующее:

1. Напечатать ключевое слово Given;

2. Напечатать уравнения в любом порядке ниже слова Given. Удостоверьтесь, что для ввода знака = используется Ctrl + =;

3. Напечатать функцию Find, соответствующую системе уравнений;

4. Нажать Ctrl + .(клавиша CTRL, сопровождаемая точкой). MathCAD отобразит символьный знак равенства →;

5. Щелкнуть мышью на функции Find.

hello_html_8fcb784.png

Дифференциальные уравнения и системы в MathCad

Для решения дифференциальных уравнений с начальными условиями система Mathcad имеет ряд встроенных функций:

rkfixed – функция для решения ОДУ и систем ОДУ методом Рунге–Кутта четвертого порядка с постоянным шагом;

Rkadapt – функция решения ОДУ и систем ОДУ методом Рунге–Кутта с переменным шагом;

Odesolve – функция, решающая ОДУ блочным методом.

rkfixed(y, x1, x2, p, D)

y – вектор начальных условий из k элементов (k – количество уравнений в системе);

x1 и x2 – левая и правая границы интервала, на котором ищется решение ОДУ или системы ОДУ;

p – число точек внутри интервала (x1, x2), в которых ищется решение;

D – вектор, состоящий из k-элементов, который содержит первую производную искомой функции или первые производные искомых функций, если речь идет о решении системы.

Результатом работы функции является матрица из p+1 строк, первый столбец которой содержит точки, в которых получено решение, а остальные столбцы – сами решения.

hello_html_7b822bef.png

При решении системы дифференциальных уравнений нужно создать вектор начальных условий из двух элементов, например, вектор v, который затем используется при формировании вектора-функции правой части дифференциального уравнения. При обращении к функции rkfixed указывается имя вектора v, и границы интервала, на котором ищется решение уравнения, например, (0 ; 5), количество точек, в которых ищется решение – 100, вектор-функция, описывающая правую часть дифференциального уравнения – D. В результате получается матрица s, в первом столбце которой содержатся значения аргумента искомых функций, во втором и третьем столбцах – значения самих функций при соответствующем значении аргумента. При построении графика можно воспользоваться первым столбцом полученной матрицы как аргументом, а вторым и третьим столбцами – как функциями.

hello_html_m5d557f8.png

Заключение

И так, перечислим основные достоинства MATHCAD`a.

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

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

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

Но, к сожалению, популярный во всем мире пакет MATHCAD фирмы MathSoft, в России распространен еще слабо, как и все программные продукты подобно рода.

Наверное, это оттого, что люди, живущие в России, ещё не привыкли к тому, что решить систему дифференциальных уравнений из пяти переменных шестого порядка можно не только с помощью карандаша и бумаги, но и с помощью компьютера и MATHCAD`a. Зачем человеку с высшим образованием, который знает и может решить эту систему, решать её на бумаге, когда можно переложить эту рутинную работу на плечи мощных вычислительных машин.

Список литературы

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

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

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

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

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

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

Foodband

Добавить комментарий

Ваш адрес email не будет опубликован. Обязательные поля помечены *