P00 · 第 2 章 · 第一编 物理量与单位

尺度分析、数量级与估算

使用十进数量级、区间上下界、幂律尺度和无量纲比进行可复核估算,以 Fermi 分解组织未知信息、识别主导项,并明确假设、误差边界与量纲方法的限制。

报告页面错误
预备知识物理量、量纲与单位制函数、复合与图像

本章目标

  1. 把正数量写成科学记数法,判断最接近的十进数量级并解释边界约定。
  2. 从 $y=Cx^\alpha$ 读取尺度指数、比值规律和常数的单位。
  3. 区分几何相似下长度、面积、体积及表面积体积比的缩放。
  4. 把 Fermi 问题分解为可估因子,为每个输入声明单位、中心值和上下界。
  5. 对乘积、商和幂律传播保守区间,并识别对答案最敏感的假设。
  6. 使用无量纲比、极限情形和独立算法检查估算,同时说明量纲分析不能给出的信息。
页面阅读位置0% · 仅保存在此浏览器
章节未开始
本册完成进度0/6 章 · 0%
本页目录

为什么要估算,而不是等到所有数据齐全

真实问题很少把所需数字整齐列在题面上。设计实验前要判断传感器量程,读到结果后要识别是否错了三个数量级,建立模型时要决定哪些效应可能主导。此时,与其给出没有依据的许多小数位,不如先构造一个可追溯估算:目标量是什么、依赖哪些因子、每个因子的合理范围是什么、结论对哪些假设最敏感。

估算不是随便猜。它是一条缩短但闭合的推理链:保持单位和量纲一致,用简单模型连接输入与目标,给出中心值和边界,再用另一种观察方式核对。一个写明“约 103kg10^3\,\mathrm{kg},在这些假设下可能相差两倍”的答案,通常比没有误差说明的 1037.428kg1037.428\,\mathrm{kg} 更诚实也更有用。

科学记数法与十进数量级

十进数量级

任意正数可唯一写成

x=a×10n,1a<10,x=a\times10^n, \qquad 1\le a<10,

其中 nn 为整数。若把“数量级”定义为最接近 xx1010 的整数幂,则当 a<103.16a<\sqrt{10}\approx3.16 时数量级为 10n10^n,当 a10a\ge\sqrt{10} 时数量级为 10n+110^{n+1}。也有资料只把指数 nn 称为数量级,因此报告时应说明约定。

例如 2.0×103m2.0\times10^3\,\mathrm m 最接近 103m10^3\,\mathrm m,而 8.0×103m8.0\times10^3\,\mathrm m 最接近 104m10^4\,\mathrm m。若只问两量相差几个数量级,使用比值最清楚:

Δn=log10(x2x1).\Delta n=\log_{10}\left(\frac{x_2}{x_1}\right).

只有同量纲量之比才可直接取对数。把 5kg5\,\mathrm{kg}2m2\,\mathrm m 比“数量级大小”没有物理含义。数量级比较适合跨越很大范围,却会抹去因子二或三的差异;当安全裕度只有 20%20\% 时,单一数量级远远不够。

幂律尺度与相似变换

幂律尺度

若在给定范围内

y=Cxα,y=Cx^\alpha,

则称 yyxx 服从指数为 α\alpha 的幂律尺度。把 xx 放大 λ\lambda 倍时,yy 放大 λα\lambda^\alpha 倍:

y(λx)y(x)=λα.\frac{y(\lambda x)}{y(x)}=\lambda^\alpha.

α\alpha 是无量纲的尺度指数;常数 CC 一般带单位,使等式量纲齐次。

对数坐标把幂律变成直线:

logyy0=log(Cx0αy0)+αlogxx0.\log\frac{y}{y_0} =\log\left(C\frac{x_0^\alpha}{y_0}\right) +\alpha\log\frac{x}{x_0}.

这里用参考量 x0x_0y0y_0 先构造无量纲比,避免对有量纲量直接取对数。图上线性并不能自动证明单一物理机制;有限范围内,其他函数也可能近似直线,仍需检查残差、范围和模型假设。

几何相似是最直接的尺度律。把物体所有长度乘 λ\lambda,相似面积乘 λ2\lambda^2,相似体积乘 λ3\lambda^3。因此表面积与体积之比按 λ1\lambda^{-1} 缩放。这个结论依赖形状相似;若厚度、姿态或内部结构随大小改变,就不能机械使用同一指数。

例 1:放大模型的材料与表面负担

一个实心相似模型的原长度为 0.40m0.40\,\mathrm m、表面积为 0.60m20.60\,\mathrm{m^2}、体积为 0.020m30.020\,\mathrm{m^3}。把所有长度放大到原来的 λ=3\lambda=3 倍,并假设材料密度保持 ρ=800kgm3\rho=800\,\mathrm{kg\,m^{-3}}

新表面积和体积为

A=λ2A=9(0.60m2)=5.4m2,A'=\lambda^2A=9(0.60\,\mathrm{m^2})=5.4\,\mathrm{m^2},
V=λ3V=27(0.020m3)=0.54m3.V'=\lambda^3V=27(0.020\,\mathrm{m^3})=0.54\,\mathrm{m^3}.

质量从 m=ρV=16kgm=\rho V=16\,\mathrm{kg} 增到 m=432kgm'=432\,\mathrm{kg}。表面积体积比由 30m130\,\mathrm{m^{-1}} 降为 10m110\,\mathrm{m^{-1}}。若散热能力大致随表面积、产热大致随体积,几何放大后两者不能保持同一比例;但这只是趋势判断,真实散热还依赖温差、流动和材料结构。

Fermi 估算:把未知目标拆成可估因子

Fermi 估算适合数据不完整但需要先验尺度的问题。一个可审计流程是:

  1. 明确目标量、时间或空间边界,并写出所需单位;
  2. 构造最简单的乘除关系,避免一开始引入难以估计的细节;
  3. 为每个输入写中心值、下界、上界、单位和依据;
  4. 计算中心估计,并用端点组合得到保守范围;
  5. 找出改变一倍就最能影响答案的输入,优先改进它;
  6. 用另一种分解、已知容量或极限情形核对;
  7. 报告结论、假设、范围和未包含的机制,不伪造有效数字。

若目标为正量的乘积

Q=Ci=1nxiai,Q=C\prod_{i=1}^{n}x_i^{a_i},

并且每个 xix_i 有正区间 [xi,,xi,+][x_{i,-},x_{i,+}],可按指数正负选择端点得到保守上下界。粗略的相对变化分析为

ΔQQiaiΔxixi.\frac{\Delta Q}{Q}\approx \sum_i a_i\frac{\Delta x_i}{x_i}.

它用于识别敏感因子,不应冒充严格概率不确定度。若输入相关、区间很宽或存在相减抵消,应直接计算边界或建立更合适模型;后续不确定度章节会区分统计传播与最坏界。

例 2:一天呼吸多少体积的空气

目标是静息与日常轻度活动混合的一天内吸入空气总体积。采用简化关系

Vday=(每分钟呼吸次数)(每次体积)(1440min/day).V_{\mathrm{day}}=(\text{每分钟呼吸次数}) (\text{每次体积})(1440\,\mathrm{min/day}).

取中心值 12min112\,\mathrm{min^{-1}}0.50L/breath0.50\,\mathrm{L/breath},得到

Vday=12breathmin×0.50Lbreath×1440minday=8640L/day=8.64m3/day.V_{\mathrm{day}} =12\,\frac{\mathrm{breath}}{\mathrm{min}} \times0.50\,\frac{\mathrm L}{\mathrm{breath}} \times1440\,\frac{\mathrm{min}}{\mathrm{day}} =8640\,\mathrm{L/day}=8.64\,\mathrm{m^3/day}.

若为估算而取呼吸频率 8820min120\,\mathrm{min^{-1}}、每次体积 0.350.350.75L0.75\,\mathrm L,端点给出约 4.04.021.6m3/day21.6\,\mathrm{m^3/day}。这个宽范围明确暴露了活动状态与个体差异;它不适合诊断任何个人。独立核对可把中心通气率 6L/min6\,\mathrm{L/min} 乘一天分钟数,得到同一中心结果。

上下界、中心值与误差边界

估算输入往往只知道区间。对全为正数的乘法,直接使用端点最透明。例如屋顶汇水体积在忽略损失时为 V=hAV=hA。若降雨深度 hh151525mm25\,\mathrm{mm},水平投影面积 AA110110130m2130\,\mathrm{m^2},先换算 h=0.015h=0.0150.025m0.025\,\mathrm m,可得

1.65m3V3.25m3.1.65\,\mathrm{m^3}\le V\le3.25\,\mathrm{m^3}.

使用中心值 h=0.020mh=0.020\,\mathrm mA=120m2A=120\,\mathrm{m^2}V0=2.4m3V_0=2.4\,\mathrm{m^3}。若沟槽、蒸发和溢流使可收集比例 η\eta 约为 0.70.70.90.9,实际模型应写 Vc=ηhAV_c=\eta hA,范围改为约 1.161.162.93m32.93\,\mathrm{m^3}。把效率遗漏后再补一句“有误差”,不如从一开始把它列为无量纲假设。

区间中点不总该用算术平均。若只知道某参数位于 10210^210410^4 且跨越的是乘法尺度,几何中心 102104=103\sqrt{10^2\cdot10^4}=10^3 更自然,因为它到上下界各相差因子十。最终报告可以写“中心约 10310^3,合理范围 10210^210410^4”,而非给出虚假的个位数精度。

例 3:带收集效率的屋顶雨水估算

某屋顶水平投影面积约 A=120m2A=120\,\mathrm{m^2},一次降雨深度约 h=20mm=0.020mh=20\,\mathrm{mm}=0.020\,\mathrm m,收集效率取 η=0.80\eta=0.80。中心估计为

Vc=ηAh=(0.80)(120m2)(0.020m)=1.92m3.V_c=\eta Ah =(0.80)(120\,\mathrm{m^2})(0.020\,\mathrm m) =1.92\,\mathrm{m^3}.

AA 的范围是 110110130m2130\,\mathrm{m^2}hh151525mm25\,\mathrm{mm}η\eta0.700.700.900.90,因三个因子均为正,保守边界为

V=(0.70)(110)(0.015)m31.16m3,V_-=(0.70)(110)(0.015)\,\mathrm{m^3}\approx1.16\,\mathrm{m^3},
V+=(0.90)(130)(0.025)m32.93m3.V_+=(0.90)(130)(0.025)\,\mathrm{m^3}\approx2.93\,\mathrm{m^3}.

因此合适报告是“约 1.9m31.9\,\mathrm{m^3},在给定输入范围下约 1.21.22.9m32.9\,\mathrm{m^3}”。这个区间是最坏端点组合,不是具有某一置信水平的概率区间。

无量纲比与主导尺度

决定能否忽略某个效应的通常不是某个有单位参数单独“大”或“小”,而是它与竞争尺度的无量纲比。例如物体特征长度 L=0.50mL=0.50\,\mathrm m,边缘影响厚度 δ=2.0mm\delta=2.0\,\mathrm{mm},比值

ε=δL=2.0×103m0.50m=4.0×103\varepsilon=\frac{\delta}{L} =\frac{2.0\times10^{-3}\,\mathrm m}{0.50\,\mathrm m} =4.0\times10^{-3}

说明边缘层厚度仅为主尺度的 0.4%0.4\%。但这不能自动证明它对所有观测量都可忽略;若目标恰是表面通量,小体积区域仍可能贡献主要效应。必须把“哪个比值小”与“目标量怎样依赖它”同时说明。

对和式 Q=A+BQ=A+B,若 B/A1|B/A|\ll1,在不发生抵消时可近似 QAQ\approx A,相对舍弃量约为 B/A|B/A|。若 AABB 符号相反且大小接近,则小小的输入变化会造成严重相对误差,不能以“每项都约为 10610^6”推断差值也约为 10610^6。主导项判断必须在目标表达式中进行。

例 4:何时可舍弃小修正项

某模型给出时间

t=Lv+τ,t=\frac{L}{v}+\tau,

其中 L=2.0kmL=2.0\,\mathrm{km}v=20ms1v=20\,\mathrm{m\,s^{-1}},响应延迟 τ=0.50s\tau=0.50\,\mathrm s。先统一单位: L/v=2000m/(20ms1)=100sL/v=2000\,\mathrm m/(20\,\mathrm{m\,s^{-1}})=100\,\mathrm s。无量纲比为

τL/v=0.50100=5.0×103.\frac{\tau}{L/v}=\frac{0.50}{100}=5.0\times10^{-3}.

若只要求 1%1\% 的总时间精度,忽略 τ\tau 带来的相对变化约 0.5%0.5\%,近似可接受;若目标是单独研究毫秒级响应延迟,当然不能舍弃它。近似是否合理由任务容差决定,而非由“0.50 看起来小”决定。

量纲方法与估算不能替代什么

量纲能排除不齐次表达式,也能在给定变量集合和幂律假设下确定指数关系,但它有清楚边界:

  • 不能决定 22π\pi 等纯数系数;
  • 不能在 sinx\sin xexe^x、多项式等不同无量纲函数之间自动选择;
  • 不能发现题目遗漏的变量,如边界粗糙度、形状或材料相变;
  • 不能把相关输入当作彼此独立,也不能凭区间给出概率置信水平;
  • 不能用数量级一致证明模型正确,错误机制也可能产生相同尺度;
  • 不能替代校准、重复测量和与真实数据的比较。

因此好的估算结尾应是可检验的下一步:哪一个输入最值得测量,哪一条假设最可能失效,哪种独立算法可发现单位或计数错误。

探索实验:用纸张叠层寻找尺度律

准备同一批纸、直尺或卡尺。分别取 N=20,40,80,160N=20,40,80,160 张,压平但不要额外施加不一致的力,测量纸叠总厚度 HH,单位为毫米。不要试图直接测一张纸,因为仪器分辨率可能与单张厚度同量级。先提出模型

H=H0+Nt,H=H_0+Nt,

其中 tt 是单张平均厚度,单位 mm/sheet\mathrm{mm/sheet}H0H_0 是零点、夹持或表面起伏造成的等效截距,单位毫米。

对每组计算 H/NH/N,再比较把 NN 加倍时 (HH0)(H-H_0) 是否近似加倍。若暂时忽略 H0H_0,可用最大叠层估计 tH/160t\approx H/160,并用相邻两组差分

tH2H1N2N1t\approx\frac{H_2-H_1}{N_2-N_1}

独立核对。差分会消去恒定的 H0H_0。记录尺的最小分度、纸张是否受压、边缘是否对齐,以及重复放置后的读数范围。若 H/NH/NNN 系统变化,可能是零点项、压缩或纸张不均匀,不能只把偏离称为“随机误差”。这个实验同时展示尺度指数 11、分辨率放大策略和模型残差的来源。

练习

练习

按“最接近的十次幂”约定,分别给出 2.8×107m2.8\times10^{-7}\,\mathrm m6.0×107m6.0\times10^{-7}\,\mathrm m 的数量级,并求二者比值。

查看提示
先把数写成 a×10na\times10^n,再与分界系数 103.16\sqrt{10}\approx3.16 比较。
查看解答

2.8<3.162.8<3.16,故第一量的数量级为 107m10^{-7}\,\mathrm m6.0>3.166.0>3.16,故第二量最接近 106m10^{-6}\,\mathrm m。二者同量纲,可以相除:

6.0×1072.8×1072.14.\frac{6.0\times10^{-7}}{2.8\times10^{-7}}\approx2.14.

它们实际只相差约因子 2.12.1,这也说明把相邻数量级误读成“必定相差十倍”会丢失信息。

练习

某相似结构的质量满足 mL3m\propto L^3。原结构 L1=0.50mL_1=0.50\,\mathrm mm1=12kgm_1=12\,\mathrm{kg}。若 L2=0.80mL_2=0.80\,\mathrm m 且材料与形状不变,估算新质量并列出关键假设。

查看提示
使用比值式 y2/y1=(x2/x1)αy_2/y_1=(x_2/x_1)^\alpha,无需先求常数 CC
查看解答
m2m1=(0.80m0.50m)3=1.63=4.096,\frac{m_2}{m_1}=\left(\frac{0.80\,\mathrm m}{0.50\,\mathrm m}\right)^3 =1.6^3=4.096,

所以 m249kgm_2\approx49\,\mathrm{kg}。假设包括三个方向按同一比例缩放、内部不是空心或重新加固、材料密度相同。若只改变长度而厚度不变,指数不会是 33

练习

估算 800800 人场所一天的饮用水需求。假设每人每天饮水 1.21.22.5L2.5\,\mathrm L,中心取 1.8L1.8\,\mathrm L。报告中心值与范围,暂不包括清洁和烹饪用水。

查看提示
用人数乘每人每日用量;先把升换成 m3\mathrm{m^3},再给低、高端点。
查看解答

中心估计为

800person×1.8Lpersonday=1440L/day=1.44m3/day.800\,\mathrm{person}\times1.8\, \frac{\mathrm L}{\mathrm{person\,day}} =1440\,\mathrm{L/day}=1.44\,\mathrm{m^3/day}.

端点范围为 0.960.962.00m3/day2.00\,\mathrm{m^3/day}。结论只覆盖直接饮用,且默认人数全天等效在场;若用于储水设计,还需加入到访波动、补给周期与安全余量。

练习

设备功率在 454555W55\,\mathrm W,每天运行 779h9\,\mathrm h。求每日能耗的保守范围,分别用 Wh\mathrm{Wh}J\mathrm J 表示,并给中心估计。

查看提示
对正量 E=PtE=Pt,最小值取两个下界,最大值取两个上界;统一使用瓦和秒。
查看解答

端点为

E=45W×7h=315Wh,E+=55W×9h=495Wh.E_-=45\,\mathrm W\times7\,\mathrm h=315\,\mathrm{Wh}, \quad E_+=55\,\mathrm W\times9\,\mathrm h=495\,\mathrm{Wh}.

中心取 P=50WP=50\,\mathrm Wt=8ht=8\,\mathrm h,得 400Wh400\,\mathrm{Wh}。因 1Wh=3600J1\,\mathrm{Wh}=3600\,\mathrm J,范围为 1.1341.1341.782MJ1.782\,\mathrm{MJ},中心为 1.44MJ1.44\,\mathrm{MJ}。区间是输入端点组合,不代表某个统计置信度。

练习

模型 Q=A+BxQ=A+Bx 中, A=200NA=200\,\mathrm NB=0.40Nm1B=0.40\,\mathrm{N\,m^{-1}}x=5.0mx=5.0\,\mathrm m。若要求 2%2\% 精度,是否可用 QAQ\approx A?若 x=50mx=50\,\mathrm m 呢?

查看提示
先计算两项各自带单位的数值,再取被舍弃项与保留项之比。
查看解答

x=5.0mx=5.0\,\mathrm m 时, Bx=2.0NBx=2.0\,\mathrm N,与 AA 的比为 1.0%1.0\%,忽略后相对完整结果的误差约 2/2020.99%2/202\approx0.99\%,满足要求。x=50mx=50\,\mathrm m 时, Bx=20NBx=20\,\mathrm N,忽略后误差约 20/2209.1%20/220\approx9.1\%,不满足要求。BB 的单位保证 BxBxAA 都是牛顿,才能相加。

练习

心跳频率中心取 70min170\,\mathrm{min^{-1}}。用两种等价分解估算一年心跳次数,并以一位有效数字报告数量级。取一年 365day365\,\mathrm{day},同时说明该模型忽略什么。

查看提示
方法一用总秒数除每次间隔;方法二用每分钟次数乘总分钟数。两种分解的单位应都约消为次数。
查看解答

方法一直接相乘:

70beatmin×60minh×24hday×365dayyear3.68×107beat/year.70\,\frac{\mathrm{beat}}{\mathrm{min}} \times60\,\frac{\mathrm{min}}{\mathrm h} \times24\,\frac{\mathrm h}{\mathrm{day}} \times365\,\frac{\mathrm{day}}{\mathrm{year}} \approx3.68\times10^7\,\mathrm{beat/year}.

方法二先求每次间隔 60/700.86s/beat60/70\approx0.86\,\mathrm{s/beat},再用一年约 3.15×107s3.15\times10^7\,\mathrm s 相除,也得约 3.7×1073.7\times10^7 次。以一位有效数字可报 4×1074\times10^7 次每年,最接近的十进数量级为 10810^8 次。模型忽略睡眠、运动、健康状态和个体随时间变化,不宜用于个人医学推断。

关系、资源与后续学习

书籍 · 2016

University Physics Volume 1

Samuel J. Ling, Jeff Sanny, William Moebs

用于核对 P00 的基础术语、量纲规则、估算步骤、测量报告和入门不确定度计算。

打开官方来源

OpenStax《University Physics Volume 1》第一章给出科学记数法、数量级、单位转换、量纲分析和 Fermi 估算的入门路线。使用其中的估算题时,应额外练习把每个假设写成区间,并给出独立核对,而不是只追求与参考数值完全相同。

下一步进入测量误差与校准:估算提供实验前的预期尺度,真实测量则会告诉我们中心值、波动、偏差和仪器限制。之后再学习不确定度传播,把本章的保守范围与有明确统计含义的报告区分开。