C++小例子——01 向量类myVector

用C++定义了一个向量类myVector,其包含三个成员变量(x,y,z),可通过三个分量来构造,也可通过另一个myVector对象来构造,重载操作符+、-、*、/来完成相关运算,重载赋值符=,重载输出流<<来实现输出,定义友元函数dotProduct、crossProduct、norm来计算点乘、叉乘、求模……
摘要由CSDN通过智能技术生成

向量(vector) a ⃗ = [ x , y , z ] \vec a=[x,y,z] a =[x,y,z]应用较广,比如节点的坐标 x ⃗ = [ x , y , z ] \vec x=[x,y,z] x =[x,y,z]、面的面积矢量 S ⃗ f = [ S f x , S f y , S f z ] \vec{S}_f=[Sf_x,Sf_y,Sf_z] S f=[Sfx,Sfy,Sfz],面空间速度矢量 u ⃗ = [ u , v , w ] \vec u=[u,v,w] u =[u,v,w]等,都离不开它。同时,向量还有点乘、叉乘、求模(也称为幅值或二范数)的运算,以及向量与向量或标量的加减操作、向量与标量的乘除运算等。那么,咱们可以用C++定义一个向量类myVector,其包含三个成员变量 x , y , z x,y,z x,y,z,可通过三个分量来构造,也可通过另一个myVector对象来构造,重载操作符+、-、*、/来完成相关运算,重载赋值符=,重载输出流<<来实现输出,定义友元函数dotProduct、crossProduct、norm来计算点乘、叉乘、求模……

注意:这里的myVector和C++里自带的库文件vector不同,这里只是含3个分量的空间矢量,而C++中的vector则更类似于一维变长度数组。

矢量类 myVector的框架如下:

成员变量

  1. 矢量在x方向的三个分量
  2. 矢量在y方向的三个分量
  3. 矢量在z方向的三个分量

成员函数

  1. 构造函数1,由三个分量构造,缺省值为0
  2. 构造函数2,由另一个矢量构造,即,拷贝构造
  3. 析构函数,由于没有指针和动态释放空间,所以啥也不用做
  4. print函数,输出(x, y, z)
  5. 提取x,y,z分量的getX,getY,getZ函数
  6. 设置x,y,z分量的setX,setY,setZ函数
  7. 同时设置x,y,z分量的setVal函数,缺省值是0
  8. 赋值函数=,两个矢量间的赋值,即,拷贝赋值
  9. 操作符+:矢量+矢量,矢量+标量
  10. 操作符-:矢量-矢量,矢量-标量
  11. 操作符+=:矢量+=矢量,矢量+=标量
  12. 操作符-=:矢量-=矢量,矢量-=标量
  13. 操作符* :矢量 * 标量,标量 * 矢量
  14. 操作符/ :矢量 / 标量
  15. 操作符 *=:矢量 *= 标量
  16. 操作符 /=:矢量 /= 标量
  17. 点乘dotProduct:,完成两向量的点积
  18. 叉乘crossProduct:,完成两向量的叉积
  19. 输出流 <<:输出向量
  20. 幅值norm:计算向量幅值/二范数/模
  21. 单位化unit:向量除以其幅值

写好的程序如下:

myVector.h头文件

// 2020-06-09
// class myVector
// definition and functions of a vector [x, y, z]
// can be used for points, velocity, etc. in CFD
#ifndef myVector_H
#
  • 4
    点赞
  • 18
    收藏
    觉得还不错? 一键收藏
  • 6
    评论
以下成员: 1. 私有变量: - `double* data`:动态分配的数组,存储向量的每个分量的值; - `int size`:向量的维数。 2. 构造函数: - `MyVector(int size, double value = 0.0)`:构造一个维数为`size`,每个分量的值都为`value`的向量。 3. 拷贝构造函数: - `MyVector(const MyVector& v)`:构造一个与向量`v`相同的新向量。 4. 析构函数: - `~MyVector()`:释放动态分配的数组空间。 5. 成员函数: - `int getSize() const`:返回向量的维数。 - `double get(int i) const`:返回向量第`i`个分量的值。 - `void set(int i, double value)`:设置向量第`i`个分量的值为`value`。 - `double norm() const`:返回向量的模长。 - `void normalize()`:将向量归一化。 - `MyVector add(const MyVector& v) const`:返回当前向量向量`v`的和。 - `double dot(const MyVector& v) const`:返回当前向量向量`v`的点积。 - `MyVector cross(const MyVector& v) const`:返回当前向量向量`v`的叉积。 示例代码: ```cpp #include <iostream> #include <cmath> class MyVector { private: double* data; int size; public: MyVector(int size, double value = 0.0) { this->size = size; data = new double[size]; for (int i = 0; i < size; i++) { data[i] = value; } } MyVector(const MyVector& v) { size = v.size; data = new double[size]; for (int i = 0; i < size; i++) { data[i] = v.data[i]; } } ~MyVector() { delete[] data; } int getSize() const { return size; } double get(int i) const { return data[i]; } void set(int i, double value) { data[i] = value; } double norm() const { double sum = 0.0; for (int i = 0; i < size; i++) { sum += data[i] * data[i]; } return std::sqrt(sum); } void normalize() { double n = norm(); if (n != 0.0) { for (int i = 0; i < size; i++) { data[i] /= n; } } } MyVector add(const MyVector& v) const { MyVector w(size); for (int i = 0; i < size; i++) { w.data[i] = data[i] + v.data[i]; } return w; } double dot(const MyVector& v) const { double sum = 0.0; for (int i = 0; i < size; i++) { sum += data[i] * v.data[i]; } return sum; } MyVector cross(const MyVector& v) const { if (size != 3 || v.size != 3) { std::cerr << "Error: cross product is defined only for 3-dimensional vectors.\n"; exit(1); } MyVector w(size); w.data[0] = data[1] * v.data[2] - data[2] * v.data[1]; w.data[1] = data[2] * v.data[0] - data[0] * v.data[2]; w.data[2] = data[0] * v.data[1] - data[1] * v.data[0]; return w; } }; int main() { MyVector v(3, 1.0); std::cout << "v = [" << v.get(0) << ", " << v.get(1) << ", " << v.get(2) << "]\n"; std::cout << "v.norm() = " << v.norm() << "\n"; v.normalize(); std::cout << "v.normalize() = [" << v.get(0) << ", " << v.get(1) << ", " << v.get(2) << "]\n"; MyVector u(3, 2.0); std::cout << "u = [" << u.get(0) << ", " << u.get(1) << ", " << u.get(2) << "]\n"; std::cout << "u.norm() = " << u.norm() << "\n"; u.normalize(); std::cout << "u.normalize() = [" << u.get(0) << ", " << u.get(1) << ", " << u.get(2) << "]\n"; MyVector w = v.add(u); std::cout << "w = v + u = [" << w.get(0) << ", " << w.get(1) << ", " << w.get(2) << "]\n"; double dot = v.dot(u); std::cout << "v.dot(u) = " << dot << "\n"; MyVector cross = v.cross(u); std::cout << "v.cross(u) = [" << cross.get(0) << ", " << cross.get(1) << ", " << cross.get(2) << "]\n"; return 0; } ```

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 6
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值