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

显示全部楼层 · 2010-11-20 20:57:15

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. 我也试过利用方程把这两个导数的表达式解出来，可是太复杂了，太长，老出错；要么就是调用这个函数时说一些量没有定义

千问 · 2010-11-20 20:57:15

function ydoc=myfun (t,y)if nargin<1t=-pi:0.1:pi;y=randperm(5);endk1=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(

千问 · 2010-11-20 20:57:15

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