Home > atmlab > randomize > iaaft > surrogate_2d_2_horizontal.m

surrogate_2d_2_horizontal

PURPOSE ^

surrogate_2d_2d_horizontal

SYNOPSIS ^

This is a script file.

DESCRIPTION ^

 surrogate_2d_2d_horizontal

 This is the main program of the Iterative Amplitude Adapted Fourier
 Transform  (IAAFT) method to make surrogate fields. This version makes 2D
 fields based on the statistical properties of 2D fields. The amplitude distribution
 is supposed to be valid for the entire field (no 'vertical' profile).

 The IAAFT method was developped by Schreiber and Schmitz (see e.g. Phys. 
 Rev Lett. 77, pp. 635-, 1996) for statistical non-linearity tests for time series.
 This method makes fields that have a specified amplitude distribution and
 power spectral coefficients. It works by iteratively adaptation the amplitude 
 distribution and the Fourier coefficients (the phases are not changed in this 
 step). Do not use this program without understanding the function
 iaaft_loop_2d_horizontal and tuning its variables to your needs.

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

DOWNLOAD ^

surrogate_2d_2_horizontal.m

SOURCE CODE ^

0001 % surrogate_2d_2d_horizontal
0002 %
0003 % This is the main program of the Iterative Amplitude Adapted Fourier
0004 % Transform  (IAAFT) method to make surrogate fields. This version makes 2D
0005 % fields based on the statistical properties of 2D fields. The amplitude distribution
0006 % is supposed to be valid for the entire field (no 'vertical' profile).
0007 %
0008 % The IAAFT method was developped by Schreiber and Schmitz (see e.g. Phys.
0009 % Rev Lett. 77, pp. 635-, 1996) for statistical non-linearity tests for time series.
0010 % This method makes fields that have a specified amplitude distribution and
0011 % power spectral coefficients. It works by iteratively adaptation the amplitude
0012 % distribution and the Fourier coefficients (the phases are not changed in this
0013 % step). Do not use this program without understanding the function
0014 % iaaft_loop_2d_horizontal and tuning its variables to your needs.
0015 
0016 % This Matlab version was written by Victor Venema,
0017 % Victor.Venema@uni-bonn.de, http:\\www.meteo.uni-bonn.de\victor, or
0018 % http:\\www.meteo.uni-bonn.de\victor\themes\surrogates\
0019 % for the generation of surrogate cloud fields.
0020 % First version: May 2003.
0021 % This version:  November 2003.
0022 
0023 % Copyright (C) 2003 Victor Venema
0024 % This program is free software; you can redistribute it and/or
0025 % modify it under the terms of the GNU General Public License
0026 % as published by the Free Software Foundation; version 2
0027 % of the License.
0028 %
0029 % This program is distributed in the hope that it will be useful,
0030 % but WITHOUT ANY WARRANTY; without even the implied warranty of
0031 % MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
0032 % GNU General Public License for more details.
0033 %
0034 % To get a copy of the GNU General Public License look on internet
0035 % or write to the Free Software Foundation, Inc., 59 Temple Place -
0036 % Suite 330, Boston, MA  02111-1307, USA.
0037 
0038 % Load data.
0039 [fourier_coeff_2d, sorted_values, x, y, template, meanValue, no_values_x, no_values_y] = load_2d_data_horizontal(1);
0040 
0041 % Main iterative loop for 2d-surrogates
0042 [surrogate, error_amplitude, error_spec] = iaaft_loop_2d_horizontal(fourier_coeff_2d, sorted_values);
0043 surrogate = surrogate + meanValue;
0044 
0045 % plot results
0046 plot_2d_surrogate(x, y, template,  'template')
0047 plot_2d_surrogate(x, y, surrogate, 'surrogate')

Generated on Mon 15-Sep-2014 13:31:28 by m2html © 2005