On 15.11.12 22:24, Arda wrote:
> I am using Symbolic Toolbox, and facing serious speed issues.
> Before describing the problem I am facing  here is my intention:
>
> I want to create a matrix A (size n by n) whose i , j entries are defined as follows:
>
> A(i,j) = (i out of j)* p^i * p ^(ji) * (1)^(ji)
I assume you mean (1p)^(ji)? Otherwise, the p factors combine to p^j.
> This A matrix has a symbolic variable p, but as the size grows the computation(high orders of p, and combination(j,i) ) becomes too slow > over 2 hours for a matrix of size 400*400.
>
> And I would like to get a vector x = A*b
> where b is an 1*n vector whose values are to be read from data. I defined b as a symbolic variable so that for different data files, I would simply plug in and compute A*b in terms of p.
That sounds numerically unstable. Have you checked that writing the
formula this way is actually a good idea on your data?
> Here is the code I am using:
>
> %%%%%%%%%%%%%%%%
> %%%fill matrix
>
> n =400;
> %n = length(rhs);
> A = sym(zeros(n,n));
> sym p;
That should be
p = sym('p');
since the sym command does not assign to a variable (and using syms
inside procedures may confuse the linter and possibly the optimizer).
> In a thread I'have read that using MatlabFunction is a better idea, but still the following operation : Z = matlabFunction(x_rhs,'vars',[p, sym('B',[1 n])])
> takes hours to perfom.
It may be worth looking at one of the entries you created. If it
contains very complicated constants, you may be better of switching to
vpa arithmetic early (i.e., when putting the constants into the
expressions in the first place), either by explicitly calling the vpa
function or by converting your doubles with one of the additional
arguments to sym, such as 'd'. If, on the other hand, your formulas look
more complicated than they need to be, add a call to simplify, probably
best after computing A*b'. Although with your particular formulas, I
don't expect that would do much.
HTH,
Christopher
