3. Сложные и неявные функции

Как дифференцировать функцию, когда её аргументы сами от чего-то зависят (цепное правило), и как найти производные функции, которую никто не выразил явно, а только задал уравнением. Это задачи №3 и №4 КР1.

на КР1: задачи 3 и 4практика 28.09практика 05.10

Что известно про КР

В составе КР1 (сообщение в чате практики от 05.10) есть «сложная функция» и «неявная функция». По словам одногруппника, на КР будут задачи как в правой части доски 05.10: найти частные производные 1 и 2 порядка функции $z=z(x,y)$, заданной уравнением $x^2+y^2+z^2=a^2$ или $z^3-3xyz=a^3$. Обе решены ниже в разделе 8. В примере КР1 из РПД курса есть ещё «найти дифференциал неявной функции в точке» и «проверить тождество со вторыми производными» — тоже разобраны.

Часть А Сложные функции

1. Цепное правило и дерево зависимостей

Сложная функция — это когда $u$ зависит от $x, y$, а $x$ и $y$ сами зависят от других переменных. Например, температура $u=f(x,y)$ в точке комнаты, а муха летит по траектории $x=x(t)$, $y=y(t)$. Тогда температура, которую чувствует муха, — функция времени $u(t)=f(x(t),y(t))$, и вопрос: чему равна $\dfrac{du}{dt}$?

Ответ складывается из двух вкладов: $u$ меняется, потому что меняется $x$, и потому что меняется $y$. Каждый вклад — произведение «скорость $u$ по $x$» на «скорость $x$ по $t$». Удобно рисовать дерево зависимостей: сверху функция, ниже промежуточные переменные, внизу независимые.

uxytt uxy stst uxyz xy ∂u/∂x∂u/∂y dx/dtdy/dt ∂u/∂z∂z/∂y∂z/∂x u = f(x, y), x(t), y(t) u = f(x, y), x(s, t), y(s, t) u = f(x, y, z), z = z(x, y)
Дерево зависимостей. Производная по нижней переменной = сумма по всем путям от u до неё, на каждом пути производные перемножаются
Цепное правило

Один параметр. Если $u=f(x,y)$ дифференцируема, а $x=x(t)$, $y=y(t)$ дифференцируемы, то

$$\frac{du}{dt}=\pd{u}{x}\,\frac{dx}{dt}+\pd{u}{y}\,\frac{dy}{dt}.$$

Два параметра. Если $x=x(s,t)$, $y=y(s,t)$, то

$$\pd{u}{s}=\pd{u}{x}\pd{x}{s}+\pd{u}{y}\pd{y}{s},\qquad \pd{u}{t}=\pd{u}{x}\pd{x}{t}+\pd{u}{y}\pd{y}{t}.$$

Общий случай. $u=f(x_1,\dots,x_n)$, каждая $x_i=x_i(t_1,\dots,t_m)$:

$$\pd{u}{t_j}=\sum_{i=1}^{n}\pd{f}{x_i}\,\pd{x_i}{t_j},\qquad j=1,\dots,m.$$

Условие: $f$ дифференцируема в соответствующей точке (достаточно непрерывных частных производных), внутренние функции имеют производные в рассматриваемой точке.

Как не запутаться
  1. Нарисуй дерево (хотя бы в уме): какие промежуточные переменные и от чего они зависят.
  2. Для $\partial u/\partial t$ найди все пути от $u$ до $t$. Сколько путей — столько слагаемых.
  3. На каждом пути перемножь производные по рёбрам.
Связь с программированием

В матричном виде общее правило — произведение матриц Якоби: $J_{u\circ x}=J_u\cdot J_x$. Ровно это делает backpropagation в нейросетях: градиент идёт по графу вычислений от выхода к входам, по каждому ребру домножается на локальную производную, а по сходящимся путям складывается.

Пример 1

$z=x^2y$, $x=\cos t$, $y=\sin t$. Найти $\dfrac{dz}{dt}$.

Решение

$z_x=2xy$, $z_y=x^2$, $x'(t)=-\sin t$, $y'(t)=\cos t$. По цепному правилу

$$\frac{dz}{dt}=2xy\cdot(-\sin t)+x^2\cos t=-2\cos t\sin^2 t+\cos^3 t.$$

Проверка подстановкой: $z=\cos^2 t\,\sin t$, $z'=-2\cos t\sin t\cdot\sin t+\cos^2t\cos t$ — то же самое.

Ответ: $\dfrac{dz}{dt}=\cos^3 t-2\cos t\sin^2 t$.

Пример 2

$z=x^2-y^2$, $x=s+t$, $y=st$. Найти $z_s$, $z_t$.

Решение

$z_x=2x$, $z_y=-2y$; $x_s=1$, $x_t=1$, $y_s=t$, $y_t=s$.

$$z_s=2x\cdot1-2y\cdot t=2(s+t)-2st^2,\qquad z_t=2x\cdot 1-2y\cdot s=2(s+t)-2s^2t.$$

В ответе выражаем всё через $s,t$, если не сказано иное.

Полная и частная производная

Самый коварный случай — третье дерево: $u=f(x,y,z)$, а $z=z(x,y)$. Переменная $x$ входит в $u$ и напрямую, и через $z$. Тогда есть две разные «производные по $x$»:

$$\frac{\partial u}{\partial x}=\pd{f}{x}+\pd{f}{z}\,\pd{z}{x}.$$

Если $y=y(x)$ и $z=z(x)$ и всё зависит от одного $x$, то $u$ — функция одной переменной, и пишут полную производную $\dfrac{du}{dx}=\pd{f}{x}+\pd{f}{y}\dfrac{dy}{dx}+\pd{f}{z}\dfrac{dz}{dx}$.

Пример 3

$u=x^2+yz$, где $z=xy$. Найти частную производную $f$ по $x$ при замороженном $z$ и производную сложной функции $u(x,y)$ по $x$.

Решение

$f(x,y,z)=x^2+yz$: $\pd{f}{x}=2x$ (здесь $z$ считается константой).

С учётом зависимости $z$ от $x$: $\pd{u}{x}=\pd{f}{x}+\pd{f}{z}\pd{z}{x}=2x+y\cdot y=2x+y^2$.

Проверка: $u=x^2+xy^2$, $u_x=2x+y^2$. Видно, что «просто частная» $2x$ — не то же самое.

Пример 4. Производная $u=\varphi(t)^{\psi(t)}$

вопрос экзамена №3 в РПД Найти $u'(t)$ для $u=\varphi(t)^{\psi(t)}$ и, в частности, для $u=(1+t^2)^{\sin t}$.

Решение

Положим $u=x^y$, где $x=\varphi(t)\gt 0$, $y=\psi(t)$. Частные производные: $\pd{u}{x}=y\,x^{y-1}$ (степенная по $x$), $\pd{u}{y}=x^y\ln x$ (показательная по $y$). По цепному правилу

$$u'=y\,x^{y-1}\varphi'+x^y\ln x\cdot\psi'=\varphi^{\psi}\left(\psi'\ln\varphi+\psi\frac{\varphi'}{\varphi}\right).$$

Это та же формула, что дало бы логарифмическое дифференцирование на первом курсе, но теперь видно, откуда два слагаемых: два пути в дереве $u\to x\to t$ и $u\to y\to t$.

Для $u=(1+t^2)^{\sin t}$: $\varphi=1+t^2$, $\psi=\sin t$:

$$u'=(1+t^2)^{\sin t}\left(\cos t\,\ln(1+t^2)+\frac{2t\sin t}{1+t^2}\right).$$

2. Задачи с абстрактной f

На практике чаще всего функция $f$ не дана, а дано только, как устроены её аргументы: «$u=f(x+y,\,xy)$, найти $u''_{xy}$». Ответ тогда выражается через производные $f$.

Обозначения

Пусть $u=f(\xi,\eta)$, где $\xi=\xi(x,y)$, $\eta=\eta(x,y)$. Обозначим

$$f'_1=\pd{f}{\xi}(\xi,\eta),\quad f'_2=\pd{f}{\eta}(\xi,\eta),\quad f''_{11}=\frac{\partial^2 f}{\partial\xi^2},\quad f''_{12}=\frac{\partial^2 f}{\partial\xi\,\partial\eta},\quad f''_{22}=\frac{\partial^2 f}{\partial\eta^2}.$$

Индекс — номер аргумента, а не переменная $x$. Все они вычисляются в точке $(\xi(x,y),\eta(x,y))$. Ещё пишут $f_u, f_v$ или $f_\xi, f_\eta$ — то же самое. Считаем $f$ дважды непрерывно дифференцируемой, поэтому $f''_{12}=f''_{21}$.

Главная ловушка второй производной

$f'_1$ — это не число и не константа. Это функция тех же аргументов $(\xi,\eta)$, а значит, сама сложная функция от $x,y$. Поэтому её производная по $x$ снова считается по цепному правилу:

$$\frac{\partial}{\partial x}\bigl(f'_1\bigr)=f''_{11}\,\xi_x+f''_{12}\,\eta_x,\qquad \frac{\partial}{\partial x}\bigl(f'_2\bigr)=f''_{21}\,\xi_x+f''_{22}\,\eta_x.$$

Частая ошибка — написать $\frac{\partial}{\partial x}f'_1=f''_{11}$ и потерять второе слагаемое.

f′₁ξηyy f′₂ξηyy f″₁₁f″₁₂1x f″₂₁f″₂₂1x
Для u = f(x+y, xy): f′₁ и f′₂ сами зависят от ξ = x+y и η = xy, поэтому (f′₁)′_y = f″₁₁·1 + f″₁₂·x, (f′₂)′_y = f″₂₁·1 + f″₂₂·x
Алгоритм для вторых производных
  1. Найди первые производные по цепному правилу: $u_x=f'_1\xi_x+f'_2\eta_x$.
  2. Дифференцируй результат как произведение: $(f'_1\cdot\xi_x)'_y=(f'_1)'_y\,\xi_x+f'_1\,\xi_{xy}$.
  3. Каждое $(f'_k)'_y$ раскрой по цепному правилу ещё раз: $(f'_k)'_y=f''_{k1}\xi_y+f''_{k2}\eta_y$.
  4. Собери подобные, используя $f''_{12}=f''_{21}$.

Задачи с доски 28.09

На фото доски видны пункты 1) и 3) первой задачи, задача «Найти» из двух пунктов и пять пунктов на дифференциалы. Запись $u=u(x^2+y^2+z^2)$ на доске означает: $u$ — функция одного аргумента $t=x^2+y^2+z^2$; ниже пишем $u=f(t)$, чтобы не путать функцию с её значением.

Пример 5. $u=f(x^2+y^2+z^2)$

с практики 28.09 Найти производные 1 и 2 порядка.

Решение

$t=x^2+y^2+z^2$, $t_x=2x$, $t_y=2y$, $t_z=2z$. Первые производные:

$$u_x=2x\,f'(t),\qquad u_y=2y\,f'(t),\qquad u_z=2z\,f'(t).$$

$u_{xx}$: дифференцируем произведение $2x\cdot f'(t)$, причём $f'(t)$ — сложная функция, $(f'(t))'_x=f''(t)\cdot t_x=2x f''(t)$:

$$u_{xx}=2f'(t)+2x\cdot 2xf''(t)=2f'+4x^2f''.$$

Смешанная: $u_{xy}=(2xf'(t))'_y=2x\cdot f''(t)\cdot 2y=4xyf''$. Остальные — по симметрии.

Ответ (всё в точке $t=x^2+y^2+z^2$):

$$u_{xx}=2f'+4x^2f'',\quad u_{yy}=2f'+4y^2f'',\quad u_{zz}=2f'+4z^2f'',$$ $$u_{xy}=4xyf'',\quad u_{xz}=4xzf'',\quad u_{yz}=4yzf''.$$
Пример 6. $u=f(x,\,xy,\,xyz)$

с практики 28.09 Найти производные 1 и 2 порядка.

Решение

Аргументы: $\xi=x$, $\eta=xy$, $\zeta=xyz$. Дерево:

uξηζ xxy xyz f′₁f′₂f′₃ 1yx yzxzxy
Пути от u до x выделены: их три, поэтому в u_x три слагаемых

Первые производные.

$$u_x=f'_1+y\,f'_2+yz\,f'_3,\qquad u_y=x\,f'_2+xz\,f'_3,\qquad u_z=xy\,f'_3.$$

Вспомогательные формулы. Для любого $k$ функция $f'_k$ зависит от $x,y,z$ так же, как $u$, поэтому

$$(f'_k)'_x=f''_{k1}+y f''_{k2}+yz f''_{k3},\quad (f'_k)'_y=x f''_{k2}+xz f''_{k3},\quad (f'_k)'_z=xy f''_{k3}.$$

$u_{xx}$. $u_{xx}=(f'_1)'_x+y(f'_2)'_x+yz(f'_3)'_x$ (множители $y$, $yz$ от $x$ не зависят):

$$u_{xx}=f''_{11}+2yf''_{12}+2yzf''_{13}+y^2f''_{22}+2y^2zf''_{23}+y^2z^2f''_{33}.$$

$u_{xy}$. Дифференцируем $u_x$ по $y$; теперь множители $y$ и $yz$ тоже зависят от $y$, отсюда слагаемые $f'_2$ и $zf'_3$:

$$u_{xy}=(f'_1)'_y+f'_2+y(f'_2)'_y+zf'_3+yz(f'_3)'_y$$ $$=xf''_{12}+xzf''_{13}+xyf''_{22}+2xyzf''_{23}+xyz^2f''_{33}+f'_2+zf'_3.$$

$u_{xz}$. $u_{xz}=(f'_1)'_z+y(f'_2)'_z+yf'_3+yz(f'_3)'_z$:

$$u_{xz}=xyf''_{13}+xy^2f''_{23}+xy^2zf''_{33}+yf'_3.$$

$u_{yy}$, $u_{yz}$, $u_{zz}$.

$$u_{yy}=x^2f''_{22}+2x^2zf''_{23}+x^2z^2f''_{33},$$ $$u_{yz}=x^2yf''_{23}+x^2yzf''_{33}+xf'_3,\qquad u_{zz}=x^2y^2f''_{33}.$$

Все производные $f$ берутся в точке $(x,\,xy,\,xyz)$.

Пример 7. $u=f(x+y,\,xy)$

с практики 28.09 Найти $\dfrac{\partial^2u}{\partial x\,\partial y}$.

Решение

$\xi=x+y$, $\eta=xy$; $\xi_x=\xi_y=1$, $\eta_x=y$, $\eta_y=x$.

$u_x=f'_1\cdot 1+f'_2\cdot y$. Дифференцируем по $y$ (картинка выше):

$$u_{xy}=(f'_1)'_y+f'_2+y\,(f'_2)'_y=(f''_{11}+xf''_{12})+f'_2+y(f''_{21}+xf''_{22}).$$

Ответ: $u_{xy}=f''_{11}+(x+y)f''_{12}+xyf''_{22}+f'_2$.

Пример 8. Лапласиан

с практики 28.09 Найти $\Delta u=u_{xx}+u_{yy}+u_{zz}$ для $u=f(x+y+z,\;x^2+y^2+z^2)$.

Решение

$\xi=x+y+z$, $\eta=x^2+y^2+z^2$; $\xi_x=1$, $\eta_x=2x$.

$u_x=f'_1+2xf'_2$. Отсюда, с $(f'_k)'_x=f''_{k1}+2xf''_{k2}$:

$$u_{xx}=f''_{11}+2xf''_{12}+2f'_2+2x(f''_{21}+2xf''_{22})=f''_{11}+4xf''_{12}+4x^2f''_{22}+2f'_2.$$

$u_{yy}$ и $u_{zz}$ — то же самое с заменой $x$ на $y$ и $z$. Складываем:

$$\Delta u=3f''_{11}+4(x+y+z)f''_{12}+4(x^2+y^2+z^2)f''_{22}+6f'_2.$$

Можно записать красиво через сами аргументы: $\Delta u=3f''_{11}+4\xi f''_{12}+4\eta f''_{22}+6f'_2$.

Дифференциалы 1 и 2 порядка (вторая половина доски 28.09)

Напоминание: если $x,y$ — независимые переменные, то $du=u_xdx+u_ydy$, $d^2u=u_{xx}dx^2+2u_{xy}dx\,dy+u_{yy}dy^2$. Для $u=f(t)$, где $t=t(x,y)$, быстрее считать по формулам (почему они такие — в разделе 3):

$$du=f'(t)\,dt,\qquad d^2u=f''(t)\,dt^2+f'(t)\,d^2t,$$

где $dt$ и $d^2t$ — дифференциалы функции $t(x,y)$ по независимым $x,y$.

Пример 9. $u=f(t)$, $t=x+y$

с практики 28.09 $x,y$ — независимые. Найти $du$, $d^2u$.

Решение

$dt=dx+dy$, $d^2t=0$ (функция линейная). Поэтому

$$du=f'(x+y)(dx+dy),\qquad d^2u=f''(x+y)(dx+dy)^2.$$

Проверка через производные: $u_{xx}=u_{xy}=u_{yy}=f''$, и $d^2u=f''(dx^2+2dxdy+dy^2)$ — совпадает.

Пример 10. $u=f(t)$, $t=y/x$

с практики 28.09 Найти $du$, $d^2u$.

Решение

$t_x=-\dfrac{y}{x^2}$, $t_y=\dfrac1x$, $dt=\dfrac{x\,dy-y\,dx}{x^2}$.

$t_{xx}=\dfrac{2y}{x^3}$, $t_{xy}=-\dfrac1{x^2}$, $t_{yy}=0$, значит $d^2t=\dfrac{2y}{x^3}dx^2-\dfrac{2}{x^2}dx\,dy=\dfrac{2(y\,dx^2-x\,dx\,dy)}{x^3}$.

$$du=f'\!\left(\tfrac yx\right)\frac{x\,dy-y\,dx}{x^2},\qquad d^2u=f''\!\left(\tfrac yx\right)\frac{(x\,dy-y\,dx)^2}{x^4}+f'\!\left(\tfrac yx\right)\frac{2(y\,dx^2-x\,dx\,dy)}{x^3}.$$
Пример 11. $u=f\bigl(\sqrt{x^2+y^2}\bigr)$

с практики 28.09 Найти $du$, $d^2u$.

Решение

Обозначим $r=\sqrt{x^2+y^2}$. $r_x=\dfrac xr$, $r_y=\dfrac yr$, $dr=\dfrac{x\,dx+y\,dy}{r}$.

$r_{xx}=\dfrac{r-x\cdot x/r}{r^2}=\dfrac{y^2}{r^3}$, $r_{xy}=-\dfrac{xy}{r^3}$, $r_{yy}=\dfrac{x^2}{r^3}$, поэтому $d^2r=\dfrac{y^2dx^2-2xy\,dx\,dy+x^2dy^2}{r^3}=\dfrac{(y\,dx-x\,dy)^2}{r^3}$.

$$du=f'(r)\frac{x\,dx+y\,dy}{r},\qquad d^2u=f''(r)\frac{(x\,dx+y\,dy)^2}{r^2}+f'(r)\frac{(y\,dx-x\,dy)^2}{r^3}.$$
Пример 12. $u=f(xyz)$

с практики 28.09 Найти $du$, $d^2u$.

Решение

$t=xyz$: $dt=yz\,dx+xz\,dy+xy\,dz$. Вторые производные $t$: $t_{xx}=t_{yy}=t_{zz}=0$, $t_{xy}=z$, $t_{xz}=y$, $t_{yz}=x$, значит $d^2t=2(z\,dx\,dy+y\,dx\,dz+x\,dy\,dz)$.

$$du=f'(xyz)(yz\,dx+xz\,dy+xy\,dz),$$ $$d^2u=f''(xyz)(yz\,dx+xz\,dy+xy\,dz)^2+2f'(xyz)(z\,dx\,dy+y\,dx\,dz+x\,dy\,dz).$$
Пример 13. $u=f(x,y,z)$, $x=t$, $y=t^2$, $z=t^3$

с практики 28.09 Найти $du$, $d^2u$. Здесь независимая переменная одна — $t$.

Решение

$dx=dt$, $dy=2t\,dt$, $dz=3t^2dt$; $d^2x=0$, $d^2y=2\,dt^2$, $d^2z=6t\,dt^2$.

Первый дифференциал по инвариантной формуле:

$$du=f'_1dx+f'_2dy+f'_3dz=(f'_1+2tf'_2+3t^2f'_3)\,dt.$$

Второй. Здесь $x,y,z$ — не независимые, поэтому к «обычной» квадратичной форме добавляются слагаемые с $d^2x,d^2y,d^2z$:

$$d^2u=\sum_{i,j}f''_{ij}\,dx_i\,dx_j+f'_1d^2x+f'_2d^2y+f'_3d^2z.$$

Подставляем:

$$d^2u=\bigl(f''_{11}+4t^2f''_{22}+9t^4f''_{33}+4tf''_{12}+6t^2f''_{13}+12t^3f''_{23}+2f'_2+6tf'_3\bigr)dt^2.$$

Проверка: это то же, что $\dfrac{d^2}{dt^2}f(t,t^2,t^3)\,dt^2$ (сверено в sympy).

3. Инвариантность первого дифференциала и неинвариантность второго

Инвариантность формы первого дифференциала

Для $u=f(x,y)$ формула

$$du=\pd{f}{x}dx+\pd{f}{y}dy$$

верна и когда $x,y$ — независимые переменные, и когда $x=x(s,t)$, $y=y(s,t)$ (тогда $dx$, $dy$ — дифференциалы этих функций).

Почему: при зависимых $x,y$ по определению $du=u_sds+u_tdt$. Подставим цепное правило $u_s=f_xx_s+f_yy_s$, $u_t=f_xx_t+f_yy_t$ и перегруппируем: $du=f_x(x_sds+x_tdt)+f_y(y_sds+y_tdt)=f_xdx+f_ydy$.

Со вторым дифференциалом так не выйдет. Для $u=f(t)$ считаем $d^2u=d(du)=d(f'(t)\,dt)$. Если $t$ — независимая переменная, $dt$ — константа, и $d^2u=f''dt^2$. Если $t=t(x,y)$, то $dt=t_xdx+t_ydy$ уже зависит от точки, и по правилу произведения

Второй дифференциал сложной функции $$d^2u=f''(t)\,dt^2+f'(t)\,d^2t.$$

Второе слагаемое исчезает, только если $d^2t=0$, например когда $t$ — независимая переменная или линейная функция независимых.

Пример 14

$u=t^2$, $t=xy$. Сравнить честный $d^2u$ с «наивным» $f''(t)dt^2$.

Решение

Честно: $u=x^2y^2$, $u_{xx}=2y^2$, $u_{xy}=4xy$, $u_{yy}=2x^2$:

$$d^2u=2y^2dx^2+8xy\,dx\,dy+2x^2dy^2.$$

Наивно: $f''=2$, $dt=y\,dx+x\,dy$, $f''dt^2=2y^2dx^2+4xy\,dx\,dy+2x^2dy^2$ — не хватает $4xy\,dx\,dy$.

Недостающее — ровно $f'(t)d^2t$: $f'=2t=2xy$, $d^2t=2\,dx\,dy$, произведение $4xy\,dx\,dy$. С ним всё сходится.

$u=f(t)$, $t=x^2+y$ ($x,y$ независимые). Чему равен $d^2u$?
$d^2t=t_{xx}dx^2+2t_{xy}dx\,dy+t_{yy}dy^2=2\,dx^2$, значит $d^2u=f''dt^2+f'd^2t=f''(2x\,dx+dy)^2+2f'dx^2$.

4. Замена переменных в выражениях

Цепное правило позволяет проверять, что функция определённого вида удовлетворяет уравнению, и переписывать выражения в других координатах. Кратко, на трёх примерах.

Пример 15. Бегущая волна

Проверить, что $u=f(x-at)$ при любой дважды дифференцируемой $f$ удовлетворяет уравнению $u_{tt}=a^2u_{xx}$.

Решение

Аргумент $\xi=x-at$, $\xi_x=1$, $\xi_t=-a$.

$u_x=f'(\xi)$, $u_{xx}=f''(\xi)$; $u_t=-af'(\xi)$, $u_{tt}=a^2f''(\xi)$. Значит $u_{tt}=a^2u_{xx}$.

Смысл: график $f$ сдвигается вправо со скоростью $a$, не меняя формы (подвигай ползунок ниже, там $a=1$).

Пример 16

Показать, что $z=f(x^2-y^2)$ удовлетворяет уравнению $y\,z_x+x\,z_y=0$.

Решение

$z_x=2xf'$, $z_y=-2yf'$. Тогда $y\cdot 2xf'+x\cdot(-2yf')=0$.

Пример 17. Полярные координаты

Пусть $x=r\cos\varphi$, $y=r\sin\varphi$, $u(x,y)=v(r,\varphi)$. Выразить $u_x^2+u_y^2$ через $v_r$, $v_\varphi$.

Решение

Удобнее идти «обратным» путём: $v$ — сложная функция от $r,\varphi$ через $x,y$:

$$v_r=u_x\cos\varphi+u_y\sin\varphi,\qquad v_\varphi=-u_x\,r\sin\varphi+u_y\,r\cos\varphi.$$

Тогда

$$v_r^2+\frac{v_\varphi^2}{r^2}=u_x^2(\cos^2\varphi+\sin^2\varphi)+u_y^2(\sin^2\varphi+\cos^2\varphi)+2u_xu_y(\cos\varphi\sin\varphi-\sin\varphi\cos\varphi)=u_x^2+u_y^2.$$

Ответ: $u_x^2+u_y^2=v_r^2+\dfrac{1}{r^2}v_\varphi^2$.

Проверка на $u=x^2+y^2=r^2$: слева $4x^2+4y^2=4r^2$, справа $(2r)^2+0=4r^2$.

5. Тренировка: сложные функции

Пример 18

$z=\ln(x^2+y^2)$, $x=uv$, $y=\dfrac uv$. Найти $z_u$, $z_v$.

Решение

$z_x=\dfrac{2x}{x^2+y^2}$, $z_y=\dfrac{2y}{x^2+y^2}$; $x_u=v$, $x_v=u$, $y_u=\dfrac1v$, $y_v=-\dfrac{u}{v^2}$. Учтём $x^2+y^2=u^2\left(v^2+v^{-2}\right)$.

$$z_u=\frac{2xv+2y/v}{x^2+y^2}=\frac{2u(v^2+v^{-2})}{u^2(v^2+v^{-2})}=\frac2u,$$ $$z_v=\frac{2xu-2yu/v^2}{x^2+y^2}=\frac{2u^2(v-v^{-3})}{u^2(v^2+v^{-2})}=\frac{2(v^4-1)}{v(v^4+1)}.$$
Пример 19

$z=\arctg\dfrac xy$, $x=u+v$, $y=u-v$. Найти $z_u$, $z_v$.

Решение

$z_x=\dfrac{1}{1+x^2/y^2}\cdot\dfrac1y=\dfrac{y}{x^2+y^2}$, $z_y=-\dfrac{x}{x^2+y^2}$. $x^2+y^2=2(u^2+v^2)$.

$$z_u=z_x+z_y=\frac{y-x}{x^2+y^2}=\frac{-2v}{2(u^2+v^2)}=-\frac{v}{u^2+v^2},\qquad z_v=z_x-z_y=\frac{x+y}{x^2+y^2}=\frac{u}{u^2+v^2}.$$
Пример 20

$u=f\left(\dfrac xy,\dfrac yz\right)$. Найти первые производные и показать, что $xu_x+yu_y+zu_z=0$.

Решение
$$u_x=\frac1yf'_1,\qquad u_y=-\frac{x}{y^2}f'_1+\frac1zf'_2,\qquad u_z=-\frac{y}{z^2}f'_2.$$ $$xu_x+yu_y+zu_z=\frac xyf'_1-\frac xyf'_1+\frac yzf'_2-\frac yzf'_2=0.$$

Это частный случай теоремы Эйлера: $u$ однородна степени 0 (не меняется при $x,y,z\to\lambda x,\lambda y,\lambda z$).

Пример 21

$z=f\left(x^2-y^2,\;e^{xy}\right)$. Найти $z_x$, $z_y$, $z_{xy}$.

Решение

$\xi=x^2-y^2$, $\eta=e^{xy}$: $\xi_x=2x$, $\xi_y=-2y$, $\eta_x=ye^{xy}$, $\eta_y=xe^{xy}$.

$$z_x=2xf'_1+ye^{xy}f'_2,\qquad z_y=-2yf'_1+xe^{xy}f'_2.$$

$z_{xy}$: $(f'_k)'_y=-2yf''_{k1}+xe^{xy}f''_{k2}$, и $(ye^{xy})'_y=e^{xy}+xye^{xy}$:

$$z_{xy}=2x\bigl(-2yf''_{11}+xe^{xy}f''_{12}\bigr)+(1+xy)e^{xy}f'_2+ye^{xy}\bigl(-2yf''_{21}+xe^{xy}f''_{22}\bigr)$$ $$=-4xyf''_{11}+2(x^2-y^2)e^{xy}f''_{12}+xye^{2xy}f''_{22}+(1+xy)e^{xy}f'_2.$$
Пример 22

$z=f\left(x,\dfrac yx\right)$. Найти $z_x$ и $z_{xx}$.

Решение

Здесь первый аргумент — сам $x$: $\xi=x$, $\eta=y/x$, $\eta_x=-y/x^2$.

$$z_x=f'_1-\frac{y}{x^2}f'_2.$$

$(f'_k)'_x=f''_{k1}-\dfrac{y}{x^2}f''_{k2}$, $\left(-\dfrac{y}{x^2}\right)'_x=\dfrac{2y}{x^3}$:

$$z_{xx}=f''_{11}-\frac{y}{x^2}f''_{12}+\frac{2y}{x^3}f'_2-\frac{y}{x^2}\left(f''_{21}-\frac{y}{x^2}f''_{22}\right)=f''_{11}-\frac{2y}{x^2}f''_{12}+\frac{y^2}{x^4}f''_{22}+\frac{2y}{x^3}f'_2.$$
Пример 23

Доказать, что $z=x^nf\left(\dfrac yx\right)$ удовлетворяет уравнению $xz_x+yz_y=nz$.

Решение

$z_x=nx^{n-1}f-x^n\cdot\dfrac{y}{x^2}f'=nx^{n-1}f-x^{n-2}yf'$, $z_y=x^n\cdot\dfrac1xf'=x^{n-1}f'$.

$xz_x+yz_y=nx^nf-x^{n-1}yf'+x^{n-1}yf'=nz$.

Пример 24

$u=e^{x-2y}$, $x=\sin t$, $y=t^3$. Найти $\dfrac{du}{dt}$.

Решение
$$\frac{du}{dt}=e^{x-2y}\cos t+e^{x-2y}\cdot(-2)\cdot3t^2=e^{\sin t-2t^3}\left(\cos t-6t^2\right).$$
Пример 25

$z=f(x,y)$, где $y=\varphi(x)$. Найти $\dfrac{dz}{dx}$ и $\dfrac{d^2z}{dx^2}$.

Решение

Всё зависит от одного $x$, так что производная полная:

$$\frac{dz}{dx}=f_x+f_y\varphi'.$$

Вторая: $f_x$ и $f_y$ сами сложные функции от $x$, $(f_x)'=f_{xx}+f_{xy}\varphi'$, $(f_y)'=f_{yx}+f_{yy}\varphi'$:

$$\frac{d^2z}{dx^2}=f_{xx}+2f_{xy}\varphi'+f_{yy}\varphi'^2+f_y\varphi''.$$

Если вдоль кривой $y=\varphi(x)$ функция $f\equiv0$, обе производные равны нулю — отсюда сразу получаются формулы для неявной функции из части Б.

Пример 26

из примера КР1 в РПД Дана функция $z=\ln(x+e^{-y})$. Показать, что $\pd{z}{x}\cdot\dfrac{\partial^2z}{\partial x\,\partial y}-\pd{z}{y}\cdot\dfrac{\partial^2z}{\partial x^2}=0$.

Решение

Это сложная функция $z=\ln D$, где $D=x+e^{-y}$, $D_x=1$, $D_y=-e^{-y}$.

$$z_x=\frac1D,\qquad z_y=\frac{-e^{-y}}{D},\qquad z_{xx}=-\frac{1}{D^2},\qquad z_{xy}=\left(D^{-1}\right)'_y=-\frac{D_y}{D^2}=\frac{e^{-y}}{D^2}.$$

Подставляем:

$$z_xz_{xy}-z_yz_{xx}=\frac1D\cdot\frac{e^{-y}}{D^2}-\frac{-e^{-y}}{D}\cdot\left(-\frac1{D^2}\right)=\frac{e^{-y}}{D^3}-\frac{e^{-y}}{D^3}=0.$$

Часть Б Неявные функции

6. Неявная функция: идея и картинка

Уравнение $F(x,y)=0$ задаёт на плоскости кривую — это линия уровня $0$ функции $F$. Например, $x^2+y^2-1=0$ — окружность. Вся окружность не является графиком $y(x)$: одному $x$ отвечают два $y$. Но если взять маленький кусочек кривой рядом с точкой, он почти всегда локально является графиком функции $y(x)$ (верхняя полуокружность $y=\sqrt{1-x^2}$ или нижняя). Такую функцию называют неявной: её задаёт уравнение, а формулы может и не быть.

Где это ломается? В точках, где касательная к кривой вертикальна. У окружности это $(\pm1,0)$: рядом с ними одному $x$ отвечают две точки кривой (или ни одной). Там $F_y=2y=0$. Это не совпадение, см. теорему ниже.

Картинка подсказывает формулу производной: градиент $(F_x,F_y)$ перпендикулярен кривой, значит направление касательной $(F_y,-F_x)$, и её угловой коэффициент $y'=-F_x/F_y$. При $F_y=0$ делить не на что — касательная вертикальна.

Пример поинтереснее — декартов лист $x^3+y^3-3axy=0$ (на картинке $a=1$). $F_y=3y^2-3ax$. Решая вместе с уравнением кривой, получаем: $y^2=ax$, подставляем $x=y^2/a$: $y^6/a^3+y^3=3y^3$, откуда $y=0$ или $y^3=2a^3$. Точка $\left(\sqrt[3]{4}\,a,\ \sqrt[3]{2}\,a\right)\approx(1{,}59;\,1{,}26)$ — вертикальная касательная. Точка $(0,0)$ особая: там $F_x=F_y=0$, кривая сама себя пересекает, и никакой единственной функции $y(x)$ рядом нет.

7. Теорема о неявной функции и производные y′, y″

Теорема о неявной функции

Пусть

  1. $F(x_0,y_0)=0$;
  2. $F$, $F_x$, $F_y$ непрерывны в окрестности точки $(x_0,y_0)$;
  3. $F_y(x_0,y_0)\ne0$.

Тогда найдётся прямоугольник $|x-x_0|\lt\delta$, $|y-y_0|\lt\varepsilon$, в котором уравнение $F(x,y)=0$ задаёт единственную функцию $y=y(x)$, причём $y(x_0)=y_0$, $y(x)$ непрерывно дифференцируема и

$$y'(x)=-\frac{F_x(x,y)}{F_y(x,y)}.$$

Откуда формула. Существование — содержание теоремы. А саму формулу легко вывести: раз $y(x)$ существует, $F(x,y(x))\equiv0$ — тождество по $x$. Дифференцируем его как сложную функцию (дерево: $F\to x$ напрямую и $F\to y\to x$):

$$F_x+F_y\cdot y'=0\quad\Longrightarrow\quad y'=-\frac{F_x}{F_y}.$$
Важно

В $F_x$ переменная $y$ заморожена — это обычная частная производная $F(x,y)$. Зависимость $y$ от $x$ учитывается отдельным слагаемым $F_yy'$. Не смешивай: либо пиши формулу $-F_x/F_y$ с частными производными, либо дифференцируй тождество, помня, что $y=y(x)$.

Два способа найти $y''$
  1. Дифференцировать тождество дважды. Продифференцировать $F(x,y(x))=0$ по $x$ (всюду, где стоит $y$, помнить про $y'$), получить равенство с $y'$; продифференцировать его ещё раз — появится $y''$; выразить $y''$ и подставить найденный $y'$.
  2. Дифференцировать формулу $y'$. Взять $y'=-F_x/F_y$ (выражение от $x$ и $y$) и продифференцировать по $x$ как сложную функцию, считая $y=y(x)$; вместо $y'$ подставить найденное выражение.

Результат можно (и часто нужно) упрощать с помощью исходного уравнения.

Общая формула (знать не обязательно) $$y''=-\frac{F_{xx}F_y^2-2F_{xy}F_xF_y+F_{yy}F_x^2}{F_y^3}.$$

Получается способом 2 в общем виде. На КР проще посчитать руками, чем вспоминать знаки.

Пример 27. Окружность

$x^2+y^2=1$. Найти $y'$ и $y''$ обоими способами.

Решение

Способ 1. Тождество $x^2+y^2(x)=1$. Дифференцируем: $2x+2yy'=0$, $y'=-\dfrac xy$. Ещё раз: $2+2y'^2+2yy''=0$, откуда

$$y''=-\frac{1+y'^2}{y}=-\frac{1+x^2/y^2}{y}=-\frac{x^2+y^2}{y^3}=-\frac1{y^3}.$$

Способ 2. $y'=-\dfrac xy$. Дифференцируем дробь, $y=y(x)$:

$$y''=-\frac{1\cdot y-x\cdot y'}{y^2}=-\frac{y+x^2/y}{y^2}=-\frac{x^2+y^2}{y^3}=-\frac1{y^3}.$$

В последнем шаге использовали уравнение кривой $x^2+y^2=1$. Проверка на верхней полуокружности: $y=\sqrt{1-x^2}$, $y''=-(1-x^2)^{-3/2}=-1/y^3$.

Пример 28. Декартов лист

$x^3+y^3-3axy=0$. Найти $y'$ и $y''$.

Решение

$F_x=3x^2-3ay$, $F_y=3y^2-3ax$:

$$y'=-\frac{x^2-ay}{y^2-ax}=\frac{ay-x^2}{y^2-ax}.$$

$y''$ способом 1. Тождество, делённое на 3, продифференцировано: $x^2+y^2y'-ay-axy'=0$. Ещё раз по $x$:

$$2x+2yy'^2+y^2y''-2ay'-axy''=0\ \Longrightarrow\ y''(y^2-ax)=2ay'-2x-2yy'^2.$$

Обозначим $P=ay-x^2$, $Q=y^2-ax$, $y'=P/Q$. Тогда $y''=\dfrac{2(aPQ-xQ^2-yP^2)}{Q^3}$. Раскрываем скобки:

$$aPQ-xQ^2-yP^2=-a^3xy+3ax^2y^2-xy^4-x^4y=-xy\left(x^3+y^3-3axy\right)-a^3xy=-a^3xy,$$

потому что точка лежит на кривой. Ответ: $y''=-\dfrac{2a^3xy}{(y^2-ax)^3}$.

Пример 29. Производные в точке

$e^y+xy=e$. Найти $y'(0)$ и $y''(0)$.

Решение

Сначала $y(0)$: $e^y=e$, $y=1$. Теорема применима: $F_y=e^y+x=e\ne0$.

Тождество: $e^yy'+y+xy'=0$. При $x=0$, $y=1$: $ey'+1=0$, $y'(0)=-\dfrac1e$.

Ещё раз: $e^yy'^2+e^yy''+y'+y'+xy''=0$. При $x=0$: $e\cdot\dfrac1{e^2}+ey''-\dfrac2e=0$, откуда $y''(0)=\dfrac{1}{e^2}$.

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

8. Функция двух переменных z(x, y), заданная неявно

Уравнение $F(x,y,z)=0$ задаёт поверхность. Если $F_z\ne0$ в точке, рядом с ней поверхность — график функции $z=z(x,y)$. Тождество $F(x,y,z(x,y))\equiv0$ дифференцируем по $x$ (при этом $y$ постоянна, а $z$ зависит от $x$):

$$F_x+F_z\,z_x=0,\qquad F_y+F_z\,z_y=0.$$
Формулы

Если $F(x_0,y_0,z_0)=0$, $F$ непрерывно дифференцируема и $F_z(x_0,y_0,z_0)\ne0$, то рядом с точкой уравнение задаёт единственную $z=z(x,y)$, и

$$z_x=-\frac{F_x}{F_z},\qquad z_y=-\frac{F_y}{F_z}.$$
Как найти все производные 1 и 2 порядка
  1. Найди $z_x$, $z_y$ по формулам (или дифференцируя тождество).
  2. $z_{xx}$: продифференцируй выражение для $z_x$ по $x$; где встретится $z_x$ — подставь найденное.
  3. $z_{xy}$: продифференцируй $z_x$ по $y$; где встретится $z_y$ — подставь.
  4. $z_{yy}$: аналогично или по симметрии, если уравнение симметрично по $x$ и $y$.
  5. Упрости, используя исходное уравнение.

Третий путь — через дифференциалы: $dF=F_xdx+F_ydy+F_zdz=0$, откуда сразу $dz=-\dfrac{F_xdx+F_ydy}{F_z}$, а коэффициенты при $dx$ и $dy$ — это $z_x$ и $z_y$.

Пример 30. $x^2+y^2+z^2=a^2$

ДЗ 05.10, будет на КР (по словам одногруппника) Найти частные производные 1 и 2 порядка, считая $z=z(x,y)$.

Решение

Первые. Дифференцируем тождество $x^2+y^2+z^2(x,y)=a^2$ по $x$: $2x+2zz_x=0$. По $y$: $2y+2zz_y=0$.

$$z_x=-\frac xz,\qquad z_y=-\frac yz\qquad(z\ne0).$$

(То же по формуле: $F_x=2x$, $F_y=2y$, $F_z=2z$.)

$z_{xx}$. Дифференцируем $z_x=-\dfrac xz$ по $x$, помня, что $z=z(x,y)$:

$$z_{xx}=-\frac{1\cdot z-x\,z_x}{z^2}=-\frac{z+x^2/z}{z^2}=-\frac{x^2+z^2}{z^3}.$$

$z_{xy}$. Дифференцируем $z_x=-\dfrac xz$ по $y$ ($x$ — константа):

$$z_{xy}=-x\cdot\left(-\frac{z_y}{z^2}\right)=\frac{x\,z_y}{z^2}=-\frac{xy}{z^3}.$$

$z_{yy}$ по симметрии: $z_{yy}=-\dfrac{y^2+z^2}{z^3}$.

Ответ:

$$z_x=-\frac xz,\quad z_y=-\frac yz,\quad z_{xx}=-\frac{x^2+z^2}{z^3}=-\frac{a^2-y^2}{z^3},\quad z_{xy}=-\frac{xy}{z^3},\quad z_{yy}=-\frac{y^2+z^2}{z^3}=-\frac{a^2-x^2}{z^3}.$$

Проверка явным выражением. На верхней полусфере $z=\sqrt{a^2-x^2-y^2}$, $z_x=\dfrac{-x}{\sqrt{a^2-x^2-y^2}}=-\dfrac xz$ — совпадает. Неявный способ хорош тем, что сразу работает и для нижней полусферы ($z\lt0$).

Через дифференциалы то же самое выглядит так: $x\,dx+y\,dy+z\,dz=0$, $dz=-\dfrac{x\,dx+y\,dy}{z}$, а

$$d^2z=-\frac{(x^2+z^2)dx^2+2xy\,dx\,dy+(y^2+z^2)dy^2}{z^3}.$$

Где не работает: при $z=0$ (экватор сферы) $F_z=0$, касательная плоскость вертикальна, $z(x,y)$ там не определяется однозначно.

Пример 31. $z^3-3xyz=a^3$

ДЗ 05.10, будет на КР (по словам одногруппника) Найти частные производные 1 и 2 порядка, считая $z=z(x,y)$.

Решение

Первые. $F=z^3-3xyz-a^3$: $F_x=-3yz$, $F_y=-3xz$, $F_z=3z^2-3xy$. Обозначим $Q=z^2-xy$ (нужно $Q\ne0$).

$$z_x=-\frac{-3yz}{3(z^2-xy)}=\frac{yz}{Q},\qquad z_y=\frac{xz}{Q}.$$

Проверка дифференцированием тождества по $x$: $3z^2z_x-3yz-3xyz_x=0$, $z_x(z^2-xy)=yz$ — то же.

$z_{xx}$. Дифференцируем $\dfrac{yz}{Q}$ по $x$; $(yz)'_x=yz_x$, $Q'_x=2zz_x-y$:

$$z_{xx}=\frac{yz_x\cdot Q-yz(2zz_x-y)}{Q^2}=\frac{yz_x(Q-2z^2)+y^2z}{Q^2}=\frac{-yz_x(z^2+xy)+y^2z}{Q^2}.$$

Подставляем $z_x=\dfrac{yz}{Q}$:

$$z_{xx}=\frac{y^2z\bigl(-(z^2+xy)+Q\bigr)}{Q^3}=\frac{y^2z(-z^2-xy+z^2-xy)}{Q^3}=-\frac{2xy^3z}{(z^2-xy)^3}.$$

$z_{xy}$. Дифференцируем $\dfrac{yz}{Q}$ по $y$; $(yz)'_y=z+yz_y$, $Q'_y=2zz_y-x$:

$$z_{xy}=\frac{(z+yz_y)Q-yz(2zz_y-x)}{Q^2}=\frac{zQ+xyz+yz_y(Q-2z^2)}{Q^2}=\frac{z^3-yz_y(z^2+xy)}{Q^2}.$$

Подставляем $z_y=\dfrac{xz}{Q}$:

$$z_{xy}=\frac{z^3Q-xyz(z^2+xy)}{Q^3}=\frac{z\bigl(z^4-xyz^2-xyz^2-x^2y^2\bigr)}{Q^3}=\frac{z\left(z^4-2xyz^2-x^2y^2\right)}{(z^2-xy)^3}.$$

$z_{yy}$: уравнение симметрично по $x\leftrightarrow y$, поэтому меняем их местами в $z_{xx}$.

Ответ:

$$z_x=\frac{yz}{z^2-xy},\quad z_y=\frac{xz}{z^2-xy},\quad z_{xx}=-\frac{2xy^3z}{(z^2-xy)^3},\quad z_{yy}=-\frac{2x^3yz}{(z^2-xy)^3},\quad z_{xy}=\frac{z(z^4-2xyz^2-x^2y^2)}{(z^2-xy)^3}.$$

Проверка в sympy (так проверялись все ответы этой главы):

from sympy import symbols, Function, diff, solve, factor
x, y, a = symbols('x y a')
z = Function('z')(x, y)
F = z**3 - 3*x*y*z - a**3        # F(x, y, z(x, y)) = 0 - тождество по x, y
zx = solve(diff(F, x), diff(z, x))[0]
zy = solve(diff(F, y), diff(z, y))[0]
zxx = solve(diff(F, x, 2), diff(z, x, 2))[0].subs(diff(z, x), zx)
zxy = solve(diff(F, x, y), diff(z, x, y))[0].subs({diff(z, x): zx, diff(z, y): zy})
print(factor(zx), factor(zxx), factor(zxy))
# -y*z/(x*y - z**2)   2*x*y**3*z/(x*y - z**2)**3   (x**2*y**2 + 2*x*y*z**2 - z**4)*z/(x*y - z**2)**3

У sympy знаменатель $xy-z^2=-(z^2-xy)$, знаки совпадают с нашими.

Пример 32. Дифференциал в точке

из примера КР1 в РПД Найти дифференциал функции $z(x,y)$, заданной неявно уравнением $e^{z-1}=\cos x\cos y+1$, в точке $M_0\left(0;\frac{\pi}{2};1\right)$.

Решение

Проверим, что точка на поверхности: $e^0=1$, $\cos0\cdot\cos\frac\pi2+1=1$. Да.

$F=e^{z-1}-\cos x\cos y-1$: $F_x=\sin x\cos y$, $F_y=\cos x\sin y$, $F_z=e^{z-1}$. В $M_0$: $F_x=0$, $F_y=1$, $F_z=1\ne0$ — теорема применима.

$$z_x(M_0)=-\frac{0}{1}=0,\qquad z_y(M_0)=-\frac11=-1.$$

Ответ: $dz(M_0)=0\cdot dx-1\cdot dy=-dy$.

$F(x,y,z)=x^2z+z^3-y=0$. Чему равно $z_x$?
$F_x=2xz$ (z заморожена), $F_z=x^2+3z^2$, $z_x=-F_x/F_z$. Вариант 1 — частая ошибка: в знаменателе взяли производную только одного слагаемого.

9. ДЗ 05.10: неявные функции одной переменной

На доске 05.10 (фото msg55, msg56) слева: «Найти $y'_x$ и $y''_{xx}$» для пунктов 2), 3), 4); пункта 1) на фото нет.

Пример 33. $\ln\sqrt{x^2+y^2}=\arctg\dfrac yx$

ДЗ 05.10 Найти $y'$ и $y''$.

Решение

Дифференцируем тождество. Слева $\ln\sqrt{x^2+y^2}=\frac12\ln(x^2+y^2)$, справа производная арктангенса:

$$\frac{x+yy'}{x^2+y^2}=\frac{1}{1+y^2/x^2}\cdot\frac{xy'-y}{x^2}=\frac{xy'-y}{x^2+y^2}.$$

Знаменатели одинаковые: $x+yy'=xy'-y$, откуда $y'(x-y)=x+y$:

$$y'=\frac{x+y}{x-y}.$$

$y''$ (способ 2, дифференцируем дробь):

$$y''=\frac{(1+y')(x-y)-(x+y)(1-y')}{(x-y)^2}=\frac{2xy'-2y}{(x-y)^2}.$$

$xy'-y=\dfrac{x(x+y)-y(x-y)}{x-y}=\dfrac{x^2+y^2}{x-y}$, поэтому

$$y''=\frac{2(x^2+y^2)}{(x-y)^3}.$$

Через формулу: $F=\frac12\ln(x^2+y^2)-\arctg\frac yx$, $F_x=\dfrac{x+y}{x^2+y^2}$, $F_y=\dfrac{y-x}{x^2+y^2}$, $y'=-F_x/F_y$ — то же.

Что это за кривая. В полярных координатах $\ln r=\arctg\frac yx=\varphi$ (при $x\gt0$), то есть $r=e^\varphi$ — логарифмическая спираль.

Пример 34. $y-\varepsilon\sin y=x$, $0\lt\varepsilon\lt1$

ДЗ 05.10 Найти $y'$ и $y''$.

Решение

Дифференцируем: $y'-\varepsilon\cos y\cdot y'=1$, то есть

$$y'=\frac{1}{1-\varepsilon\cos y}.$$

Зачем условие $0\lt\varepsilon\lt1$: $F_y=1-\varepsilon\cos y\ge1-\varepsilon\gt0$ при всех $y$. Значит, теорема о неявной функции работает в каждой точке, и знаменатель никогда не обнуляется. Более того, $x=y-\varepsilon\sin y$ строго возрастает по $y$, поэтому $y(x)$ определена на всей прямой.

$y''$. $y'=(1-\varepsilon\cos y)^{-1}$, дифференцируем как сложную функцию от $y(x)$:

$$y''=-(1-\varepsilon\cos y)^{-2}\cdot\varepsilon\sin y\cdot y'=-\frac{\varepsilon\sin y}{(1-\varepsilon\cos y)^3}.$$
Откуда это уравнение

$E-e\sin E=M$ — уравнение Кеплера из небесной механики: по «средней аномалии» $M$ (линейна по времени) находят положение планеты на эллиптической орбите с эксцентриситетом $0\lt e\lt1$. Явно $E(M)$ не выражается, а производные — пожалуйста.

Пример 35. $y=\arctg\dfrac yx$

ДЗ 05.10 Найти $y'$ и $y''$.

Решение

На обоих фото (msg55 и msg56) условие записано именно так: $y=\arctg\frac yx$. С ним есть странности, разберём честно.

1. Тривиальное решение. $y\equiv0$ удовлетворяет уравнению при любом $x\ne0$: $0=\arctg0$. Для него $y'=y''=0$, и считать нечего.

2. Нетривиальные решения есть только при $0\lt x\lt1$. Так как $|\arctg|\lt\frac\pi2$, то $|y|\lt\frac\pi2$. При $y\ne0$ уравнение равносильно $\tg y=\frac yx$, то есть

$$x=y\ctg y,\qquad 0\lt|y|\lt\frac\pi2.$$

Функция $y\ctg y$ на $\left(0,\frac\pi2\right)$ убывает от $1$ до $0$, поэтому нетривиальные решения есть только при $0\lt x\lt1$: две симметричные ветки $\pm y(x)$. При $x\to1$ они сливаются с $y=0$.

3. Формальное решение. $F=y-\arctg\frac yx$:

$$F_x=-\frac{1}{1+y^2/x^2}\cdot\left(-\frac{y}{x^2}\right)=\frac{y}{x^2+y^2},\qquad F_y=1-\frac{1}{1+y^2/x^2}\cdot\frac1x=\frac{x^2+y^2-x}{x^2+y^2}.$$ $$y'=-\frac{F_x}{F_y}=-\frac{y}{x^2+y^2-x}.$$

$y''$: обозначим $Q=x^2+y^2-x$, $Q'=2x-1+2yy'$, и заметим $y'Q=-y$:

$$y''=-\frac{y'Q-yQ'}{Q^2}=\frac{y+y(2x-1+2yy')}{Q^2}=\frac{2y(x+yy')}{Q^2}.$$

$x+yy'=x-\dfrac{y^2}{Q}=\dfrac{x^3+xy^2-x^2-y^2}{Q}=\dfrac{(x-1)(x^2+y^2)}{Q}$, итого

$$y''=\frac{2y(x-1)(x^2+y^2)}{(x^2+y^2-x)^3}.$$

На тривиальной ветке $y=0$ формулы дают $y'=y''=0$ (при $x\ne1$), всё согласовано.

4. Особая точка $(1,0)$. Там $F_y=\frac{1+0-1}{1}=0$: теорема о неявной функции не работает, и действительно, через эту точку проходят и $y=0$, и обе нетривиальные ветки.

5. Проверка через обратную функцию. На нетривиальной ветке $x=y\ctg y$ — явная функция от $y$, поэтому $y'=\dfrac{1}{x'(y)}=\dfrac{\sin^2y}{\sin y\cos y-y}$. Если подставить $x=y\ctg y$ в $-\dfrac{y}{x^2+y^2-x}$, получится то же (сверено в sympy численно и символьно).

Вывод. Ответ: $y'=-\dfrac{y}{x^2+y^2-x}$, $y''=\dfrac{2y(x-1)(x^2+y^2)}{(x^2+y^2-x)^3}$, с оговорками выше. Возможно, на доске опечатка, и имелось в виду похожее и более типичное уравнение $y=x+\arctg y$. Сверь с одногруппниками; на всякий случай решение и для него:

$y'=1+\dfrac{y'}{1+y^2}$ $\Rightarrow$ $y'\cdot\dfrac{y^2}{1+y^2}=1$ $\Rightarrow$ $y'=\dfrac{1+y^2}{y^2}=1+y^{-2}$, $\quad y''=-2y^{-3}y'=-\dfrac{2(1+y^2)}{y^5}$.

10. Системы неявных функций и якобиан

Иногда уравнений несколько, и они задают сразу несколько функций. Например, система

$$\begin{cases}F(x,y,u,v)=0,\\ G(x,y,u,v)=0\end{cases}$$

может задавать $u=u(x,y)$, $v=v(x,y)$. Условие, аналогичное $F_y\ne0$, — ненулевой якобиан по тем переменным, которые выражаем:

$$J=\frac{\partial(F,G)}{\partial(u,v)}=\begin{vmatrix}F_u&F_v\\ G_u&G_v\end{vmatrix}\ne0.$$
Как считать производные
  1. Продифференцировать оба уравнения по $x$, считая $u,v$ функциями: $F_x+F_uu_x+F_vv_x=0$, $G_x+G_uu_x+G_vv_x=0$.
  2. Это линейная система относительно $u_x$, $v_x$ с определителем $J$. Решить (по Крамеру).
  3. То же по $y$.
Пример 36

$x=u\cos v$, $y=u\sin v$. Найти $u_x$, $u_y$, $v_x$, $v_y$.

Решение

Дифференцируем по $x$:

$$\begin{cases}1=u_x\cos v-u\sin v\cdot v_x,\\ 0=u_x\sin v+u\cos v\cdot v_x.\end{cases}$$

Определитель $\begin{vmatrix}\cos v&-u\sin v\\ \sin v&u\cos v\end{vmatrix}=u$, нужно $u\ne0$. По Крамеру:

$$u_x=\frac{\begin{vmatrix}1&-u\sin v\\0&u\cos v\end{vmatrix}}{u}=\cos v,\qquad v_x=\frac{\begin{vmatrix}\cos v&1\\ \sin v&0\end{vmatrix}}{u}=-\frac{\sin v}{u}.$$

Аналогично по $y$ (правая часть $(0,1)$): $u_y=\sin v$, $v_y=\dfrac{\cos v}{u}$.

Смысл: это полярные координаты, $u=r$, $v=\varphi$. Проверка: $r=\sqrt{x^2+y^2}$, $r_x=\frac xr=\cos\varphi$ — совпадает.

11. Тренировка: неявные функции

Пример 37

$x^2-xy+2y^2+x-y-1=0$. Найти $y'$ и $y''$ в точке $(0;1)$.

Решение

Точка на кривой: $0-0+2+0-1-1=0$. Дифференцируем:

$$2x-y-xy'+4yy'+1-y'=0.$$

В $(0,1)$: $-1+4y'+1-y'=0$, $y'(0)=0$.

Ещё раз: $2-y'-y'-xy''+4y'^2+4yy''-y''=0$. В точке ($y'=0$): $2+4y''-y''=0$, $y''(0)=-\dfrac23$.

Пример 38

$y=1+xe^y$. Найти $y'$ и $y''$.

Решение

$y'=e^y+xe^yy'$, $y'=\dfrac{e^y}{1-xe^y}$. Из уравнения $xe^y=y-1$, поэтому $y'=\dfrac{e^y}{2-y}$.

$y''=\dfrac{e^yy'(2-y)+e^yy'}{(2-y)^2}=\dfrac{e^yy'(3-y)}{(2-y)^2}=\dfrac{e^{2y}(3-y)}{(2-y)^3}$.

Пример 39

$x^y=y^x$ ($x,y\gt0$). Найти $y'$.

Решение

Логарифмируем: $y\ln x=x\ln y$. Дифференцируем: $y'\ln x+\dfrac yx=\ln y+\dfrac{x}{y}y'$.

$$y'\left(\ln x-\frac xy\right)=\ln y-\frac yx\ \Longrightarrow\ y'=\frac{y(x\ln y-y)}{x(y\ln x-x)}.$$
Пример 40

$e^z=xyz$. Найти $z_x$, $z_y$, $z_{xx}$, $z_{xy}$.

Решение

$F=e^z-xyz$: $F_x=-yz$, $F_y=-xz$, $F_z=e^z-xy$. Используя $e^z=xyz$:

$$z_x=\frac{yz}{e^z-xy}=\frac{yz}{xyz-xy}=\frac{z}{x(z-1)},\qquad z_y=\frac{z}{y(z-1)}.$$

$z_{xy}$: $x$ — константа, $\left(\dfrac{z}{z-1}\right)'_y=\dfrac{z_y(z-1)-zz_y}{(z-1)^2}=-\dfrac{z_y}{(z-1)^2}$:

$$z_{xy}=-\frac{1}{x}\cdot\frac{z_y}{(z-1)^2}=-\frac{z}{xy(z-1)^3}.$$

$z_{xx}$: дифференцируем $\dfrac{z}{x(z-1)}$ по $x$:

$$z_{xx}=\frac{z_x\cdot x(z-1)-z\bigl((z-1)+xz_x\bigr)}{x^2(z-1)^2}=\frac{-xz_x-z(z-1)}{x^2(z-1)^2}=-\frac{z(z^2-2z+2)}{x^2(z-1)^3}.$$

(В предпоследнем шаге $xz_x(z-1)-xzz_x=-xz_x$, затем подставили $xz_x=\frac{z}{z-1}$.)

Пример 41

$\dfrac xz=\ln\dfrac zy$. Найти $z_x$, $z_y$, $z_{xx}$.

Решение

$F=\dfrac xz-\ln z+\ln y$: $F_x=\dfrac1z$, $F_y=\dfrac1y$, $F_z=-\dfrac{x}{z^2}-\dfrac1z=-\dfrac{x+z}{z^2}$.

$$z_x=\frac{1/z}{(x+z)/z^2}=\frac{z}{x+z},\qquad z_y=\frac{z^2}{y(x+z)}.$$ $$z_{xx}=\frac{z_x(x+z)-z(1+z_x)}{(x+z)^2}=\frac{xz_x-z}{(x+z)^2}=\frac{\frac{xz}{x+z}-z}{(x+z)^2}=-\frac{z^2}{(x+z)^3}.$$
Пример 42. $dz$ и $d^2z$ в точке

$x^2+2y^2+3z^2+xy-z-9=0$, точка $(1;-2;1)$. Найти $dz$ и $d^2z$.

Решение

Точка на поверхности: $1+8+3-2-1-9=0$.

По $x$: $2x+6zz_x+y-z_x=0$. В точке: $2+6z_x-2-z_x=0$, $z_x=0$.

По $y$: $4y+6zz_y+x-z_y=0$. В точке: $-8+5z_y+1=0$, $z_y=\frac75$.

Вторые: дифференцируем первое тождество по $x$: $2+6z_x^2+6zz_{xx}-z_{xx}=0$ $\Rightarrow$ $2+5z_{xx}=0$, $z_{xx}=-\frac25$.

Первое по $y$: $6z_yz_x+6zz_{xy}+1-z_{xy}=0$ $\Rightarrow$ $5z_{xy}+1=0$, $z_{xy}=-\frac15$.

Второе по $y$: $4+6z_y^2+6zz_{yy}-z_{yy}=0$ $\Rightarrow$ $4+\frac{294}{25}+5z_{yy}=0$, $z_{yy}=-\frac{394}{125}$.

$$dz=\frac75\,dy,\qquad d^2z=-\frac25dx^2-\frac25dx\,dy-\frac{394}{125}dy^2.$$
Пример 43

$x\cos y+y\cos z+z\cos x=1$. Найти $z_x$, $z_y$.

Решение

$F_x=\cos y-z\sin x$, $F_y=-x\sin y+\cos z$, $F_z=-y\sin z+\cos x$.

$$z_x=-\frac{\cos y-z\sin x}{\cos x-y\sin z},\qquad z_y=-\frac{\cos z-x\sin y}{\cos x-y\sin z}.$$

Типичная ошибка здесь — в $F_x$ забыть слагаемое от $z\cos x$: $z$ заморожена, но $\cos x$ дифференцируется.

Дифференцируем тождество $x^2+y^2+z^2=a^2$ по $x$ во второй раз. Что получится из слагаемого $z\cdot z_x$?
$z\cdot z_x$ — произведение двух функций от $x$, обе меняются: $(z\,z_x)'_x=z_x\cdot z_x+z\cdot z_{xx}$.

12. Типичные ошибки

Ошибки, которые чаще всего стоят баллов
  1. Забыть, что $z$ зависит от $x$. При дифференцировании тождества $F(x,y,z)=0$ каждое $z$ даёт множитель $z_x$: $(z^3)'_x=3z^2z_x$, а не $3z^2$ и не $0$.
  2. Обратная ошибка в формуле $-F_x/F_z$. В $F_x$ переменная $z$ заморожена, никаких $z_x$ там нет. Либо формула, либо дифференцирование тождества — не смешивать.
  3. Путать $f'_1$ и $\partial f/\partial x$. $f'_1$ — производная по первому аргументу $f$, а не по $x$. Например, для $u=f(x+y,xy)$: $u_x=f'_1+yf'_2$, а «$\partial f/\partial x$» без уточнения вообще неоднозначно.
  4. Считать $f'_1$ константой при второй производной: $(f'_1)'_x=f''_{11}\xi_x+f''_{12}\eta_x$, оба слагаемых.
  5. Потерять слагаемое от произведения: $(y\,f'_2)'_y=f'_2+y(f'_2)'_y$. Первое слагаемое пропускают постоянно.
  6. $d^2u=f''dt^2$ для зависимого $t$. Нужно $f''dt^2+f'd^2t$.
  7. Знак в $y'=-F_x/F_y$. Минус обязателен. Проверка на окружности: $y'=-x/y$, в первой четверти кривая убывает.
  8. Не подставить $y'$ в выражение для $y''$. В ответе $y''$ должны остаться только $x,y$ (и параметры), без $y'$.
  9. Не проверить $F_y\ne0$ (или $F_z\ne0$) в точке, где просят производную. Если ноль — теорема не применима.
  10. Неявная функция в точке: сначала найти недостающую координату из уравнения (например, $y(0)$), потом подставлять.

Мини-шпаргалка главы

Цепное правило

$$\frac{du}{dt}=f_x\,x'+f_y\,y',\qquad \pd{u}{s}=f_x\pd{x}{s}+f_y\pd{y}{s}.$$

Слагаемых столько, сколько путей в дереве от $u$ до переменной.

Абстрактная $f(\xi,\eta)$

$$u_x=f'_1\xi_x+f'_2\eta_x,\qquad (f'_k)'_x=f''_{k1}\xi_x+f''_{k2}\eta_x.$$

$f''_{12}=f''_{21}$. Произведения дифференцировать по правилу Лейбница.

Дифференциалы $u=f(t)$

$$du=f'dt\ \text{(инвариантно)},\qquad d^2u=f''dt^2+f'd^2t.$$

$\varphi^\psi$

$$\left(\varphi^\psi\right)'=\varphi^\psi\left(\psi'\ln\varphi+\psi\frac{\varphi'}{\varphi}\right).$$

Неявная $y(x)$: $F(x,y)=0$

$$y'=-\frac{F_x}{F_y},\quad F_y\ne0.$$

$y''$: дифференцировать тождество дважды или формулу $y'$ (помня $y=y(x)$), подставить $y'$, упростить уравнением.

Неявная $z(x,y)$: $F(x,y,z)=0$

$$z_x=-\frac{F_x}{F_z},\quad z_y=-\frac{F_y}{F_z},\quad F_z\ne0.$$

$dz=-\dfrac{F_xdx+F_ydy}{F_z}$.

Сфера $x^2+y^2+z^2=a^2$

$$z_x=-\frac xz,\ z_{xx}=-\frac{x^2+z^2}{z^3},\ z_{xy}=-\frac{xy}{z^3}.$$

$z^3-3xyz=a^3$

$$z_x=\frac{yz}{z^2-xy},\ z_{xx}=-\frac{2xy^3z}{(z^2-xy)^3},$$ $$z_{xy}=\frac{z(z^4-2xyz^2-x^2y^2)}{(z^2-xy)^3}.$$

Система: продифференцировать все уравнения, решить линейную систему; нужно $\dfrac{\partial(F,G)}{\partial(u,v)}\ne0$.