稿Matlab在麦克斯韦速率分布律教学中应用.docVIP

MATLAB在麦克斯韦速率分布律教学中的应用 王明美 （合肥师范学院物理与电子工程系，安徽，合肥，230061） 摘要：利用MATLAB软件对麦克斯韦速率分布律的速率分布曲线和不同温度下的速率分布曲线进行了模拟，并对分子速率的三个统计值和一定速率区间的分子百分率的进行了编程计算。 关键词：MATLAB; 麦克斯韦; 速率分布曲线; 速率分布律;分子速率 MATLAB In the teaching of?Maxwell speed distribution law WANG Ming-mei （Department of Physics and Electronic Engineering, Hefei Normal University, Hefei 230061 ,China） Abstract: Using MATLAB to simulate the speed distribution curve of the Maxwell speed Distribution law and the speed distribution curve of different temperature, to program and calculate the three statistics and the molecular percentage rate of certain speed range of molecular speed. Key words: MATLAB; Maxwell; distribution curve; speed distribution law; molecular speed 1 引言 Matlab是当今科学界最具影响力的软件，它起源于矩阵运算，已发展成一种高度集成的计算机语言，可以提供强大的科学运算灵活的程序设计流程，高质量的图形可视化与界面设计功能，在工业研究与开发，理科教学，数值分析和科学计算方面有广泛的应用。[1][4] 2 速率分布曲线 麦克斯韦经过理论研究，指出在平衡状态中气体分子速率分布函数的具体形式是 ，其中的叫做麦克斯韦速率分布函数，表示速率分布函数的曲线叫做麦克斯韦速率分布曲线。[2] 用Matlab模拟氢分子T=273K时的速率分布曲线的程序如下： m0=3.35e-27; T=273;k=1.38e-23; v=0:100:5000; f=4*pi*(v.^2).*(m0/(2*pi*k*T))^1.5.*exp(-m0*(v.^2)/(2*k*T)); %分布函数表达式 plot(v,f) %绘图 xlabel(v/(m/s));ylabel(f(v)); 所得图形如图1所示。 图1 氢分子T=273K时的速率分布曲线 2 速率分布律与温度的关系 用Matlab模拟氮分子在两种温度下的速率分布.已知氦分子的质量，，， 。程序如下： m0=4.648e-25; T1=300;T2=1200; k=1.38e-23; v=0:10:1000; %速率区间取0到1000m/s f1=4*pi*(v.^2).*(m0/(2*pi*k*T1))^1.5.*exp(-m0*(v.^2)/(2*k*T1)); f2=4*pi*(v.^2).*(m0/(2*pi*k*T2))^1.5.*exp(-m0*(v.^2)/(2*k*T2)); plot(v,f1) hold on %保持图形[4] plot(v,f2) %画图 xlabel(v/(m/s));ylabel(f(v)); 所得图形如图2所示[3]。 图2 氮分子在两种温度下的速率分布 3 计算分子速率的三个统计值 气体分子速率的三个统计值分别是分子速率的平均速率、方均根速率和最概然速率。 计算分子的平均速率 计算公式为为方便起见，令，，。编程如下： syms x b pi%建立函数关系 f=4/(sqrt(pi))*b^(-1/2)*exp(-x^2)*x^3;%被积函数表达式 A=int(f,0,inf) % 运行结果 A =2/pi^(1/2)/b^(1/2)，这样从运行结果得 （2）计算分子的方均根速率 分子的方均根速率的计算公式为，编程如下： syms x b pi%建立函数关系 f=4/(sqrt(pi))*b^(-1)*exp(-x^2)*x^4; %被积函数表达式 B=int(f,0,inf) % 运行结果 B =3/2/b，这样从运行结果得 (3)计算最概然速率 最概然速率的求法是，令，求出所有驻点；对每个驻点求二阶导数的值；如果某个驻点有则可以确定在此处处分布函数有极大值[5]。以下编程求解，为方便起见，仍令。 第