07-Jan-2022 18:12:15

colored_noise_test():
  MATLAB/Octave version 9.8.0.1380330 (R2020a) Update 2.
  Test colored_noise().

r8vec_sftf_test():
  MATLAB version
  r8vec_sftf() computes the "slow" Fourier transform (forward)
  of a vector of real data.
  The original data can be recovered using R8VEC_SFTB.

  Fourier coefficients:

    0.640639
    0.078468   -0.152306
    0.058940    0.085121
   -0.099885   -0.087606
    0.260392    0.098680
    0.145092   -0.025968
   -0.184128    0.072947
   -0.084795    0.078096

  Compare data R and recovered data R2:

    0.814724    0.814724
    0.905792    0.905792
    0.126987    0.126987
    0.913376    0.913376
    0.632359    0.632359
    0.097540    0.097540
    0.278498    0.278498
    0.546882    0.546882
    0.957507    0.957507
    0.964889    0.964889
    0.157613    0.157613
    0.970593    0.970593
    0.957167    0.957167
    0.485376    0.485376
    0.800280    0.800280

r8vec_sftf_test():
  Normal end of execution.

colored_noise_test01():
  Generating 128 sample points.
  1/F^ALPHA noise has ALPHA = 0.000000
  Variance is 1.000000

  Noise sample:

         1:    -0.204966
         2:    -0.124144
         3:       1.4897
         4:      1.40903
         5:      1.41719
         6:     0.671497
         7:     -1.20749
         8:     0.717239
         9:      1.63024
        10:     0.488894
  Data written to file "alpha_0.00.txt".

colored_noise_test01():
  Generating 128 sample points.
  1/F^ALPHA noise has ALPHA = 0.250000
  Variance is 1.000000

  Noise sample:

         1:    -0.285651
         2:    -0.867073
         3:     -1.10321
         4:     -1.35148
         5:    -0.799478
         6:     -2.24093
         7:     0.546169
         8:     0.367127
         9:   -0.0977445
        10:     -0.12621
  Data written to file "alpha_0.25.txt".

colored_noise_test01():
  Generating 128 sample points.
  1/F^ALPHA noise has ALPHA = 0.500000
  Variance is 1.000000

  Noise sample:

         1:      1.08263
         2:      1.27674
         3:    -0.230227
         4:       0.3784
         5:    -0.760836
         6:     -1.41058
         7:      0.57306
         8:    0.0255533
         9:     -0.05486
        10:     0.895062
  Data written to file "alpha_0.50.txt".

colored_noise_test01():
  Generating 128 sample points.
  1/F^ALPHA noise has ALPHA = 0.750000
  Variance is 1.000000

  Noise sample:

         1:      0.69916
         2:     0.531834
         3:     0.775658
         4:     -1.08554
         5:     -1.27356
         6:    -0.959288
         7:    -0.403577
         8:     0.787803
         9:    -0.149723
        10:      1.20034
  Data written to file "alpha_0.75.txt".

colored_noise_test01():
  Generating 128 sample points.
  1/F^ALPHA noise has ALPHA = 1.000000
  Variance is 1.000000

  Noise sample:

         1:   -0.0692143
         2:     -0.54193
         3:   -0.0438072
         4:     0.151834
         5:     0.103912
         6:    -0.563445
         7:     -1.49979
         8:    -0.504825
         9:      -1.8203
        10:     -2.12378
  Data written to file "alpha_1.00.txt".

colored_noise_test01():
  Generating 128 sample points.
  1/F^ALPHA noise has ALPHA = 1.250000
  Variance is 1.000000

  Noise sample:

         1:     -1.62732
         2:    -0.850729
         3:    -0.346142
         4:    -0.630389
         5:     -1.68106
         6:     0.819045
         7:     -2.13337
         8:     -1.89816
         9:     -2.81249
        10:     -2.77607
  Data written to file "alpha_1.25.txt".

colored_noise_test01():
  Generating 128 sample points.
  1/F^ALPHA noise has ALPHA = 1.500000
  Variance is 1.000000

  Noise sample:

         1:     -1.26909
         2:    -0.453834
         3:      2.32973
         4:      2.38275
         5:       1.1868
         6:      2.01302
         7:     0.617643
         8:    -0.654749
         9:     0.435629
        10:    -0.394227
  Data written to file "alpha_1.50.txt".

colored_noise_test01():
  Generating 128 sample points.
  1/F^ALPHA noise has ALPHA = 1.750000
  Variance is 1.000000

  Noise sample:

         1:      1.26653
         2:     0.857043
         3:     0.614606
         4:      -1.5909
         5:     -2.10156
         6:     -3.28858
         7:     -3.39283
         8:     -2.67088
         9:     -1.09133
        10:     0.577339
  Data written to file "alpha_1.75.txt".

colored_noise_test01():
  Generating 128 sample points.
  1/F^ALPHA noise has ALPHA = 2.000000
  Variance is 1.000000

  Noise sample:

         1:      2.75256
         2:      2.89088
         3:     0.983809
         4:     0.618816
         5:    -0.229294
         6:    -0.994047
         7:     -2.12174
         8:     -2.04355
         9:    0.0630768
        10:    -0.652771
  Data written to file "alpha_2.00.txt".

colored_noise_test02():
  Generating 200 realizations
  Generating 128 sample points.
  1/F^ALPHA noise has ALPHA = 0.000000
  Variance is 1.000000

  Plot saved as "alpha_0.00_paths.png"

colored_noise_test02():
  Generating 200 realizations
  Generating 128 sample points.
  1/F^ALPHA noise has ALPHA = 0.500000
  Variance is 1.000000

  Plot saved as "alpha_0.50_paths.png"

colored_noise_test02():
  Generating 200 realizations
  Generating 128 sample points.
  1/F^ALPHA noise has ALPHA = 1.000000
  Variance is 1.000000

  Plot saved as "alpha_1.00_paths.png"

colored_noise_test02():
  Generating 200 realizations
  Generating 128 sample points.
  1/F^ALPHA noise has ALPHA = 1.500000
  Variance is 1.000000

  Plot saved as "alpha_1.50_paths.png"

colored_noise_test02():
  Generating 200 realizations
  Generating 128 sample points.
  1/F^ALPHA noise has ALPHA = 2.000000
  Variance is 1.000000

  Plot saved as "alpha_2.00_paths.png"

make_plots()
  We assume colored_noise_test() has just been run,
  creating files alpha_0.00.txt through alpha_2.00.txt.
  This program computes a common axis size for these files
  and plots each data set as a separate PNG file.

  Number of entries in sequence is 128

  Y-range of plots is -9.974505 to 11.963239

make_plots():
  Normal end of execution.

colored_noise_test():
  Normal end of execution.

07-Jan-2022 18:12:45