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

Открыть Live Script

Этот пример показывает, как аналитически найти и вычислить производные с помощью Symbolic Math 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.

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

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

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

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

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

Поддержка