Частные производные

Частная производная — это обычная производная, если смотреть на функцию вдоль одной оси, заморозив остальные.

Частная производная $\dfrac{\partial f}{\partial x}$ — скорость изменения $f$ при изменении только $x$, когда остальные переменные держат постоянными.

Идея заморозки

У функции $f(x,y)$ две «ручки настройки». Чтобы понять влияние одной, фиксируем другую. Производная по $x$ при фиксированном $y$:

$$\frac{\partial f}{\partial x} = \lim_{h \to 0} \frac{f(x+h,\,y) - f(x,\,y)}{h}$$

На практике дифференцируют по обычным правилам, считая $y$ константой. Для $f = x^2 y + \sin x$:

$$\frac{\partial f}{\partial x} = 2xy + \cos x, \qquad \frac{\partial f}{\partial y} = x^2$$

Геометрия

$\partial f/\partial x$ — это наклон сечения поверхности плоскостью $y = \text{const}$. Вы идёте по рельефу строго на восток и измеряете крутизну подъёма.

Численная проверка (центральная разность)

import math

def f(x, y):
    return x*x*y + math.sin(x)

def partial_x(f, x, y, h=1e-6):
    return (f(x+h, y) - f(x-h, y)) / (2*h)

def partial_y(f, x, y, h=1e-6):
    return (f(x, y+h) - f(x, y-h)) / (2*h)

x0, y0 = 1.0, 2.0
print("df/dx числ. =", round(partial_x(f, x0, y0), 6))
print("df/dx точн. =", round(2*x0*y0 + math.cos(x0), 6))
print("df/dy числ. =", round(partial_y(f, x0, y0), 6))
print("df/dy точн. =", round(x0*x0, 6))

Вывод:

df/dx числ. = 4.540302
df/dx точн. = 4.540302
df/dy числ. = 1.0
df/dy точн. = 1.0

Как работает под капотом

Численная центральная разность $\frac{f(x+h)-f(x-h)}{2h}$ точнее односторонней: её ошибка пропорциональна $h^2$, а не $h$. Это потому, что симметричная схема гасит линейный член ошибки в разложении Тейлора. Поэтому в коде мы и используем именно её — совпадение с аналитикой до шестого знака не случайно.

Частые ошибки

  • Дифференцировать по $x$, забыв заморозить $y$ (например, дифференцировать $y$ как переменную).
  • Брать слишком малое $h$ в численной схеме — машинная погрешность съест точность.
  • Путать $\partial$ (частная) и $d$ (полная) производную — это разные операции.

Итог

  • Частная производная по $x$ — производная при замороженных остальных переменных.
  • Геометрически это наклон сечения поверхности.
  • Численно её надёжно даёт центральная разность с ошибкой $\sim h^2$.
Проверьте себя
1. Как вычислить ∂f/∂x для f(x,y)?
AДифференцировать по x, считая y константой
BДифференцировать одновременно по x и y
CДифференцировать по y, считая x константой
DВзять корень из f
2. Для f(x,y)=x²y чему равна ∂f/∂y?
A2xy
B
C2x
D
3. Почему центральная разность точнее односторонней?
AОна использует меньше точек
BЕё ошибка ~h² вместо ~h
CОна не требует производной
DОна работает только для полиномов