©2010-2026 宁波自信网络信息技术有限公司 版权所有
客服电话：0574-28810668 投诉电话：18658249818
浙ICP备2021020529号-1 | 浙B2-20240490
关注我们 ：
1、 东北石油大学 课 程 设 计 课 程 计算物理 题 目 MC方法模拟理想气体分子的麦克斯韦速率分布 院 系 电子科学学院 专业班级 应用物理10-2 学生姓名 学生学号 指导教师
2、 2014年 3月 14日 东北石油大学课程设计任务书 课程 计算物理和MATLAB课程设计 题目 MC方法模拟理想气体分子的麦克斯韦速率分布 专业 应物 姓名 学号 主要内容、基本要求、主要参考资料等 主要内容： 麦克斯韦速率分布函数为 通过利用M-C的方法模拟麦克斯韦速率分布图。 基本要求： 1、给出M-C研究麦克斯韦速率分布的方法及步骤。 2、画出M-C的模拟图，且与理论值相比较。 主要参考资料： [1]马文蔚，解希顺等. 物理
3、学[M],高等教育出版社，2007. [2]张志涌. 精通MATLAB6.5[M]. 北京：北京航空航天大学出版社，2003. [3]彭芳麟.计算物理基础[M]，高等教育出版社，2009. 完成期限 2014.2.28-2014.3.14 指导教师 专业负责人 2014年 2 月 28 日 目 录 目 录 1 第1章 概述 2 第2章 MATLAB的基础知识 3 2.1 Matlab基本命令 3 2.2 Matlab与绘图有关的基本命令 3 第3章 麦克斯韦速率分布
4、的数学描述 5 3.1 速率分布和分布函数 5 3.2 理想气体分子的麦克斯韦速率分布律 6 第4章 Matlab仿真实验及结果讨论 7 第5章 结论 8 参考文献 9 第1章 概述 引言麦克斯韦速率分布律是分子动理论的重要结论之一,它是研究气体分子碰撞、大量分子热运动服从统计规律性等问题的重要理论依据,正确理解它对学习热学非常有用。但是由于推导困难、公式复杂、曲线难画,麦克斯韦速率分布律学习起来比较费力,成了热学理论教学中的一个难点。Matlab是由美国的MathWorks公司于20世纪80年代中期出品的一款商业数学软件,由于优秀的数值计算能力和卓越的数据可视能力,使其
5、很快在数学软件中大放异彩,与Mathematica、Maple、MathCAD一起并称为四大数学软件。通过Matlab的图形用户界面(Graphical User Interfaces,GUI)可以设计出一个界面友好、方便人机交互的图形界面程序。在物理教学中,适当利用Matlab GUI建立直观的物理图像,有助于学生更好地理解物理概念,提高学习兴趣,调动学习积极性 在分子物理学和热学中可知，由大量分子组成的气体，因分子间的频繁碰撞，各个分子的速度大小和方向瞬息万变。任一时刻，某个分子具有多大的运动速率完全是偶然的，可以是零到无穷大之间的任何值。1859年，英国物理学家麦克斯韦（J.C.Max
6、well）从理论上导出了气体分子速率分布律——麦克斯韦速率分布律。 利用MATLAB软件灵活的数值计算与符号计算、简单的语句表达、简洁完善的图形绘制、丰富的工具箱函数和简易的扩展功能，将复杂的数学公式绘制成曲线，并研究单个参数的影响。 第2章 Matlab的基础知识 2.1 Matlab基本命令 表2-1 MATLAB基本命令 主题词 含义 主题词 含义 format 设置数据显示格式 feval 函数求值 who 显示变量名 input 提示输入 whos 显示变量信息 disp 输出 clear 清除内存变
7、量 tic 启动秒表 save 保存工作变量到文件 toc 时间读数（秒） load 从文件装载变量 help 帮助 linspace 区间等分 lookfor 查找 length 获取数组长度 type 列程序清单 size 矩阵大小 which 查找文件目录 max 最大值 double 双精度 min 最小值 str2num 字符串转化为数值 sum 求和 num2str 数值转化为字符串 find 条件检索 2.2 Matlab与绘图有关的基本命令 表2-2 常用作图命令和函数 主题词 含 义
8、 主题词 含 义 plot 基本二维图形 clabel 等高线高度标志 fplot 一元函数图像 grid 格栅 ezplot 画二维曲线的符号命令 hold 图形保持 plot3 空间曲线 axis 定制坐标轴 meshgrid 网格数据生成 view 改变视点 mesh 网面图 subplot 子图 surf 曲面图 figure 新图形窗口 contour 等高线图 clf 清除图形 contour3 三维等高线图 close 关闭图形窗口 title 标题 ylabel y轴说明 在线条多于一
9、条时，若用户没有指定使用颜色，则plot循环使用由当前坐标轴颜色顺序属性（current axes ColorOrder property）定义的颜色，以区别不同的线条。在用完上述属性值后，plot又循环使用由坐标轴线型顺序属性（axes LineStyleOrder property）定义的线型，以区别不同的线条。用法 plot(X,Y) 当X,Y均为实数向量，且为同维向量（可以不是同型向量），X=[x(i)]，Y=[y(i)]，则plot(X,Y)先描出点(x(i)，y(i))，然后用直线依次相连；若X，Y为复数向量，则不考虑虚数部分。若X，Y均为同维同型实数矩阵，X = [X(i)]，Y
10、 = [Y(i)]，其中X(i),Y(i)为列向量，则plot(X,Y)依次画出plot(X(i),Y(i))，矩阵有几列就有几条线；若X，Y中一个为向量，另一个为矩阵，且向量的维数等于矩阵的行数或者列数，则矩阵按向量的方向分解成几个向量，再与向量配对分别画出，矩阵可分解成几个向量就有几条线；在上述的几种使用形式中，若有复数出现，则复数的虚数部分将不被考虑。 plot(Y) 若Y为实数向量，Y的维数为m，则plot(Y)等价于plot(X,Y)，其中x=1：m；若y为实数矩阵，则把y按列的方向分解成几个列向量，而y 的行数为n，则plot(Y)等价于plot(X,Y)其中x=[1;2;…;n]
11、在上述的几种使用形式中，若有复数出现，则复数的虚数部分将不被考虑。plot(X1,Y1,X2,Y2,…)，其中Xi与Yi成对出现，plot(X1,Y1,X2,Y2,…)将分别按顺序取两数据Xi与Yi进行画图。若其中仅仅有Xi或Yi是矩阵，其余的为向量，向量维数与矩阵的维数匹配，则按匹配的方向来分解矩阵，再分别将配对的向量画出。plot(X1,Y1,LineSpec1,X2,Y2,LineSpec2…) 将按顺序分别画出由三参数定义Xi,Yi,LineSpeci的线条。其中参数LineSpeci指明了线条的类型，标记符号，和画线用的颜色。在plot 命令中我们可以混合使用三参数和二参数的形式。
12、 第3章 麦克斯韦速率分布的数学描述 气体中个别分子的速度具有怎样的数值和方向完全是偶然的，但就大量分子的整体来看，在一定的条件下，气体分子的速度分布也遵从一定的统计规律。这个规律也叫麦克斯韦速率分布律。 3.1 速率分布和分布函数 为了描述平衡态下气体分子的速率分布，先将分子速率范围分成许多相等的速率区间，然后通过试验或理论推导找出分布在各个速率区间内的分子数与总分子数的比率。这些比率便给出了分子的速率分布。下表给出了时空气分子的速率分布: 速率区间 分子数比率 速率区间 分子数比率 <100 1.4 400
13、~500 20.5 100~200 8.4 500~600 15.1 200~300 16.2 600~700 9.2 300~400 21.5 >700 7.7 由表可知，300~500的分子数占总数的比率最大。其余的分子数占总分子数的比率都比较小。为了精确地描述分子速率分布，应将速率区间取得足够小，使。这时可将表示为微分，以表示分布在的分子数，比率是速率的函数，而且可以认为与成正比，因而可表示为 速率分布函数物理意义：速率在v附近单位速率区间内的分子数占总分子数的百分比，或
14、者说为某一分子的速率在v附近单位速率区间内的概率。 分布在有限速率区间内的分子数为 = 分布在整个速率区间的分子数显然为分子总数，所以 归一化条件 3.2 理想分子气体的麦克斯韦速率分布律 麦克斯韦速率分布律指出，在平衡状态下，理想气体分子速率分布在区间内的分子数占总分子数的百分比为 T为气体的热力学温度，m为气体分子的质量，k为玻尔兹曼常量。可得麦克斯韦速率分布函数为 第4章
15、MATLAB仿真实验及结果讨论 先把麦克斯韦速率分布律列成一个子程序，以便经常调用，并把一些常用的常数也放在其中，这样主程序就简单了许多。在主程序中，定义变量的数值（如：T,mu和v等的数值大小或取值范围），然后调用maxwell（T,mu,v）函数进行数值计算，并利用plot（v,y）函数将计算结果绘制成图形。如此三遍，便可以绘制出如上所示的三条不同条件下的麦克斯韦速率分布律曲线。 由图可以知道，速率分布曲线的形状与气体温度T和分子质量u有关。温度升高时，分子热运动加剧，即速率较大的分子数占总分子数的百分率增大。另外，分布曲线还需满足归一化条件，所以温度升高时，分布曲线向右移动，的极大值
16、减小，曲线变得较为平坦。 M文件： function f=maxwell(T,mu,v) %mu--分子量，千克.摩尔~1 %v--分子速度 %T--气体的绝对温度 R=8.31; %气体常数 k=1.381*10^(-23); %玻尔兹曼常数 NA=6.022*10^23; %阿伏伽德罗常数 m=mu/NA; %分子质量 f=4*pi*(m/(2*pi*k*T))^(3/2).*exp(-m*v.^2./(2*k*T)).*v.*v;
17、 %麦克斯韦分布律 [*]程序主体： clear T=200;mu=28e-3; %给出T,mv v=eps:1500; %给出自变量数组 y=maxwell(T,mu,v); %调用M函数文件 plot(v,y),hold on %画出分布曲线 v1=400:600; %给定速度范围 y1=maxwell(T,mu,v1); %该范围的分布 fill([v1,600,400],[y1,0,0]
18、'b') trapz(y1) %求该范围概率积分 T=300;mu=28e-3;y=maxwell(T,mu,v);plot(v,y) %改变T，画曲线 T=200;mu=2e-3;y=maxwell(T,mu,v);plot(v,y) %改变mu，画曲线 gtext('T=300,\mu=28*10^{-3}'), %此语句可以生成希 gtext('T=200,\mu=28*10^{-3}'), 腊字母及上标（指
19、 gtext('T=200,\mu=2*10^{-3}'), 数） hold off ans =0.3326 麦克斯韦速率分布曲线 （f(v)为纵轴，v为横轴） 第5章 结论 由上可以看出，MATLAB软件在分子物理学和热学方面有很大用处，它可以方便、快捷、清晰形象地解决分子物理学和热学中的一些复杂问题，将数值计算和绘图这类问题变成数字化，只需更改几个数字就可以得到其他的结果。在实际应用中，更多复杂的问题都可以用MATLAB来进行方便快捷
20、的解决，根据实际情况，灵活应用MATLAB。 参考文献 1] 刘勍，温志贤.MATLAB基础及应用［M］南京：东南大学出版社，2011. [2] 王少杰，顾牡.新编基础物理学［M］北京：科学出版社，2009. [3] 孙冻梅.MATLAB在力学热学教学中的应用［J］《教育教学论坛》，2013年06期 [4] 刘卫国.MATLAB程序设计与应用（第二版）［M］北京：高等教育出版社，2006. [5] 张世功.MATLAB用于物理实验常见的数据处理［J］大学物理实验,2007,20（4）85—86.
21、[6] 陈扬.MATLAB6.x图形编程与图像处理.西安：西安电子科技大学出版社，2002. [7] 吴利华.MATLAB在大学物理实验教学中的应用［J］实验科学与技术,2009,7(2):101- 103. [8] 蔡旭辉等.MATLAB 基础与应用教程[M].北京:人民邮电出版社,2011 [9] 王沫然. MATLAB与科学计算［M］.北京:电子工业出版社，2006. 东北石油大学课程设计成绩评价表 课程名称 计算物理 题目名称 两点电荷间电场线的MATLAB模拟 学生姓名 王兆祥 学号 090901340106 指导教师姓名 职称 序号 评价
22、项目 指 标 满分 评分 1 工作量、工作态度和出勤率 按期圆满的完成了规定的任务，难易程度和工作量符合教学要求，工作努力，遵守纪律，出勤率高，工作作风严谨，善于与他人合作。 20 2 课程设计质量 课程设计选题合理，计算过程简练准确，分析问题思路清晰，结构严谨，文理通顺，撰写规范，图表完备正确。 45 3 创新 工作中有创新意识，对前人工作有一些改进或有一定应用价值。 5 4 答辩 能正确回答指导教师所提出的问题。 30 总分 评语： 指导教师： 年 月 日
©2010-2026 宁波自信网络信息技术有限公司 版权所有
客服电话：0574-28810668 投诉电话：18658249818
浙ICP备2021020529号-1 | 浙B2-20240490
关注我们 ：
