matrix multiplication dimension issues

1 view (last 30 days)
Douglas Brenner
Douglas Brenner on 6 Dec 2018
Answered: Douglas Brenner on 7 Dec 2018
I have an article which I'm trying to follow. It calculates two matricies that I called synch and async.
I read two files into two matricies; find the starting locations I want; and select n points. The I transpose the matricies (spect1 and spect2) each ending up being 200 x 2.
I then do matrix multiplication (synch = spect1(2:n, : )*spect2(2:n, :)'/(n-2);) on line 43 and get what I am pretty sure is the right output. However, this is not the formula in the article. The formula in the article is synch = spect1(:, 2:n)*spect2(:,2:n)'/(n-2); This gives me an Index exceeds matrix dimensions error. Not transposing the matricies after reading them in doesn't help. I get the wrong answer.
The calculation for async is on line 63. The article says async = spect1(:, 2:n) * N *spect2(:, 2:n)' /(n-2) but both ordering of the dimensions gives the Index exceeds matrix dimensions error.
Any help would be appreicated.

Answers (3)

James Tursa
James Tursa on 6 Dec 2018
Edited: James Tursa on 6 Dec 2018
If the formula in the article has the calculation as async = spect1(:, 2:n) * N * spect2(:, 2:n)'/(n-2) , and N is an (n-1)x(n-1) size matrix, it seems pretty clear that your spect1 and spect2 should be 2x200 instead of 200x2. So don't transpose the data in the beginning. I.e., let spect1 and spect2 be 2x200 and then use the formulas as-is from the article. Then you should be able to do this:
sync = spect1(:, 2:n) * spect2(:, 2:n)'/(n-2);
async = spect1(:, 2:n) * N * spect2(:, 2:n)'/(n-2);
Is the data real or complex? Are you sure the data in spect1 and spect2 is in the correct order? (maybe your spect1 data should be in spect2 and vice-versa)
  2 Comments
James Tursa
James Tursa on 6 Dec 2018
Edited: James Tursa on 6 Dec 2018
In your code, the N matrix is (n-1)x(n-1). How can you possibly use the async formula to end up with an (n-1)x(n-1) matrix? Why are you expecting this size as a result? spect1 and spect2 would have to be 199x200 size in order for this to agree with your expectations.

Sign in to comment.


Douglas Brenner
Douglas Brenner on 6 Dec 2018
The basic science is to take two spectra and perform a 2d cross correlation resulting in a real and a complex matrix. They should be nxn or maybe in the case of the imaginary matrix (sync) 199 x 199. So I know what I should get from what I have just haven't been able together there.

Douglas Brenner
Douglas Brenner on 7 Dec 2018
Found the answer.
n is not the number of points in a spectrum but the number of spectra.
Thanks for trying

Community Treasure Hunt

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

Start Hunting!

Translated by