Question

matlab求一阶，二阶导数，并画图

for t=0:0.005:5
x=4*cos(3*t);
y=4*sin(3*t);
z=2*t;
a=24*3^1/2;
b=12*3^1/2;
l1=((-3^1/2/12*b+3^1/2/12*a+x).^2+(-b/4+a/4+y).^2+z.^2)^1/2;
v11=diff(l1);
v1=diff(v11)
plot(v11)
plot(v1)
hold on
end
这样求出来的只的l1的数据，没有一阶，二阶导数的数据，并出的图是一条直线，有大侠指导一下不。谢谢你。

Answer 1

t=0:0.005:5;

x=4*cos(3*t);

y=4*sin(3*t);

z=2*t;

c=3^(1/2)/12;

a=288*c;

b=a/2;

l=((-c*b+c*a+x).^2+(-b/4+a/4+y).^2+z.^2).^(1/2);

n=length(t);

s1=diag(ones(1,n-1),-1)+diag(ones(1,n));

s1(:,end)=[];

s2=s1(1:end-1,1:end-1);

t1=t*s1/2;

t2=t1*s2/2;

v11=diff(l);

v12=diff(l,2);

subplot(2,1,1)

plot(t1,v11)

subplot(2,1,2)

plot(t2,v12)