09-Jan-2023 21:34:10

colored_noise_test():
  MATLAB/Octave version 4.2.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.527266
   -0.077655   -0.014646
   -0.100172   -0.097967
   -0.071440   -0.055828
    0.076659   -0.004159
   -0.157591   -0.060167
    0.062417   -0.046256
    0.102857    0.132317

  Compare data R and recovered data R2:

    0.362342    0.362342
    0.099120    0.099120
    0.562651    0.562651
    0.443318    0.443318
    0.712404    0.712404
    0.831448    0.831448
    0.179620    0.179620
    0.838691    0.838691
    0.577579    0.577579
    0.417855    0.417855
    0.360943    0.360943
    0.844488    0.844488
    0.445239    0.445239
    0.758622    0.758622
    0.474674    0.474674

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.505813
         2:     -1.48758
         3:      1.37948
         4:      1.27688
         5:     -0.50486
         6:    0.0222772
         7:    -0.707908
         8:     -1.36703
         9:     0.413306
        10:    -0.640873
  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.824753
         2:     -1.49078
         3:     -3.09391
         4:     0.218635
         5:      1.34089
         6:    -0.408596
         7:     0.514974
         8:     0.363432
         9:      1.00345
        10:     0.534807
  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.56977
         2:     0.749213
         3:     0.541593
         4:      1.05776
         5:     0.872976
         6:    -0.198261
         7:      -1.3169
         8:    -0.876586
         9:    -0.940121
        10:     -0.55159
  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.843926
         2:     0.991899
         3:     0.123969
         4:     0.503984
         5:    0.0568735
         6:    -0.189097
         7:     0.419531
         8:      1.07463
         9:      1.53761
        10:      1.00215
  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:      1.27923
         2:     0.941118
         3:    -0.162349
         4:     -0.81666
         5:     -1.57748
         6:     0.980924
         7:      1.95567
         8:    -0.837928
         9:    -0.690124
        10:     0.479334
  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:     -0.28416
         2:     -1.75218
         3:     -1.92329
         4:     -2.13841
         5:     -1.42739
         6:      -1.0713
         7:     -1.62551
         8:    -0.832008
         9:     0.201192
        10:   -0.0326442
  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.05441
         2:      1.67089
         3:       1.4444
         4:    -0.508296
         5:      1.15097
         6:      2.55119
         7:      1.28035
         8:      3.37457
         9:      2.96199
        10:      2.43726
  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:     0.373159
         2:    -0.824277
         3:    -0.836017
         4:   -0.0255315
         5:     -1.44653
         6:     -2.14563
         7:     -1.84033
         8:    -0.417851
         9:    -0.858121
        10:    -0.522014
  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:    -0.910519
         2:    -0.208858
         3:    -0.578378
         4:     -1.96381
         5:     -1.60179
         6:     -2.20157
         7:     -2.81794
         8:     -1.69824
         9:     -1.61194
        10:     -2.96288
  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

  Graphics 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

  Graphics 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

  Graphics 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

  Graphics 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

  Graphics 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 -4.585078 to 16.329756

make_plots():
  Normal end of execution.

colored_noise_test():
  Normal end of execution.

09-Jan-2023 22:04:22