M02初等模型量纲分析和无量纲化
- 格式:pdf
- 大小:439.71 KB
- 文档页数:34
无量纲化法公式无量纲化法公式是一种在科学和工程领域中广泛应用的工具,它能够将具有不同量纲的物理量转化为无量纲的形式,从而方便我们进行分析和比较。
咱先来说说啥叫无量纲化。
比如说,有两个物理量,一个是长度,单位是米;另一个是时间,单位是秒。
这俩家伙的量纲完全不同,直接比较或者运算就会很麻烦。
但通过无量纲化,就可以把它们变成能够在同一尺度上进行讨论和处理的形式。
那无量纲化法公式具体是咋操作的呢?常见的方法有很多种,像标准化、归一化等等。
标准化就是把数据减去平均值再除以标准差,这样得到的结果均值为 0 ,标准差为 1 。
归一化呢,就是把数据映射到 0 到 1 的区间内。
我记得有一次在给学生们讲这个知识点的时候,有个小家伙特别有意思。
我刚在黑板上写下无量纲化法公式,他就瞪大眼睛问我:“老师,这一堆符号看着好复杂,到底有啥用啊?”我笑着跟他说:“别着急,等会儿你就知道它的厉害了。
”然后我就开始举例子,假设我们要比较不同汽车的速度和油耗。
速度的单位是千米每小时,油耗的单位是升每百公里。
这两个量要是直接比,那根本没法比。
但是用无量纲化法公式处理一下,就能很清楚地看出哪辆车更经济实惠。
那孩子听完,恍然大悟地点点头,那表情可认真了。
无量纲化法公式在实际应用中可太有用啦!比如在流体力学中,雷诺数就是一个无量纲量,它能帮助我们判断流体的流动状态是层流还是湍流。
在传热学里,努塞尔数能告诉我们热传递的效率。
再比如说,在研究经济数据的时候,不同地区的 GDP 数值差异巨大,直接比较没有太大意义。
但通过无量纲化处理,就能更客观地比较不同地区经济发展的相对水平。
总之,无量纲化法公式虽然看起来有点复杂,但一旦掌握了它,就能在处理各种数据和物理量的时候更加得心应手,让我们能更清晰地看到事物的本质和规律。
希望通过我的这番讲解,能让您对无量纲化法公式有了更清楚的认识。
不管是在学习还是工作中,遇到需要处理不同量纲的数据时,可别忘了这个神奇的工具哟!。
无量纲化处理方法
无量纲化处理方法是指将不同单位或量纲的数据转化为无单位的纯数值,使得不同量级的数据可以进行比较和统一处理。
常用的无量纲化处理方法有:
1. 最大最小归一化:将数据按照最大值和最小值进行线性变换,使得数据的取值范围在0到1之间。
公式为:
$$X_{new} = \frac{X-X_{min}}{X_{max}-X_{min}}$$
这种方法适用于对数据的绝对值范围不关心,只关心数据在
特定区间内分布情况的情况。
2. 标准化:将数据按照均值和标准差进行线性变换,使得数据的均值为0,标准差为1。
公式为:
$$X_{new} = \frac{X-\mu}{\sigma}$$
这种方法适用于数据的分布符合高斯分布的情况。
3. 小数定标规范化:将数据除以一个固定的基数,通常选择
10的某个次幂,使得数据的绝对值都小于1。
公式为:
$$X_{new} = \frac{X}{10^m}$$
其中,m取决于数据集中的最大绝对值。
4. 非线性变换:通过某种函数对数据进行变换,将其转化为无量纲的纯数值。
常见的非线性变换方法有对数变换、指数变换等。
这种方法适用于数据分布存在偏态或不符合线性关系的情况。
无量纲化处理方法的选择要根据具体的数据特点和所需的分析
目的来确定,合适的无量纲化方法可以提升数据处理和分析的效果。
§5 量 纲 分 析 法 建 模量纲分析(Dimensional Analysis)是20世纪初提出的在物理领域中建立数学模型的一种方法,它在经验和实验的基础上利用物理定律的量纲齐次原则,确定各物理量之间的关系.本节在一个例子的引导下先介绍量纲齐次原则和著名的BuckinghamPi 定理,然后用这个定理讨论一个力学问题的建模方法,并介绍量纲分析在物理模拟中的应用.最后给出一种简化模型的方法——无量纲化.一、量纲齐次原则许多物理量是有量纲的,有些物理量的量纲是基本的,另一些物理量的量纲则可 以由基本量纲根据其定义或某些物理定律推导出来.例如在研究动力学问题时常把长度l 、质量m 和时间t 的量纲作为基本量纲,记以相应的大写字母L ,M 和T .于 是速度v 、加速度a 的量纲可以按照其定义分别用1-LT 和2-LT表示,力f 的量纲则应根据牛顿第二定律用质量和加速度量纲的乘积2-LMT 表示.有些物理常数也有量纲,如万有引力定律221r m m k f =中的引力常数k ,由 221m m fr k =可知其量纲应从力f 、距离r 和质量m 的量纲求出,为2-LMT ·2L ·2-M =213--T M L .通常,一个物理量q 的量纲记作[q],于是上述各物理量的量纲为[l]=L ,[m]=M ,[t]=T ,[v]=LT -1,[a ]=LT -2,[f] =LMT -2,[k]= 213--T M L .对于无量纲量α,我们记[α]=1(因为可视为[α]=000T M L ).用数学公式表示一个物理定律时,等号两端必须保持量纲的一致,或称量纲齐次性(Dimensional Homogeneity).量纲分析就是利用量纲齐次原则来寻求理量之间的关系[6,20].在叙述主要定理之前先看一个例子.单摆运动 这是一个熟知的物理现象,质量为m 的小球系在长度为l 的线的一端,稍偏离平衡位置后小球在重力mg 作用下(g 为重力加速度)做往复摆动,忽略阻力.求摆动周期t 的表达式.在这个问题中出现的物理量有t ,m ,l ,g ,设它们之间有关系式其中1α,2α,3α是待定常数,λ是无量纲的比例系数.取(1)式的量纲表达式即[][][][]321αααg l m t =将[t]=T ,[m]=M ,[l]=L ,[g]=LT -2代入得按照量纲齐次原则应有(3)的解为1α=0,2α=1/2,3α=-1/2,代人(1)式得g l t λ= (4) (4)式与用力学规律得到的结果是一致的.为了导出量纲分析建模的一般方法,将这个例子中各个变量之间的关系写作进而假设(5)式形如 π=4321y y y y g l m t (6)其中1y ~4y 是待定常数,π是无量纲常数.将t ,m ,l ,g 的量纲用基本量纲L ,M ,T表示为100][T M L t =,010][T M L m =,001][T M L l =,201][-=T M L l ,则(6)的量纲表达式可写作(注意到000][T M L =π)即 000241243T M L T M L y y y y y =-+ (7)此方程组有一个基本解T T y y y y y )1,1,0,2(),,,(4321-== (9)代回(6)式得 π=-g l t 12 (10)而(5)式等价于0)(=πF (11)(10),(11)两式就是用量纲齐次原则从(5)式得到的结果.前面给出的(4)式只是它的特殊表达形式.把从(5)式到(11)式的推导过程一般化,就是著名的Pi 定理.定理 设有m 个物理量m q q q ,,,21 ,是与量纲单位的选取无关的物理定律*,n X X X ,,,21 是基本量纲,n ≤m . m q q q ,,,21 的量纲可表为m j X q n i ai i ij ,...,2,1,][1==∏= (13)矩阵m n ij a A ⨯=}{称量纲矩阵.若A 的秩r RankA = (14)设线性齐次方程组(y 是m 维向量) 0=Ay (15)的m-r 个基本解为r m s y y y y T sm s s s -==,,2,1,),,,(21 (16)则∏==m j y j i sj q1π为m-r 个相互独立的无量纲量.且与(12)式等价.F 表示一个未定的函数关系.[航船的阻力] 长l 、吃水深度h 的船以速度v 航行,若不考虑风的影响,那么航船受的阻力f除依(8)赖于船的诸变量l ,h ,v 以外,还与水的参数——密度ρ、粘性系数μ,以及重力加速度g 有关.下面用量纲分析方法确定阻力f 和这些物理量之间的关系.我们按照Pi 定理中(12)~(18)式的步骤进行.1.航船问题中涉及的物理量有:阻力f ,船长l ,吃水深度h ,速度v ,水的密度ρ,水的粘性系数μ,重力加速度g ,要寻求的关系式记作2.这是一个力学问题,基本量纲选为L ,M ,T .上述各物理量的量纲表为其中μ的量纲由基本关系xv p ∂∂=μ得到.其中p 是压强(单位面积受的力),所以2][-=LMT p 212---=⋅MT L L ;v 是流速,x 是尺度,所以111---=⋅=⎥⎦⎤⎢⎣⎡∂∂T L LT x v . 并且有n=3<m=7.3.由(20)立即可写出量纲矩阵并且计算 )(3r RankA == (22)4.解齐次方程0=Ay (23)方程(23)有m-r=7—3=4个基本解,可取为5.(24)式给出4个相互独立的无量纲量而(19)式与 等价,Φ是未定的函数,(25)、(26)两式表达了航船问题中各物理量间的全部关系. 6,为得到阻力f 的显示表达式,由(25)及(26)中4π的式子可写出其中ψ表示一个未定函数.在流体力学中无量纲量)(lg 2/12-=πv称Froude 数,)(3πμρ=lv 称Reynold 数(雷诺数),分别记作μρlv v Fr ==Re ,lg (28) 则(27)式又表示为 Re),,1(22Fr hl f ρψυ= (29)这就是用量纲分析方法确定的航船阻力与各物理量之间的关系,这个结果用通常的机理分析是难以得到的.虽然这里函数ψ的形式无从知道,但是在下面将会看到这个表达式在物理模拟中的用途.评注 从上面的例子可以看出,量纲分析方法在建立物理问题的数学模型中能够得到一些重要的、有用的结果,但是也有较大的局限性.在应用和评价这个方法时以下几点值得注意.1.正确确定各物理量 面对一个实际问题将哪些物理量包括在量纲分析的基本关系式f(·)=0中,对所得结果的合理性是至关重要的.对于航船问题,如果在(19)式中忽略了水的密度ρ或粘性系数μ,则得到的结果就会不同.各物理量的确定主要靠经验和物理知识,无法绝对保证所得结果是正确或有用的.2.合理选取基本量纲 基本量纲选少了,无法表示各物理量,当然不行;选多了也会使问题复杂化.在一般情况下力学问题选取L ,M ,T 即可,热学问题加上温度量纲Θ,电学问题加上电量量纲Q .3.恰当构造基本解 线性齐次方程组的基本解可以有许多不同的构造方法,虽然基本解组能够相互线性表出,但是为了特定的建模目的恰当地构造基本解,能够更直接地得到我们所期望的结果.4.结果的效用和局限性 量纲齐次原则和n 定理是具有普遍意义的又是相当初 等的方法,它不需要非常专门的物理知识及高等的数学方法,就可以得到用其他方法 难以得到的结果,如(29)式.一般地说,从未知定律f(m q q q ,,,21 )=0到用量纲分析方法得到的等价形式F(r m -πππ,,,21 )=0,不仅物理量个数减少了r 个,而且原始物理量m q q q ,,,21 ,组合成了一些有用的无量纲量r m -πππ,,,21 ,下面将进一步讨论它们的用途.另一方面,用这个方法得到的结果是有局限的,“不彻底”的.F(·)=0中仍然包含着一些未定函数和常数 (无量纲量),诸如物理定律中经常 出现的三角函数sin(·)、指数函数exp(·)不可能用量纲分析法得到,因为这些函 数的自变量和函数值都是无量纲的.二、量纲分析的应用——物理模拟中的比例模型我们在1.1节曾介绍过物理模型,它是在实验室条件下按照缩小了的比例尺寸构造的,目的是根据相应的比例来研究原型的某些性质.量纲分析的结果可以指导这种比例关系的确定.以本节提到的单摆运动为例.已经得到模型中摆动周期t 与摆长l 的关系为若记原型中相应的各个物理量为t ',l ',g ',因为λ是无量纲量,在模型与原型中不变,又显然有g=g ,,所以由(30)式立即得到这样,如果模型摆的尺寸按照摆长比例l: l ' =1:4设计制造,那么测定了模型摆的周期t 以后,就可以知道原型摆的周期为t '=2t .可以看出,这里主要用了无量纲量在模型和原型中保持不变的性质.下面利用航船问题的结果讨论怎样构造航船模型,以确定原型航船在海洋中受的阻力,并且当速度不大时可以忽略雷诺数Re 的影响.以g v h l f ,,,,,ρ和g v h l f '''''',,,,,ρ,分别记模型和原型中的各物理量,由(28)、(29)式(略去Re)得注意(32),(33)两式中的函数ψ是一样的.当无量纲量成立时,由(32)、(33)式可得只要模型船和原型船的形状相似,就可以保证(34)的第1式成立.而注意到g=g ',(34)的第2式给出如果在模拟中用与海水有相同密度的水,即ρρ'=,则由(35),(36)式可得于是确定了模型船和原型船的比例l l ':,并测得了模型船的阻力f 后,就能够确定原型航船的阻力f 了.三、无量纲化我们不拟对无量纲化方法作一般阐述,而是通过一个例子介绍这种方法如何用来对模型进行简化.抛射问题 在星球表面以初速v 竖直向上发射火箭*,记星球半径为r ,星球表面重力加速度为g ,忽略阻力,讨论发射高度x 随时间t 的变化规律.设J 轴竖直向上,在发射时刻f=0火箭高度x=O(星球表面).火箭和星球的质量分别记作1m 和2m ,则由牛顿第二定律和万有引力定律可得以x=O 时x=-g 代入(38)式,并注意到初始条件,抛射问题满足如下方程(39)的解可以表示为即发射高度x 是以r ,v ,g 为参数的时间f 的函数.这里的目的不是研究这个函数的具体形式(虽然可以通过求解方程(39)直接得到),而是讨论用无量纲化方法简化它的途径.(40)式包含3个独立参数r ,v ,g ,由(40)式得到的进一步的结果,如火箭到达最高点的时间0==x M t t 。
无量纲化处理的方法我折腾了好久无量纲化处理的方法,总算找到点门道。
说实话,这事我一开始也是瞎摸索。
最初,我试过极差化法。
就是把数据中的最大值和最小值找出来,然后用每个数据减去最小值再除以最大值与最小值的差。
我当时处理一批数据嘛,想着这个看起来好像挺简单的。
但是我犯了个错,我没先检查数据里有没有异常值。
结果后来才发现有些数据偏差特别大是因为里面混进了错误数据,这就导致我的无量纲化结果不太靠谱。
那对于没有错误数据的正常情况呢,极差化法确实对一些数值分布比较均匀的数据处理起来比较顺手。
就好比是把一群个头高低不同的人,让最矮的站在同一个平台上(减最小值),然后按照最高和最矮的人的身高差来调整比例(除以极差),这样就都在一个相对可比的范围内了。
然后我还试过标准化方法,就是把数据减去均值再除以标准差。
这个标准差还算起来有点复杂呢。
我记得那次我处理的数据量特别大,计算均值和标准差的时候差点把我搞晕了。
算的时候吧,就有点儿像从一群羊里找出平均体重,然后再看每只羊相对这个平均体重是胖了还是瘦了多少的感觉。
这里标准差啊,就像是这群羊体重的波动情况。
这个方法对于那些分布呈正态的数据效果特别好。
还有一种归一化方法,直接把数据除以数据的总和。
我在做一个统计每个城市不同行业占经济比重的研究时用了这个。
就是假如一个城市所有行业的经济总和看成是一个蛋糕,每个行业就是这个蛋糕的一块,我们把每块蛋糕的大小除以整个蛋糕的大小,就得到了无量纲化后的相对比重。
不过这个方法数据要是有负数就不太好使了,这也是我试出来的一个经验教训。
我也听说还有其他方法,但目前我还没试过。
不过就我试过的这几种方法,我觉得一定要根据数据的特征来选择合适的无量纲化方法。
要先看看数据有没有异常,数值是怎么分布的,如果数值基本上是正态分布呢,标准化方法是个不错的选择;要是想简单处理且不存在负数的情况下,归一化方法能用;如果数据里没异常且比较均匀,极差化法也行。
反正我感觉每次做无量纲化处理的时候,都得多想想数据的情况,这样才能选到合适的方法。
无量纲化标准化处理在数据分析和建模过程中,经常会遇到不同量纲的数据,这给模型的训练和预测带来了困难。
为了解决这个问题,我们需要对数据进行无量纲化标准化处理。
本文将介绍无量纲化和标准化的概念、方法和应用,帮助读者更好地理解和运用这些技术。
无量纲化是指将数据转换为没有特定单位的形式,通常是将数据缩放到一个特定的范围。
常见的无量纲化方法包括最小-最大缩放和Z-score标准化。
最小-最大缩放将数据缩放到一个特定的范围,通常是[0, 1]或[-1, 1],公式为:\[X_{norm} = \frac{X X_{min}}{X_{max} X_{min}}\]其中,\(X_{min}\)和\(X_{max}\)分别是数据的最小值和最大值。
Z-score标准化将数据缩放成均值为0,标准差为1的分布,公式为:\[X_{std} = \frac{X \mu}{\sigma}\]其中,\(\mu\)和\(\sigma\)分别是数据的均值和标准差。
标准化是指将数据转换为均值为0,标准差为1的分布。
标准化可以使数据更易于比较和分析,同时有助于提高模型的训练和预测性能。
常见的标准化方法包括Z-score标准化和小数定标标准化。
小数定标标准化将数据除以一个固定的数值,使得数据的绝对值都在0到1之间,公式为:\[X_{std} = \frac{X}{10^k}\]其中,\(k\)是使得数据的绝对值都在0到1之间的最小整数。
无量纲化和标准化的应用非常广泛,例如在聚类分析、主成分分析、回归分析和神经网络等领域都有重要的作用。
在聚类分析中,无量纲化和标准化可以使得不同特征对聚类结果的影响更加均衡。
在主成分分析中,无量纲化和标准化可以使得各个特征对主成分的贡献更加平等。
在回归分析和神经网络中,无量纲化和标准化可以加快模型的收敛速度,提高模型的训练和预测性能。
总之,无量纲化和标准化是数据分析和建模过程中非常重要的一环。
通过合适的无量纲化和标准化方法,可以使得数据更易于比较和分析,提高模型的训练和预测性能。
《数学建模(一)》课程教学大纲【课程基本情况】一、课程代码:000373二、课程类别及性质:公共选修课三、课程学时学分:54学时(教学:24 实践:30)2学分四、教学对象:12、13级学生五、课程教材:《数学模型》、姜启源谢金星叶俊等、高等教育出版社六、开设系(部):信科系七、先修课:高等数学、线性代数【教学目的】通过本课程的学习,使学生能够较好地理解数学模型、数学建模的含义,了解数学建模的重要性。
通过示例的学习使同学们基本掌握建立数学模型的方法和步骤,并能通过数学方法、数学软件求解模型,而且能够对模型的精准性进行分析。
通过学习,培养了同学们的把实际问题表述成数学问题的能力,从而提高了他们的抽象思维能力。
并且通过MATLAB、LINGO 数学软件的应用,提高了他们的计算机应用水平。
【教学内容、基本要求及学时分配】第一章建立数学模型教学时数:2学时第一节从现实对象到数学模型基本要求:掌握数学模型、数学建模的含义。
第二节数学建模的重要意义基本要求:了解数学建模的重要性。
第三节数学建模的示例(不讲授)基本要求:掌握三个示例的建模过程;重点:模型的建立、模型的求解。
第四节数学建模的基本方法和步骤基本要求:掌握数学建模的基本方法和步骤;重点:建模的基本方法和步骤。
第五节数学模型的特点和分类基本要求:了解数学模型的特点和分类。
第六节数学建模能力的培养(不讲授)基本要求:了解建立数学模型所需要的能力。
第二章初等模型教学时数:4学时第一节公平的席位分配基本要求:掌握公平席位的建模方法;重点:建立数量指标。
第二节录像机计数器的用途基本要求:掌握录像机计数器的建模方法;重点:模型的假设及模型的构成。
难点:建立模型的过程。
第三节双层玻璃的功效基本要求:掌握双层玻璃的功效的建模方法及模型应用;重点:模型的构成。
第四节汽车刹车距离基本要求:掌握t秒准则的建立方法;重点:模型建立的过程。
第五节划艇比赛的成绩(不讲授)第六节动物的身长和体重(不讲授)第七节实物交换(不讲授)第八节核军备竞赛(不讲授)第九节扬帆远航(不讲授)第十节量纲分析与无量纲化(不讲授)第三章简单的优化模型教学时数:4学时第一节存贮模型基本要求:掌握存贮模型在两种情况下的建模方法;重点:模型假设。
数据的无量纲化处理数据的无量纲化处理是指将具有不同量纲的数据转化为统一的标准化数据,以消除不同量纲对数据分析和建模的影响。
无量纲化处理可以有效地提高数据的可比性和可解释性,使得数据分析和建模更加准确和可靠。
本文将介绍数据的无量纲化处理的几种常见方法,并结合实例进行详细说明。
1. 标准化(Standardization)标准化是将数据转化为均值为0,标准差为1的分布,常用的标准化方法有Z-Score标准化和MinMax标准化。
Z-Score标准化公式如下:$$Z = \frac{X - \mu}{\sigma}$$其中,Z为标准化后的值,X为原始值,μ为原始数据的均值,σ为原始数据的标准差。
MinMax标准化公式如下:$$X_{\text{new}} = \frac{X - X_{\text{min}}}{X_{\text{max}} - X_{\text{min}}} $$其中,X为原始值,X_new为标准化后的值,X_min为原始数据的最小值,X_max为原始数据的最大值。
2. 归一化(Normalization)归一化是将数据缩放到0和1之间的范围,常用的归一化方法有Min-Max归一化和Decimal Scaling归一化。
Min-Max归一化公式如下:$$X_{\text{new}} = \frac{X - X_{\text{min}}}{X_{\text{max}} - X_{\text{min}}} $$其中,X为原始值,X_new为归一化后的值,X_min为原始数据的最小值,X_max为原始数据的最大值。
Decimal Scaling归一化公式如下:$$X_{\text{new}} = \frac{X}{10^j}$$其中,X为原始值,X_new为归一化后的值,j为使得归一化后的值在0和1之间的最小整数。
3. 对数转换(Log Transformation)对数转换是将数据进行对数运算,常用的对数转换方法有自然对数转换和底数为10的对数转换。
§5 量 纲 分 析 法 建 模量纲分析(Dimensional Analysis)是20世纪初提出的在物理领域中建立数学模型的一种方法,它在经验和实验的基础上利用物理定律的量纲齐次原则,确定各物理量之间的关系.本节在一个例子的引导下先介绍量纲齐次原则和著名的BuckinghamPi 定理,然后用这个定理讨论一个力学问题的建模方法,并介绍量纲分析在物理模拟中的应用.最后给出一种简化模型的方法——无量纲化.一、量纲齐次原则许多物理量是有量纲的,有些物理量的量纲是基本的,另一些物理量的量纲则可 以由基本量纲根据其定义或某些物理定律推导出来.例如在研究动力学问题时常把长度l 、质量m 和时间t 的量纲作为基本量纲,记以相应的大写字母L ,M 和T .于 是速度v 、加速度a 的量纲可以按照其定义分别用1-LT 和2-LT表示,力f 的量纲则应根据牛顿第二定律用质量和加速度量纲的乘积2-LMT 表示.有些物理常数也有量纲,如万有引力定律221r m m k f =中的引力常数k ,由 221m m fr k =可知其量纲应从力f 、距离r 和质量m 的量纲求出,为2-LMT ·2L ·2-M =213--T M L .通常,一个物理量q 的量纲记作[q],于是上述各物理量的量纲为[l]=L ,[m]=M ,[t]=T ,[v]=LT -1,[a ]=LT -2,[f] =LMT -2,[k]= 213--T M L .对于无量纲量α,我们记[α]=1(因为可视为[α]=000T M L ).用数学公式表示一个物理定律时,等号两端必须保持量纲的一致,或称量纲齐次性(Dimensional Homogeneity).量纲分析就是利用量纲齐次原则来寻求理量之间的关系[6,20].在叙述主要定理之前先看一个例子.单摆运动 这是一个熟知的物理现象,质量为m 的小球系在长度为l 的线的一端,稍偏离平衡位置后小球在重力mg 作用下(g 为重力加速度)做往复摆动,忽略阻力.求摆动周期t 的表达式.在这个问题中出现的物理量有t ,m ,l ,g ,设它们之间有关系式其中1α,2α,3α是待定常数,λ是无量纲的比例系数.取(1)式的量纲表达式即[][][][]321αααg l m t =将[t]=T ,[m]=M ,[l]=L ,[g]=LT -2代入得按照量纲齐次原则应有(3)的解为1α=0,2α=1/2,3α=-1/2,代人(1)式得g l t λ= (4) (4)式与用力学规律得到的结果是一致的.为了导出量纲分析建模的一般方法,将这个例子中各个变量之间的关系写作进而假设(5)式形如 π=4321y y y y g l m t (6)其中1y ~4y 是待定常数,π是无量纲常数.将t ,m ,l ,g 的量纲用基本量纲L ,M ,T表示为100][T M L t =,010][T M L m =,001][T M L l =,201][-=T M L l ,则(6)的量纲表达式可写作(注意到000][T M L =π)即 000241243T M L T M L y y y y y =-+ (7)此方程组有一个基本解T T y y y y y )1,1,0,2(),,,(4321-== (9)代回(6)式得 π=-g l t 12 (10)而(5)式等价于0)(=πF (11)(10),(11)两式就是用量纲齐次原则从(5)式得到的结果.前面给出的(4)式只是它的特殊表达形式.把从(5)式到(11)式的推导过程一般化,就是著名的Pi 定理.定理 设有m 个物理量m q q q ,,,21 ,是与量纲单位的选取无关的物理定律*,n X X X ,,,21 是基本量纲,n ≤m . m q q q ,,,21 的量纲可表为m j X q n i ai i ij ,...,2,1,][1==∏= (13)矩阵m n ij a A ⨯=}{称量纲矩阵.若A 的秩r RankA = (14)设线性齐次方程组(y 是m 维向量) 0=Ay (15)的m-r 个基本解为r m s y y y y T sm s s s -==,,2,1,),,,(21 (16)则∏==m j y j i sj q1π为m-r 个相互独立的无量纲量.且与(12)式等价.F 表示一个未定的函数关系.[航船的阻力] 长l 、吃水深度h 的船以速度v 航行,若不考虑风的影响,那么航船受的阻力f除依(8)赖于船的诸变量l ,h ,v 以外,还与水的参数——密度ρ、粘性系数μ,以及重力加速度g 有关.下面用量纲分析方法确定阻力f 和这些物理量之间的关系.我们按照Pi 定理中(12)~(18)式的步骤进行.1.航船问题中涉及的物理量有:阻力f ,船长l ,吃水深度h ,速度v ,水的密度ρ,水的粘性系数μ,重力加速度g ,要寻求的关系式记作2.这是一个力学问题,基本量纲选为L ,M ,T .上述各物理量的量纲表为其中μ的量纲由基本关系xv p ∂∂=μ得到.其中p 是压强(单位面积受的力),所以2][-=LMT p 212---=⋅MT L L ;v 是流速,x 是尺度,所以111---=⋅=⎥⎦⎤⎢⎣⎡∂∂T L LT x v . 并且有n=3<m=7.3.由(20)立即可写出量纲矩阵并且计算 )(3r RankA == (22)4.解齐次方程0=Ay (23)方程(23)有m-r=7—3=4个基本解,可取为5.(24)式给出4个相互独立的无量纲量而(19)式与 等价,Φ是未定的函数,(25)、(26)两式表达了航船问题中各物理量间的全部关系. 6,为得到阻力f 的显示表达式,由(25)及(26)中4π的式子可写出其中ψ表示一个未定函数.在流体力学中无量纲量)(lg 2/12-=πv称Froude 数,)(3πμρ=lv 称Reynold 数(雷诺数),分别记作μρlv v Fr ==Re ,lg (28) 则(27)式又表示为 Re),,1(22Fr hl f ρψυ= (29)这就是用量纲分析方法确定的航船阻力与各物理量之间的关系,这个结果用通常的机理分析是难以得到的.虽然这里函数ψ的形式无从知道,但是在下面将会看到这个表达式在物理模拟中的用途.评注 从上面的例子可以看出,量纲分析方法在建立物理问题的数学模型中能够得到一些重要的、有用的结果,但是也有较大的局限性.在应用和评价这个方法时以下几点值得注意.1.正确确定各物理量 面对一个实际问题将哪些物理量包括在量纲分析的基本关系式f(·)=0中,对所得结果的合理性是至关重要的.对于航船问题,如果在(19)式中忽略了水的密度ρ或粘性系数μ,则得到的结果就会不同.各物理量的确定主要靠经验和物理知识,无法绝对保证所得结果是正确或有用的.2.合理选取基本量纲 基本量纲选少了,无法表示各物理量,当然不行;选多了也会使问题复杂化.在一般情况下力学问题选取L ,M ,T 即可,热学问题加上温度量纲Θ,电学问题加上电量量纲Q .3.恰当构造基本解 线性齐次方程组的基本解可以有许多不同的构造方法,虽然基本解组能够相互线性表出,但是为了特定的建模目的恰当地构造基本解,能够更直接地得到我们所期望的结果.4.结果的效用和局限性 量纲齐次原则和n 定理是具有普遍意义的又是相当初 等的方法,它不需要非常专门的物理知识及高等的数学方法,就可以得到用其他方法 难以得到的结果,如(29)式.一般地说,从未知定律f(m q q q ,,,21 )=0到用量纲分析方法得到的等价形式F(r m -πππ,,,21 )=0,不仅物理量个数减少了r 个,而且原始物理量m q q q ,,,21 ,组合成了一些有用的无量纲量r m -πππ,,,21 ,下面将进一步讨论它们的用途.另一方面,用这个方法得到的结果是有局限的,“不彻底”的.F(·)=0中仍然包含着一些未定函数和常数 (无量纲量),诸如物理定律中经常 出现的三角函数sin(·)、指数函数exp(·)不可能用量纲分析法得到,因为这些函 数的自变量和函数值都是无量纲的.二、量纲分析的应用——物理模拟中的比例模型我们在1.1节曾介绍过物理模型,它是在实验室条件下按照缩小了的比例尺寸构造的,目的是根据相应的比例来研究原型的某些性质.量纲分析的结果可以指导这种比例关系的确定.以本节提到的单摆运动为例.已经得到模型中摆动周期t 与摆长l 的关系为若记原型中相应的各个物理量为t ',l ',g ',因为λ是无量纲量,在模型与原型中不变,又显然有g=g ,,所以由(30)式立即得到这样,如果模型摆的尺寸按照摆长比例l: l ' =1:4设计制造,那么测定了模型摆的周期t 以后,就可以知道原型摆的周期为t '=2t .可以看出,这里主要用了无量纲量在模型和原型中保持不变的性质.下面利用航船问题的结果讨论怎样构造航船模型,以确定原型航船在海洋中受的阻力,并且当速度不大时可以忽略雷诺数Re 的影响.以g v h l f ,,,,,ρ和g v h l f '''''',,,,,ρ,分别记模型和原型中的各物理量,由(28)、(29)式(略去Re)得注意(32),(33)两式中的函数ψ是一样的.当无量纲量成立时,由(32)、(33)式可得只要模型船和原型船的形状相似,就可以保证(34)的第1式成立.而注意到g=g ',(34)的第2式给出如果在模拟中用与海水有相同密度的水,即ρρ'=,则由(35),(36)式可得于是确定了模型船和原型船的比例l l ':,并测得了模型船的阻力f 后,就能够确定原型航船的阻力f 了.三、无量纲化我们不拟对无量纲化方法作一般阐述,而是通过一个例子介绍这种方法如何用来对模型进行简化.抛射问题 在星球表面以初速v 竖直向上发射火箭*,记星球半径为r ,星球表面重力加速度为g ,忽略阻力,讨论发射高度x 随时间t 的变化规律.设J 轴竖直向上,在发射时刻f=0火箭高度x=O(星球表面).火箭和星球的质量分别记作1m 和2m ,则由牛顿第二定律和万有引力定律可得以x=O 时x=-g 代入(38)式,并注意到初始条件,抛射问题满足如下方程(39)的解可以表示为即发射高度x 是以r ,v ,g 为参数的时间f 的函数.这里的目的不是研究这个函数的具体形式(虽然可以通过求解方程(39)直接得到),而是讨论用无量纲化方法简化它的途径.(40)式包含3个独立参数r ,v ,g ,由(40)式得到的进一步的结果,如火箭到达最高点的时间0==x M t t 。