MC方法模拟理想气体分子麦克斯韦速率分布律.doc

东北石油大学

课 程 设 计

课 程 计算物理
题 目 MC方法模拟理想气体分子的麦克斯韦速率分布
院 系 电子科学学院
专业班级 应用物理10-2
学生姓名
学生学号
指导教师

2014年 3月 14日

东北石油大学课程设计任务书

课程 计算物理和MATLAB课程设计
题目 MC方法模拟理想气体分子的麦克斯韦速率分布
专业 应物 姓名 学号
主要内容、基本要求、主要参考资料等
主要内容：

麦克斯韦速率分布函数为

通过利用M-C的方法模拟麦克斯韦速率分布图。

基本要求：

1、给出M-C研究麦克斯韦速率分布的方法及步骤。
2、画出M-C的模拟图，且与理论值相比较。

主要参考资料：

[1]马文蔚，解希顺等. 物理学[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章 麦克斯韦速率分布的数学描述 5
3.1 速率分布和分布函数 5
3.2 理想气体分子的麦克斯韦速率分布律 6
第4章 Matlab仿真实验及结果讨论 7
第5章 结论 8
参考文献 9

第1章 概述
引言麦克斯韦速率分布律是分子动理论的重要结论之一,它是研究气体分子碰撞、大量分子热运动服从统计规律性等问题的重要理论依据,正确理解它对学习热学非常有用。但是由于推导困难、公式复杂、曲线难画,麦克斯韦速率分布律学习起来比较费力,成了热学理论教学中的一个难点。Matlab是由美国的MathWorks公司于20世纪80年代中期出品的一款商业数学软件,由于优秀的数值计算能力和卓越的数据可视能力,使其很快在数学软件中大放异彩,与Mathematica、Maple、MathCAD一起并称为四大数学软件。通过Matlab的图形用户界面(Graphical User Interfaces,GUI)可以设计出一个界面友好、方便人机交互的图形界面程序。在物理教学中,适当利用Matlab GUI建立直观的物理图像,有助于学生更好地理解物理概念,提高学习兴趣,调动学习积极性
在分子物理学和热学中可知，由大量分子组成的气体，因分子间的频繁碰撞，各个分子的速度大小和方向瞬息万变。任一时刻，某个分子具有多大的运动速率完全是偶然的，可以是零到无穷大之间的任何值。1859年，英国物理学家麦克斯韦（J.C.Maxwell）从理论上导出了气体分子速率分布律——麦克斯韦速率分布律。
利用MATLAB软件灵活的数值计算与符号计算、简单的语句表达、简洁完善的图形绘制、丰富的工具箱函数和简易的扩展功能，将复杂的数学公式绘制成曲线，并研究单个参数的影响。

第2章 Matlab的基础知识
2.1 Matlab基本命令
表2-1 MATLAB基本命令
主题词
含义
主题词
含义
format
设置数据显示格式
feval
函数求值
who
显示变量名
input
提示输入
whos
显示变量信息
disp
输出
clear
清除内存变量
tic
启动秒表
save
保存工作变量到文件
toc
时间读数（秒）
load
从文件装载变量
help
帮助
linspace
区间等分
lookfor
查找
length
获取数组长度
type
列程序清单
size
矩阵大小
which
查找文件目录
max
最大值
double
双精度
min
最小值
str2num
字符串转化为数值
sum
求和
num2str
数值转化为字符串
find
条件检索

2.2 Matlab与绘图有关的基本命令
表2-2 常用作图命令和函数
主题词
含 义
主题词
含 义
plot
基本二维图形
clabel
等高线高度标志
fplot
一元函数图像
grid
格栅
ezplot
画二维曲线的符号命令
hold
图形保持
plot3
空间曲线
axis
定制坐标轴
meshgrid
网格数据生成
view
改变视点
mesh
网面图
subplot
子图
surf
曲面图
figure
新图形窗口
contour
等高线图
clf
清除图形
contour3
三维等高线图
close
关闭图形窗口
title
标题
ylabel
y轴说明

在线条多于一条时，若用户没有指定使用颜色，则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 = [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]；在上述的几种使用形式中，若有复数出现，则复数的虚数部分将不被考虑。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 命令中我们可以混合使用三参数和二参数的形式。

第3章 麦克斯韦速率分布的数学描述
气体中个别分子的速度具有怎样的数值和方向完全是偶然的，但就大量分子的整体来看，在一定的条件下，气体分子的速度分布也遵从一定的统计规律。这个规律也叫麦克斯韦速率分布律。

3.1 速率分布和分布函数
为了描述平衡态下气体分子的速率分布，先将分子速率范围分成许多相等的速率区间，然后通过试验或理论推导找出分布在各个速率区间内的分子数与总分子数的比率。这些比率便给出了分子的速率分布。下表给出了时空气分子的速率分布:
速率区间
分子数比率
速率区间
分子数比率
<100
1.4
400~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附近单位速率区间内的分子数占总分子数的百分比，或者说为某一分子的速率在v附近单位速率区间内的概率。

分布在有限速率区间内的分子数为
=
分布在整个速率区间的分子数显然为分子总数，所以

归一化条件

3.2 理想分子气体的麦克斯韦速率分布律
麦克斯韦速率分布律指出，在平衡状态下，理想气体分子速率分布在区间内的分子数占总分子数的百分比为

T为气体的热力学温度，m为气体分子的质量，k为玻尔兹曼常量。可得麦克斯韦速率分布函数为

第4章 MATLAB仿真实验及结果讨论
先把麦克斯韦速率分布律列成一个子程序，以便经常调用，并把一些常用的常数也放在其中，这样主程序就简单了许多。在主程序中，定义变量的数值（如：T,mu和v等的数值大小或取值范围），然后调用maxwell（T,mu,v）函数进行数值计算，并利用plot（v,y）函数将计算结果绘制成图形。如此三遍，便可以绘制出如上所示的三条不同条件下的麦克斯韦速率分布律曲线。
由图可以知道，速率分布曲线的形状与气体温度T和分子质量u有关。温度升高时，分子热运动加剧，即速率较大的分子数占总分子数的百分率增大。另外，分布曲线还需满足归一化条件，所以温度升高时，分布曲线向右移动，的极大值减小，曲线变得较为平坦。

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;
%麦克斯韦分布律

[*]程序主体：
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],'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}'), 腊字母及上标（指
gtext('T=200,\mu=2*10^{-3}'), 数）
hold off
ans =0.3326

麦克斯韦速率分布曲线

（f(v)为纵轴，v为横轴）

第5章 结论
由上可以看出，MATLAB软件在分子物理学和热学方面有很大用处，它可以方便、快捷、清晰形象地解决分子物理学和热学中的一些复杂问题，将数值计算和绘图这类问题变成数字化，只需更改几个数字就可以得到其他的结果。在实际应用中，更多复杂的问题都可以用MATLAB来进行方便快捷的解决，根据实际情况，灵活应用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.
[6] 陈扬.MATLAB6.x图形编程与图像处理.西安：西安电子科技大学出版社，2002.
[7] 吴利华.MATLAB在大学物理实验教学中的应用［J］实验科学与技术,2009,7(2):101- 103.
[8] 蔡旭辉等.MATLAB 基础与应用教程[M].北京:人民邮电出版社,2011
[9] 王沫然. MATLAB与科学计算［M］.北京:电子工业出版社，2006.
东北石油大学课程设计成绩评价表

课程名称
计算物理
题目名称
两点电荷间电场线的MATLAB模拟
学生姓名
王兆祥
学号
090901340106
指导教师姓名

职称

序号
评价项目
指 标
满分
评分
1
工作量、工作态度和出勤率
按期圆满的完成了规定的任务，难易程度和工作量符合教学要求，工作努力，遵守纪律，出勤率高，工作作风严谨，善于与他人合作。
20

2
课程设计质量
课程设计选题合理，计算过程简练准确，分析问题思路清晰，结构严谨，文理通顺，撰写规范，图表完备正确。
45

3
创新
工作中有创新意识，对前人工作有一些改进或有一定应用价值。
5

4
答辩
能正确回答指导教师所提出的问题。
30

总分

评语：

指导教师： 年 月 日