# How to Estimate the exponent a?

2 views (last 30 days)
Ege Tunç on 28 Dec 2018
Commented: Walter Roberson on 7 Jan 2019
P=e^-aT in matlab?
Jan on 28 Dec 2018
"P~e^-aT where x is the size of the shield" - x does not occur in this equation. Later we find: "where T is the size of the shield".
Walter Roberson on 7 Jan 2019

Jan on 28 Dec 2018
You have the equation P ~ exp(-aT). You determine the value of P by the simulation and have T as input parameter. Then it seems to be trivial to solve this equation to obtain a:
a = -log(P) / T
Instead of collecting the results in a diagram only, you need to store the final positions in an array. E.g.:
% x_t(1) = 0; % Remove this
% y_t(1) = 0; % Remove this
x_final = zeros(1, M);
y_final = zeros(1, M);
for m = 1:M
x_t = zeros(1, N);
x_t = zeros(1, N)
for n = 1:N-1 % Not until N, but N-1!
...
end
x_final(m) = x_t(N);
y_final(m) = y_t(N);
end
Now you can determine P by comparing the final positions with the goal.
Jan on 28 Dec 2018
Edited: Jan on 28 Dec 2018
While "A=++" is no valid Matlab syntax, I assume you need something like this:
P = sum(x_final > T) / M;
a = -log(P) / T
I do not understand, why you ask for "x(N+1)=x(N)+1". I showed you a way to store the final positions and the method to calculate the parameter a based on the results of the simulation. Keep the rest as you like.

### Categories

Find more on Schedule Model Components in Help Center and File Exchange

### Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!