paraSet.m 1.2 KB

12345678910111213141516171819202122232425262728293031323334353637383940414243444546474849
  1. function [Mission_Param, Environment_Param, Normalized_Param] = paraSet()
  2. % global Mission_Param Environment_Param Normalized_Param
  3. %% 环境参数
  4. Re = 6371004; % Earth radius (m)
  5. g0 = 9.801; % gravitational acceleration at Re (m/s^2)
  6. M0 = 1500; % initial mass (Kg)
  7. Isp = 2000; % engine specific impulse (s)
  8. m0 = 1; % normalized initial mass
  9. T = 1; % thrust magnitude (N)
  10. %% 归一化
  11. Lscale = Re;
  12. Vscale = sqrt(Re*g0);
  13. tscale = sqrt(Re/g0);
  14. mscale = M0;
  15. Fscale = M0*g0;
  16. %% 约束
  17. Mission_Param.P0 = 11623000/Re;
  18. Mission_Param.ex0 = 0.75;
  19. Mission_Param.ey0 = 0;
  20. Mission_Param.hx0 = 0.0612;
  21. Mission_Param.hy0 = 0;
  22. Mission_Param.L0 = pi;
  23. Mission_Param.Pf = 42165000/Re;
  24. Mission_Param.exf = 0;
  25. Mission_Param.eyf = 0;
  26. Mission_Param.hxf = 0;
  27. Mission_Param.hyf = 0;
  28. %% 环境
  29. Environment_Param.Re = Re;
  30. Environment_Param.g0 = g0;
  31. Environment_Param.M0 = M0;
  32. Environment_Param.Isp = Isp;
  33. Environment_Param.m0 = m0;
  34. Environment_Param.T = T/Fscale;
  35. Environment_Param.c = Isp*g0/sqrt(Re*g0); % Isp*g0/Vscale
  36. %% 归一化
  37. Normalized_Param.Lscale = Lscale;
  38. Normalized_Param.Vscale = Vscale;
  39. Normalized_Param.tscale = tscale;
  40. Normalized_Param.mscale = mscale;
  41. Normalized_Param.Fscale = Fscale;
  42. end