0

我正在尝试编写一个matlab代码来模拟炮弹的弹丸运动，包括空气阻力和空气密度随温度变化的影响，但是我到目前为止的代码只计算炮弹轨迹的直线是不正确的。谁能指出我哪里出了问题并指出正确的方向，包括空气密度和温度的影响？谢谢。

clear;

%input parameters
v0=input('Enter the muzzle velocity (m/s) ');
theta=input('Enter the quadrant elevation (degrees) ');
%T0=input('Enter the value for the ground temperature in degreees ');
%T0=T0+275.16;

b2bym0=4e-5;
g=9.8;
dt=1e-2;

%define initial conditions
x0=0;
y0=0;
vx0=v0*cosd(theta);
vy0=v0*sind(theta);
fdragx=-b2bym0*v0*vx0;
fdragy=-b2bym0*v0*vy0;

n=1000; %iterations

%Tratio=(T0/300)^(2.5);

%define data array
%t=zeros(1000);

x=zeros(1000); %x-position

y=zeros(1000); %y-position

vx=zeros(1000); %x-velocity

vy=zeros(1000); %y-velocity



for i=1:n

    t(i)=i*dt;

    vx(i)=vx0+fdragx*dt;
    vy(i)=vy0+fdragy*dt;

    x(i)=x0+vx(i)*dt;
    y(i)=y0+vy(i)*dt;

    x0=x(i);
    y0=y(i);

    vx0=vx(i);
    vy0=vy(i);



end
plot(x,y,'g+')
4

2 回答 2

1

它看起来不像你正在为你的 y 速度建模 g 力的向下加速度

于 2013-10-09T00:27:46.973 回答
0

好像有三个问题。

首先，您在更新时错过了重力vy。固定的

第二阻力没有随着速度而更新。固定的

第三，您使用初始值而不是以前的值来计算新的位置/速度。在循环中尝试更改这些行。2:n如果 matlab 索引为 1，您可能必须更新您的 for 循环。

vx(i)=vx(i-1)+fdragx*dt;
vy(i)=vy(i-1)+(-g+fdragy)*dt;
x(i)=x(i-1)+vx(i)*dt;
y(i)=y(i-1)+vy(i)*dt;

编辑：没有看到初始条件的更新，忽略第三条评论。

于 2013-10-09T01:25:10.587 回答