Лекции.Орг


Поиск:




Категории:

Астрономия
Биология
География
Другие языки
Интернет
Информатика
История
Культура
Литература
Логика
Математика
Медицина
Механика
Охрана труда
Педагогика
Политика
Право
Психология
Религия
Риторика
Социология
Спорт
Строительство
Технология
Транспорт
Физика
Философия
Финансы
Химия
Экология
Экономика
Электроника

 

 

 

 


Решение с помощью функции lsolve().




1. Задаем матрицу системы А и матрицу свободных членов В. Для этого пишем А, вводим двоеточие (которое автоматически заменится на знак присваивания «:=»), затем на панели инструментов для работы с матрицами и векторами нажимаем кнопку создания матрицы или вектора и вводим матрицу системы. Аналогично вводим столбец свободных членов.

 

2. Вводим lsolve(A,B)=

Пример:

 

Решение с помощь функции Given…Find

1. Вводим начальные приближения: x1:=1, x2:=1, x3:=1.

2. Вводим слово Given

3. Ниже вводим систему уравнений, используя при этом вместо обычного знака равно знак булева равенства (вводится нажатием Ctrl+=).

4. Ниже пишем: Find(x1,x2,x3)=

 

Пример:

Решение СЛАУ методом обратной матрицы

1. Вводим матрицу системы и матрицу свободных членов, как и при решении с помощью lsolve().

2. Вводим x:=A-1B. Для ввода обозначения обратной матрицы нужно воспользоваться соответствующей кнопкой панели инструментов «Матрицы».

3. Вводим x=

Пример:

 

 

Решение нелинейных уравнений с помощью MathCad.

Локализация корней.

Предположим нам нужно уточнить корни уравнения f(x):=x cos(x).

Для локализации корней удобно построить график функции. Для этого нужно:

1. Определить функцию, т.е. ввести f(x):=x cos(x)

2. Выбрать Вставить-Графики-Зависимость XY.

3. Под осью абсцисс ввести х, слева от оси ординат ввести f(x).

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

Основные свойства графика можно настроить выбрав пункт Формат в контекстном меню, вызываемом щелчком правой кнопкой мыши на графике.

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

Уточнение корней.

Уточнение корней производится с помощью функции root():

1. Устанавливаем точность вычислений, изменяя значения пункта Порог сложности на закладке Толерантность окна, вызываемого командой Формат-Результат. На закладке Формат номера того же окна выбираем способ отображения результата и количество отображаемых десятичных знаков.

2. Вводим root(f(x),x)=

Пример:

 

 

 

 

Интерполирование функций средствами MathСad

В MathCad имеются средства для сплайновой интерполяции.

Интерполирование производится с помощью функции interp (vs,vx,vy,x), которая возвращает вектор значений функции y в интересующих точках х.

Здесь и далее vx – вектор исходных значений аргумента, vy – вектор исходных значений функции, vs – вектор вторых производных, возвращаемых одной из функций cspline (vx,vy), pspline (vx,vy), lspline (vx,vy), которые возвращают вектор vs при использовании кубических, параболических и линейных сплайнов соответственно.

Для интерполирования линиями можно воспользоваться функцией linterp(vx, vy, x).

Конечно, используя вычислительные возможности MathCad, всегда можно воспользоваться другим методом интерполяции (например, с помощью полиномов Ньютона и Лагранжа).

Пример.

Имеются экспериментально полученные точки:

 

Значение аргумента x Значение функции y
-1.06 -0.837 -0.684 -0.315 -0.117 -0.0 0.115 0.5 1.22 0.854 0.513 0.271 0.217 0.198 0.218 0.277

Построить график функции y(x) используя точки, отличные от экспериментально полученных.

Решение.

 

 

 

Аппроксимация зависимостей с помощью MathCad

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

regress(vx,vy, k) – возвращает вектор, который использует функция interp для нахождения полинома степени k, который наилучшим образом приближает значения x и y данных, хранящихся в векторах vx, vy.

interp (vs,vx,vy, x) – возвращает приближенное значение y, соответствующее значению x. Где vs – вектор, получаемый с помощью regress.

Пример.

Имеются экспериментально полученные точки:

 

Значение аргумента x Значение функции y
  1,5 3,7

Определить вид аппроксимирующего полинома 4-й степени и с его помощью построить график функции y(x).

Решение.

Решим задачу методом наименьших квадратов, а затем с помощью встроенных функций MathCad.

 

Введем исходные данные

 

Вычислим коэффициенты системы уравнений

 

 

Определим матрицу коэффициентов и матрицу свободных членов системы
Найдем коэффициенты полинома, решив систему

 

Определим вид полинома

 

 

 

 

Построим график Р(х)
Во всех точках х для сравнения вычислим значения функции с помощью встроенных средств аппроксимации

 

 

Строим графики

 

Как видно, графики совпадают.

 

Вычисление определенного интеграла
и производной средствами MathCad

Предположим нужно вычислить интеграл функции на отрезке [0; 2] и производную этой функции в точке х=1. Решение такой задачи в MathCad имеет вид:

 

 

 

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

Для вычисления производных более высоких порядков используется специальная конструкция:


Решение ОДУ средствами MathCad

Для решения линейных обыкновенных дифференциальных уравнений в MathCad функция odesolve(x,b,[step]), которая возвращает значение функции, зависящей от х и являющейся решением линейного ОДУ. Здесь х – аргумент, b – правый конец отрезка, [step] – необязательный параметр, количество шагов для нахождения решения. Количество начальных условий должно равняться порядку уравнения.

Пример. Решим уравнение на отрезке [0, 4π]. Известно, что y(0)=0, а y’(0)=1.

Решение:

Примечание: В блоке Given используются знаки булева равенства, вставляемые нажатием Ctrl+= и знак производной (штрих), вставляемый нажатием Ctrl+F7.

Решение ОДУ первого порядка вида может быть получено с помощью функций rkfixed(y0,a,b,n,D), которая возвращает матрицу, состоящую из двух столбцов. В первом столбце хранятся значения аргумента, во втором – функции (результаты решения). Здесь y0 – начальное значение функции y, а – начало отрезка, b – конец отрезка, n – количество отрезков разбиения, D – первая производная от у.

Пример. Решим уравнение на отрезке [0, π]. Известно, что y(0)=1.

Решение:

 

 

Для решения ОДУ первого порядка так же можно использовать функцию rkadapt(), аналогичную рассмотренной выше rkfixed(), за исключением того, что решение находится не с фиксированным шагом, а с автоматическим его подбором.

 

 

 





Поделиться с друзьями:


Дата добавления: 2017-02-25; Мы поможем в написании ваших работ!; просмотров: 1414 | Нарушение авторских прав


Поиск на сайте:

Лучшие изречения:

Велико ли, мало ли дело, его надо делать. © Неизвестно
==> читать все изречения...

2489 - | 2155 -


© 2015-2024 lektsii.org - Контакты - Последнее добавление

Ген: 0.008 с.