Home    Search    Reference Manuals    Return    Track Group Documentation


#!/usr/bin/octave -q

#   PSD          FREQ
#  [rad²/(1/m)]  [1/m]

 PSD_FI= [
  2.63768E-05  0.00000E+00
  2.63768E-05  1.00000E-02
  2.63569E-05  1.00334E-02
  2.15010E-05  1.33779E-02
  1.84403E-05  1.67224E-02
  1.55614E-05  2.00669E-02
  1.32894E-05  2.34114E-02
  1.02048E-05  2.67559E-02
  8.73987E-06  3.01003E-02
  9.50913E-06  3.34448E-02
  9.14280E-06  3.67893E-02
  8.31665E-06  4.01338E-02
  6.64736E-06  4.34783E-02
  5.99966E-06  4.68227E-02
  6.56259E-06  5.01672E-02
  6.06927E-06  5.35117E-02
  5.20516E-06  5.68562E-02
  4.93624E-06  6.02007E-02
  4.71950E-06  6.35452E-02
  4.32023E-06  6.68896E-02
  4.30656E-06  7.02341E-02
  4.62060E-06  7.35786E-02
  4.19490E-06  7.69231E-02
  3.69699E-06  8.02676E-02
  3.75002E-06  8.36120E-02
  3.66234E-06  8.69565E-02
  3.36124E-06  9.03010E-02
  3.23160E-06  9.36455E-02
  3.38242E-06  9.69900E-02
  3.21842E-06  1.00334E-01
  2.80510E-06  1.03679E-01
  2.42021E-06  1.07023E-01
  2.19268E-06  1.10368E-01
  1.86382E-06  1.13712E-01
  1.40720E-06  1.17057E-01
  1.25010E-06  1.20401E-01
  1.37530E-06  1.23746E-01
  1.50606E-06  1.27090E-01
  1.30958E-06  1.30435E-01
  9.25072E-07  1.33779E-01
  7.64852E-07  1.37124E-01
  7.94406E-07  1.40468E-01
  7.23987E-07  1.43813E-01
  6.36554E-07  1.47157E-01
  6.33253E-07  1.50502E-01
  5.76703E-07  1.53846E-01
  4.66796E-07  1.57191E-01
  3.47556E-07  1.60535E-01
  3.61715E-07  1.63880E-01
  4.45766E-07  1.67224E-01
  4.14932E-07  1.70569E-01
  3.40368E-07  1.73913E-01
  3.21766E-07  1.77258E-01
  3.47054E-07  1.80602E-01
  3.06490E-07  1.83946E-01
  2.69969E-07  1.87291E-01
  2.70445E-07  1.90635E-01
  2.21689E-07  1.93980E-01
  2.21024E-07  1.97324E-01
  2.32538E-07  2.00669E-01
  1.90394E-07  2.04013E-01
  1.60738E-07  2.07358E-01
  1.22614E-07  2.10702E-01
  1.30019E-07  2.14047E-01
  1.59650E-07  2.17391E-01
  1.36417E-07  2.20736E-01
  1.19795E-07  2.24080E-01
  1.20126E-07  2.27425E-01
  1.11714E-07  2.30769E-01
  1.05495E-07  2.34114E-01
  1.01202E-07  2.37458E-01
  8.91224E-08  2.40803E-01
  8.77216E-08  2.44147E-01
  9.08385E-08  2.47492E-01
  7.93177E-08  2.50836E-01
  7.54270E-08  2.54181E-01
  7.70132E-08  2.57525E-01
  6.66692E-08  2.60870E-01
  6.31542E-08  2.64214E-01
  6.14598E-08  2.67559E-01
  5.49481E-08  2.70903E-01
  5.39017E-08  2.74247E-01
  5.28510E-08  2.77592E-01
  4.98851E-08  2.80936E-01
  4.89032E-08  2.84281E-01
  4.54240E-08  2.87625E-01
  3.97151E-08  2.90970E-01
  4.14080E-08  2.94314E-01
  3.90612E-08  2.97659E-01
  3.33230E-08  3.01003E-01
  3.65546E-08  3.04348E-01
  3.85405E-08  3.07692E-01
  3.32903E-08  3.11037E-01
  2.93885E-08  3.14381E-01
  3.28373E-08  3.17726E-01
  3.80065E-08  3.21070E-01
  3.31444E-08  3.24415E-01
  2.43206E-08  3.27759E-01
  2.40125E-08  3.31104E-01
  2.65691E-08  3.34448E-01
  2.49603E-08  3.37793E-01
  1.99400E-08  3.41137E-01
  1.98136E-08  3.44482E-01
  2.33135E-08  3.47826E-01
  2.16325E-08  3.51171E-01
  1.84490E-08  3.54515E-01
  1.92630E-08  3.57860E-01
  1.97011E-08  3.61204E-01
  1.85709E-08  3.64549E-01
  1.85309E-08  3.67893E-01
  1.81962E-08  3.71237E-01
  1.62436E-08  3.74582E-01
  1.52348E-08  3.77926E-01
  1.51315E-08  3.81271E-01
  1.44317E-08  3.84615E-01
  1.32916E-08  3.87960E-01
  1.13762E-08  3.91304E-01
  1.29366E-08  3.94649E-01
  1.63088E-08  3.97993E-01
  1.53394E-08  4.01338E-01
  1.20677E-08  4.04682E-01
  9.75464E-09  4.08027E-01
  9.75149E-09  4.11371E-01
  1.12084E-08  4.14716E-01
  1.13665E-08  4.18060E-01
  1.04555E-08  4.21405E-01
  1.02855E-08  4.24749E-01
  9.64717E-09  4.28094E-01
  8.88024E-09  4.31438E-01
  9.11645E-09  4.34783E-01
  9.25824E-09  4.38127E-01
  8.90130E-09  4.41472E-01
  8.58492E-09  4.44816E-01
  8.36368E-09  4.48161E-01
  8.02500E-09  4.51505E-01
  8.01601E-09  4.54850E-01
  8.07724E-09  4.58194E-01
  7.65946E-09  4.61538E-01
  7.16543E-09  4.64883E-01
  6.92779E-09  4.68227E-01
  6.85231E-09  4.71572E-01
  6.60071E-09  4.74916E-01
  6.32462E-09  4.78261E-01
  6.38942E-09  4.81605E-01
  6.61269E-09  4.84950E-01
  6.56185E-09  4.88294E-01
  6.39560E-09  4.91639E-01
  6.56899E-09  4.94983E-01
  6.58089E-09  4.98328E-01
  5.98404E-09  5.01672E-01
  5.22680E-09  5.05017E-01
  4.69636E-09  5.08361E-01
  4.60851E-09  5.11706E-01
  4.75687E-09  5.15050E-01
  4.75557E-09  5.18395E-01
  4.62919E-09  5.21739E-01
  4.48060E-09  5.25084E-01
  4.37359E-09  5.28428E-01
  4.33536E-09  5.31773E-01
  4.27571E-09  5.35117E-01
  4.15379E-09  5.38462E-01
  4.08087E-09  5.41806E-01
  4.06223E-09  5.45151E-01
  4.01843E-09  5.48495E-01
  3.88362E-09  5.51839E-01
  3.72652E-09  5.55184E-01
  3.66953E-09  5.58528E-01
  3.62516E-09  5.61873E-01
  3.62766E-09  5.65217E-01
  3.66158E-09  5.68562E-01
  3.59737E-09  5.71906E-01
  3.51440E-09  5.75251E-01
  3.41583E-09  5.78595E-01
  3.30261E-09  5.81940E-01
  3.28244E-09  5.85284E-01
  3.29861E-09  5.88629E-01
  3.26429E-09  5.91973E-01
  3.21172E-09  5.95318E-01
  3.15191E-09  5.98662E-01
  3.06323E-09  6.02007E-01
  2.97888E-09  6.05351E-01
  2.94136E-09  6.08696E-01
  2.88387E-09  6.12040E-01
  2.77993E-09  6.15385E-01
  2.71635E-09  6.18729E-01
  2.70600E-09  6.22074E-01
  2.70075E-09  6.25418E-01
  2.65092E-09  6.28763E-01
  2.54508E-09  6.32107E-01
  2.44665E-09  6.35451E-01
  2.36649E-09  6.38796E-01
  2.25599E-09  6.42140E-01
  2.13231E-09  6.45485E-01
  2.07117E-09  6.48829E-01
  2.08841E-09  6.52174E-01
  2.10865E-09  6.55518E-01
  2.08908E-09  6.58863E-01
  2.06033E-09  6.62207E-01
  2.02474E-09  6.65552E-01
  1.98534E-09  6.68896E-01
  1.99110E-09  6.72241E-01
  2.04997E-09  6.75585E-01
  2.07206E-09  6.78930E-01
  2.01616E-09  6.82274E-01
  1.96124E-09  6.85619E-01
  1.91528E-09  6.88963E-01
  1.85756E-09  6.92308E-01
  1.79129E-09  6.95652E-01
  1.71952E-09  6.98997E-01
  1.67284E-09  7.02341E-01
  1.64718E-09  7.05686E-01
  1.63305E-09  7.09030E-01
  1.63428E-09  7.12375E-01
  1.65118E-09  7.15719E-01
  1.65616E-09  7.19064E-01
  1.62274E-09  7.22408E-01
  1.59527E-09  7.25753E-01
  1.58002E-09  7.29097E-01
  1.55736E-09  7.32441E-01
  1.54029E-09  7.35786E-01
  1.48557E-09  7.39130E-01
  1.39573E-09  7.42475E-01
  1.33000E-09  7.45819E-01
  1.28459E-09  7.49164E-01
  1.24143E-09  7.52508E-01
  1.20370E-09  7.55853E-01
  1.16208E-09  7.59197E-01
  1.12565E-09  7.62542E-01
  1.11139E-09  7.65886E-01
  1.08640E-09  7.69231E-01
  1.04799E-09  7.72575E-01
  1.02762E-09  7.75920E-01
  1.01450E-09  7.79264E-01
  9.97493E-10  7.82609E-01
  9.88004E-10  7.85953E-01
  9.78198E-10  7.89298E-01
  9.58170E-10  7.92642E-01
  9.37215E-10  7.95987E-01
  9.08387E-10  7.99331E-01
  8.75255E-10  8.02676E-01
  8.53574E-10  8.06020E-01
  8.41054E-10  8.09365E-01
  8.30973E-10  8.12709E-01
  8.30307E-10  8.16054E-01
  8.35337E-10  8.19398E-01
  8.31325E-10  8.22742E-01
  8.30269E-10  8.26087E-01
  8.28890E-10  8.29431E-01
  8.14120E-10  8.32776E-01
  7.90122E-10  8.36120E-01
  7.69397E-10  8.39465E-01
  7.64227E-10  8.42809E-01
  7.60543E-10  8.46154E-01
  7.46117E-10  8.49498E-01
  7.34317E-10  8.52843E-01
  7.36061E-10  8.56187E-01
  7.33472E-10  8.59532E-01
  7.14946E-10  8.62876E-01
  6.99156E-10  8.66221E-01
  6.90911E-10  8.69565E-01
  6.81893E-10  8.72910E-01
  6.71204E-10  8.76254E-01
  6.66192E-10  8.79599E-01
  6.61005E-10  8.82943E-01
  6.49232E-10  8.86288E-01
  6.47116E-10  8.89632E-01
  6.44607E-10  8.92977E-01
  6.21641E-10  8.96321E-01
  5.94418E-10  8.99666E-01
  5.73764E-10  9.03010E-01
  5.65065E-10  9.06354E-01
  5.77914E-10  9.09699E-01
  5.89499E-10  9.13043E-01
  5.79627E-10  9.16388E-01
  5.62481E-10  9.19732E-01
  5.50532E-10  9.23077E-01
  5.37454E-10  9.26421E-01
  5.22336E-10  9.29766E-01
  5.14899E-10  9.33110E-01
  5.11498E-10  9.36455E-01
  4.96072E-10  9.39799E-01
  4.77770E-10  9.43144E-01
  4.73048E-10  9.46488E-01
  4.70505E-10  9.49833E-01
  4.58641E-10  9.53177E-01
  4.43324E-10  9.56522E-01
  4.41722E-10  9.59866E-01
  4.50646E-10  9.63211E-01
  4.53580E-10  9.66555E-01
  4.56658E-10  9.69900E-01
  4.52104E-10  9.73244E-01
  4.34013E-10  9.76589E-01
  4.22303E-10  9.79933E-01
  4.12194E-10  9.83278E-01
  3.95471E-10  9.86622E-01
  3.84526E-10  9.89967E-01
  3.84333E-10  9.93311E-01
  3.85786E-10  9.96656E-01
  3.84519E-10  1.00000E+00
  3.84519E-10  1.00000E+02
  ];

 fid= fopen ("PSD_FI.psd", "w");
 if (fid == -1) error ("Can not open file PSD_FI.psd"); end

 fprintf(fid,"\n PSD_FI=\n");

  for i= 1:length(PSD_FI)
   if (PSD_FI(i,2) < 1./25.)
    PSD_FI(i,1)= 10.*PSD_FI(i,1);
   else
    PSD_FI(i,1)=     PSD_FI(i,1);
   end
   fprintf(fid," %e %e\n", PSD_FI(i,1), PSD_FI(i,2));
  end

 fclose (fid);

 fprintf("File PSD_FI.psd has been written\n");