MATLAB: эффективное вычисление производных полиномов

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

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

Производные полиномов

Полином — это математическое выражение, состоящее из переменных, возведенных в разные степени, объединенных с коэффициентами. Например, полином P(x) = 3x² + 2x + 5 — это полином второй степени.

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

Например, рассмотрим полином:

P(x) = 3x³ + 4x² + 2x + 1

Производная P'(x) вычисляется как:

• Для члена 3x³: умножьте 3 на 3 (показатель степени), в результате 9x².
• Для члена 4x²: умножьте 4 на 2, в результате 8x.
• Для члена 2x: умножьте 2 на 1, в результате 2.
• Константный член (1) имеет производную 0.

Таким образом, производная:

P'(x) = 9x² + 8x + 2

Производные в MATLAB

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

Для нахождения производной полинома MATLAB предоставляет функцию polyder.

Синтаксис

k = polyder(p)
k = polyder(a,b)
[q,d] = polyder(a,b)

Пояснение синтаксиса

k = polyder(p) вычисляет производную полинома, заданного коэффициентами в p, в результате новый полином k(x), представляющий производную d/dx p(x).

k = polyder(a,b) вычисляет производную произведения двух полиномов a и b, в результате новый полином k(x), представляющий d/dx[a(x)·b(x)].

[q, d] = polyder(a, b) вычисляет производную частного двух полиномов a и b, возвращая два полинома: q(x) (числитель) и d(x) (знаменатель), представляющие производную a(x)/b(x).

Пример 1: Вычисление производной с помощью polyder(p)

Рассмотрим полином:

P(x) = 4x³ + 3x² + 2x + 1

Этот полином можно представить вектором его коэффициентов в MATLAB:

p = [4 3 2 1];

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

k = polyder(p);

При выполнении кода в командном окне MATLAB вывод:

>> p = [4 3 2 1]; k = polyder(p)
k = 12 6 2 >>

Для члена 4x³ производная 12x² (умножьте коэффициент 4 на показатель 3 и уменьшите показатель на 1).

Для члена 3x² производная 6x.

Для члена 2x производная 2.

Константный член 1 имеет производную 0.

Таким образом, производный полином:

k(x) = 12x² + 6x + 2

В MATLAB результат k будет: [12 6 2]

Пример 2: Еще один пример нахождения производных полинома

Рассмотрим следующий полином:

p(x) = 5x⁴ - 2x³ + 7x² - 3x + 8

Этот полином можно представить вектором его коэффициентов в MATLAB:

p = [5 -2 7 -3 8]

Для нахождения производной этого полинома используем функцию polyder в MATLAB:

k = polyder(p)

Эта команда вернет коэффициенты производной полинома p.

При выполнении кода в командном окне MATLAB вывод:

>> p = [5 -2 7 -3 8]; k = polyder(p)
k = 20 -6 14 -3 >>

Вектор k = [20 -6 14 -3] представляет полином:

k(x) = 20x³ - 6x² + 14x - 3

Пример 3: Производная произведения двух полиномов с помощью polyder(a, b)

Рассмотрим два полинома:

a(x) = 2x² + 3x + 1
b(x) = 4x + 5

Эти полиномы можно представить векторами их коэффициентов в MATLAB:

a = [2 3 1]
b = [4 5]

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

k = polyder(a, b);

Это вернет коэффициенты производной произведения a(x) и b(x).

При выполнении кода в командном окне MATLAB вывод:

>> a = [2 3 1]; b = [4 5]; k = polyder(a, b)
k = 24 44 19 >>

Таким образом, производный полином: k(x) = 24x² + 44x + 19

Пример 4: Производная двух заданных полиномов

Рассмотрим два разных полинома:

a(x) = 3x³ + 2x² + x + 4
b(x) = x² - 5x + 6

Эти полиномы можно представить следующими векторами коэффициентов в MATLAB:

a = [3 2 1 4];
b = [1 -5 6];

Для вычисления производной произведения этих двух полиномов используем функцию polyder с векторами a и b в качестве входов:

k = polyder(a, b);

Эта команда вернет коэффициенты производной произведения a(x) и b(x).

При выполнении кода в командном окне MATLAB вывод:

>> a = [3 2 1 4]; b = [1 -5 6]; k = polyder(a, b)
k = 15 -52 27 22 -14 >>

Таким образом, производный полином:

k(x) = 15x⁴ - 52x³ + 27x² + 22x - 14

Пример 5: Производная частного двух полиномов с помощью [q, d] = polyder(a, b)

Рассмотрим два полинома:

a(x) = 4x² + 3x + 2
b(x) = x² - 2x + 1

Эти полиномы можно представить векторами их коэффициентов в MATLAB:

a = [4 3 2];
b = [1 -2 1];

Для вычисления производной частного a(x) / b(x) используем функцию polyder с двумя выходными аргументами q и d:

[q,d] = polyder(a,b)

Это вернет два полинома: q(x) (числитель) и d(x) (знаменатель) производной a(x) / b(x).

При выполнении кода в командном окне MATLAB вывод:

>> a = [4 3 2]; b = [1 -2 1]; [q,d] = polyder(a,b)
q = -11 4 7
d = 1 -4 6 -4 1 >>