微分方程的求解MATLAB

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

 

评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值