matlab用龙格库塔法求解微分方程组题目是这样的:将微分方程组重写为一阶以后是这样的初始矢量不知道在matlab中应该怎么写,分数不多,但是是我有的全部了 TAT
来源:学生作业帮助网 编辑:作业帮 时间:2024/06/27 22:08:19
![matlab用龙格库塔法求解微分方程组题目是这样的:将微分方程组重写为一阶以后是这样的初始矢量不知道在matlab中应该怎么写,分数不多,但是是我有的全部了 TAT](/uploads/image/z/7830813-21-3.jpg?t=matlab%E7%94%A8%E9%BE%99%E6%A0%BC%E5%BA%93%E5%A1%94%E6%B3%95%E6%B1%82%E8%A7%A3%E5%BE%AE%E5%88%86%E6%96%B9%E7%A8%8B%E7%BB%84%E9%A2%98%E7%9B%AE%E6%98%AF%E8%BF%99%E6%A0%B7%E7%9A%84%EF%BC%9A%E5%B0%86%E5%BE%AE%E5%88%86%E6%96%B9%E7%A8%8B%E7%BB%84%E9%87%8D%E5%86%99%E4%B8%BA%E4%B8%80%E9%98%B6%E4%BB%A5%E5%90%8E%E6%98%AF%E8%BF%99%E6%A0%B7%E7%9A%84%E5%88%9D%E5%A7%8B%E7%9F%A2%E9%87%8F%E4%B8%8D%E7%9F%A5%E9%81%93%E5%9C%A8matlab%E4%B8%AD%E5%BA%94%E8%AF%A5%E6%80%8E%E4%B9%88%E5%86%99%2C%E5%88%86%E6%95%B0%E4%B8%8D%E5%A4%9A%2C%E4%BD%86%E6%98%AF%E6%98%AF%E6%88%91%E6%9C%89%E7%9A%84%E5%85%A8%E9%83%A8%E4%BA%86+TAT)
matlab用龙格库塔法求解微分方程组题目是这样的:将微分方程组重写为一阶以后是这样的初始矢量不知道在matlab中应该怎么写,分数不多,但是是我有的全部了 TAT
matlab用龙格库塔法求解微分方程组
题目是这样的:
将微分方程组重写为一阶以后是这样的
初始矢量
不知道在matlab中应该怎么写,
分数不多,但是是我有的全部了 TAT
matlab用龙格库塔法求解微分方程组题目是这样的:将微分方程组重写为一阶以后是这样的初始矢量不知道在matlab中应该怎么写,分数不多,但是是我有的全部了 TAT
解方程很简单,调用ode45即可
关键是求运动时间,假设z方向坐标变为0即为落至地面运动结束,可求得运动时间的近似值及落地点坐标
代码如下:
clear all
clc
c=0.07;
u0=25;
a=pi/6;
f=@(t,x)([x(2);
-c*sqrt(x(2)^2+(x(4)-7-.35*x(5))^2+x(6)^2)*x(2);
x(4);
-c*sqrt(x(2)^2+(x(4)-7-.35*x(5))^2+x(6)^2)*(x(4)-7-.35*x(5));
x(6);
-9.81-c*sqrt(x(2)^2+(x(4)-7-.35*x(5))^2+x(6)^2)*x(6)]);
[t,Y]=ode45(f,0:.01:2,[0 u0*cos(a) 0 0 1.4 u0*sin(a)]);
n=max(find(Y(:,5)>0));
t_end=t(n)
s=[Y(n,1),Y(n,3),0]
plot3(Y(1:n,1),Y(1:n,3),Y(1:n,5))
其中t_end为运动时间
s为落地点坐标
看了半天,怎么都觉得方程不太对
再检查一下方程吧,尤其是第二个方程等号右边到底有没有负号
最后一个方程等号右边也有问题,请仔细检查