plotting Y1(2t) to get the correct figure instead of Y1(t)

Question

Guillaume SIBILLE 2020-9-17

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/595648-plotting-y1-2t-to-get-the-correct-figure-instead-of-y1-t

回答： Ayush Gupta 2020-10-9

I've been trying to solve the following system of differentiel equations and to plot one of these functions.

The solver works correctly but I have an unusual problem.

I get the correct graph when I plot y1(2*t) (like I did on this code) instead of y1(t), and I don't know why this is the case. I'm pretty sure that the other functions I made (u, ha and h) are correctly set up.

If anyone has a clue, I'd be glad to hear it!

Thanks in advance

%parameter

We = 780

Re = 1145

t_e = 1.05*We^(-2/3)

s_e = sqrt(3*t_e)

v_e = 0.5* sqrt(3/t_e)

b_e = ((sqrt(12)*t_e^1.5)/pi)

%time interval and initial conditions

t_interval = linspace(t_e,10,1000) ;

init_cond = [s_e b_e v_e];

%solution

[t,y] = ode113(@(t,y) odefcn(t,y, We) , t_interval , init_cond);

%plots y1(2t)

plot(t(1:1:end/2),y(2:2:end,1),'b');

loglog(t(1:1:end/2),y(2:2:end,1),'b')

axis([0.1 10 0.5 5])

function dydt = odefcn(t,y, We)

dydt = [ y(3);

((u(y(1),t)-y(3))*(h(y(1),t)))*4/(y(2)*pi);

(((u(y(1),t)-y(3))^2)*h(y(1),t)-2/We)*8/(y(2)^2*pi);];

end

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

请先登录，再进行评论。

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

Answer 1

Ayush Gupta 2020-10-9

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/595648-plotting-y1-2t-to-get-the-correct-figure-instead-of-y1-t#answer_509411

The code works fine and there is no problem in it. One possible reason for this could be that these equations or functions are being mis-matched or wrong somewhere which is because of the context of the problem might not be understood fully.