Extracting elements from matrix to do summation

Now I have 3 matrixes

Oi is a 3x1 matrix

Dj is a 1x3 matrix

cij is a 3x3 matrix

And then I have 2 equations to solve:

Ai and Bj are unknowns and Ai should be a 3x1 matrix and Bj should be a 1x3 matrix

and it is actually only a 6 equations 6 unknown question.

But I found difficulties typing the equations in matlab.

Because I can't simply type something like

syms= i j

A(i)=sum(B(j)*D(j)*c(i,j),i=1,3)

Because things like A(i) is not gonna work

So how should I type the equations correctly?

I can't do it manually as I have equation like

where n equal to like 25

Answer 1

[Image Analyst]

What about

j = 1 : 25; % Define indexes over which to sum
% Now do the multiplication.  i is some number that's predefined.
A25 = B(j) .* D(j) .* c(i,j); % Array of 25 terms.
% Now sum all the terms.
A(i) = sum(A25);

Answer 2

[Roger Stafford]

Matlab can indeed save you some writing. Let these arrays already be defined:

 c = [c11,c12,...,c1n;
      c21,c22,...,c2n;
      .....
      cn1,cn2,...,cnn]

 O = [O1;O2;...;On]  (Column vector)

 D = [D1,D2,...,Dn]  (Row vector)

Now do the following in matlab:

 M = [eye(n),-bsxfun(@times,c,D);
     -bsxfun(@times,c,O).',eye(n)];
 [U,S,V] = svd(M);

Because your equations are homogeneous, you will either have as the only solution that all A's and B's are zero, or there will be infinitely many solutions. You can determine the nature of these solutions from the singular values in the diagonal of array S. If all singular values there are non-zero, then the only solution is for all A's and B's to be zero. If all but the rightmost singular value are non-zero and that is zero, then all multiples of V(:,2*n) are solutions with A's first and then the B's next. That is, A = t*V(1:n,2*n), B =t*V(n+1:2*n).' for any scalar t. If the rightmost two singular values are zero, then any linear combination of V(:,n) and V(:,n-1) are also solutions. If the rightmost k of the singular values are zero, then any linear combination of V(:,n), V(:,n-1),...,V(:,n-k+1) is a solution. That is the nature of homogeneous linear equations.