📚 IB Mathematics: Differential Equations in Physical Models | IB数学：微分方程的物理模型
Differential equations are the language of change in physics. They describe how quantities evolve over time or space, from the cooling of a cup of coffee to the orbit of a planet. In IB Mathematics, particularly at Higher Level, students learn to construct, solve, and interpret differential equations that model real-world physical systems. This article will guide you through the most important models, step by step, with both mathematical detail and physical intuition.
微分方程是物理学中描述变化的语言。它们刻画了量随时间或空间如何演化，从一杯咖啡的冷却到行星的轨道，无一例外。在IB数学中，尤其是在高级水平课程中，学生需要学会构建、求解并解释用于模拟真实物理系统的微分方程。本文将带你逐步了解最重要的物理模型，兼顾数学细节与物理直觉。
1. Why Differential Equations Matter in Physics | 物理中微分方程的重要性
Newton’s second law \( F = ma \) is already a differential equation, because acceleration is the second derivative of position. When the force depends on position, velocity, or time, the equation becomes a differential equation that must be solved to find the motion of an object. Similarly, the rate of radioactive decay, the charge in a circuit, and the growth of a population all obey differential equations. Understanding this connection allows physicists to predict the future from the present state of a system.
牛顿第二定律 \( F = ma \) 本身就是一个微分方程，因为加速度是位置对时间的二阶导数。当力依赖于位置、速度或时间时，方程就变成必须求解以确定物体运动的微分方程。类似地，放射性衰变速率、电路中的电荷以及种群增长都遵循微分方程。理解这种联系使物理学家能够根据系统的当前状态预测未来。
2. Steps for Constructing a Physical Model | 构建物理模型的步骤
To model any physical situation with a differential equation, follow these steps: (1) identify the dependent and independent variables; (2) express the rate of change using derivatives; (3) write down the physical law that relates the rate to the state variables; (4) define the initial or boundary conditions; (5) solve the equation analytically or numerically, and interpret the result in the context of the problem.
要用微分方程对任何物理情境建模，可以遵循以下步骤：(1) 确定自变量和因变量；(2) 用导数表示变化率；(3) 写出物理定律，将变化率与状态变量联系起来；(4) 给出初始条件或边界条件；(5) 解析或数值求解方程，并结合具体问题解释结果。
3. Newton’s Law of Cooling | 牛顿冷却定律
Newton’s law of cooling states that the rate of change of the temperature \( T \) of an object is proportional to the difference between its temperature and the ambient temperature \( T_a \). The model is \( \frac{dT}{dt} = -k(T – T_a) \), where \( k > 0 \) is a constant that depends on the surface area and heat transfer coefficient. The negative sign indicates that the object cools when it is hotter than the surroundings.
牛顿冷却定律指出，物体温度 \( T \) 的变化率与物体和周围环境温度 \( T_a \) 之差成正比。其模型为 \( \frac{dT}{dt} = -k(T – T_a) \)，其中 \( k > 0 \) 是由表面积和传热系数决定的常数。负号表示当物体比周围环境温度高时，它会冷却。
This is a first-order linear differential equation. Its general solution is \( T(t) = T_a + (T_0 – T_a)e^{-kt} \), where \( T_0 \) is the initial temperature. As \( t \to \infty \), \( T \) approaches \( T_a \), but never exactly reaches it in finite time. This exponential approach is a common feature of many physical relaxation processes.
这是一阶线性微分方程。其通解为 \( T(t) = T_a + (T_0 – T_a)e^{-kt} \)，其中 \( T_0 \) 是初始温度。当 \( t \to \infty \) 时，\( T \) 趋近于 \( T_a \)，但在有限时间内永远不会精确达到。这种指数逼近是许多物理弛豫过程的共同特征。
4. Radioactive Decay and Exponential Decay | 放射性衰变与指数衰减
Radioactive decay is governed by the simple differential equation \( \frac{dN}{dt} = -\lambda N \), where \( N(t) \) is the number of undecayed nuclei and \( \lambda \) is the decay constant. The solution is \( N(t) = N_0 e^{-\lambda t} \). The half-life \( T_{1/2} \) satisfies \( N(T_{1/2}) = N_0/2 \), giving \( T_{1/2} = \frac{\ln 2}{\lambda} \).
放射性衰变遵循简单的微分方程 \( \frac{dN}{dt} = -\lambda N \)，其中 \( N(t) \) 是未衰变原子核的数目，\( \lambda \) 为衰变常数。其解为 \( N(t) = N_0 e^{-\lambda t} \)。半衰期 \( T_{1/2} \) 满足 \( N(T_{1/2}) = N_0/2 \)，从而得到 \( T_{1/2} = \frac{\ln 2}{\lambda} \)。
This model also applies to capacitor discharge in an RC circuit, absorption of light in a medium, and elimination of drugs from the bloodstream. The key idea is that the rate of decrease is proportional to the amount present, which leads to a constant proportional decay per unit time.
该模型同样适用于RC电路中的电容放电、介质中的光吸收以及药物从血液中的消除。核心思想是减少速率与当前量成正比，从而使得单位时间内的衰减比例恒定。
5. Falling Bodies with Air Resistance | 受空气阻力的落体运动
When an object falls under gravity and experiences air resistance proportional to velocity, its equation of motion is \( m\frac{dv}{dt} = mg – kv \), where \( m \) is mass, \( g \) is gravitational acceleration, and \( k \) is the drag coefficient. Dividing by \( m \), let \( \alpha = k/m \), then \( \frac{dv}{dt} = g – \alpha v \). This is a linear first-order equation.
当物体在重力作用下下落并受到正比于速度的空气阻力时，其运动方程为 \( m\frac{dv}{dt} = mg – kv \)，其中 \( m \) 是质量，\( g \) 是重力加速度，\( k \) 是阻力系数。两边除以 \( m \)，令 \( \alpha = k/m \)，得 \( \frac{dv}{dt} = g – \alpha v \)。这是一阶线性方程。
The solution is \( v(t) = \frac{g}{\alpha}(1 – e^{-\alpha t}) \), starting from rest. As \( t \to \infty \), the velocity approaches the terminal velocity \( v_T = \frac{g}{\alpha} = \frac{mg}{k} \). This occurs when the drag force balances the weight, so acceleration becomes zero. In IB problems, you may need to sketch the velocity-time graph showing exponential approach to the terminal speed.
从静止开始的解为 \( v(t) = \frac{g}{\alpha}(1 – e^{-\alpha t}) \)。当 \( t \to \infty \) 时，速度趋于极限速度 \( v_T = \frac{g}{\alpha} = \frac{mg}{k} \)。此时阻力与重力平衡，加速度为零。在IB习题中，你可能需要绘制速度-时间图像，展示其指数趋近于极限速度的过程。
6. Logistic Growth Model | 逻辑斯蒂增长模型
For populations with limited resources, the growth rate decreases as the population \( P \) approaches the carrying capacity \( K \). The logistic differential equation is \( \frac{dP}{dt} = rP\left(1 – \frac{P}{K}\right) \), where \( r \) is the intrinsic growth rate. This nonlinear equation can be solved by separation of variables, yielding the logistic function.
对于资源有限的种群，其增长率会随着种群数量 \( P \) 逼近环境承载量 \( K \) 而下降。逻辑斯蒂微分方程为 \( \frac{dP}{dt} = rP\left(1 – \frac{P}{K}\right) \)，其中 \( r \) 为内禀增长率。这一非线性方程可通过分离变量求解，得到逻辑斯蒂函数。
The solution is \( P(t) = \frac{K}{1 + A e^{-rt}} \), where \( A = \frac{K – P_0}{P_0} \). The graph is an S-shaped curve: exponential growth initially, then slowing, and finally leveling off at \( P = K \). In physics, similar equations model the spread of a virus or the charging of a capacitor with a nonlinear element.
其解为 \( P(t) = \frac{K}{1 + A e^{-rt}} \)，其中 \( A = \frac{K – P_0}{P_0} \)。图像是S形曲线：初期指数增长，随后减缓，最后趋平于 \( P = K \)。在物理学中，类似的方程也用于模拟病毒传播或含有非线性元件的电容充电过程。
7. Simple Harmonic Motion | 简谐振动
Simple harmonic motion (SHM) occurs when the restoring force is proportional to displacement and opposite in direction: \( F = -kx \). Newton’s law gives \( m\frac{d^2x}{dt^2} = -kx \), or \( \frac{d^2x}{dt^2} + \omega^2 x = 0 \), where \( \omega = \sqrt{k/m} \) is the angular frequency. This is a second-order linear differential equation with constant coefficients.
当回复力正比于位移且方向相反时，即 \( F = -kx \)，就会发生简谐振动。牛顿定律给出 \( m\frac{d^2x}{dt^2} = -kx \)，即 \( \frac{d^2x}{dt^2} + \omega^2 x = 0 \)，其中 \( \omega = \sqrt{k/m} \) 为角频率。这是常系数二阶线性微分方程。
The general solution is \( x(t) = A\cos(\omega t + \phi) \), where \( A \) is the amplitude and \( \phi \) is the phase constant. The values of \( A \) and \( \phi \) are determined by the initial displacement and velocity. In IB, you may be asked to derive the equation from energy considerations or to use \( x = A\sin(\omega t + \phi) \) as an equivalent form. The period is \( T = \frac{2\pi}{\omega} \).
通解为 \( x(t) = A\cos(\omega t + \phi) \)，其中 \( A \) 是振幅，\( \phi \) 是初相位。\( A \) 和 \( \phi \) 由初始位移和初始速度决定。在IB中，你可能需要从能量角度推导该方程，或者使用等价形式 \( x = A\sin(\omega t + \phi) \)。周期为 \( T = \frac{2\pi}{\omega} \)。
8. Damped Harmonic Motion | 阻尼振动
Real oscillators experience damping, often proportional to velocity. The equation becomes \( m\frac{d^2x}{dt^2} + c\frac{dx}{dt} + kx = 0 \), where \( c \) is the damping coefficient. Dividing by \( m \), we get \( \frac{d^2x}{dt^2} + 2\beta\frac{dx}{dt} + \omega_0^2 x = 0 \), with \( \beta = \frac{c}{2m} \) and \( \omega_0 = \sqrt{k/m} \).
真实振子会受到阻尼，通常与速度成正比。方程为 \( m\frac{d^2x}{dt^2} + c\frac{dx}{dt} + kx = 0 \)，其中 \( c \) 是阻尼系数。除以 \( m \)，得到 \( \frac{d^2x}{dt^2} + 2\beta\frac{dx}{dt} + \omega_0^2 x = 0 \)，其中 \( \beta = \frac{c}{2m} \)，\( \omega_0 = \sqrt{k/m} \)。
The solution depends on the discriminant. If \( \beta < \omega_0 \), the system is underdamped and oscillates with exponentially decaying amplitude: \( x(t) = A e^{-\beta t}\cos(\omega_d t + \phi) \), where \( \omega_d = \sqrt{\omega_0^2 - \beta^2} \). If \( \beta = \omega_0 \), it is critically damped, and if \( \beta > \omega_0 \), it is overdamped. Critical damping is important in suspension systems and precision instruments.
解取决于判别式。若 \( \beta < \omega_0 \)，系统为欠阻尼，作振幅指数衰减的振动：\( x(t) = A e^{-\beta t}\cos(\omega_d t + \phi) \)，其中 \( \omega_d = \sqrt{\omega_0^2 - \beta^2} \)。若 \( \beta = \omega_0 \)，则为临界阻尼；若 \( \beta > \omega_0 \)，则为过阻尼。临界阻尼在悬挂系统和精密仪器中非常重要。
9. Forced Oscillations and Resonance | 受迫振动与共振
When an external periodic force \( F_0\cos(\omega t) \) acts on an oscillator, the equation becomes \( m\frac{d^2x}{dt^2} + c\frac{dx}{dt} + kx = F_0\cos(\omega t) \). The solution consists of a transient term that dies out and a steady-state term with the same frequency as the driving force. The amplitude of the steady-state response depends strongly on the driving frequency \( \omega \).
当外部周期力 \( F_0\cos(\omega t) \) 作用于振子时，方程为 \( m\frac{d^2x}{dt^2} + c\frac{dx}{dt} + kx = F_0\cos(\omega t) \)。解包含随时间衰减的瞬态项和与驱动力同频的稳态项。稳态响应的振幅强烈依赖于驱动频率 \( \omega \)。
Resonance occurs when \( \omega \) is close to the natural frequency \( \omega_0 \). In the absence of damping, the amplitude grows without bound if \( \omega = \omega_0 \). With damping, the amplitude peaks at a frequency slightly below \( \omega_0 \). Understanding resonance is crucial in engineering, from bridge design to speaker construction. In IB, you might be asked to explain why soldiers break step when crossing a bridge.
当 \( \omega \) 接近固有频率 \( \omega_0 \) 时会发生共振。若无阻尼，当 \( \omega = \omega_0 \) 时振幅无限增长。有阻尼时，振幅在略低于 \( \omega_0 \) 的频率处达到峰值。理解共振在工程中至关重要，从桥梁设计到扬声器制造。在IB中，你可能会被要求解释为什么士兵过桥时要碎步。
10. Electrical Circuits: RL and RC Models | 电路模型：RL与RC
In an RL circuit containing a resistor and an inductor, Kirchhoff’s voltage law gives \( L\frac{di}{dt} + Ri = V(t) \), where \( i(t) \) is the current. For a constant voltage \( V_0 \), the solution is \( i(t) = \frac{V_0}{R}(1 – e^{-Rt/L}) \). The time constant is \( \tau = L/R \), which is the time for the current to reach about 63% of its final value.
在含有电阻和电感的RL电路中，基尔霍夫电压定律给出 \( L\frac{di}{dt} + Ri = V(t) \)，其中 \( i(t) \) 是电流。对于恒定电压 \( V_0 \)，其解为 \( i(t) = \frac{V_0}{R}(1 – e^{-Rt/L}) \)。时间常数为 \( \tau = L/R \)，即电流达到终值约63%所需的时间。
Similarly, an RC circuit with a capacitor obeys \( R\frac{dq}{dt} + \frac{q}{C} = V(t) \), where \( q \) is charge. The solution for charging is \( q(t) = CV_0(1 – e^{-t/(RC)}) \), and for discharging, \( q(t) = q_0 e^{-t/(RC)} \). These equations are mathematically identical to the Newton’s cooling and radioactive decay models, showing the unity of physical laws.
类似地，含有电容的RC电路满足 \( R\frac{dq}{dt} + \frac{q}{C} = V(t) \)，其中 \( q \) 是电荷。充电过程解为 \( q(t) = CV_0(1 – e^{-t/(RC)}) \)，放电过程为 \( q(t) = q_0 e^{-t/(RC)} \)。这些方程在数学上与牛顿冷却和放射性衰变模型完全一致，体现了物理规律的统一性。
11. Coupled Differential Equations in Physics | 物理中的耦合微分方程组
Many physical systems involve two or more interacting quantities, leading to coupled differential equations. For example, the Lotka-Volterra predator-prey model uses \( \frac{dx}{dt} = \alpha x – \beta xy \) and \( \frac{dy}{dt} = \delta xy – \gamma y \), where \( x \) is prey and \( y \) is predator. In physics, coupled oscillators and two-body problems require solving systems of differential equations.
许多物理系统涉及两个或更多相互作用的量，从而产生耦合微分方程。例如，Lotka-Volterra捕食者-猎物模型使用 \( \frac{dx}{dt} = \alpha x – \beta xy \) 和 \( \frac{dy}{dt} = \delta xy – \gamma y \)，其中 \( x \) 是猎物，\( y \) 是捕食者。在物理学中，耦合振子和两体问题都需要求解微分方程组。
In IB Mathematics, you are not required to solve coupled nonlinear systems analytically, but you may use phase plane analysis or Euler’s method to approximate solutions. Understanding how derivatives link multiple variables is essential for university-level physics and engineering.
在IB数学中，并不要求解析求解耦合非线性系统，但你可以使用相平面分析或欧拉方法来近似求解。理解导数如何将多个变量联系起来，对于大学水平的物理和工程学至关重要。
12. Numerical Methods: Euler’s Method | 数值方法：欧拉法
When a differential equation cannot be solved analytically, we use numerical methods. Euler’s method approximates the solution by taking small steps: \( y_{n+1} = y_n + h \cdot f(x_n, y_n) \), where \( h \) is the step size. The error decreases as \( h \) decreases, though very small steps require more computation. This method is often used in IB coursework to model physical systems with variable coefficients or nonlinear terms.
当微分方程无法解析求解时，我们使用数值方法。欧拉法通过小步进近似解：\( y_{n+1} = y_n + h \cdot f(x_n, y_n) \)，其中 \( h \) 是步长。误差随 \( h \) 减小而减小，但过小的步长会增加计算量。在IB课程中，这种方法常用于模拟具有变系数或非线性项的物理系统。
For example, to model the motion of a pendulum with large amplitude, the equation \( \frac{d^2\theta}{dt^2} = -\frac{g}{L}\sin\theta \) cannot be solved with elementary functions. Euler’s method, or better, a Runge-Kutta method, provides a numerical solution that can be plotted and analyzed. In IB, you may be asked to apply Euler’s method to a first-order differential equation and comment on the accuracy.
例如，要模拟大角度单摆运动，方程 \( \frac{d^2\theta}{dt^2} = -\frac{g}{L}\sin\theta \) 无法用初等函数求解。欧拉法或更优的龙格-库塔法可以提供数值解，供绘图和分析。在IB中，你可能会被要求将欧拉法应用于一阶微分方程并讨论其精度。
Differential equations are more than abstract symbols; they are the bridge between mathematical theory and physical reality. From cooling to oscillating, from circuits to populations, the same mathematical structures appear again and again. By mastering these models, IB students gain powerful tools for understanding the universe and are well prepared for further studies in mathematics, physics, and engineering.
微分方程不仅仅是抽象的符号，它们是数学理论与物理现实之间的桥梁。从冷却到振动，从电路到生物种群，相同的数学结构反复出现。通过掌握这些模型，IB学生将获得理解宇宙的强大工具，并为大学阶段进一步学习数学、物理和工程学做好充分准备。
Published by TutorHao | IB Mathematics Revision Series | aleveler.com
更多咨询请联系16621398022（同微信）
屏轩国际教育cambridge primary/secondary checkpoint, cat4, ukiset,ukcat,igcse,alevel,PAT,STEP,MAT, ibdp,ap,ssat,sat,sat2课程辅导，国外大学本科硕士研究生博士课程论文辅导