Sum of matrix times scalars

1 view (last 30 days)
Mohammed Kagalwala
Mohammed Kagalwala on 21 Oct 2019
Edited: Mohammed Kagalwala on 21 Oct 2019
Hi,
I have a series of 3x3 matrices T1, T2, T3 ... TN. In addition, I have a series of scalars a1, a2, ... aN. I wish to perform the following sum WITHOUT the use of a for-loop. T1*a1 + T2*a2 + .... TN*aN. I have solved this using a for loop currently, but wish to take advantage of Matlabs superior data structures i.e. 3D arrays, cells, etc.
Thank you for your help !
  3 Comments
Mohammed Kagalwala
Mohammed Kagalwala on 21 Oct 2019
They are in different individual variable names however no reason they can't be grouped in a cell or 3D matrix

Sign in to comment.

Accepted Answer

Matt J
Matt J on 21 Oct 2019
Edited: Matt J on 21 Oct 2019
Hold your 3x3 matrices in a 3x3xN array caled T and your scalars in an Nx1 vector called a and do,
result=reshape(T,9,[])*a(:);
result=reshape(result,3,3);
  2 Comments
Guillaume
Guillaume on 21 Oct 2019
Note that if a is a Nx1 vector, then a(:) and a are the same.
Another option is:
result = sum(T .* reshape(a, 1, 1, []), 3);
Mohammed Kagalwala
Mohammed Kagalwala on 21 Oct 2019
Edited: Mohammed Kagalwala on 21 Oct 2019
Thank you for your help ! The reshape solutions both work, however on average Matt's solution is faster ( I tested with N = 10000). My other method was coming up with an S1 and S2 such that T*S1*q*S2 = to the 3x3 sum I want. Here T is 3x3N, and q is Nx1. I'm going to try and extend the reshape idea to the following double sum, , where = 3x3 matrix.

Sign in to comment.

More Answers (0)

Categories

Find more on Linear Algebra in Help Center and File Exchange

Community Treasure Hunt

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

Start Hunting!