0001 function [I_true,f_true,g_true,h_true,par] = generate_reference_image(par)
0002
0003
0004
0005
0006
0007 I_true = LSM_62_sym_generate_image(par,1);
0008
0009 par.init_rand_samples = init_rand_seed(par.init_rand_samples);
0010
0011
0012
0013 [f_true,g_true,h_true,par] = ...
0014 LSM_62_sym_generate_true_observing_images(I_true,par);
0015
0016 fs_true = conv2(f_true,[1,0,-1]/2,'valid');
0017 empirical_bandwidth = std(fs_true(:))/std(f_true(:))/(2*pi);
0018 disp(['estimated bandwidth : ', num2str(empirical_bandwidth)])
0019
0020
0021 end
0022