matrix multiplication dimension issues
1 view (last 30 days)
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.
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)