exceeds the number of array

Question

Jocelyn 2020-11-16

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/650058-exceeds-the-number-of-array

评论： Walter Roberson 2020-11-16

Hi,

I keep getting error messages about line 32 (Vy (i) = line). I double checked the brackets and they seem to match up. Below is the error message I am recieving.

Index exceeds the number of array elements (1).

Error in CC551_A5_P3 (line 32)

Vy(i) = Vx(i)*(tand(phi_o(i))-(g*t(i)/Vx_o)*(1+(0.5*Vx_o*k1*t(i)))); % y direction velocity

Thank you!

% Given Information In Problem
Vx_o = 2060; % muzzle velocity (ft/s)
k1 = 0.00082; % (1/ft) note: converted from 0.00025 (1/m)
p = 0.0751; % (lbm/ft^3) row - standard sea level met data 
a = 1120; % (ft/s) standard sea level met data
% initial conditions
t = 0; % (s)
g = 32.174; %gravitational constant (ft/sec^2) [at sea level]
%setting up inital conditions at i = 1  
   t(1) = 0;
   Vx(1) = 2060;
   x(1) = 0;
   phi_o(1) = 0;
   phi(1) = 0;
for i = 1:6
    
    x(i)= (i-1)*200*3;   % x3 = converting range from yards to feet
    
    Vx(i) = Vx_o*exp(-k1*x(i)); % striking velocity
        
    t(i) = ((x(i)/Vx_o)*((Vx_o/Vx(i))-1))/log(Vx_o/Vx(i)); % time of flight
    
    Vy(i) = Vx(i)*(tand(phi_o(i))-(g*t(i)/Vx_o)*(1+(0.5*Vx_o*k1*t(i)))); % y direction velocity
    
    phi(i) = (atand(Vy(i)/Vx(i)));                % impact angle  
    
    phi_o(i) = atand(tand(phi(i))+((g*t(i)/Vx_o)*(0.5*(1+(Vx_o/Vx(i))))));
end
Range_Table = table(x(:)/3, Vx(:), t(:), phi_o(:), phi(:),...
    'VariableNames',{'yards', 'impact velocity (ft/s)', 'time of flight (s)', 'inital QE angle (degrees)', 'impact angle (degrees)'})

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Walter Roberson 2020-11-16

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/650058-exceeds-the-number-of-array#answer_546453

编辑：Walter Roberson 2020-11-16

在 MATLAB Online 中打开

    Vy(i) = Vx(i)*(tand(phi_o(i))-(g*t(i)/Vx_o)*(1+(0.5*Vx_o*k1*t(i)))); % y direction velocity

That needs phi_o(i) to exist. But

    phi(i) = (atand(Vy(i)/Vx(i)));                % impact angle  
    phi_o(i) = atand(tand(phi(i))+((g*t(i)/Vx_o)*(0.5*(1+(Vx_o/Vx(i))))));
    

it does not exist until a few lines later.

Consider the possibility that you should be iterating i=2:6 or i=2:7 and that you should be indexing at (i-1) in places.

2 个评论
显示无隐藏无

Jocelyn 2020-11-16

在 MATLAB Online 中打开

I know that the phi and phi_o values should change value each time through the loop.

I was trying to set phi_o = 0 for the first time in the loop so that the Vy(i) = could be calculated, and then the phi and phi_o values could be calculated and these values could be used during the second iteration of the loop and so forth until the loop stops.

When I subbed in i -1 into the (i) and changed the loop for i = 2:6 the table produced NaN, and zeros for the phi and phi o columns

% Given Information In Problem
Vx_o = 2060; % muzzle velocity (ft/s)
k1 = 0.00082; % (1/ft) note: converted from 0.00025 (1/m)
p = 0.0751; % (lbm/ft^3) row - standard sea level met data 
a = 1120; % (ft/s) standard sea level met data
% initial conditions
t = 0; % (s)
g = 32.174; %gravitational constant (ft/sec^2) [at sea level]
%setting up inital conditions at i = 1  
   t(1) = 0;
   Vx(1) = 2060;
   x(1) = 0;
   phi_o(1) = 0;
for i = 2:6
    
    x(i)= (i-1)*200*3;   % x3 = converting range from yards to feet
    
    Vx(i) = Vx_o*exp(-k1*x(i-1)); % striking velocity
        
    t(i) = ((x(i-1)/Vx_o)*((Vx_o/Vx(i-1))-1))/log(Vx_o/Vx(i-1)); % time of flight
    
    Vy(i) = Vx(i-1)*(tand(phi_o(i-1))-(g*t(i-1)/Vx_o)*(1+(0.5*Vx_o*k1*t(i-1)))); % y direction velocity
    
    phi(i) = (atand(Vy(i-1)/Vx(i-1)));                % impact angle  ---> x60 to convert from degrees to minutes
    
    phi_o(i) = atand(tand(phi(i-1))+((g*t(i-1)/Vx_o)*(0.5*(1+(Vx_o/Vx(i-1))))));
end
Range_Table = table(x(:)/3, Vx(:), t(:), phi_o(:), phi(:),...
    'VariableNames',{'yards', 'impact velocity (ft/s)', 'time of flight (s)', 'inital QE angle (degrees)', 'impact angle (degrees)'})