why does my code take so long to run?

1 次查看（过去 30 天）
Omar B.
Omar B. 2020-10-7
评论： Rena Berman 2020-10-12
I am trying to compute v' * f(A) * v, where f is a given function, A matrix, and v is a unite vector.
My code takes along time doing run and did not get the result. Could please help me to fix this problem?
% want to compute v'*f(A)*v
syms f(x)
N=1000;
Ns=1:N;
R=1./(Ns);
A=sym(toeplitz(R));% the input matrix
v=ones(N,1);
v=v/norm(v);
f(x) = atan(sqrt(x)); % the function f(x)=arctan(sqrt(x))/x
B = funm(A,f); % the resulting matrix
y=A\v;
exact=v'*B*y
  5 个评论
显示 3更早的评论
Rik
Rik 2020-10-10
编辑：Rena Berman 2020-10-12
@Omar, Please explain why you want your question removed.
I also made a capture of the Google cache from before you removed your question here.

请先登录，再进行评论。

请先登录，再回答此问题。

回答（1 个）

Hiro Yoshino
Hiro Yoshino 2020-10-8
It looks that your problem is not that complicated. Why don't you write your code without using symbolic expressions? I bet it is faster.
  3 个评论
显示 1更早的评论
Walter Roberson
Walter Roberson 2020-10-8
There is a challenge that the numeric funm is different than the symbolic funm .
For the symbolic funm, you pass in a symbolic function or symbolic expression, and funm automatically takes derivatives in order to build a taylor series (or however it is represented). No control is offered over the number of terms used, and all the derivatives are likely to be exact derivatives.
If we look at even just the sqrt() part, as a matrix function, sqrtm(), that mathematically involves a closed form svd, sqrt of the diagonal, and reconstruct -- all as closed form expressions. For a 1000 x 1000 matrix, that is going to take a very long time.
For the numeric funm, you pass in a function handle to a function that must accept a vector and must also accept a derivative number to calculate and evaluate the input at. But you do not have diff() available to help because it is a numeric function handle you are working on, not a symbolic function. So you need to have a formula to calculate the N'th derivative.... Or you need to slip into the symbolic world, probably memoizing as you go.
Walter Roberson
Walter Roberson 2020-10-8
Even for N=10 the symbolic version takes much much much too long.

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Linear Algebra 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by