Multi dimensional vector operation

1 次查看（过去 30 天）
MarshallSc
MarshallSc 2021-11-14
评论： Matt J 2021-11-15
Can someone please help me with figuring out the code for the operation below:
I have two sets of matrices, for example:
r_x=rand(3,3);r_y=rand(3,3);r_z=rand(3,3);
h_x=rand(3,3);h_y=rand(3,3);h_z=rand(3,3);
That is shown as, for x for example:
What I want to do is to take each component of the matrices in the 3 dimensions and form a column vector such that, for example for component (1,1):
r_11=[r_x(1,1);r_y(1,1);r_z(1,1)]; h_11=[h_x(1,1);h_y(1,1);h_z(1,1)]
And solve three sets of simultaneous equations:
Which essentially do this calculation:
I want to solve for the L matrix (r and h are known). I'd appreciate it if someone can please help me. Thank you.
  2 个评论
Matt J
Matt J 2021-11-14
编辑：Matt J 2021-11-14
Does the superscript 2 in L^2 mean the elements are squared? Because the equations are linear in L^2, it is unlikely that the solutions for L^2 will all be non-negative, and therefore unlikely that their square roots will be real-valued, if that's what you were hoping for.
MarshallSc
MarshallSc 2021-11-15
Yes, they are quadratic. My actual values are in complex conjugate format so it is not a problem if they're not real values.

请先登录，再进行评论。

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

采纳的回答

Matt J
Matt J 2021-11-14
编辑：Matt J 2021-11-14
L_squared=([hx(:),hy(:),hz(:)]\[rx(:),ry(:),rz(:)]).';
  2 个评论
MarshallSc
MarshallSc 2021-11-15
Thanks Matt for your answer. Shouldn't we identify each components of the matrix first (convert them into vectors containing x,y,z for each point) and then perform the operation on them, maybe in a for loop. The syntax you wrote appear to render a 3*3 matrix, while the actual size at the end of the calculation should be a 3*3 array each has a 3*1 vector. I converted them into 3*3 cell first and then want to do the subsequent operation. Something like this:
for i=1:3
for j=1:3
R{i,j}=[rx(i,j);ry(i,j);rz(i,j)];
H{i,j}=[hx(i,j);hy(i,j);hz(i,j)];
end
end
For each component of H or R in 3D, I need to have the 3*3 matrix for each point. For example:
for i=1:10
for j=1:10
L{i,j}=sqrt(R{i,j}/H{i,j});
end
end
But this one gives me 3*3 matrix with 6 elements as zero. I know the division of a 3*1 by 3*1 would be a 3*1, but based on the equation I noted in my original script, the Ls are obtained by simultaneous equations and would not be zero. Sorry if it's a little bit confusing.
Matt J
Matt J 2021-11-15
For each component of H or R in 3D, I need to have the 3*3 matrix for each point.
You can't derive a 3x3 matrix L from two 3x1 vectors. That gives you only 3 equations for 9 unknowns.

请先登录，再进行评论。

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Mathematics 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by