Вычисление плоскости касательной, чтобы появиться

В этом примере показано, как аппроксимировать градиенты функции конечными разностями. Это затем показывает, как построить плоскость касательной к точке на поверхности при помощи этих аппроксимированных градиентов.

Создайте функцию $$f(x,y)=x^2+y^2$$ использование указателя на функцию.

f = @(x,y) x.^2 + y.^2;

Аппроксимируйте частные производные $$f(x,y)$$ относительно $$x$$ и $$y$$ при помощи gradient функция. Выберите длину конечной разности, которая совпадает с размером mesh.

[xx,yy] = meshgrid(-5:0.25:5);
[fx,fy] = gradient(f(xx,yy),0.25);

Плоскость касательной к точке на поверхности, $$P=(x_0,y_0,f(x_0,y_0))$$, дают

fx и fy матрицы являются приближениями к частным производным $$\frac{\partial f}{\partial x}$$ и $$\frac{\partial f}{\partial y}$$. Интересным местом в этом примере, где плоскость касательной соответствует функциональной поверхности, является (x0,y0) = (1,2). Значением функции в этом интересном месте является f(1,2) = 5.

Аппроксимировать плоскость касательной z необходимо найти значение производных при интересе. Получите индекс той точки и найдите аппроксимативные производные там.

x0 = 1;
y0 = 2;
t = (xx == x0) & (yy == y0);
indt = find(t);
fx0 = fx(indt);
fy0 = fy(indt);

Создайте указатель на функцию уравнением плоскости касательной z.

z = @(x,y) f(x0,y0) + fx0*(x-x0) + fy0*(y-y0);

Постройте исходную функцию $$f(x,y)$$, точка P, и часть плоского z это - касательная к функции в P.

surf(xx,yy,f(xx,yy),'EdgeAlpha',0.7,'FaceAlpha',0.9)
hold on
surf(xx,yy,z(xx,yy))
plot3(1,2,f(1,2),'r*')

Посмотрите профиль стороны.

view(-135,9)

Смотрите также

gradient

Похожие темы

• Создания указателя на функцию