标签:
问题及代码:
/* *Copyright (c)2015,烟台大学计算机与控制工程学院 *All rights reserved. *文件名称:project.cpp *作 者:陈文青 *完成日期:2015年5月31日 *版 本 号:v1.0 * *问题描述: (6)与圆心相连的直线:给定一点p,其与圆心相连成的直线,会和圆有两个交点,如图。在上面定义的Point(点)类和Circle(圆)类基础上, 设计一种方案,输出这两点的坐标。 *程序输入: *程序输出: */ #include <iostream> #include<cmath> using namespace std; class Circle; //在求交点的函数中用到了之后声明的圆类,所以要提前声明 //定义点类 class Point { public: Point(double xx=0,double yy=0):x(xx),y(yy) {} //点的构造函数 ~Point() {} friend ostream &operator<<(ostream &,const Point &); //重载输出运算符 friend void get_crossover_point(Circle c,Point p); //直线与圆的交点 protected: double x; double y; }; ostream &operator<<(ostream & out,const Point &p) //重载输出运算符 { out<<"("<<p.x<<","<<p.y<<")"; return out; } //由Point类派生出一个Circle(圆)类 class Circle:public Point { public: Circle(double xx=0,double yy=0,double rr=0):Point(xx,yy),r(rr) {} //圆的构造函数 ~Circle() {} friend ostream &operator<<(ostream &,const Circle &); //重载圆的输出运算符 friend void get_crossover_point(Circle c,Point p); //求直线与圆的交点 private: double r; }; ostream &operator<<(ostream &out,const Circle &c) //重载圆的输出运算符 { out<<"原点O:("<<c.x<<", "<<c.y<<"),半径r="<<c.r; return out; } void get_crossover_point(Circle c,Point p) //直线与圆的交点 { Point p1,p2; p1.x = (c.x + sqrt(c.r*c.r/(1+((c.y-p.y)/(c.x-p.x))*((c.y-p.y)/(c.x-p.x))))); p2.x = (c.x - sqrt(c.r*c.r/(1+((c.y-p.y)/(c.x-p.x))*((c.y-p.y)/(c.x-p.x))))); p1.y = (p.y + (p1.x -p.x)*(c.y-p.y)/(c.x-p.x)); p2.y = (p.y + (p2.x -p.x)*(c.y-p.y)/(c.x-p.x)); cout<<"直线与圆的交点为:"; cout<<"p1:"<<p1<<'\t'<<"p2"<<p2<<endl; } int main( ) { Circle c1(0,0,1),c2(1,0,3); Point p1(1,0),p2(5,0); cout<<"点p1:"<<p1<<'\t'<<"圆c1:"<<c1<<endl; get_crossover_point(c1,p1); cout<<"点p2:"<<p2<<'\t'<<"圆c2:"<<c2<<endl; get_crossover_point(c2,p2); return 0; }
运行结果:
知识点总结:
在求交点坐标里,get_crossover_point作为友元函数,既引用了Point类,又引用了Circle类,所以声明的时候,在两个类内都要声明一下,这样才不用担心Point内protected的函数不能正常调用。而这样做,也必须把后定义的Circle类在Point类之前声明一下。不然也会导致Circle类内的数据成员无法正常引用。
学习心得:
这个程序实际上是有bug存在的:
在测试数据的时候,我偶然发现运行结果输出了以下结果。
因为在老师给的公式中,如果分母为0的话,再强大的计算机也是算不出来的,因为它违背了数学原理。所以运行结果给我的输出结果是 难(nan)。
所以,仍然有待修改。标签:
原文地址:http://blog.csdn.net/destiny_forever/article/details/46288847