0

我正在尝试使用 ode45 函数求解一个微分方程，该函数涉及两个自由度，即 x 和 y 方向。在下面的 oem 函数中，我试图将需要的所有四个变量存储到矩阵中以插入 ode45，但我无法让 oem 函数识别 ds(3) 或 ds(4)。我试过解决它，但我碰壁了。

function odeprob2 
%ode45 function [T, S] = ode45(@oem, linspace(0,0.1,100), [0,2]);
fprintf('Final Position in the X Direction = %4.2f \n', S(end,1));
fprintf('Final Velocity in the X Direction= %4.2f \n', S(end,2));
fprintf('Final Position in the Y Direction = %4.2f \n', S(end,3));
fprintf('Final Velocity in the Y Direction= %4.2f \n', S(end,4));
plot(T, S(:,1), '-k')
xlabel('Time, sec')
ylabel('Position, m')

function ds = oem(t, s)
  %Constants: 
  rad = 15; %mm
  dens = 8000; %kg/m^3
  deno = 980; %kg/m^3
  g = 9.81 ; %m/s^2
ds = zeros(4, 1);
ds(1) = s(2);
ds(2) = ((-1.6) * s(2)) / (((4/3)*pi*(rad^3))*dens);
ds(3) = s(4);
ds(4) = ((((4/3) * pi * (rad^3)) * (dens - deno) * g) - (1.6 * s(4))) / (((4/3) *pi*(rad^3))*dens);
4

0 回答 0