MATLAB在熱物理学中的应用.docVIP

MATLAB在热物理学中的应用　thermal physics Abstract: Based on MATLAB’s features of numerical calculation, visualization of graphics processing,opening and scalable architecture,introduced the applications of language of Matlab with high-performance in thermal physics of university physics,include in the Maxwell speed distribution、specific heat at constant volume regression analysis of ideal gas,and the chemical industry thermodynamics. Key Words: MATLAB; Maxwell speed distribution;ideal gas; thermodynamic 引言 大学物理学是工科学生的一门必修课, 由于大学物理数学处理比较复杂,恰当地使用可视化以展现数学公式的物理图像, 使其变得直观、形象。MATLAB 是一套高性能的数值计算和可视化软件,下面我们从一些典型的实例出发, 介绍 MATLAB 在热物理学方面的具体应用。 1 MATLAB-PDEtool介绍 MATLAB-PDEtool提供了一个功能强大使用灵活的二维有限元偏微分方程求解环境，其图形用户界面更是使用十分方便、直观一般来说，MATLAB-PDEtool包括3个步骤： （1） 定义一个PDE问题，它包括确定二维求解区域、边界条件和PDE系数。MATLAB-PDEtool能够求解的PDE型式有：椭圆型、抛物线型、双曲线型、特征值型。当使用GUI时，可以在画图模式下确定求解区域；在边界模式下选择方程形式和设置方程系数。 （2） 数值求解，它包括剖分、离散方程和得到一个数值解。在GUI中，在剖分模式下形成满意的网格；在求解模式下通过选择数值计算方法求解。 （3） 图形化显示结果。通常用于表现有限元计算结果的图形有：变形网格图、云图、等值线图、矢量图、网格图、表面图、流线图等。 2 MATLAB在麦克斯韦速率分布中的应用 气体动力学理论中麦克斯韦速率分布律是大学物理讲授与学习中的一个难点和重点。这是因为公式比较复杂抽象，数学推导证明比较繁琐。如果借助 Matlab 就可以比较方便地解决这些问题[1]。首先, 推导三种速率和归一化条件。已知分布函数表达式为： （1） 最大概然速率分布可由下式求出： （2） 平均速率的定义是： （3） 方均根速率为： （4） 归一化条件是： （5） 后三项求解比较复杂,其中用到Gamma函数,传统方法是查数学用表得到结果。如果应用 Matlab的符号计算功能, 只需要简单几行语句就可以解决这些问题。用到相关函数有：符号变量创建函数 syms，求微分函数diff， 求积分函数 int，符号化简函数 simple，字符串转化函数 eval。 具体 Matlab 语句如下： Syms mktvanp g=exp(- (m*v^2)/(2*k*t))*v^2； a=m/(2*k*t);f=int(exp(- a*v^2)*v^n,v,0,inf)； b=f*4*pi*(a/pi)^(3/2)； pretty(solve(diff(g,v)))； %求解最概然速率 n=2;eval(simple(eval(b)))，%证明归一化条件 n=3;pretty(simple(eval(b)))，%求解平均速率 n=4;pretty(eval(simple(eval(b^0.5))))， %求解方均根速率 程序中第五行语句给出最大概然速率： （6） 其物理意义是“若把整个速率范围分成许多相等的小区间，则 vp 所在的区间的分子数占分子总数的百分比最大”[2]。所以可知在( 0，+∞) 速率区间的分子数占分子总数的百分比(对应着曲线与X轴所围的面积) 恒等于1，即满足归一化条件，第六行语句证明了这个结论。 第七行语句给出平均速率：