Filter löschen
Filter löschen

How to find gradient of a vector field in matlab symbolic

19 Ansichten (letzte 30 Tage)
Chandan
Chandan am 12 Dez. 2023
Beantwortet: Chandan am 13 Dez. 2023
I am trying to find gradient of a vector field in matlab symbolic , whose output will be matrix but it am getting error
  2 Kommentare
KSSV
KSSV am 12 Dez. 2023
Copy your code here. Don't attach it as an screen shot/ image.
Chandan
Chandan am 12 Dez. 2023
clc ; clear all
syms x y z
syms u(x,y,z) v(x,y,z) w(x,y,z)
V = [u v w];
S = [x y z];
gradient(u,S)
ans(x, y, z) = 
gradient(V,S) % How to take gradient of vector field
Error using sym/gradient
Invalid argument at position 1. Argument must be scalar or a symbolic function with scalar formula.

Melden Sie sich an, um zu kommentieren.

Akzeptierte Antwort

Walter Roberson
Walter Roberson am 12 Dez. 2023
Verschoben: Walter Roberson am 12 Dez. 2023
The fundamental problem you are having is that gradient does not accept a vector the first parameter.
syms x y z
syms u(x,y,z) v(x,y,z) w(x,y,z)
V = [u(x,y,z) v(x,y,z) w(x,y,z)];
S = [x y z];
temp = arrayfun(@(EXPR) gradient(EXPR,S), V, 'uniform', 0);
result(x,y,z) = [temp{:}]
result(x, y, z) = 

Weitere Antworten (2)

Sulaymon Eshkabilov
Sulaymon Eshkabilov am 12 Dez. 2023
If you assign an expression for V, you will get this:
clc ; clearvars
syms x y z
syms u(x,y,z) v(x,y,z) w(x,y,z)
V(x,y,z) = 2*u+3*v-w % Some e.g. expression
V(x, y, z) = 
S = [x y z];
du = gradient(u,S)
du(x, y, z) = 
dV = gradient(V,S)
dV(x, y, z) = 
% OR simply
dV= gradient(V,[x,y,z])
dV(x, y, z) = 
  1 Kommentar
Dyuman Joshi
Dyuman Joshi am 12 Dez. 2023
V is not a combination of u, v and w, but an array with u, v and w as elements.

Melden Sie sich an, um zu kommentieren.


Chandan
Chandan am 13 Dez. 2023
% thanks for answer. It worked with the Jacobian too.
clc; clearvars;
syms x y z
syms u(x, y, z) v(x, y,z) w(x,y,z)
V = [u v w];
S = [x y z];
jacobian(V,S)
ans(x, y, z) = 

Kategorien

Mehr zu Symbolic Math Toolbox finden Sie in Help Center und File Exchange

Produkte


Version

R2023a

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by