Дифференциальное исчисление функций многих переменных Частные производные >diff(f,x1$n1,x2$n2,…, xm$nm), где x1,…, xm – переменные, по которым производится. - презентация

1 Дифференциальное исчисление функций многих переменных Частные производные >diff(f,x1$n1,x2$n2,…, xm$nm), где x1,…, xm – переменные, по которым производится дифференцирование, а после знака $ указаны соответствующие порядки дифференцирования записывается в виде: diff(f,x,y).

2 Примеры Найти и функции > f:=arctan(x/y): >Diff(f,x)=simplify(diff(f,x)); > Diff(f,y)=simplify(diff(f,y));

3 >restart; f:=(x-y)/(x+y): > Diff(f,x$2)=simplify(diff(f,x$2)); > Diff(f,x,y)=diff(f,x,y);

4 Локальные и условные экстремумы функций многих переменных. Для исследования функции на локальный и условный экстремум используется команда из стандартной библиотеки >extrema(f,{cond},{x,y,…},'s'), где cond – ограничения для поиска условного экстремума, которые записываются в виде равенств. После ограничений в фигурных скобках указываются все переменные, от которых зависит функция f, а затем в кавычках записывается s – имя переменной, которой будут присвоены координаты точек экстремума. Если ограничений не указывать, то будет производиться поиск локального экстремума.

5 Наибольшее и наименьшее значения >maximize(f,{x1,…,xn},range), >minimize(f,{x1,…,xn}, range), range-интервалы для каждой переменной, указывающие область поиска наибольшего и наименьшего значений.

6 Примеры Найти экстремумы функции >restart: readlib(extrema): > f:=2*x^4+y^4-x^2-2*y^2: > extrema(f,{},{x,y},'s');s; >subs([x=1/2,y=1],f);

7 Примеры > restart: readlib(maximize):readlib(minimize): > f:=x^2+2*x*y-4*x+8*y: > maximize(f,{x,y},{x=0..1,y=0..2}); 17 > minimize(f,{x,y},{x=0..1,y=0..2}); -4

8 Примеры Найти условные экстремумы функции f(х,у)=xy+yz при условиях x2+y2=2, y+z=2, x>0, y>0, z>0 >restart: readlib(extrema): f:=x*y+y*z: >assume(x>0);assume(y>0);assume(z> 0); >simplify(extrema(f,{x^2+y^2=2,y+z= 2},{x,y,z},'s'));

9 {min(3/2RootOf(_Z2+4_Z+1)+1/2, 0), max(3/2RootOf(_Z2+4_Z+1)+1/2, 2)} >convert(%,radical); convert(s,radical); {min max }

10 Примеры на условный экстремум >restart: with(simplex): f:=-x+2*y+3*z: > cond:={x+2*y-3*z

11 Векторный анализ. Библиотека linalg. В Maple grad вычисляется одноименной командой >grad(f,[x,y,z],c), где f – функция, [x,y,z] – набор переменных, от которых она зависит. c- параметр позволяет вычислять данную дифференциальную операцию в различных криволинейных координатах (по умолчанию используется прямоугольная декартова система координат). Для вычисления дифференциальной операции в цилиндрических координатах следует записать coords=cylindrical,в сферических координатах – coords=spherical.

12 Пример Дана функция Найти Найти производную функции u(x,y) по направлению вектора q=[1,1]. >restart: with(linalg): > u:=arctan(y/x): g:=simplify(grad(u, [x, y]));

13 Пример >q:=vector([1,1]);e:=normalize(q); q:=[1, 1] е:= duq:=simplify(dotprod(g,e)); duq=

14 Ряды и произведения >sum(expr, n=a..b), где expr – выражение, зависящее от индекса суммирования, a..b – пределы индекса суммирования, указывающие, что суммировать следует от n=a до n=b (n=infinity). >product(P(n),n=a..b)

15 Пример Найти полную и N-частичную суммы ряда, общий член которого равен: an= >restart:a[n]:=1/((3*n-2)*(3*n+1)); >S[N]:=Sum(a[n], n=1..N)=sum(a[n], n=1..N); > S:=limit(rhs(S[N]), N=+infinity);

16 Пример Вычислить бесконечное произведение > Product((n^3-1)/(n^3+1),n=2..infinity) =simplify(product((n^3-1)/(n^3+1), n=2..infinity));

17 Разложение функции в степенной ряд и ряд Тейлора Разложение функции f(x) в степенной ряд в окрестности точки а >series(f(x), x=a, n), где а – точка, в окрестности которой производится разложение, n – число членов ряда. >taylor(f(x), x=a, n) раскладывает функции f(x) в окрестности точки x=a до порядка n-1 по формуле Тейлора. f(x1,…,xn) >readlib(mtaylor):>mtaylor(f(x), [x1,…,xn], n)

18 Пример Разложить в степенной ряд в окрестности х0=0, порядок O(x^5) f(x)=series(exp(-x)*sqrt(x+1), x=0, 5); >convert(%,polynom)

19 Пример Разложить в ряд Тейлора в окрестности точки (0, 0) до 6-го порядка >readlib(mtaylor): > f=mtaylor(sin(x^2+y^2), [x=0,y=0], 7);

20 Пример >taylor(cos(x),x,8): p:=convert(%,polynom); > plot([cos(x),p],x=- 2*Pi..2*Pi,color=[blue,green],thicknes s=[3,3],linestyle=[1,3]);

21