Issue with running Multiple gravity assists code

4 次查看（过去 30 天）

Suraj Parasuram 2018-12-25

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/437209-issue-with-running-multiple-gravity-assists-code

评论： Suraj Parasuram 2018-12-26

mga.m

For an interplanetary mission design involving multiple gravity assists, I have been trying to use the code. It requires to create a structure array and everytime I create and pass it to the function it throws an error saying "Too many input arguments". I am unable to understand how to circumvent the error. I would appreaciate if somebody could show me how to pass the parameters. I have also attached the MATLAB code titled as mga.m. This is my attempt:

t(1)=2000;
t(2)=112;
t(3)=376;
MGAproblem.sequence=[3 2 3];
MGAproblem.DVLaunch=2.98;
MGAproblem.objective.type='Orbit insertion';
MGAproblem.objective.rp=1E+5;
MGAproblem.objective.e=0; 
yplot=1;
[J,rp,DVvec,Itime,FBalt,alph,Vrelin,Vrelout] = mga(t,MGAproblem,yplot)

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

请先登录，再进行评论。

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

回答（1 个）

Walter Roberson 2018-12-26

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/437209-issue-with-running-multiple-gravity-assists-code#answer_353903

At one time the code accepted yplot as aa third parameter but that was changed. Now set MGAproblem.yplot = 1

9 个评论
显示 7更早的评论隐藏 7更早的评论

Walter Roberson 2018-12-26

在 MATLAB Online 中打开

The code has

%Reading the values for the minimum allowed pericentre radius of each planet
%
rpmin=zeros(6,1);
rpmin(2)=problem.rpmin.venus;   %(Venus)
rpmin(3)=problem.rpmin.earth;   %(Earth)
rpmin(4)=problem.rpmin.mars;     %(Mars)
rpmin(5)=problem.rpmin.jupiter; %(Jupiter)
rpmin(6)=problem.rpmin.saturn;  %(Saturn)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%Reading the penalty factors for low pericenter fly-bys. These can be
%interpreted as km/sec penalties for each km of violation
%
rpfact=zeros(6,1);
rpfact(2)=problem.rpfact.venus;
rpfact(3)=problem.rpfact.earth;
rpfact(4)=problem.rpfact.mars;
rpfact(5)=problem.rpfact.jupiter;
rpfact(6)=problem.rpfact.saturn;