Need some help finding finding the gradient and creating a vector space to use contour and quiver plots.

Question

Hunter Manning 2022-10-19

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1831048-need-some-help-finding-finding-the-gradient-and-creating-a-vector-space-to-use-contour-and-quiver-pl

回答： Nithin Kumar 2023-9-8

The potential of an electric field due to a point charge is given (in spherical coordinates) by V= k/r where k is a constant. Assuming k = 1, find V in cartesian coordinates. Use Matlab to determine the gradient in cartesian coordinates via the symbolic toolbox , and then create a vector space and use the contour and quiver plots to illustrate the gradient of the function in the z = 0 plane.

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Nithin Kumar 2023-9-8

在 MATLAB Online 中打开

Hi Hunter,

I understand that you want to determine the gradient in cartesian coordinates using Symbolic toolbox. To perform this analysis in MATLAB, kindly refer to the following steps:

1. Define the symbolic variables and constants:

syms x y z r k;
k = 1; % Assuming k = 1
r = sqrt(x^2 + y^2 + z^2); % Define r in Cartesian coordinates

2. Calculate the electric potential 'V':

V = k / r;

3. Calculate the gradient of 'V' in Cartesian coordinates:

gradV = gradient(V, [x, y, z]);

4. Create a vector space in the xy-plane (z = 0) and evaluate the gradient:

[x_vals, y_vals] = meshgrid(-5:0.5:5, -5:0.5:5);
z_val = 0;
gradV_eval = subs(gradV, [x, y, z], [x_vals, y_vals, z_val]);

5. Create "contour" and "quiver" plots to illustrate the gradient in the "z = 0" plane:

figure;
% Contour plot of the electric potential
subplot(1, 2, 1);
contour(x_vals, y_vals, double(subs(V, [x, y, z], [x_vals, y_vals, z_val])));
title('Contour Plot of Electric Potential (V)');
% Quiver plot of the gradient
subplot(1, 2, 2);
quiver(x_vals, y_vals, double(gradV_eval(:,:,1)), double(gradV_eval(:,:,2)));
title('Quiver Plot of Gradient of V');
xlabel('x');
ylabel('y');
% Adjust plot parameters as required
axis equal;