What am I missing here?

Zifeng Qiu

2020 7 11

0 个回答

5 次查看（30 天）

0 个投票

This is a Matlab grader problem for calculating the jumper's distance, I follow the code that my professor was provided, what am I missing here to get the final disnace?

function Dist = BJump
    z0 = [0;0;pi/8;10];
    dt = 0.1;
    
    T = zeros(1,1000);
    T(1) = 0;
    Z = zeros(4,1000);
    Z(:,1) = z0;
    
    for j = 1:10000-1
        K1 = physics(T(j),Z(:,j));
        K2 = physics(T(j) + dt/2,Z(:,j) + dt/2*K1);
        K3 = physics(T(j) + dt/2,Z(:,j) + dt/2*K2);
        K4 = physics(T(j) + dt,Z(:,j) + dt*K3);
        Z(:,j+1) = Z(:,j) + dt/6*(K1 + 2*K2 + 2*K3 + K4);
        T(j+1) = T(j) + dt;
    end
    
    function dzdt=physics(t,z)
        dzdt = 0*z;
        dzdt(1) = z(4)*cos(z(3));
        dzdt(2) = z(4)*sin(z(3));
        dzdt(3) = -9.81/z(4)*cos(z(3));
        D = (0.72)*(0.94)*(0.5)/2*(dzdt(1)^2 + dzdt(2)^2);
        dzdt(4) = -D/80-9.81*sin(z(3));
    end
end