基于matlab分子物理学和热学分析与设计 结业论文

基于matlab分子物理学和热学分析与设计 结业论文

ID:324209

大小:120.00 KB

页数:5页

时间:2017-07-23

基于matlab分子物理学和热学分析与设计  结业论文_第1页
基于matlab分子物理学和热学分析与设计  结业论文_第2页
基于matlab分子物理学和热学分析与设计  结业论文_第3页
基于matlab分子物理学和热学分析与设计  结业论文_第4页
基于matlab分子物理学和热学分析与设计  结业论文_第5页
资源描述:

《基于matlab分子物理学和热学分析与设计 结业论文》由会员上传分享,免费在线阅读,更多相关内容在学术论文-天天文库

1、基于MATLAB分子物理学和热学分析与设计摘要:MATLAB是一种以矩阵运算为基础的交互式程序语言,是专门针对科学和工程中计算和绘图的需求开发的。与其他计算机语言相比,具有简洁和智能化的特点,人机性能好。本文利用MATLAB的数值计算、可视化图形处理和程序设计等特点,设计了一个有关于大学物理中分子物理学和热学的应用,旨在将复杂的麦克斯韦速度分布律的数值计算和绘图简单化。关键词:分子物理学;热学;理想气体;麦克斯韦速率分布律;分析设计;MATLAB引言在分子物理学和热学中可知,由大量分子组成的气体,因分子间的频繁

2、碰撞,各个分子的速度大小和方向瞬息万变。任一时刻,某个分子具有多大的运动速率完全是偶然的,可以是零到无穷大之间的任何值。1859年,英国物理学家麦克斯韦(J.C.Maxwell)从理论上导出了气体分子速率分布律——麦克斯韦速率分布律。利用MATLAB软件灵活的数值计算与符号计算、简单的语句表达、简洁完善的图形绘制、丰富的工具箱函数和简易的扩展功能,将复杂的数学公式绘制成曲线,并研究单个参数的影响。1、麦克斯韦速率分布律[*]2.1速率分布和分布函数为了描述平衡态下气体分子的速率分布,先将分子速率范围分成许多相等

3、的速率区间,然后通过试验或理论推导找出分布在各个速率区间内的分子数与总分子数的比率。这些比率便给出了分子的速率分布。下表给出了时空气分子的速率分布:速率区间分子数比率速率区间分子数比率<1001.4400~50020.5100~2008.4500~60015.1200~30016.2600~7009.2300~40021.5>7007.7由表可知,300~500的分子数占总数的比率最大。其余的分子数占总分子数的比率都比较小。为了精确地描述分子速率分布,应将速率区间取得足够小,使。这时可将表示为微分,以表示分布在

4、的分子数,比率是速率的函数,而且可以认为与成正比,因而可表示为速率分布函数物理意义:速率在v附近单位速率区间内的分子数占总分子数的百分比,或者说为某一分子的速率在v附近单位速率区间内的概率。分布在有限速率区间内的分子数为=分布在整个速率区间的分子数显然为分子总数,所以归一化条件2.2理想气体分子的麦克斯韦速率分布律麦克斯韦速率分布律指出,在平衡状态下,理想气体分子速率分布在区间内的分子数占总分子数的百分比为T为气体的热力学温度,u为气体分子的质量,k为玻尔兹曼常量。可得麦克斯韦速率分布函数为麦克斯韦速率分布曲线

5、(f(v)为纵轴,v为横轴)1、MATLAB仿真与分析先把麦克斯韦速率分布律列成一个子程序,以便经常调用,并把一些常用的常数也放在其中,这样主程序就简单了许多。在主程序中,定义变量的数值(如:T,mu和v等的数值大小或取值范围),然后调用maxwell(T,mu,v)函数进行数值计算,并利用plot(v,y)函数将计算结果绘制成图形。如此三遍,便可以绘制出如上所示的三条不同条件下的麦克斯韦速率分布律曲线。由图可以知道,速率分布曲线的形状与气体温度T和分子质量u有关。温度升高时,分子热运动加剧,即速率较大的分子数

6、占总分子数的百分率增大。另外,分布曲线还需满足归一化条件,所以温度升高时,分布曲线向右移动,的极大值减小,曲线变得较为平坦。M文件:functionf=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;%麦克

7、斯韦分布律[*]程序主体:clearT=200;mu=28e-3;%给出T,mvv=eps:1500;%给出自变量数组y=maxwell(T,mu,v);%调用M函数文件plot(v,y),holdon%画出分布曲线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;

8、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}'),数)holdoffans=0.33261、结论由上可以看出,MATLAB软件在分子物理学

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。