如何利用matlab求微分方程空间曲线的密切平面、法平面，从切平面？

Question

问题相关代码，请勿粘贴截图

 syms a b t;
r=[a(cos(t)),a(sin(t)),b*t];
r1=diff(r,t);%求r的一阶导
r2=diff(r1,t);%求r的二阶导
r3=diff(r2,t);%求r的三阶导
A1=r1/mo(r1);
A=simplify(A1);%求单位切向量
T1=cross(r1,r2)/mo(cross(r1,r2));
T=simplify(T1);%求副法向量
B1=cross(T,A);
B=simplify(B1);%求主法向量
R=(X,Y,Z);
cross(R-r,A,B)=0;%求密切平面方程
d=cross(R-r,A);
cross(d,B)=0;
dot((R-r1),A)=0;%求法平面
dot((R-r1),B)=0;%求副法平面

运行结果及报错内容

Answer 1

你好，修改如下：

 syms a b t;
r=[a*(cos(t)),a*(sin(t)),b*t];
r1=diff(r,t);%求r的一阶导
r2=diff(r1,t);%求r的二阶导
r3=diff(r2,t);%求r的三阶导
A1=r1/norm(r1);
A=simplify(A1);%求单位切向量
T1=cross(r1,r2)/norm(cross(r1,r2));
T=simplify(T1);%求副法向量
B1=cross(T,A);
B=simplify(B1);%求主法向量
syms X Y Z
R=[X,Y,Z];
SF = cross(A,B)*(R-r)'%求密切平面方程
d = cross(R-r,A);
SF0 = cross(d,B)
SF2 = dot((R-r1),A)%求法平面
SF3 = dot((R-r1),B)%求副法平面