台球问题分析.pptVIP

张乃孝：算法与数据结构——C语言描述 * * 台球运动过程模拟 2006 算法与数据结构作业题 目的：分析物理问题的方法、 熟悉变量、数组、指针以及输入输出函数等C语言内容 问题：完全相同的16个硬球，直径为r，质量为m。不同球之间距离为r/2，15个排列成为等边三角型。另外一个以速度V运动。如果球与壁碰撞的地点与顶点距离小于a，就认为落袋。 要求：请计算各个球落袋的时间次序，在屏幕和文件中输出。 2米 3米 V a 问题的描述 基本思路 每个球受到的作用力是可以计算的，那么每次时间增加量为?t后球的新位置和速度也是可以唯一确定的（牛顿定理）。 在新位置，可以判断球是否已经落袋。 随时间的逐渐增加，所有的球的运动轨迹是可以模拟的。 1、 在桌面上运动的球可能发生的情况： 当球之间的距离小于一定的距离，就会有排斥作用（碰撞）； 当球碰到球台的边缘的时候，就会被反射回来； 当它们到达球袋区域，就会落袋，离开桌面； 2、 物理理论和模型 ? 牛顿运动方程 球之间相互作用力的描述 球与球台边缘碰撞过程力的描述 一、问题的分析 物理过程分析清楚之后，下一步内容就是数据结构和算法的工作了。。。。。。 t1 t2 t3 t4 tn 2-1 数据结构 二、程序设计 1、一些常数是必须提前知道的： #define TOTALBALLNUMBER 17 #define ORIGINX 0 //球阵列的中心位置 #define ORIGINY 0 #define HEIGHT 3.0 //高，Y轴 #define WIDTH 2.0 //宽，X轴 #define CORNER 0.1 //球袋的尺寸，大小有什么影响？ #define INTERVAL 0.1 //球之间的间隔 #define MAXDISTANCE 0.3 //这样所有的球都有一定的初速度，为什么？ 2米 3米 V a 需要一些变量来在程序中如何描述这些球的状态 BilliardsType BallsOntable[TOTALBALLNUMBER ]; BilliardsType FalledBall[TOTALBALLNUMBER ]; 描述球的集合 BilliardsType Ball; 描述单个球 2、定义球的结构 typedef struct { int No; float x,y, vx,vy, ax,ay, mass; } BilliardsType; 有许多算法可以采用：分子动力学方法、几何计算方法等。 我们采用分子动力学的方法，比较容易描述时间尺度上的过程 init(); //各个球的初始位置，运动球的初始速度、位置等 while(桌面上球集合为空) { 计算?t后各个球的位置； 判断是否有球落入袋中。如果有，把它从桌面球集合中除去，加入到落袋球的序列中； ? 计算桌面上各个球在新位置上所受到的作用力； 计算在这些作用力作用下各个球的新的速度； } ? 输出结果（各个球落袋的次序、时间、速度、角度等信息） 2-2、 算法 初试化函数init() 作用力计算ForceCal() 运动位置BallNewState() RemoveFromBallList（） ResultOutput 程序的流程 判断一个集合是否为空: isEmpty(Balllist); 置一个集合为空: emptySet(Balllist); 从集合中去掉一个元素: removeFromSet(Balllist, Ball); 向集合里增加一个元素: addToSet(Balllist, Ball); 计算下一刻的球的新位置：newposition(Balllist); 计算下一刻的球的新速度：newvolecity(Balllist); 检查落入袋中的球： FallinCheck(Balllist, Balllist); 需要以下的一些操作： BilliardsType BallsOntable[TOTALBALLNUMBER ],ball; BilliardsType FalledBall[TOTALBALLNUMBER ]; //描述球的集合 Float t=0,delta_t=0.001; Int counter,CurrentBallNo= TOTALBALLNUMBER; Init(Ontable); //初始化所有台球的位置、速度 EmptySet(FalledBall