Question

请问这个问题如何用matlab解决

function ydoc=myfun (t,y)
k1=1;k2=1;k3=0.5;k4=1.5;k5=5;
l1=0.01;l2=0.6;l3=0.2;
a1=0.75;a2=0.75;
vd=-k1*tanh(l1*y(2));
wd=-k2*tanh(l2*y(3))+k3*tanh(l3*y(1))*sin(t);
ydoc=[(y(5)+wd)*y(2);
-(y(5)+wd)*y(1)+y(4)+vd;
y(5)+wd;
-k4*sign(y(4))*(abs(y(4)))^a1;
-k5*sign(y(5))*(abs(y(5)))^a2];
问题是：1.画图（doc(vd), doc(wd)）-----vd的导数和wd 的导数
的图像怎么画？？？
2. 我也试过利用方程把这两个导数的表达式解出来，可是太复杂了，太长，老出错；要么就是调用这个函数时说一些量没有定义

Answer 1

function ydoc=myfun (t,y)
if nargin<1
t=-pi:0.1:pi;
y=randperm(5);
end
k1=1;k2=1;k3=0.5;k4=1.5;k5=5;
l1=0.01;l2=0.6;l3=0.2;
a1=0.75;a2=0.75;
vd=-k1*tanh(l1*y(2));
wd=-k2*tanh(l2*y(3))+k3*tanh(l3*y(1))*sin(t);
p1=(y(5)+wd)*y(2);
p2=-(y(5)+wd)*y(1)+y(4)+vd;
p3=y(5)+wd;
p4=-k4*sign(y(4))*(abs(y(4)))^a1;
p5=-k5*sign(y(5))*(abs(y(5)))^a2;
ydoc=cat(2,p1,p2,p3,p4,p5);
L=1:length(ydoc);
plot(L,ydoc,'g-')

Answer 2

用曲线拟合工具箱。还可以自动产生代码