Matlab_BackUp/lowerTest.m at master · mhaneveer/Matlab_BackUp · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
function datawrite = lowerTest(NORAD_ID)

verificationPath = 'D:\Documents\TU_Delft\MSc\Thesis\02_Program\Data\Test\';

TLE_path = [verificationPath num2str(NORAD_ID) 'test.txt'];

stepsize = 1         ; %minuten

mfe = [0.0 1.0 1.0];
typerun   = 'm';
typeinput = 'm';
mfe = [0.0 1.0 1.0];

TLEsteps =[0 0];
whichconst = 84;

[tsince,jd,ro,ro,ro,vo,vo,vo,...
    p,a,ecc,incl,node,argp,nu,m,arglat,truelon,lonper] ...
    = SGP4_Matlab(TLE_path,typeinput,typerun,whichconst,mfe,TLEsteps);

clear TLEsteps tsince

minJD = min(jd);
timeframe = abs(max(jd)-min(jd));
steps = floor(timeframe * (1440/stepsize));

for i = 1: steps
    JD_extra(i) = minJD + (stepsize/1440) * i - (stepsize/1440);
end

flipJD = fliplr(jd);
jdo = 1;

for ii = 1:steps
       jd1 = flipJD(jdo);
       jd2 = flipJD(jdo+1);
       if JD_extra(ii) > jd1 + (abs(jd1-jd2)/2)
                  jdo = jdo + 1;
                  jd1 = flipJD(jdo);

                  if jdo == length(jd)
                      jd2 = flipJD(jdo) + 9999;

                      for j = ii:steps
                            epoch(j) = jd1;
                            tsince(j)= JD_extra(j) - jd1;
                      end
                      break
                  end

                jd2 = flipJD(jdo+1);

                epoch(ii) = jd1;
                tsince(ii)= JD_extra(ii) - jd1;
       end

       if JD_extra(ii) <= jd1 + (abs(jd1 - jd2)/2)
                epoch(ii) = jd1;
                tsince(ii)= JD_extra(ii) - jd1;
       end
end

TLEsteps = [epoch'  tsince']

       [tsince,jd, ro1,ro2,ro3,vo1,vo2,vo3, ...
        p,a,ecc,incl,node,argp,nu,m,arglat,truelon,lonper]...
        = SGP4_Matlab(TLE_path,typeinput,typerun,whichconst,mfe,TLEsteps);

jd      = fliplr(jd);
ro1     = fliplr(ro1);
ro2     = fliplr(ro2);
ro3     = fliplr(ro3);
vo1     = fliplr(vo1);
vo2     = fliplr(vo2);
vo3     = fliplr(vo3);
incl    = fliplr(incl);

r = [ro1; ro2; ro3]';
v = [vo1; vo2; vo3]';

datawrite = [jd' r v,incl'];