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

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

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

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

Пример.

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

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

Определить вид аппроксимирующего полинома 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(y₀,a,b,n,D), которая возвращает матрицу, состоящую из двух столбцов. В первом столбце хранятся значения аргумента, во втором – функции (результаты решения). Здесь y₀ – начальное значение функции y, а – начало отрезка, b – конец отрезка, n – количество отрезков разбиения, D – первая производная от у.

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

Решение:

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

43