function dx=cili1(e,x)
dx=zeros(4,1);
syms m n e f p;
m=0.015;
p=0.005;
n=0.005;
f=sqrt(1+p-0.1*(1+3/2*0.1));
e=[-0.2 0.2];
R1=-10*p*p/(3*f*f)+3*p;
R2=-44*p*p/(3*f*f)+2*p;
R3=-2*p*p/(f*f)+p;
R4=-30*p*p/(f*f)+3*p;
dx(1)=-0.5*m*x(1)-R3/(8*f)*x(1)*x(3)*x(3)*sin(2*(x(4)-x(2)))-0.5*n*f*sin(x(2));
dx(2)=-e+R1/(8*f)*x(1)*x(1)+R2/(8*f)*x(3)*x(3)+R3/(8*f)*x(3)*x(3)*cos(2*(x(4)-x(2)))-0.5/(x(1))*n*f*cos(x(2));
dx(3)=-0.5*m*x(3)-R3/(8*f)*x(3)*x(1)*x(1)*sin(2*(x(2)-x(4)))-0.5*n*f*cos(x(4));
dx(4)=-e+R4/(8*f)*x(3)*x(3)+R2/(8*f)*x(1)*x(1)+R3/(8*f)*x(1)*x(1)*cos(2*(x(2)-x(4)))+0.5/(x(2))*n*f*cos(x(4));
end
最后
以上就是粗暴唇膏最近收集整理的关于微分方程的求解MATLAB的全部内容，更多相关微分方程内容请搜索靠谱客的其他文章。
本图文内容来源于网友提供，作为学习参考使用，或来自网络收集整理，版权属于原作者所有。
发表评论 取消回复