How to decompose a 4th order tensor

10 views (last 30 days)
Derek Cooper
Derek Cooper on 3 Feb 2021
Edited: Derek Cooper on 4 Feb 2021
Hello
Is there any matlab function to decompose a 4th order tensor?

Answers (2)

Christine Tobler
Christine Tobler on 3 Feb 2021
The value of K1111 is overwritten with the next value in every loop iteration - is this your intention, or did you maybe mean to add them up? That is, for example,
K1111 = K1111 + inv(J) * F(a,A) * F(b,B) * K(A,B,K,D) * F(K,K) * F(d,D);
I'm also getting the impression you're using the same variable name K for both the iterator in the for-loop and the tensor K you're operating on. The formula you wrote up was using F(c, C) instead of F(K, K) in the code.
  2 Comments
Derek Cooper
Derek Cooper on 3 Feb 2021
Thanks for you answer.
The sought is components of matrix M and yes i want to add them upp. and K(A,B,C,D) is taken with respect of positions of K in an old M matrix. I hope you answerstand. The positions are the same as in M.
Christine Tobler
Christine Tobler on 3 Feb 2021
But your code is using K(A, B, K, D), not K(A, B, C, D). Can you confirm if this is correct? I don't think the code matches up with the definition you've written above, but I'm not sure which of the two fits what you want to compute.

Sign in to comment.


Derek Cooper
Derek Cooper on 3 Feb 2021
Edited: Derek Cooper on 3 Feb 2021
Sorry you are right, it should be K(A,B,C,D), and is components of D. so maybe its appropriate to write it as D(A,B,C,D).
  2 Comments
Christine Tobler
Christine Tobler on 3 Feb 2021
Is D a scalar index from 1 to 3, or is it a 4-dimensional array? It seems to be used for both, maybe there should be a different variable name for one of them, too.
Derek Cooper
Derek Cooper on 3 Feb 2021
D is a 3x3 matrix, not a 4 dimentional array. The reason why ABCD is used is because if K1111 is sought, then the value of D(1,1) should be taken since K1111 belongs to the same position as D(1,1)

Sign in to comment.

Categories

Find more on Creating and Concatenating Matrices in Help Center and File Exchange

Tags

Community Treasure Hunt

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

Start Hunting!