0001 function par = simulated_prepare_true_values(par)
0002
0003
0004
0005
0006
0007 range_g = 0:255;
0008 par.vg = var_g(range_g(:),par.stg);
0009 par.vh = var_h(range_g(:),par.sth);
0010
0011 par.vg = (par.vg+1/12*ones(256,1));
0012 par.vh = (par.vh+1/12*ones(256,1));
0013
0014 par.sigma_max = sqrt(max([par.vg(:);par.vh(:)]));
0015
0016
0017 par.A_true = [par.A_t; 0 0 1];
0018
0019 par.B_true = sqrtm(par.A_true);
0020 par.Bi_true = inv(par.B_true);
0021
0022
0023 par.a7_true = sqrt(par.R_true(1));
0024 par.a8_true = par.R_true(2)/(1+par.a7_true);
0025 par.S_true = [par.a7_true,par.a8_true];
0026
0027
0028 par.x_true = [reshape(par.A_true(1:2,:),6,1);par.R_true'];
0029
0030
0031 par.A_a = par.A_true;
0032 par.Ai_a = inv(par.A_a);
0033 par.R_a = par.R_true;
0034 par.Ri_a = 1/(par.R_a(1,1))*[1,-par.R_a(2)];
0035
0036 end
0037