matrix multiplication dimension issues
1 次查看(过去 30 天)
显示 更早的评论
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.
0 个评论
回答(3 个)
James Tursa
2018-12-6
编辑:James Tursa
2018-12-6
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 个评论
James Tursa
2018-12-6
编辑:James Tursa
2018-12-6
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.
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Frequently-used Algorithms 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!