How can I vectorize this code that uses the fzero function?
显示 更早的评论
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.
采纳的回答
Star Strider
2015-8-3
1 个投票
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 个评论
Martin Schwesinger
2015-8-3
Star Strider
2015-8-3
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.
更多回答（0 个）
类别
在 帮助中心 和 File Exchange 中查找有关 Interpolation 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
