用ode45函数可以求解这个微分方程组。其自定义函数为
dx(1)=x(2);
dx(2)=-2*beta*x(2)-x(1)-mu_s*sign(x(2)-v0)+3*(mu_s-mu_m)/(2*vm)*(x(2)-v0)-(mu_s-mu_m)/(2*vm^3)*(x(2)-v0)^3;
求解结果
τ X1(τ) X2(τ)
0 0 0
0.1 0.00049648 0.00017502
0.2 0.00049883 -0.0019544
0.3 0.00049962 -0.00023919
0.4 0.00047744 0.0022456
0.5 0.00045653 -0.0016108
0.6 0.00043671 -0.0001398
0.7 0.00043339 0.0021207
0.8 0.0004207 0.0011055
0.9 0.00041579 -0.0024293
1 0.00041282 0.0002398
。。。。。。。。。。。
https://pan.baidu.com/s/1iu7UlW4rRn5P5Re1DYIt6Q,