reducing number of operations in matrix multiplication

Question

Douglas Brenner 2018-10-10

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/423350-reducing-number-of-operations-in-matrix-multiplication

回答： Douglas Brenner 2018-10-11

I want to perform the follow operation: synch = spect1(1:num_pts ,2)*spect2(1:num_pts,2)'; The problem is that num_pts > 300,000. Would take too long. Fortunately, I don't need the entire output. I only need the diagonal elements plus elements from the say 10 adjacent diagonals on either side. Anybody know how to do this or want to take a guess? Thanks

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Douglas Brenner 2018-10-11

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/423350-reducing-number-of-operations-in-matrix-multiplication#answer_340886

编辑：Torsten 2018-10-11

在 MATLAB Online 中打开

    testm = [1,2,3,4,5;2,4,5,7,10]
    num_pts = 5;
    testm = transpose(testm)
    corr = testm(:,2)*testm(:,2)'
    diag(corr)
    main_diag = testm(:,2).*testm(:,2)
    for i=1:3
      diags1{i} = testm(i+1:num_pts ,2).*testm(1:num_pts-i-1,2);
      diags2{i} = testm(i+1:num_pts ,2).*testm(1:num_pts-i-1,2);
    end

failed at diags1{i} = testm(i+1:num_pts ,2).*testm(1:num_pts-i-1,2); Error using .* Matrix dimensions must agree. Error in test (line 8) diags1{i} = testm(i+1:num_pts ,2).*testm(1:num_pts-i-1,2);

Got to run.

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Torsten 2018-10-11

The number of elements of testm(i+1:num_pts,2) and testm(1:num_pts-i-1,2) are not the same ...

请先登录，再进行评论。

Answer 2

Douglas Brenner 2018-10-11

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/423350-reducing-number-of-operations-in-matrix-multiplication#answer_340902

Clearly but that was the suggested solution. Trying for i=1:2 diags1{i} = testm(i+1:num_pts-i-1 ,2).*testm(1:num_pts-i-1,2); diags2{i} = testm(i+1:num_pts-i-1 ,2).*testm(1:num_pts-i-1,2); end

gives the same error as does other attempts to fix it