For loop to find roots

3 次查看（过去 30 天）
Giovanni Curiazio
评论： Askic V 2023-2-6
I am trying to find the roots of a function, with a parameter that changes. The function is that of a parabola
V=rand(10,1)
g=9.81
z0=500
I should use the following expression in for loop right?
t_roots=roots([-.5*g,V,z0])

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

回答（1 个）

Askic V
Askic V 2023-2-6
Ran in:
Yes, if V contains paramters, then the could would look like this:
V = rand(10,1);
g = 9.81;
z0 = 500;
for vi = 1: numel(V)
t_roots = roots([-0.5*g,vi,z0])
end
t_roots = 2×1
10.1988 -9.9950
t_roots = 2×1
10.3023 -9.8946
t_roots = 2×1
10.4068 -9.7952
t_roots = 2×1
10.5124 -9.6969
t_roots = 2×1
10.6189 -9.5995
t_roots = 2×1
10.7265 -9.5033
t_roots = 2×1
10.8351 -9.4080
t_roots = 2×1
10.9448 -9.3138
t_roots = 2×1
11.0554 -9.2205
t_roots = 2×1
11.1671 -9.1283
  5 个评论
显示 3更早的评论
Giovanni Curiazio
Giovanni Curiazio 2023-2-6
编辑：Giovanni Curiazio 2023-2-6
Now I would like to achieve the same thing but considering more debris: so the starting height remains the same but at the initial velocity I would like to add up all the elements of DeltaV_c (for each direction). And finally plot the different parabolas. In fact as you can see, in the code I used only the first line of DeltaV_c, now I would like to use all the lines so that I have 10 different speeds, so 10 parabolas
Askic V
Askic V 2023-2-6
Ran in:
I guess you are trying to achive something like this?
g=9.80665; %agravity
z0= 14325; %starting height
x0=0;
y0=0;
V_0=200*.3048; %starting velocity
gamma=0.78;
DeltaV_c = randn(10,3); %Speed increments in the three directions x, y, z
% Debris velocity components after the explosion
[r,c] = size(DeltaV_c);
for i = 1:r % each row contains increments in directions
V_0x=V_0*cos(gamma)+DeltaV_c(i,1);
V_0y=0+DeltaV_c(i,2);
V_0z=V_0*sin(gamma)+DeltaV_c(i,3);
%debris flight time
t_radici = roots([-.5*g,V_0z,z0]);
t = linspace(0,t_radici(t_radici>0),1000);
%distanza percorsa nelle tre direzioni
x=x0+V_0x*t;
y=y0+V_0y*t;
z=z0+V_0z*t-.5*g*t.^2;
plot3(x,y,z)
hold on
plot3(x0,y0,z0,'ro','MarkerSize',5)
plot3(x(end),y(end),z(end),'ro','MarkerSize',5)
xlabel('Posizione x (m)')
ylabel('Posizione y (m)')
zlabel('Posizione z (m)')
grid on
end

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Operators and Elementary Operations 的更多信息

标签

产品


版本

R2022a

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by