How can I vectorize this code that uses the fzero function?
2 views (last 30 days)
Show older comments
Martin Schwesinger
on 3 Aug 2015
Commented: Star Strider
on 3 Aug 2015
Hello,
I have the following code:
precision = 1000;
E = zeros(1,precision);
ecc = 0.3;
for k = 1:precision
M = 2*pi*k/precision;
E(k) = fzero(@(x) anomaly(x,M,ecc),M)
end
-----
function y = anomaly(x,M,ecc)
y = x - ecc*sin(x) - M;
end
Is there a way to vectorize this code and get rid of the k-loop? The fzero function doesn't accept a vector argument.
0 Comments
Accepted Answer
Star Strider
on 3 Aug 2015
That is likely as efficient as you can write code to calculate several zeros using fzero. You have to specify an initial estimate of the zero, and that apparently changes in your loop. I plotted ‘E’, and it is not linear, so your code is likely optimal.
2 Comments
Star Strider
on 3 Aug 2015
My pleasure!
I apologise for the delay — life intrudes.
You would have to estimate the zero-crossing in each interval of ‘M’ regardless. The fzero call is likely as efficient as you can get, because your function is nonlinear and a reasonably precise estimate of ‘E’ with it is likely as good as you can get.
You might even calculate fewer values of ‘E’ in your loop, and then use interp1 for reasonably precise interpolation of ‘E’ at specific values of ‘M’.
It all depends on the rest of your code, and what you want to do.
More Answers (0)
See Also
Categories
Find more on Optimization 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!