Need help! I keep getting the error "not enough input arguments" how do I fix this?
조회 수: 3 (최근 30일)
이전 댓글 표시
function drv = untitiled18(t,rv);
% small r is position
% big r is distance
mu = 3.986004*(10^5); %%unit is in km^3/s^2, and is Earth's Gravitational parameter
drv = zeros(6,1);
rx = rv(1);
vx = rv(2);
ry = rv(3);
vy = rv(4);
rz = rv(5);
vz = rv(6);
R = sqrt((rx^2)+(ry^2)+(rz^2));
% X component/ i direction
drv(1) = rv(2);
drv(2) = (-mu.*rx)/(R^3); % (vx) velocity in the x direction
% Y component/ j direction
drv(3) = rv(4);
drv(4) = (-mu.*ry)/(R^3); % (vy) velocity in the y direction
% Z component/ k direction
drv(5) = rv(6);
drv(6) = (-mu.*rz)/(R^3); % (vz) velocity in the z direction
end
댓글 수: 7
Star Strider
2018년 9월 10일
I agree. I doubt there are any significant differences between those versions with respect to ode45.
What line is throwing the error?
답변 (2개)
Star Strider
2018년 9월 10일
The only thing I can think of is that you have a function called ‘rv’ somewhere. To find it, run this:
q = which('rv','-all')
댓글 수: 2
Star Strider
2018년 9월 10일
You could run your ‘propagate2BP’ function with the solved values to get some of those (I do not understand what variables those would be).
It would be easier to begin with the solved position values and use numerical differentiation one time to then get the velocity values, and two times to get acceleration.
Use the gradient (link) function to do the numerical differentiation. This will be more accurate if you define ‘t’ as:
t = linspace(0, 172800, 21681);
or something similar, to get equally-spaced solved position values.
In differentiating a matrix, gradient differentiates across both the rows and columns, and the first output is the derivative down the columns. I would do the differentiation as:
drv = gradient(rv, mean(diff(t)), 1);
Note that this differentiates all the columns, so choose the original position columns from the ‘drv’ output.
Walter Roberson
2018년 9월 10일
편집: Walter Roberson
2018년 9월 10일
The error line you post is
[t,rv] = ode45(propagate_2BP22, t, [rinputx, vinputx, rinputy, vinputy, rinputz, vinputz], options);
which is not the same as what you posted before that, which was
[t,rv] = ode45(@propagate2BP, t, [rinputx, vinputx, rinputy, vinputy, rinputz, vinputz], options);
Notice the difference between propagate2BP and @propagate2BP . You need the @ version so your propagate_2BP22 would have to be @propagate_2BP22
댓글 수: 0
참고 항목
카테고리
Help Center 및 File Exchange에서 Counter and Timer Input and Output에 대해 자세히 알아보기
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!