Question

请求高手帮忙看下这个方程组应该怎么解求大神帮助

m1*D2x1-F+2*k1*(x1-x2+h1)+2*A*N1*f+A*N3*f=0m2*D2x2-2*k1*(x1-x2+h1)+k2*x2-2*A*N2*f-2*A*N3*f+FM=0N2=m1*g*(l1+x1-x2)/l2N1=m1*g-N2N3=(F*l6-2*F*(l6-l5)-2*K1(x1-x2+h1)*l6+2*K1*(x1-x2+h1)*D-A*f*N1*l6-A*N2*f*l6)/(2*l2-2*f*l6+2*A*f*l6)当F>=2*k1*(x1-x2+h1)时：D=l5当F<2*k1*(x1-x2+h1)时：D=l6-l5 （l为小写L)当Dx1-Dx2>=0时：A=1当Dx1-Dx2<0时：A=-1以上是方程组，求的是x1和x2，其它的为系数。请各位大哥，高手帮帮忙，小弟是初学者，还望各位高手赐教

Answer 1

你说的使用那些字母直接表示，意思是说使用Matlab的符号工具箱吧，呵呵，Matlab还没有那么智能，这么复杂的问题不可能直接代劳还有一点需要说明的，你解的是微分方程，在数值求解的时候必须给出初值条件x0和微分区间tspan，这次我直接任意给一个了，你根据需要自己修改function matlabskyclctspan=[0 10];x0=rand(4,1);[t,x]=ode45(@odefun,tspan,x0);plot(t,x)function dx=odefun(t,x)%x=[x1 x2 x1' x2']F=150000;K1=5000;K2=8000;h1=0.05;l1=0.2;l2=0.5;l3=0.4;l4=0.2;l5=0.3;l6=0.35;m1=12;m2=10;M=50;f=0.16;g=9.8;x1=x(1);x2=x(2);Dx1=x(3);Dx2=x(4);if F>=2*K1*(x1-x2+h1) D=l5;else D=l6-l5;endif Dx1-Dx2>=0 A=1;else A=-1;endN2=m1*g*(l1+x1-x2)/l2;N1=m1*g-N2;N3=(F*l6-2*F*(l6-l5)-2*K1*(x1-x2+h1)*l6+2*K1*(x1-x2+h1)*D-A*f*N1*l6-A*N2*f*l6)/(2*l2-2*f*l6+2*A*f*l6);dx=[x(3) x(4) -(-F+2*K1*(x1-x2+h1)+2*A*N1*f+A*N3*f)/m1 -(-2*K1*(x1-x2+h1)+K2*x2-2*A*N2*f-2*A*N3*f+F*M)/m2]; 查看原帖>>