2 %Brian
Staber (brian.staber@gmail.com)
9 modelParameters.mu = 1e3*[1.7710, 0.0658, 0.0680, 1.4152, 0.0718];
10 modelParameters.beta = [25.4185, 0.0432];
12 optimParameters.station = 44;
13 optimParameters.np = 28;
15 load(
'/home/s/staber/Trilinos_results/nrl/data/eij.mat');
22 Exx = [exx{ID(k)}{1}, exx{ID(k)}{2}];
23 Eyy = [eyy{ID(k)}{1}, eyy{ID(k)}{2}];
24 Exy = [exy{ID(k)}{1}, exy{ID(k)}{2}];
29 fd = fopen(strcat(
'/home/s/staber/Trilinos_results/nrl/random_generator_for_pca_likelihood/station',num2str(optimParameters.station),
'/output.txt'),
'w');
30 [ln,lt] = meshgrid(1
e-2*(5:5:25)*sqrt(50^2+25^2), 1
e-2*(4:4:20)*sqrt(50^2+25^2));
35 modelParameters.lc = [ln(k), lt(k)];
36 modelParameters.delta = repmat(0.1,1,4);
38 optimParameters.tol = 1
e-6;
39 optimParameters.nmc = 100;
41 output{
k} = costFunction(modelParameters,optimParameters,
Yexpi);
47 save(strcat(
'result_station',num2str(optimParameters.station),
'_11_07_2018.mat'),
'output',
'-v7.3');
load('/home/s/staber/Trilinos_results/nrl/data/eij.mat')
fprintf(fp, '< ParameterList >\n\n')
Brian Staber(brian.staber @gmail.com) % clc clearvars close all modelParameters.mu