Python程序设计期末作业.doc

Python程序设计期末作业.doc

ID:55632088

大小:84.00 KB

页数:6页

时间:2020-05-21

Python程序设计期末作业.doc_第1页
Python程序设计期末作业.doc_第2页
Python程序设计期末作业.doc_第3页
Python程序设计期末作业.doc_第4页
Python程序设计期末作业.doc_第5页
资源描述:

《Python程序设计期末作业.doc》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库

1、球-弹簧物理模型模拟程序一、编写目的1.模拟一个非线性动力系统,体会系统参数的混沌现象。2.熟悉VPython的使用。二、物理模型描述模型空间的四周及顶面为刚性壁,底为一有质量的刚性平板,平板由一根轻弹簧支撑。模型空间内有数个刚性球作自由运动。球的初始位置在空间内均匀分布,初始速度满足正态分布。所有碰撞均为完全弹性碰撞。并假设平板只能作Z方向一维运动。三、设计思想由物理学中的能量守恒和动量守恒推导碰撞前后物体的速度变化。每次迭代先使物体状态按时间小量dt发展,再判断是否有碰撞发生,若有则由前面推导的公式改变相关物体的速度。以适当的速率重复上述步骤,再用visual库绘制出

2、对应图形,则得到这个物理模型的模拟演示动画。四、使用方法直接运行程序,则开始模拟并实时绘制。也可修改程序开始处的几个全局变量,改变模型的参数。五、参考资料VPython官方文档六、说明在Python2.7下调试通过。本程序中没有一行代码直接来自网络、书籍等处,参考资料仅为VPython官方文档。附截图及源代码。程序运行效果截图fromvisualimport*fromrandomimportrandomfrommathimportsin,cos,sqrt,log,pi#物理模型的参数以全局变量形式定义g=9.81#重力加速度k=10#弹簧胡克系数m1=1#球质量m2=10

3、#平板质量N=10#球数R=1#球半径#物体基类,提供设置、返回位置、速度信息,及迭代一个时间小量dt的方法classobject:defsetVel(self,v):self.v=vdefsetPos(self,p):self.p=pdefgetPos(self):returnself.pdefgetVel(self):returnself.vdefgo(self,dt):self.v+=self.a*dtself.p+=self.v*dt#球,继承于物体类classC_Ball(object):def__init__(self,p,v):self.p=pself.v=

4、vself.a=vector(0,0,-g)#平板,继承于物体类,有自己的迭代方法,即加速度与位置有关classC_Board(object):def__init__(self):self.p=vector(0,0,0)self.v=vector(0,0,0)self.a=vector(0,0,0)defgo(self,dt):self.a=self.p*-kself.v+=self.a*dtself.p+=self.v*dt#生成随机位置向量(模型空间内均匀分布)defrandPos():x=(random()-0.5)*2*(15-R-0.5)y=(random()-

5、0.5)*2*(15-R-0.5)z=random()*(30-R*2-0.5)+(R+0.5)returnvector(x,y,z)#生成三维标准正态分布向量,使用Box-Muller公式defrandVel():u1,u2=random(),random()x=sqrt(-2*log(u1))*cos(2*pi*u2)y=sqrt(-2*log(u1))*sin(2*pi*u2)u1,u2=random(),random()z=sqrt(-2*log(u1))*cos(2*pi*u2)returnvector(x,y,z)#物理模型,含有一个平板和N个球classC_

6、System:#放置球,并检查冲突def__init__(self):self.ball=[]foriinrange(N):self.ball.append(C_Ball(randPos(),randVel()*3))j=0whilej

7、:foriinrange(N):self.ball[i].go(dt)self.board.go(dt)foriinrange(N):#检查与平板的碰撞ifself.ball[i].getPos().z-R<=self.board.getPos().z+0.5:self.ball[i].setPos(self.ball[i].getPos()+(self.board.getPos().z-self.ball[i].getPos().z+1.5)*vector(0,0,1))v1=self.ball[i].getVel()v2=

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。