MATLAB График зависимости переменной в оде
У меня проблема в MATLAB. Я пытаюсь построить переменную внутри функции Оде. Я называю свою функцию с этим кодом:
%Above is the code for the IC and constants
%Initial conditions [x0, vx0, y0, vy0]
IC = [x_0; w_x; y_0; w_y];
%Simulationszeitvektor
tspan = 0:dt:Tsim;
%Anonymous Function
a =@(t,z) Parabel(z,m,g,u_x,u_y,r_p,A,rho_L,nv);
%Options
options = odeset('Events',@Ereignis,'Refine',4,'RelTol', 1e-12);
%Solver
[t,z] = ode45(a,tspan, IC, options);
Проблема в том, что число Рейнольдса Re. Число Рейнольдса зависит от скорости. Скорость уменьшается во время сопротивления. Перетаскивание зависит от числа Рейнольдса. Скорость рассчитывается с помощью решателя ode45 из Matlab. Файл функции выглядит следующим образом:
%My Function File
function dz = Parabel(z,m,g,u_x,u_y,r_p,A,rho_L,nv)
dz = zeros (4,1);
% relativ velocity
% z(2) is the velocity in x-direction
% z(4) is the velocity in y-direction
v_x=z(2)-u_x;
v_y=z(4)-u_y;
v=sqrt(v_x^2+v_y^2);
% Calculation Reynoldszahl
Re = (2*r_p*v)/nv;
%Drag Force
F_k = Re*constant
dz(2) = -F_k*constant;
dz(4) = -F_k*constant-g;
%The other code part is not necessary
То, что я хочу сделать, это построить Re-число Рейнольдса, потому что мне действительно нужно знать число. Я не знаю, как я могу это сделать. Мне нужно сохранить Re в рабочей области, как переменная z.