Home > src > demos > reentry_demo > reentry_cond.m

reentry_cond

PURPOSE ^

Compute condition numbers of transition matrices

SYNOPSIS ^

This is a script file.

DESCRIPTION ^

 Compute condition numbers of transition matrices
 for re-entry simulation with no noise.

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

SOURCE CODE ^

0001   %
0002   % Compute condition numbers of transition matrices
0003   % for re-entry simulation with no noise.
0004   %
0005   
0006   % Copyright (C) 2005-2006 Simo S�kk�  %
0007   % This software is distributed under the GNU General Public
0008   % Licence (version 2 or later); please refer to the file
0009   % Licence.txt, included with the software, for details.
0010   fprintf('[Generating data...]\n');
0011   
0012   nsteps = round(200/dt);
0013   x = true_m0;
0014   X = zeros(size(x,1),nsteps);
0015   Y = zeros(2,nsteps);
0016   T = zeros(1,nsteps);
0017 
0018   C = zeros(1,nsteps);
0019   t = 0;
0020   for k=1:nsteps
0021     ddt = dt / sim_iter;
0022     for i=1:sim_iter
0023       A = reentry_df_dx(x,{ddt,b0,H0,Gm0,R0});
0024       x = reentry_f(x,{ddt,b0,H0,Gm0,R0});
0025     end
0026     c = cond(A)
0027     
0028     t = t + dt;
0029     C(:,k) = c;
0030     T(k) = t;
0031   end

Generated on Fri 12-Aug-2011 15:15:16 by m2html © 2005