Reputation: 57
Im tring to substract the diagonal values with eigval
and store the new value in the matrix Diagonal
:
CovarianceMatrix=[8 -3 1;2 1 0;3 4 5];
Col=3;
Row=3;
store=1;
syms eigval;
for loop1= Col:-1:1
Rw=1;
syms eigval;
for loop2= 1:Row
if Rw==loop1
Diagonal= (CovarianceMatrix(Rw,loop1)-eigval);
Fix_Diagonal_2(loop2,store)=sym(Diagonal);
else
Diagonal= CovarianceMatrix(Rw,loop1);
Fix_Diagonal_2(loop2,store)=Diagonal;
end
Rw=Rw+1;
loop1=loop1-1;
if loop1==0
loop1=3;
end
end
store=store+1;
end
But because I'm using a symbolic variable it gives an error:
The following error occurred converting from sym to double:
Error using mupadmex
Error in MuPAD command: DOUBLE cannot convert the input expression into a double
array.
If the input expression contains a symbolic variable, use the VPA function
instead.
How can I solve this? I would like to copy the new substracted value into the diagonal matrix.
Upvotes: 0
Views: 12641
Reputation: 18484
Here's a simple bit of code that exhibits your same error so maybe it will help clarify the issue:
syms x; % Create symbolic variable
a1 = rand(2); % Floating point array 1
a2 = rand(2); % Floating point array 2
d = a1(1)-x; % This is now a symbolic expression
a2(1) = d; % Error: you can't store a symbolic expression in a double array
which in R2013b (and R2015b) returns
The following error occurred converting from sym to double:
Error using mupadmex
Error in MuPAD command: DOUBLE cannot convert the input expression into a double
array.
If the input expression contains a symbolic variable, use the VPA function
instead.
vpa
can't be used here because x
is a symbolic variable that hasn't been defined rather than a symbolic value (vpa(d)
has an inconsequential effect).
For your code, the error likely occurs on this line:
Fix_Diagonal_2(loop2,store)=Diagonal;
You can't use vpa
because eigval
is a symbolic variable with no value. You can possibly solve your problem by casting Fix_Diagonal_2
to sym
:
Fix_Diagonal_2 = sym(Fix_Diagonal_2);
You'll probably want to do that outside of the for
loops. I also don't see why you're redefining eigval
on every iteration of the outer loop.
Upvotes: 4