matlab ode45求解微分方程组
来源:学生作业帮 编辑:大师作文网作业帮 分类:数学作业 时间:2024/11/10 11:02:29
matlab ode45求解微分方程组
syms Dz H st V
m=300;p=2500;v0=250;x0=0;H0=7000;S=0.45;L=2.5;R0=1.2495;T0=288.15;
微分方程组:
dV=(p*cos(0.24*Dz)-(0.2+0.005*(0.24*Dz)^2)*0.5*R0*((T0-0.0065*H)/T0)^4.2588*V^2*S-m*0.5*R0*((T0-0.0065*H)/T0)^4.2588*V^2*sin(st))/m;
dst=(p*sin(0.24*Dz)+(0.25*0.24*Dz+0.05)*0.5*R0*((T0-0.0065*H)/T0)^4.2588*V^2*S-m*0.5*R0*((T0-0.0065*H)/T0)^4.2588*V^2*cos(st))/(m*V);
dx=V*cos(st);
dH=V*sin(st);
dH=(Dz-k*(H-(2000*cos(0.000314*1.1*x+5000))+k1*(-(796438175382409875*sin((6371505403059279*x)/18446744073709551616))/1152921504606846976))/k;
syms Dz H st V
m=300;p=2500;v0=250;x0=0;H0=7000;S=0.45;L=2.5;R0=1.2495;T0=288.15;
微分方程组:
dV=(p*cos(0.24*Dz)-(0.2+0.005*(0.24*Dz)^2)*0.5*R0*((T0-0.0065*H)/T0)^4.2588*V^2*S-m*0.5*R0*((T0-0.0065*H)/T0)^4.2588*V^2*sin(st))/m;
dst=(p*sin(0.24*Dz)+(0.25*0.24*Dz+0.05)*0.5*R0*((T0-0.0065*H)/T0)^4.2588*V^2*S-m*0.5*R0*((T0-0.0065*H)/T0)^4.2588*V^2*cos(st))/(m*V);
dx=V*cos(st);
dH=V*sin(st);
dH=(Dz-k*(H-(2000*cos(0.000314*1.1*x+5000))+k1*(-(796438175382409875*sin((6371505403059279*x)/18446744073709551616))/1152921504606846976))/k;
1、后两个方程都是dH?
2、最后一个dH的表达式有误(括号不匹配).
3、st的初值?
2、最后一个dH的表达式有误(括号不匹配).
3、st的初值?