向(矢)量叉乘法编程实现
在一个朋友的邀请下,我特意用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 }
补充:软件开发 , C语言 ,