向（矢）量叉乘法编程实现

在一个朋友的邀请下，我特意用c语言实现了三维向量的叉乘法，因为他在网上找不到关于向量乘法的相关代码。既然写出来了，我就想和大家一起分享，在写的时候对数据结构考虑的不多，重在实现运算过程，有好的建议的朋友希望踊跃留言。

 

  ///main.c/

#include <stdio.h>
#include "vector_multiply.h"
int main(void)
{
 float a[3] = {1.0, 2.0, 3.0}, b[3] = {7.0, 8.0, 9.0}, c[3] = {0};
 
 vector_multiply( a,  b, c);
 printf("%.1f:%.1f:%.1f\n", c[0], c[1], c[2]);
 
 return(0);
}

 

///vector_multiply.c///

void vector_multiply(float a[], float b[], float c[]) {
 int i, j, k, n = 0;

/*according the the formula we can get these marks*/
 int d[3][4] = { { 0, 1, -1, 0 }, { 0, -1, 1, 0 }, { 0, 1, -1, 0 } };

 

 for (k = 0; k < 3; k++) {
  for (i = 0; i < 3; i++) {
   for (j = 0; j < 3; j++)
    if (k != i && k != j)
     c[k] += a[i] * b[j] * d[k][n++];
  }
  n = 0;
 }

 return;
}

 

vector_multiply.h//

#ifndef VECTOR_MULTIPLY_H_
#define VECTOR_MULTIPLY_H_

void vector_multiply(float a[], float b[], float c[]);

#endif /*VECTOR_MULTIPLY_H_*/

///

gcc main.c vector_multiply.c -o  VectorOperations

 

 

最后再留下向量叉乘的公式：

向量：u=(u1,u2,u3) v=(v1,v2,v3)

叉积公式：u x v = | i    j     k|
                 |a1    b1   c1|
                 |a2    b2   c2|

               = { u2v3-v2u3 , u3v1-v3u1 , u1v2-u2v1 }