c++求两点的距离利用友元_C++ 求解球内某点沿某向量到球面的距离

最新推荐文章于 2022-07-13 22:07:41 发布
最新推荐文章于 2022-07-13 22:07:41 发布
阅读量301 收藏
点赞数
文章标签： c++求两点的距离利用友元
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/weixin_32304579/article/details/112404425
版权
本文介绍了如何在C++中计算三维空间中球体内一点沿给定向量到球面的距离，并详细阐述了数学公式推导、代码实现以及处理特殊情况（如点位于球面上）的方法。内容包括直线参数方程、球面方程的结合，以及针对相切和不相切情况的解决方案。
摘要由CSDN通过智能技术生成

00 本文内容

  • 问题描述
    • 图像
    • 输入输出
  • 数学公式推导
    • 三维坐标下直线的参数方程
    • 直线与球面交点求解方法
    • 判断距离是否与方向向量相同
  • C++代码实现
    • 函数化封装
    • 可根据需要只修改函数参数
  • 特殊情况
    • 实际项目中遇到返回无效值
    • 当该点在球面上，讨论相切与不相切

01 问题描述

问题描述：在三维坐标内，已知一球体sphere内有一点focus，该点沿某方向向量n0延伸，与球面交于一点f。求解该点focus沿该方向向量到球面的距离，即focus和f两点间距离。

输入：该点focus的坐标、方向向量n0的坐标、sphere球心坐标与半径R

输入：距离ans

0dda0cd51177563ba6ab9537c1556a45.png
淡蓝色球体中有一点focus，其沿着黑色箭头方向与球面交于点f2，求出两点间距离

02 数学公式推导

刚开始遇到这个问题，大致觉得可以编程实现，但是很多高中几何知识都忘了。经查阅，将数学描述放置于下方。

已知focus点的坐标

和方向向量n0
，可以得到两者所在直线的参数方程

即

而已知球面方程

将参数方程组带入球面桌面，可以得到关于参数t的一元二次方程组

其中A、B、C可用

等参数表示

求解可得

带回直线的参数方程，可得两交点坐标，可分别与focus求两点间距离。

最后判断是否沿方向向量方向，可以通过两向量对应部分符号是否相同来判断。

03 C++代码实现

这是最初版本，可根据实际需要修改。

#include <math.h>
#include<iostream>
using namespace std;

typedef struct Points
{
    double x,y,z;
}Point;

double fun(Point focus,Point normalVector,Point sphere,double R)
{
    //直接到直线参数方程带入球面式子，得到了关于参数t的一元二次方程
    double A,B,C;//ABC为关于t的一元二次方程的参数
    A=pow(normalVector.x,2) + pow(normalVector.y,2) + pow(normalVector.z,2);
    B=2*( (focus.x-sphere.x)*normalVector.x + (focus.y-sphere.y)*normalVector.y + (focus.z-sphere.z)*normalVector.z );
    C=pow(focus.x-sphere.x,2) + pow(focus.y-sphere.y,2) + pow(focus.z-sphere.z,2) - pow(R,2);
    //cout<<"A B C t1 t2: "<<A<<" "<<B<<" "<<C<<" ";

    //现在可得 A*t*t+B*t+C=0,解一元二次方程
    double t1,t2;
    t1 = (-B+sqrt(pow(B,2)-4*A*C)) / (2*A);
    t2 = (-B-sqrt(pow(B,2)-4*A*C)) / (2*A);
    //cout<<t1<<" "<<t2<<endl;

    //将t带入带原直线参数方程中，可以得到交点坐标
    Point intersection1,intersection2;//交点
    intersection1.x=normalVector.x*t1+focus.x;
    intersection1.y=normalVector.y*t1+focus.y;
    intersection1.z=normalVector.z*t1+focus.z;

    intersection2.x=normalVector.x*t2+focus.x;
    intersection2.y=normalVector.y*t2+focus.y;
    intersection2.z=normalVector.z*t2+focus.z;

    //现在知道求内点focus，两交点intersection，求法向量方向的距离
    if((intersection1.x-focus.x)*normalVector.x>0||(intersection1.y-focus.y)*normalVector.y>0||(intersection1.z-focus.z)*normalVector.z>0)
    {
        //这时候说明 intersection1 是和法向量同向的点
        double ans=sqrt(pow(intersection1.x-focus.x,2)+pow(intersection1.y-focus.y,2)+pow(intersection1.z-focus.z,2));
        return ans;
    }
    else
    {
        //这时候说明 intersection2 是和法向量同向的点
        double ans=sqrt(pow(intersection2.x-focus.x,2)+pow(intersection2.y-focus.y,2)+pow(intersection2.z-focus.z,2));
        return ans;
    }
}
int main()
{
    Point focus;//圆内该点
    Point normalVector;//法向量
    Point sphere;//园的圆心
    double R;//圆的半径

    cin>>focus.x>>focus.y>>focus.z;
    cin>>normalVector.x>>normalVector.y>>normalVector.z;
    cin>>sphere.x>>sphere.y>>sphere.z>>R;
    double ans=fun(focus,normalVector,sphere,R);
    cout<<"ans="<<ans<<endl;

}

04 特殊情况

但是很多实际情况中（给老师用后有错），该点focus很可能会出现在球面上，此时使用上述代码会出错。

因此接下来将讨论focus在球面上的情况：

  1. 方向向量与球面相切
  2. 方向向量与球面不相切，且该方向向量指向球体外部
  3. 方向向量与球面不相切，且该方向向量指向球体内部

情况1：只有一个交点，也就是focus本身，这时直接返回0

情况2：指向外部，则此时方向向量方向是0，返回0

情况3：指向内部，返回值应该为两交点间距离

修改后的代码如下：

#include <math.h>
#include<iostream>
using namespace std;

typedef struct Points
{
    double x,y,z;
}Point;

double fun(Point focus,Point normalVector,Point sphere,double R)
{
    //直接到直线参数方程带入球面式子，得到了关于参数t的一元二次方程
    double A,B,C;//ABC为关于t的一元二次方程的参数
    A=pow(normalVector.x,2) + pow(normalVector.y,2) + pow(normalVector.z,2);
    B=2*( (focus.x-sphere.x)*normalVector.x + (focus.y-sphere.y)*normalVector.y + (focus.z-sphere.z)*normalVector.z );
    C=pow(focus.x-sphere.x,2) + pow(focus.y-sphere.y,2) + pow(focus.z-sphere.z,2) - pow(R,2);
    // cout<<"A B C t1 t2: "<<A<<" "<<B<<" "<<C<<" ";

    //现在可得 A*t*t+B*t+C=0,解一元二次方程
    double t1,t2;
    t1 = (-B+sqrt(pow(B,2)-4*A*C)) / (2*A);
    t2 = (-B-sqrt(pow(B,2)-4*A*C)) / (2*A);
    // cout<<t1<<" "<<t2<<endl;

    if(t1==0&&t2==0)//该点在球面上，且方向向量与球面相切 ，无论哪个方向都是0 
        return 0;
    
    //将t带入带原直线参数方程中，可以得到交点坐标
    Point intersection1,intersection2;//交点
    intersection1.x=normalVector.x*t1+focus.x;
    intersection1.y=normalVector.y*t1+focus.y;
    intersection1.z=normalVector.z*t1+focus.z;
    intersection2.x=normalVector.x*t2+focus.x;
    intersection2.y=normalVector.y*t2+focus.y;
    intersection2.z=normalVector.z*t2+focus.z;

    double f1=(intersection1.x-focus.x)*normalVector.x+(intersection1.y-focus.y)*normalVector.y+(intersection1.z-focus.z)*normalVector.z;
    double f2=(intersection2.x-focus.x)*normalVector.x+(intersection2.y-focus.y)*normalVector.y+(intersection2.z-focus.z)*normalVector.z;

    //现在知道求内点focus，两交点intersection，求法向量方向的距离
    if(f1>0||f2<0)
    {
        //这时候说明 intersection1 是和法向量同向的点
        double ans=sqrt(pow(intersection1.x-focus.x,2)+pow(intersection1.y-focus.y,2)+pow(intersection1.z-focus.z,2));
        return ans;
    }
    else if(f1<0||f2>0)
    {
        //这时候说明 intersection2 是和法向量同向的点
        double ans=sqrt(pow(intersection2.x-focus.x,2)+pow(intersection2.y-focus.y,2)+pow(intersection2.z-focus.z,2));
        return ans;
    }

    return -1;
}
int main()
{
    Point focus;//圆内该点
    Point normalVector;//法向量
    Point sphere;//园的圆心
    double R;//圆的半径

    cin>>focus.x>>focus.y>>focus.z;
    cin>>normalVector.x>>normalVector.y>>normalVector.z;
    cin>>sphere.x>>sphere.y>>sphere.z>>R;
    double ans=fun(focus,normalVector,sphere,R);
    cout<<"ans="<<ans<<endl;

}
真的是单大宝
关注
Open3D(C++) Ransac拟合球体(详细过程版)
点云侠的博客
03-10 1149
Ransac拟合球体的C++版本Open3D详细过程实现。
友元机制计算两点之间的距离
11-27
设计一个点（Point）类，其中包括一对坐标数据成员、一个求两点之间距离的友元接口和显示坐标点的成员函数，并用数据进行测试。 （2）程序设计说明 ① 在头文件中声明类，在源文件中定义类的成员函数，采用构造函数...
空间直线与球面的相交算法
qq_21291397的博客
07-10 5526
文章目录1. 原理推导1.1. 直线公式1.2. 求交2、代码实现：2.1 关键代码：2.2 main（）函数的调用 1. 原理推导 1.1. 直线公式 空间直线的数学定义是，已知直线L上一点M0(x0,y0,c0)，以及直线L的方向向量s(m,n,p)，那么空间直线L的方程为： 以上是空间直线的标准式方程（点向式方程）。令上面式子的比值为t，那么直线的参数式方程为： 对于知道线段的起点O和终点E，显然方向向量为D=E−O。这时，根据射线的向量方程，线段上某一点P为 转化为参数方程就是： 并且，采取
C++求解圆心、半径、直线与圆交点
weixin_44529257的博客
07-13 3101
C++求解圆心、半径、直线与圆的交点
已知圆心和半径获取圆与折线的交点C++算法实现
u014629983的博客
08-23 963
算法逻辑很简单循环遍历折线线段，获取圆与线段的交点。 typedef CArray&lt;CPoint, CPoint&gt; POINT3DLIST; // 获取圆与折线交点 // 输入参数 CenterPt 圆心，r 半径，line 折线 // 输出参数 CrossPs 圆与折线交点 // 返回交点个数 int GetCircLineCross(const CPoint&amp; Ce...
输入球心三维坐标及球半径求算球体的球面坐标数据(C/C++)
verystory的博客
01-11 7238
本博文 实现的功能为：输入球心三维坐标以及球半径 得到球体的球面坐标数据这样一个功能(当然球面的点数据密集程度可控) 使用C语言写的 算法很简单 代码也很简单..
Exe7_0友元_C++友元函数编程（两种有友元方式）_cloudzfv_
09-29
本示例"Exe7_0友元_C++友元函数编程（两种有友元方式）_cloudzfv_"旨在帮助初学者理解并掌握C++中的友元函数的两种主要用法。现在，我们将详细讨论这两种友元函数的实现方式。 1. **友元函数** 友元函数不是类的...
C++计算两点之间的距离（两种实现方式）
04-10
计算两点之间距离（友元函数）
C++ 工程师校招面试考点汇总_面试_c/C++_C++校招知识点_c++面试考点_
10-04
在C++工程师的校招面试中，面试官通常会考察一系列关键知识点，这些知识点涵盖了语言基础、编程技巧、数据结构、算法以及系统设计等多个方面。以下是对这些面试考点的详细阐述： 1. **C++语言基础**：面试者需要...
C++运算符重载 成员函数与友元函数详解
12-31
代码如下:#include<iostream>using namespace std;class A{ int x,y; public: A(int xx,int yy):x(xx),y(yy){} A(){x=0;y=0;} A operator+(const A&b) //不加const限定，也可以 { return A(x+b.x,y+b.y);...
C++求两点之间距离
10-03
用C++写的一个球两点之间距离的小程序。
C++实现三维空间中点到点、点到直线、点到平面的距离计算
12-15
C++实现三维空间中点到点、点到直线、点到平面的距离计算。
计算球面上任意两点间的球面距离（C++实现）
孙悟空
10-23 5618
本文给出了求解球面上两点球面距离的原理与C++实现代码，可以计算球面两点的欧式距离（弦长）、圆心角、球面距离。
计算几何学 | 距离 | Distance | C/C++实现
qq_40998706的博客
02-16 2179
问题描述 输出线段s1、s2之间的距离。设s1的端点为p0、p1，s2的端点为p2、p3。 输入： 第1行输入问题数q。接下来q行给出q个问题。各问题线段s1、s2的坐标按照以下格式给出： xp0x_{p0}xp0​ yp0y_{p0}yp0​ xp1x_{p1}xp1​ yp1y_{p1}yp1​ xp2x_{p2}xp2​ yp2y_{p2}yp2​ xp3x_{p3}xp3​ yp3y_{p...
三维点沿指定向量方向到平面的距离计算方法及C++代码实现
zeeq的博客
12-29 3372
  设平面外一点为P(p1,p2,p3)，指定的方向向量为d=(d1,d2,d3)，平面Q方程为Ax+By+Cz=D，其法向量为n=(A,B,C)。则计算点P沿方向向量d到平面Q的距离，可以通过下式计算         (P+t* d)* n=D   [1]   其中t为待求的常数。下面将[1]式分两个步骤解释一下该式子的意义（理解了以后其实就很简单）：         P’ = P+t* d  [2]         P’ * n=D  [3]   [2]式表示三维点P沿着方向向量d移动t倍距离，到达点P
c++ 求解球面距离
Arana的博客
03-21 962
文章目录【需求规格说明】1、题目2、问题描述3、问题分析【算法设计】1、设计思想2、设计表示：【调试报告】调试问题1调试问题2调试问题3调试问题4【附录】1、源代码2、测试数据3、运行结果 【需求规格说明】 1、题目 求解球面距离： 1）地球的平均半径为6371千米，已知地球上两个城市A、B的经度和纬度，编程序求出这两个城市之间的球面距离（一般来说，球面上任意两点A和B都可以与球心确定唯一的大圆，...
向量距离(Distance)
qq_42898629的博客
10-31 9250
向量是一个有向线段，计算两个向量之间的距离，其实就是计算两个点之间的距 离。 在 Unity 引擎内，我们可以直接通过 Vector3.Distance(v1, v2)来得到 两个向量之间的距离。 数学公式: 例: 公式解析: 1.求两个向量之间的距离，本质就是求两个点之间的距离； 2.用其中的一个点与另外一个点进行相减，本质是两个向量中每个分量的相减运算，最终得到一个新的...
C++-点到点、点到线、共线等常用计算
ゃ﹏淸凨ゞ
12-11 2527
C++常用数学计算1. 一点到另外两点组成直线的距离2.获取平面两点间的距离3计算两个向量的夹角4. 判断三点共线 1. 一点到另外两点组成直线的距离 已知点p /p1/p2,求p到p1 p2的距离 根据点到直线的距离公式：d = (Ax0 + By0 + C) / sqrt(A * A + B * B) A = y2 - y1 B = x1 - x2； C = x2 * y1 - x1 * y2; 程序实现 //定义Point2f结构体 struct Point2f { float x; floa
点到线段的距离（向量）
PL博客
04-27 4684
一，很自然的想法就是判断点和线段的位置关系已经点的垂点是否在线段上，然后每种情况分别计算。 这种做法显然可行，只是觉得比较繁琐，想要更简单一点的计算方式； 二， 向量做法， 推荐看这篇博客https://blog.csdn.net/angelazy/article/details/38489293 我觉说得挺好的，看了很容易理解 我自己写的代码 ...
c++求两点的距离利用友元_【C++类和对象】计算两点间距离的两种方法
最新发布
06-06
好的，我来为您解答。 首先，我们需要定义一个 `Point` 类，表示一个点的坐标。这个类可以有两个私有成员变量 `x` 和 `y`，分别表示横坐标和纵坐标。 然后，我们可以在 `Point` 类中定义一个公有成员函数 `distanceTo`，用于计算当前点和另一个点之间的距离。这个函数接受一个 `Point` 类型的参数，表示另一个点。函数的实现可以使用勾股定理，计算两点之间的距离。 接下来，我们可以使用友元函数来计算两个点之间的距离。友元函数是在类外部定义的函数，但可以访问类的私有成员。我们可以在 `Point` 类中声明一个友元函数 `distanceBetween`，用于计算两个点之间的距离。这个函数接受两个 `Point` 类型的参数，分别表示两个点。函数的实现可以直接调用 `distanceTo` 函数，计算出两个点之间的距离。 最后，在主函数中，我们可以创建两个 `Point` 对象，分别表示两个点，并调用 `distanceTo` 函数和 `distanceBetween` 函数，计算出两个点之间的距离。 下面是代码示例： ```cpp #include <iostream> #include <cmath> using namespace std; class Point { private: double x, y; // 点的坐标 public: Point(double _x, double _y): x(_x), y(_y) {} double distanceTo(Point p) { return sqrt(pow(x - p.x, 2) + pow(y - p.y, 2)); } friend double distanceBetween(Point p1, Point p2); }; double distanceBetween(Point p1, Point p2) { return p1.distanceTo(p2); } int main() { Point p1(0, 0); Point p2(3, 4); cout << "The distance between (" << p1.distanceTo(p2) << ")" << endl; cout << "The distance between (" << distanceBetween(p1, p2) << ")" << endl; return 0; } ``` 运行结果： ``` The distance between (5) The distance between (5) ``` 另外，我们也可以将 `distanceBetween` 函数定义为 `Point` 类的静态成员函数，这样就不需要使用友元函数了。下面是代码示例： ```cpp #include <iostream> #include <cmath> using namespace std; class Point { private: double x, y; // 点的坐标 public: Point(double _x, double _y): x(_x), y(_y) {} double distanceTo(Point p) { return sqrt(pow(x - p.x, 2) + pow(y - p.y, 2)); } static double distanceBetween(Point p1, Point p2) { return p1.distanceTo(p2); } }; int main() { Point p1(0, 0); Point p2(3, 4); cout << "The distance between (" << p1.distanceTo(p2) << ")" << endl; cout << "The distance between (" << Point::distanceBetween(p1, p2) << ")" << endl; return 0; } ``` 运行结果： ``` The distance between (5) The distance between (5) ```
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值