How to Assemble the matrix shown in the image on Matlab for any number of nodes?

Question

amar 2016-1-27

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/265485-how-to-assemble-the-matrix-shown-in-the-image-on-matlab-for-any-number-of-nodes

评论： Rena Berman 2017-1-24

采纳的回答： Walter Roberson

2 个评论
显示无隐藏无

Stephen23 2016-1-30

OP continued the same question here:

http://www.mathworks.com/matlabcentral/answers/265649-how-to-assemble-this-matrix-below-with-this-pattern-for-any-number-of-nodes-n

Rena Berman 2017-1-24

(Answers dev) Restored question.

请先登录，再进行评论。

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

Answer 1

Walter Roberson 2016-1-30

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/265485-how-to-assemble-the-matrix-shown-in-the-image-on-matlab-for-any-number-of-nodes#answer_207926

在 MATLAB Online 中打开

n = 4;    %for example
main_diagonal = [ [-4, repmat(-3,1,n-2), -4], repmat([-5, repmat(-4,1,n-2), -5], 1, 2), [-6, repmat(-5,1,n-2), -6] ];
N = length(main_diagonal);
m = diag(ones(1, N-4), -4) +
    diag(ones(1, N-1), -1) +
    diag(main_diagonal) +
    diag(ones(1, N-1), 1) +
    diag(ones(1, N-4), -4)

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

请先登录，再进行评论。

Answer 2

Guillaume 2016-1-27

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/265485-how-to-assemble-the-matrix-shown-in-the-image-on-matlab-for-any-number-of-nodes#answer_207657

在 MATLAB Online 中打开

There is not enough information to know how to construct the main diagonal, but you can just sum a bunch of diag:

m = diag(ones(1, 12), -4) + ...
    diag(ones(1, 15), -1) + ...
    diag([-4 -3 -3 -4 -5 -4 -4 -5 -5 -4 -4 -5 -6 -5 -5 -6], 0) + ... %no idea what the pattern is
    diag(ones(1, 15), 1) + ...
    diag(ones(1, 12), 4)

2 个评论
显示无隐藏无

amar 2016-1-27

I want to assemble this matrix for any number of nodes 'n'

Walter Roberson 2016-1-27

在 MATLAB Online 中打开

main_diagonal = [-4 -3 -3 -4 -5 -4 -4 -5 -5 -4 -4 -5 -6 -5 -5 -6];  %we do not know the rule for this
n = length(main_diagonal);
m = diag(ones(1, n-4), -4) +
    diag(ones(1, n-1), -1) +
    diag(main_diagonal) +
    diag(ones(1, n-1), 1) +
    diag(ones(1, n-4), -4)