How is this supposed to work? I don't see where anything is done with the intermediate parameter values. The function simply perturbs the parameter values until they walk outside the boundaries. The loop only exits when the parameter set is outside of the specified region, not within.
As an example:
dEP2 = 1;
dEP3 = 1;
dSIG2 = 0.01;
dSIG3 = 0.1;
dPo = 2;
model.Po = 190;
model.Th = 500-model.Po;
model.EP2 = 3;
model.EP3 = 23;
model.SIG2 = 0.015;
model.SIG3 = 0.3;
n = 0;
t = testparams(model);
while all(t)
%model.EP1 = model.EP1; % this does nothing
%model.SIG1 = model.SIG1; % this does nothing
model.Po = model.Po + round(dPo*(randn/2));
model.Th = (500-model.Po) ;
model.EP2 = model.EP2 + dEP2*(randn/2);
model.EP3 = model.EP3 + dEP3*(randn);
model.SIG2 = model.SIG2 + dSIG2*(randn/2);
model.SIG3 = model.SIG3 + dSIG3*(randn/2);
n = n+1;
t = testparams(model);
end
n % loop executed twice
t % two parameters are out of bounds (SIG2 and EP3)
model
function tvec = testparams(model)
tvec = [model.Po+model.Th==500 model.Po>188 model.Po<210 model.Th>290 model.Th<312 ...
model.EP2>2.8 model.EP2<4.5 model.EP3>22 model.EP3<26 ...
model.SIG2>0.01 model.SIG2<0.022 model.SIG3>0.2 model.SIG3<0.6];
end
The loop exits after a variable number of passes, not just one. The loop only exits when the set of parameters contains a value outside the specified limits.
My guess is that this loop shouldn't even be part of this function, but part of the calling scope instead, but that's just a wild guess.
