练习42碰撞和热力学.DOC

练习42 碰撞和热力学 知识背景 碰撞是在力学中经常遇到的一类问题，求解这类问题一般要用到动量守恒、能量守恒定律以及其他物理规律。相信很多的读者在中学时就经常与小球碰撞之类的问题打交道，对于这类问题的模型应该还比较熟悉的。热力学中分子的运动规律也是研究物理的人所乐意探讨的，我们在本练习中也将涉及到这方面的内容。 主要内容 【本练习讲述的知识点】 本练习首先介绍小球的弹性碰撞问题，训练读者将实际问题转化为数学模型并进而用MATLAB语句加以实现的能力。我们还将介绍热力学中分子速度的求解问题。本练习中用到了roots命令，在编程中使用了限制性语句。请读者仔细体会这些内容。 练习过程 小球碰撞 我们先来看一个中学时很熟悉的题目。 有一个质量为m的小球以速度u0。正面撞击质量为M的静止小球，若碰撞是完全弹性碰撞，求碰撞后两球的速度，并讨论速度与两球质量比M/m之间的函数关系。 我们来分析这个问题。完全弹性碰撞即无能量损失的碰撞，这大大方便了问题的讨论。我们假设碰撞后两求的速度与u0同向，球m的速度为u ，M的速度为v，设质量比为k=M/m，取v r=v/u0，u r=u/u0。列出动量守恒和能量守恒方程如下： mu0=mu+Mv 动量守恒 动能守恒 对这两格式子进行变化，化为下面的形式： 整理后，得到 m球的能量损失为 整理后，得到 这是一个一元二次方程，我们可以很容易地求出u r。编制程序，并输入命令区： K=logspace(-1,1,11) for i=1:length(K) u1=roots([(1+1/K(i)),-2/K(i),(1/K(i)-1)]); u(i)=u1(abs(u1-1).001); end v=(1-u)./K; de=1-u.*u u,v semilogx(K,[u;v;de]),grid gtext(de) gtext(v) gtext(u) 得到的图形如图42-1所示（u r即图中的u，v r即图中的v）。 我们来分析图线的意义。对于m球相对动能损失 de来说，随着K值增大，de先增大，后变小。当K值为1时，de 达到最大值，即全部动能都传给了M球。对于m球相对速度u来说。随K值增大，一直减小。 K=1时，u=0；K1时，u0。对于v，则随着K值增大，v持续减小，K=1时，即M=m时，M球速度等于u。 我们再来看一下得到的计算结果： K=Columns 1 through 7 0.1000 0.1585 0.2512 0.3981 0.6310 1.0000 1.5849 Columns 8 through 11 2.5119 3.9811 6.3096 10.0000 de=Columns 1 through 7 0.3306 0.1724 0.6418 0.8147 0.9488 1.0000 0.9488 Columns 8 through 11 0.8147 0.6418 0.4724 0.3306 u=Columns 1 through 7 0.8182 0.7264 0.5985 0.4305 0.2263 0 -0.2263 Columns 8 through 11 -0.4305 -0.5985 -0.7264 -0.8182 v=Columns 1 through 7 1.8182 1.7264 1.5985 1.4305 1.2263 1.0000 0.7737 Columns 8 through 11 0.5695 0.4015 0.2736 0.1818 从计算结果中，我们可以 很清楚地看出变化的规律。利用所得数据进行插值计算，我们可以求出在小球在其他质量比时的速度。 （2）热力学 在热力学中，往往要求分子运动的速度。将涉及到指数运算，用人工方法剑 显的比较繁杂。我们通过一个经典的例子来学习用MATLAB解热力学问题。 试求解300k温度下氮气分子运动速度分布曲线，并求速度在（300，500）m/s范围的分子所占的比例。 我们来分析这个问题。对于求分子运动速度的问题，要用到麦克斯韦速度分布率。形势如下： 这是一个较为复杂的数学公式，而且是经常用到的公式。我们考虑编制一个函数M文件，这样我们便能够方便的实现调用。 编制M文件如图42－2所示。 将文件命名为mksw.m,