电子产业一站式赋能平台

PCB联盟网

搜索
查看: 944|回复: 1
收起左侧

ODE45解算不符合预期

[复制链接]

531

主题

913

帖子

4563

积分

四级会员

Rank: 4

积分
4563
发表于 2022-10-3 12:37:30 | 显示全部楼层 |阅读模式
背景:为喷水射流轨迹建模,x(1)初始速度,x(2)俯仰角度,x(3)射程,x(4)射高;
预期如下图:



ODE45解算不符合预期-1.jpg
现状:除了45°等寥寥无几的角度,其他射流轨迹都在二、三、四象限,解还带复数。

具体代码如下:Any thoughts would be appreciated!
定义函数
function dx = myfunc2226(t,x)
     dx = zeros(4,1);
     p=1.225;g=9.8;A0=0.00099;a=3.8;m=0.0503;
     F=0.5*p*x(1)^2*A0*(1+a*log(1+x(3)))*0.15;
     dx(1)=-F/m-g*sin(x(2));
     dx(2)=-g*cos(x(2))/x(1);
     dx(3)=x(1)*cos(x(2));
     dx(4)=x(1)*sin(x(2));
end


运行ode45

x0=[50.81,30,0,0]; %这里的第二个参数是θ 俯仰角度,设置为45°,则可以正常得到一个第一象限的曲线,但是30°的轨迹就在第四象限
[t,x]=ode45('myfunc2226',[0,8],x0);
n=size(x,1)
plot(x(1:n,3),x(1:n,4),'-');
回复

使用道具 举报

553

主题

1695

帖子

5576

积分

四级会员

Rank: 4

积分
5576
发表于 2022-10-3 12:38:10 | 显示全部楼层
ODE45解算不符合预期-1.jpg
补充一下解的现状
回复 支持 反对

使用道具 举报

发表回复

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

关闭

站长推荐上一条 /1 下一条


联系客服 关注微信 下载APP 返回顶部 返回列表