runreentry.m

来自「大气层内航天器飞行动力学matlab仿真模型」· M 代码 · 共 20 行

% Copyright Ashish Tewari (c) 2006
global dtr; dtr = pi/180;
global mu; mu = 3.986004e14;
global S; S = 4;
global c; c=0.5;
global m; m =  350;
global rm; rm = 6378140;
global omega; omega = 2*pi/(23*3600+56*60+4.0905);
global Gamma; Gamma=1.41; 
global f8; f8 = fopen('data8.mat', 'a'); 
long = -10*dtr;            
lat = -79.8489182889*dtr; 
rad= 6579.89967e3;  
vel= 7589.30433867; 
fpa= 0.54681217*dtr;
chi= 99.955734*dtr; 
options = odeset('RelTol', 1e-8);    
orbinit = [long; lat; rad; vel; fpa; chi]; 
[t, o] = ode45('reentry',[0, 1750],orbinit,options);
fclose('all');