基于Mathematica和Matlab的典型数学物理方程解分布图像制作.doc

第 29 卷 第 7 期2012 年 7 月吉林化工学院学报Vol． 29 No． 7JOURNAL OF JILIN INSTITUT 第 29 卷 第 7 期 2012 年 7 月 吉 林 化 工 学 院 学 报 Vol． 29 No． 7 JOURNAL OF JILIN INSTITUTE OF CHEMICAL TECHNOLOGY Jul． 2012 文章编号:1007-2853( 2012) 07-0041-04 基于 和 的典型数学物理 Mathematica Matlab 方程解分布图像制作 陈殿伟1 ，杨海英2 ( 1． 吉林化工学院 理学院，吉林 吉林 132022; 2． 吉林市第五中学 物理组，吉林 吉林 132011) 摘要: 通过计算机辅助软件 Mathematica 和 Matlab，把典型数学物理方程解的空间分布制作成三维图像． 关 键 词: Mathematica; Matlab; 数学物理方程 中图分类号: O 411． 1 文献标志码: A 典型数学物理方程包括: 波动方程; 热传导方 程; 泊松方程，Δu = 0 为拉普拉斯方程，二者都是 稳定分布方程． 这三个方程构成了“数理方程”的 主要内容． 它们都是二阶线性方程，刻画了很多物 理现象的规律［1-3］． 波动方程分布图像 1 波动方程或称波方程是一种重要的偏微分方 程，主要描述自然界中的各种的波动现象． 其表达 式为: 图 1 位移在时间和空间中的分布图像 2 ( 1) utt = a Δu + f ． 热传导方程 2 如弦的 x = 0 端固定，x = l 端受迫作谐振动 F = Asinωt，弦的初始速度为零，求弦的振动． 这个 定解问题是: 当一个物体内部各点的温度不一样时，则热 量就会从温度高的地方向温度低的地方流动，这 2 utt － a uxx = 0 ， u | x = 0 = 0，u | x = l = Asinωt ， u | t = 0 = 0，ut | t = 0 = Aω ． x = l 端为非齐次边界条件． 种现象就是热传导问题． 由于热传导过程总是表 ( 2) 现为温度随时间和位置的变化，所以，解决热传导 问题归结为求物体内部温度的分布问题． 如果研 究物质的扩散问题，物质的扩散是由浓度高的地 这是一个不受外力作用的振动，但它 在 弦 x = l的一端，有一个谐振源 F = Asinωt，而在初始 方扩散到低的地方，且服从傅立叶热传导方程相 类似的能斯特扩散定理，即物质的扩散与浓度的 时刻，速 度 为 Aω． 用 行［4］，得到如下图形: Mathematica 直 接 编 程 运 变化成正比． 与热传导方程类似，可以得到浓度 u 满足热传导方程，所以热传导方程也称扩散方程 其表达式为 图 1 中，a = 2，A = 5，l = 1，ω = 1，x ∈［0，1 ］， t∈［0，10］． 从图中看出，在弦 x = 0 的一端，由于没 有谐振源，因此，位移为 0; 而在 x = l 的一端，谐振 源 F = Asinωt，位移 u( x，t) 随着 F 的作用而变化． = a2 Δu ． u ( 3) t 对于细杆导热问题． 初始时刻杆的一端温度 为零度，另端温度为 u ，杆上温度梯度均匀，零度 0 收稿日期:2012-04-14 作者简介:陈殿伟( 1968-) ，男，吉林省吉林市人，吉林化工学院副教授，主要从事凝聚态物理方面的研究． 的一端保持温度不变，另一端跟外界绝热．杆上温度 u ( x，t) 满足下列泛定方程和定解 条件化，因为杆上的每个领域温度已经达到平衡．稳定分布方程分布图像32 的一端保持温度不变，另一端跟外界绝热． 杆上温度 u ( x，t) 满足下列泛定方程和定解 条件 化，因为杆上的每个领域温度已经达到平衡． 稳定分布方程分布图像 3 2 ut － a uxx = 0 ， u | x = 0 = 0，ux | x = l = 0 ， 在热传导问题中，如果温度分布稳定，热源强 度 f( x，y，z) 不随时间变化，热传导持续下去，最终 将达到稳定状态，空间中各点的温度不再随时间 ( 4) 2 u | t = 0 = u0 x / l ( a = k / cρ) ． 泛定方程和边界条件都是齐次的，用分离变 数法 求 解． 在 这 里 可 以 用 Mathematica 直 接 编 程［4］，就可以画出杆上温度 u ( x，t) 的分布图像， 如图 2 ～ 3 所示． ［3］ 变化 ，即 ut = 0，得到方程 Δu + f( x，y，z) = 0 ， ( 5) 此方程称为泊松方程，如果没有热源，即 f = 0 则得到拉普拉斯方程． 在扩散问题中，浓度处于稳 定的状态考虑