Дифференцирование

Открыть сценарий в реальном времени

В этом примере показано, как аналитически находить и оценивать производные с помощью символьных математических Toolbox™. В примере вы найдете 1-ю и 2-ю производные f (x) и используйте эти производные, чтобы найти локальные максимумы, минимумы и точки перегиба.

Первые производные: поиск локальных минимумов и максимумов

Вычисление первой производной выражения помогает найти локальные минимумы и максимумы этого выражения. Перед созданием символьного выражения объявите символьные переменные:

syms x

По умолчанию в результаты включаются решения, включающие мнимые компоненты. Здесь рассмотрим только реальные значения x установив предположение, что x является реальным:

assume(x, 'real')

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

f = (3*x^3 + 17*x^2 + 6*x + 1)/(2*x^3 - x + 3)

f = 
 $\frac{3 x^{3} + 17 x^{2} + 6 x + 1}{2 x^{3} - x + 3}$

Печать этого выражения показывает, что выражение имеет горизонтальные и вертикальные асимптоты, локальный минимум от -1 до 0 и локальный максимум от 1 до 2:

fplot(f)
grid

Figure contains an axes. The axes contains an object of type functionline.

Чтобы найти горизонтальную асимптоту, вычислите пределы f для x приближение положительных и отрицательных бесконечностей. Горизонтальная асимптота y = 3/2:

lim_left = limit(f, x, -inf)

lim_left = 
 $\frac{3}{2}$

lim_right = limit(f, x, inf)

lim_right = 
 $\frac{3}{2}$

Добавьте горизонтальную асимптоту к графику:

hold on
plot(xlim, [lim_right lim_right], 'LineStyle', '-.', 'Color', [0.25 0.25 0.25])

Figure contains an axes. The axes contains 2 objects of type functionline, line.

Поиск вертикальной асимптоты f, найти полюса f:

pole_pos = poles(f, x)

pole_pos = 
 $- \frac{1}{6 {(\frac{3}{4} - \frac{\sqrt{241} \sqrt{432}}{432})}^{1 / 3}} - {(\frac{3}{4} - \frac{\sqrt{241} \sqrt{432}}{432})}^{1 / 3}$

Аппроксимировать точное решение численно с помощью double функция:

double(pole_pos)

ans = -1.2896

Теперь найдите локальный минимум и максимум f. Если точка является локальным экстремумом (минимальным или максимальным), первая производная выражения в этой точке равна нулю. Вычислить производную f использование diff:

g = diff(f, x)

g = 
 $\frac{9 x^{2} + 34 x + 6}{2 x^{3} - x + 3} - \frac{(6 x^{2} - 1) (3 x^{3} + 17 x^{2} + 6 x + 1)}{{(2 x^{3} - x + 3)}^{2}}$

Найти местные крайности f, решить уравнение g == 0:

g0 = solve(g, x)

g0 = 
 $\begin{array}{l} (\begin{array}{c} \frac{\sqrt{σ_{2}}}{6 {σ_{3}}^{1 / 6}} - σ_{1} - \frac{15}{68} \\ \frac{\sqrt{σ_{2}}}{6 {σ_{3}}^{1 / 6}} + σ_{1} - \frac{15}{68} \end{array}) \\ where \\ σ_{1} = \frac{\sqrt{\frac{337491 \sqrt{6} \sqrt{\frac{3 \sqrt{3} \sqrt{178939632355}}{9826} + \frac{2198209}{9826}}}{39304} + \frac{2841 {σ_{3}}^{1 / 3} \sqrt{σ_{2}}}{578} - 9 {σ_{3}}^{2 / 3} \sqrt{σ_{2}} - \frac{361 \sqrt{σ_{2}}}{289}}}{6 {σ_{3}}^{1 / 6} {σ_{2}}^{1 / 4}} \\ σ_{2} = \frac{2841 {σ_{3}}^{1 / 3}}{1156} + 9 {σ_{3}}^{2 / 3} + \frac{361}{289} \\ σ_{3} = \frac{\sqrt{3} \sqrt{178939632355}}{176868} + \frac{2198209}{530604} \end{array}$

Аппроксимировать точное решение численно с помощью double функция:

double(g0)

ans = 2×1

   -0.1892
    1.2860

Выражение f имеет локальный максимум при x = 1.286 и локальный минимум в x = -0.189. Получение значений функций в этих точках с помощью subs:

f0 = subs(f,x,g0)

f0 = 
 $\begin{array}{l} (\begin{array}{c} \frac{3 σ_{2} - 17 {(σ_{5} - σ_{6} + \frac{15}{68})}^{2} - σ_{4} + σ_{1} + \frac{11}{34}}{σ_{6} + 2 σ_{2} - σ_{5} - \frac{219}{68}} \\ - \frac{σ_{4} + 17 {(σ_{6} + σ_{5} - \frac{15}{68})}^{2} + 3 σ_{3} + σ_{1} - \frac{11}{34}}{σ_{6} - 2 σ_{3} + σ_{5} - \frac{219}{68}} \end{array}) \\ where \\ σ_{1} = \frac{σ_{7}}{{σ_{9}}^{1 / 6} {σ_{8}}^{1 / 4}} \\ σ_{2} = {(σ_{5} - σ_{6} + \frac{15}{68})}^{3} \\ σ_{3} = {(σ_{6} + σ_{5} - \frac{15}{68})}^{3} \\ σ_{4} = \frac{\sqrt{σ_{8}}}{{σ_{9}}^{1 / 6}} \\ σ_{5} = \frac{σ_{7}}{6 {σ_{9}}^{1 / 6} {σ_{8}}^{1 / 4}} \\ σ_{6} = \frac{\sqrt{σ_{8}}}{6 {σ_{9}}^{1 / 6}} \\ σ_{7} = \sqrt{\frac{337491 \sqrt{6} \sqrt{\frac{3 \sqrt{3} \sqrt{178939632355}}{9826} + \frac{2198209}{9826}}}{39304} + \frac{2841 {σ_{9}}^{1 / 3} \sqrt{σ_{8}}}{578} - 9 {σ_{9}}^{2 / 3} \sqrt{σ_{8}} - \frac{361 \sqrt{σ_{8}}}{289}} \\ σ_{8} = \frac{2841 {σ_{9}}^{1 / 3}}{1156} + 9 {σ_{9}}^{2 / 3} + \frac{361}{289} \\ σ_{9} = \frac{\sqrt{3} \sqrt{178939632355}}{176868} + \frac{2198209}{530604} \end{array}$

Аппроксимировать точное решение численно с помощью double функция на переменной f0:

double(f0)

ans = 2×1

    0.1427
    7.2410

Добавьте маркеры точек на график в крайнем случае:

plot(g0, f0, 'ok')

Figure contains an axes. The axes contains 3 objects of type functionline, line.

Вторые производные: поиск точек перегиба

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

h = diff(f, x, 2)

h = 
 $\begin{array}{l} \frac{18 x + 34}{σ_{1}} - \frac{2 (6 x^{2} - 1) (9 x^{2} + 34 x + 6)}{{σ_{1}}^{2}} - \frac{12 x σ_{2}}{{σ_{1}}^{2}} + \frac{2 {(6 x^{2} - 1)}^{2} σ_{2}}{{σ_{1}}^{3}} \\ where \\ σ_{1} = 2 x^{3} - x + 3 \\ σ_{2} = 3 x^{3} + 17 x^{2} + 6 x + 1 \end{array}$

Теперь упростите этот результат:

h = simplify(h)

h = 
 $\frac{2 (68 x^{6} + 90 x^{5} + 18 x^{4} - 699 x^{3} - 249 x^{2} + 63 x + 172)}{{(2 x^{3} - x + 3)}^{3}}$

Поиск точек перегиба f, решить уравнение h = 0. Здесь используйте числовой решатель vpasolve для вычисления аппроксимаций решений с плавающей запятой:

h0 = vpasolve(h, x)

h0 = 
 $(\begin{array}{c} 0.57871842655441748319601085860196 \\ 1.8651543689917122385037075917613 \\ - 1.4228127856020972275345064554049 - 1.8180342567480118987898749770461 i \\ - 1.4228127856020972275345064554049 + 1.8180342567480118987898749770461 i \\ - 0.46088831805332057449182335801198 + 0.47672261854520359440077796751805 i \\ - 0.46088831805332057449182335801198 - 0.47672261854520359440077796751805 i \end{array})$

Выражение f имеет две точки перегиба: x = 1.865 и x = 0.579. Обратите внимание, что vpasolve также возвращает сложные решения. Отменить следующие действия:

h0(imag(h0)~=0) = []

h0 = 
 $(\begin{array}{c} 0.57871842655441748319601085860196 \\ 1.8651543689917122385037075917613 \end{array})$

Добавьте маркеры на график, показывающий точки перегиба:

plot(h0, subs(f,x,h0), '*k')
hold off

Figure contains an axes. The axes contains 4 objects of type functionline, line.

Документация

Дифференцирование

Первые производные: поиск локальных минимумов и максимумов

Вторые производные: поиск точек перегиба

Документация по инструментам символьной математики

Поддержка