How can find t?

Question

Gimmy Morales 2015-2-15

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/178578-how-can-find-t

编辑： Star Strider 2015-2-15

I did my homework to solve a problem with Heun's method. Thing is, I want to show the value of t(my x-axis) when v=0, with "index= find(v==0)" doesn't work... How can i find my value of t? in this program?

% Método de Heun
clear all
clc
a= 0;  b= 20;
n= 40;   v0= 300;   h = (b - a) / n;
v(1,:) = v0; t(1) = a;
f= @(t,v) -0.0035*v^2 -3
for i = 1 : n
    t(i+1) = t(i) + h
    g = f(t(i),v(i,:))
    z = v(i,:) + h * g
    v(i+1,:) = v(i,:) + h/2 * ( g + f(t(i+1),z) )
end;
plot(t, v)
whitebg('k')
title('V vs T,  dy/dx= -0.0035v²-3', 'Fontsize', 18)
xlabel('Tiempo', 'Fontsize', 16) 
ylabel('Velocidad', 'Fontsize', 16)
find(y==10)

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Gimmy Morales 2015-2-15

在 MATLAB Online 中打开

    % % Método de Heun
clear all
clc
a= 0;  b= 20;
n= 40;   v0= 300;   h = (b - a) / n;
v(1,:) = v0; t(1) = a;
f= @(t,v) -0.0035*v^2 -3;
for i = 1 : n
    t(i+1) = t(i) + h;
    g = f(t(i),v(i,:));
    z = v(i,:) + h * g;
    v(i+1,:) = v(i,:) + h/2 * ( g + f(t(i+1),z) );
end;
plot(t, v);
whitebg('k');
title('V vs T,  dy/dx= -0.0035v²-3', 'Fontsize', 18);
xlabel('Tiempo', 'Fontsize', 16) ;
ylabel('Velocidad', 'Fontsize', 16);
t=t';
twv0=[t     v]

请先登录，再进行评论。

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

Answer 1

Star Strider 2015-2-15

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/178578-how-can-find-t#answer_168084

编辑：Star Strider 2015-2-15

在 MATLAB Online 中打开

I did not run your code, but your way of finding the value of the index of ‘t’ where ‘v’ crosses zero is almost correct.

See if:

index1 = find(v<=0, 1, 'first');

or:

index2 = find(v>=0, 1, 'last');

brings you closer to what you want. Those indices should be below and above the indices of the value of ‘t’ for which ‘v’ crosses zero.

You have to deal with inexact numeric representations that are part of floating-point arithmetic, as well as your velocity variable not being exactly zero in your calculations.

Note: your function has to have at least one zero-crossing for this to work.