MATLAB: Ode45 – Does matlab keep track of derivatives which can be recalled

ode45

I've run into an issue with a system of coupled 2nd order nonlinear ODEs. I've converted the system to a system of 1st order ODEs like normal, but its not in the format that ode45 requires. ode45 requires

x_dot = f(t,x), where x is the state vector

My form is

x_dot = f(t,x,x_dot)

where the function also depends on the derivative. I searched all around to see if anyone else has had this issue, and only found one post on a forum. It just said to proceed normally with ode45 and use the old values (previous step) of x_dot to calculate the rhs.

This makes sense, but ode45 doesn't work that way (not that I know of). I only have access to the old state values and not their derivatives.

Any ideas?

Best Answer

Re-formulate your equations from

 x_dot = f(t,x,x_dot)

 g(t,x,x_dot) = f(t,x,x_dot) - x_dot = 0

and then use

 ode15i

to integrate the set of equations in g().

Related Solutions

MATLAB: Storing data from a function into and array to be recalled

a = arduino('com3', 'uno');
s = servo(a, 'D10','MaxPulseDuration', .0025, 'MinPulseDuration', .0005);
n = 60;
x=zeros(1,n);
for k = 1:n
z=readVoltage(a,'A5');
p=((1023 / 5) * z);
x(k)=p/1023;
writePosition(s,x(n));
end

# Now you can recall the array x, when in required (same workspace)

for k=1:n
    writePosition(s,x(k))
end

MATLAB: Precision/Recall (perfcurve)

perfcurve() is meant for computing performance curves for your model - such as PR curves or ROC curves. In the literature, you may find these curves used for a more rigorous examination of a model's performance than can be given by a single score.

I recommend looking up the doc linked below. The Algorithms > Pointwise Confidence Bounds section explains what goes on under-the-hood after you call the function with your particular arguments. It samples from your input data multiple times, with different thresholds, to come up with a range of P/R scores, which are used to plot the curve.

perfcurve() doc: https://www.mathworks.com/help/stats/perfcurve.html

To calculate a single value, you may extract counts of TP, FP, TN & FN from your labels and predictions; then manually compute -

Precision = TP / (TP + FP)

Recall = TP / (TP + FN)

Should be done in 2 lines of code.

Alternatively, you may use confmat() to produce a confusion matrix, to readily give you TP/FP/TN/FN counts. And then do the same computation.

%Using this method, here's an anonymous function to calculate precision.
% << for an existing Confision Matrix 'ConfusionMat' >>
precision = @(confusionMat) diag(confusionMat)./sum(confusionMat,2);
% And another for recall
recall = @(confusionMat) diag(confusionMat)./sum(confusionMat,1)';
% Computes a vector of P/R values - each value corresponds to selecting a different label as +ve.

Hope it Helps!

Best Answer

Related Solutions

MATLAB: Storing data from a function into and array to be recalled

MATLAB: Precision/Recall (perfcurve)

Related Question