clear all;close all;
syms x y z ax ay az ; %定义x,y,z为符号变量;
disp(*请将您输入的函数转换为f=Ax*ax+Ay*ay+Az*az的标准格式格式(其中Ax,Ay,Az代表函数式)*); %转化为笛卡尔坐标系中矢量的表达式;
s=input(*请输入Ax部分函数 Ax=*,*s*); %输入ax方向上的表达式Ax;
Ax=s;
dx=diff(s,x); %求Ax对x的偏导数dAx/dx;
a1=Ax;
b1=dx/Ax;
s=input(*请输入Ay部分函数 Ay=*,*s*); %输入ay方向上的表达式Ay;
Ay=s;
dy=diff(s,y); %求Ay对y的偏导数dAy/dy;
a2=Ay;
b2=dy/a2;
s=input(*请输入Az部分函数 Az=*,*s*) %输入az方向上的表达式Az;
Az=s;
dz=diff(s,z); %求Az对z的偏导数dAz/dz;
a3=Az;
b3=dz/Az;
divA=dx+dy+dz %定义函数divA等于dAx/dx加dAy/dy加dAz/dz;
curlA=[ax,ay,az;b1,b2,b3;a1,a2,a3]
syms x y z ax ay az ; %定义x,y,z为符号变量;
disp(*请将您输入的函数转换为f=Ax*ax+Ay*ay+Az*az的标准格式格式(其中Ax,Ay,Az代表函数式)*); %转化为笛卡尔坐标系中矢量的表达式;
s=input(*请输入Ax部分函数 Ax=*,*s*); %输入ax方向上的表达式Ax;
Ax=s;
dx=diff(s,x); %求Ax对x的偏导数dAx/dx;
a1=Ax;
b1=dx/Ax;
s=input(*请输入Ay部分函数 Ay=*,*s*); %输入ay方向上的表达式Ay;
Ay=s;
dy=diff(s,y); %求Ay对y的偏导数dAy/dy;
a2=Ay;
b2=dy/a2;
s=input(*请输入Az部分函数 Az=*,*s*) %输入az方向上的表达式Az;
Az=s;
dz=diff(s,z); %求Az对z的偏导数dAz/dz;
a3=Az;
b3=dz/Az;
divA=dx+dy+dz %定义函数divA等于dAx/dx加dAy/dy加dAz/dz;
curlA=[ax,ay,az;b1,b2,b3;a1,a2,a3]