Cross Product of an Array

24 次查看（过去 30 天）
Hollis Williams
Hollis Williams 2019-4-22
评论： James Tursa 2019-4-22
I have an array of size 3 x 100 and I basically want to create a new array where I take each column vector in the array and compute the cross product with the same vector each time to create another array of size 3 x 100, so in this case I take every vector and form the cross product with [0 0 1]'. What would be the easiest way of doing this?

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

采纳的回答

James Tursa
James Tursa 2019-4-22
编辑：James Tursa 2019-4-22
M = your 3xN matrix
v = your 3x1 vector
result = cross(M,repmat(v,1,size(M,2)));
  10 个评论
显示 8更早的评论
Hollis Williams
Hollis Williams 2019-4-22
Yes, I was using the size() command but made a typo, it is working now. I have created three 1 x 100 arrays, is it possible to put these together into one 3 x 100 array? So if A,B,C are all 1 x 100 arrays, I would need
A
B
C

请先登录，再进行评论。

更多回答（1 个）

Matt J
Matt J 2019-4-22
编辑：Matt J 2019-4-22
This way avoids repmatting, which may be desirable when N is large.
M = your 3xN matrix
v = your 3x1 vector
result=xprodmat(v)*M./vecnorm(M);
where
function A=xprodmat(a)
%Matrix representation of a cross product
%
% A=xprodmat(a)
%
%in:
%
% a: 3D vector
%
%out:
%
% A: a matrix such that A*b=cross(a,b)
if length(a)<3, error 'Input must be a vector of length 3'; end
ax=a(1);
ay=a(2);
az=a(3);
A=zeros(3);
A(2,1)=az; A(1,2)=-az;
A(3,1)=-ay; A(1,3)=ay;
A(3,2)=ax; A(2,3)=-ax;
end
  1 个评论
Hollis Williams
Hollis Williams 2019-4-22
I don't think repmatting is the thing which is most computationally costly in my code but I will bear this in mind to see if I need to save time later on, thanks a lot.

请先登录，再进行评论。

类别

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

标签

产品


版本

R2018b

Community Treasure Hunt

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

Start Hunting!

Translated by