用泰勒级数来描述微小振动和波动的行为 - 公式详解版

用泰勒级数来描述微小振动和波动的行为 - 公式详解版

flyfish

1. 简谐振动(Harmonic Oscillator)

势能函数(Potential Energy Function) 势能函数 V ( x ) V(x) V(x) 描述了系统在位移 x x x 处的势能。对于一个简谐振子(如弹簧振子),势能函数为: V ( x ) = 1 2 k x 2 V(x) = \frac{1}{2} k x^2 V(x)=21kx2
其中:

  • V ( x ) V(x) V(x): 势能 (Potential Energy),单位为焦耳 (Joules, J)

  • k k k: 弹簧常数 (Spring Constant),单位为牛顿每米 (N/m)

  • x x x: 位移 (Displacement),单位为米 (Meters, m)
    运动方程(Equation of Motion) 通过牛顿第二定律 F = m a F = ma F=ma 和胡克定律 F = − k x F = -kx F=kx,得到系统的运动方程: m d 2 x d t 2 = − k x m \frac{d^2 x}{dt^2} = -kx mdt2d2x=kx
    其中:

  • m m m: 质量 (Mass),单位为千克 (Kilograms, kg)

  • d 2 x d t 2 \frac{d^2 x}{dt^2} dt2d2x: 加速度 (Acceleration),单位为米每二次方秒 (m/s^2)

  • k k k: 弹簧常数 (Spring Constant),单位为牛顿每米 (N/m)

  • x x x: 位移 (Displacement),单位为米 (Meters, m)
    简谐解(Harmonic Solution)
    对于简谐振动,运动方程的解为:
    x ( t ) = A cos ⁡ ( ω t + ϕ ) x(t) = A \cos(\omega t + \phi) x(t)=Acos(ωt+ϕ)
    其中:

  • x ( t ) x(t) x(t): 时间 t t t 时的位移 (Displacement at time t t t),单位为米 (Meters, m)

  • A A A: 振幅 (Amplitude),单位为米 (Meters, m)

  • ω \omega ω: 角频率 (Angular Frequency),单位为弧度每秒 (Radians per second, rad/s)

  • t t t: 时间 (Time),单位为秒 (Seconds, s)

  • ϕ \phi ϕ: 相位角 (Phase Angle),单位为弧度 (Radians, rad)

2. 波动方程(Wave Equation)

一维波动方程(One-Dimensional Wave Equation)
描述一维波动的基本方程为:
∂ 2 u ∂ t 2 = c 2 ∂ 2 u ∂ x 2 \frac{\partial^2 u}{\partial t^2} = c^2 \frac{\partial^2 u}{\partial x^2} t22u=c2x22u
其中:

  • u ( x , t ) u(x, t) u(x,t): 在位置 x x x 和时间 t t t 的位移 (Displacement at position x x x and time t t t),单位为米 (Meters, m)

  • ∂ 2 u ∂ t 2 \frac{\partial^2 u}{\partial t^2} t22u: 对时间的二阶偏导数 (Second partial derivative with respect to time),单位为米每二次方秒 (m/s^2)

  • c c c: 波速 (Wave Speed),单位为米每秒 (Meters per second, m/s)

  • ∂ 2 u ∂ x 2 \frac{\partial^2 u}{\partial x^2} x22u: 对位置的二阶偏导数 (Second partial derivative with respect to position),单位为米每二次方米 (1/m)
    波的解(Solution of the Wave)
    波动方程的通解通常可以表示为简谐波的叠加:
    u ( x , t ) = ∑ n A n cos ⁡ ( k n x − ω n t + ϕ n ) u(x, t) = \sum_{n} A_n \cos(k_n x - \omega_n t + \phi_n) u(x,t)=nAncos(knxωnt+ϕn)
    其中:

  • u ( x , t ) u(x, t) u(x,t): 在位置 x x x 和时间 t t t 的位移 (Displacement at position x x x and time t t t),单位为米 (Meters, m)

  • A n A_n An: 第 n n n 个简谐波的振幅 (Amplitude of the n n n-th harmonic),单位为米 (Meters, m)

  • k n k_n kn: 波数 (Wave Number) 为 2 π λ n \frac{2\pi}{\lambda_n} λn2π,单位为弧度每米 (Radians per meter, rad/m)

  • λ n \lambda_n λn: 波长 (Wavelength) 为 2 π k n \frac{2\pi}{k_n} kn2π,单位为米 (Meters, m)

  • ω n \omega_n ωn: 角频率 (Angular Frequency) 为 2 π f n 2\pi f_n 2πfn,单位为弧度每秒 (Radians per second, rad/s)

  • f n f_n fn: 频率 (Frequency),单位为赫兹 (Hertz, Hz)

  • ϕ n \phi_n ϕn: 相位角 (Phase Angle),单位为弧度 (Radians, rad)

  • t t t: 时间 (Time),单位为秒 (Seconds, s)

  • x x x: 位置 (Position),单位为米 (Meters, m)

3. 微小振动的稳定性分析(Stability Analysis of Small Oscillations)

势能函数(Potential Energy Function) 对于一个一般系统,势能函数 V ( x ) V(x) V(x) 在平衡位置 x 0 x_0 x0 附近展开: V ( x ) ≈ V ( x 0 ) + 1 2 V ′ ′ ( x 0 ) ( x − x 0 ) 2 V(x) \approx V(x_0) + \frac{1}{2} V''(x_0) (x - x_0)^2 V(x)V(x0)+21V′′(x0)(xx0)2
其中:

  • V ( x ) V(x) V(x): 势能 (Potential Energy),单位为焦耳 (Joules, J)

  • x x x: 位移 (Displacement),单位为米 (Meters, m)

  • x 0 x_0 x0: 平衡位置 (Equilibrium Position),单位为米 (Meters, m)

  • V ′ ′ ( x 0 ) V''(x_0) V′′(x0): 势能函数在平衡位置的二阶导数 (Second derivative of potential energy function at equilibrium position),单位为焦耳每二次方米 (J/m²)
    运动方程(Equation of Motion)
    线性化的运动方程为:
    m d 2 ( δ x ) d t 2 ≈ − V ′ ′ ( x 0 ) δ x m \frac{d^2 (\delta x)}{dt^2} \approx -V''(x_0) \delta x mdt2d2(δx)V′′(x0)δx
    其中:

  • δ x \delta x δx: 小位移 (Small displacement),单位为米 (Meters, m)

  • m m m: 质量 (Mass),单位为千克 (Kilograms, kg)

  • d 2 ( δ x ) d t 2 \frac{d^2 (\delta x)}{dt^2} dt2d2(δx): 小位移的加速度 (Acceleration of small displacement),单位为米每二次方秒 (m/s²)

  • V ′ ′ ( x 0 ) V''(x_0) V′′(x0): 势能函数在平衡位置的二阶导数 (Second derivative of potential energy function at equilibrium position),单位为焦耳每二次方米 (J/m²)

4. 具体应用例子(Specific Application Examples)

例子 1: 单摆的微小振动(Small Oscillations of a Pendulum)

势能函数(Potential Energy Function)
单摆的势能函数为:
V ( θ ) = m g L ( 1 − cos ⁡ ( θ ) ) V(\theta) = mgL(1 - \cos(\theta)) V(θ)=mgL(1cos(θ))
其中:

  • V ( θ ) V(\theta) V(θ): 势能 (Potential Energy),单位为焦耳 (Joules, J)

  • m m m: 摆锤的质量 (Mass of the pendulum),单位为千克 (Kilograms, kg)

  • g g g: 重力加速度 (Gravitational acceleration),单位为米每二次方秒 (m/s²)

  • L L L: 摆长 (Length of the pendulum),单位为米 (Meters, m)

  • θ \theta θ: 偏角 (Angle of displacement),单位为弧度 (Radians, rad)
    θ ≈ 0 \theta \approx 0 θ0 附近展开: V ( θ ) ≈ m g L ( 1 − ( 1 − θ 2 2 ) ) = 1 2 m g L θ 2 V(\theta) \approx mgL \left( 1 - \left( 1 - \frac{\theta^2}{2} \right) \right) = \frac{1}{2} mgL \theta^2 V(θ)mgL(1(12θ2))=21mgLθ2运动方程(Equation of Motion)
    运动方程变为:
    m L d 2 θ d t 2 ≈ − m g L θ mL \frac{d^2 \theta}{dt^2} \approx -mgL \theta mLdt2d2θmgLθ
    解为简谐振动:
    θ ( t ) = θ 0 cos ⁡ ( g L t + ϕ ) \theta(t) = \theta_0 \cos\left(\sqrt{\frac{g}{L}} t + \phi \right) θ(t)=θ0cos(Lg t+ϕ)

例子 2: 弦的微小振动(Small Vibrations of a String)

波动方程(Wave Equation)
对于固定在两端的弦,波动方程为:
∂ 2 y ∂ t 2 = c 2 ∂ 2 y ∂ x 2 \frac{\partial^2 y}{\partial t^2} = c^2 \frac{\partial^2 y}{\partial x^2} t22y=c2x22y
其中:

  • y ( x , t ) y(x, t) y(x,t): 弦在位置 x x x 和时间 t t t 的位移 (Displacement of the string at position x x x and time t t t),单位为米 (Meters, m)

  • ∂ 2 y ∂ t 2 \frac{\partial^2 y}{\partial t^2} t22y: 对时间的二阶偏导数 (Second partial derivative with respect to time),单位为米每二次方秒 (m/s²)

  • c c c: 波速 (Wave Speed),单位为米每秒 (Meters per second, m/s)

  • ∂ 2 y ∂ x 2 \frac{\partial^2 y}{\partial x^2} x22y: 对位置的二阶偏导数 (Second partial derivative with respect to position),单位为米每二次方米 (1/m)
    波的解(Solution of the Wave)
    采用分离变量法,解为:
    y ( x , t ) = ∑ n = 1 ∞ ( A n cos ⁡ ( n π c L t ) + B n sin ⁡ ( n π c L t ) ) sin ⁡ ( n π x L ) y(x, t) = \sum_{n=1}^{\infty} \left( A_n \cos\left(\frac{n\pi c}{L} t\right) + B_n \sin\left(\frac{n\pi c}{L} t\right) \right) \sin\left(\frac{n\pi x}{L}\right) y(x,t)=n=1(Ancos(Lct)+Bnsin(Lct))sin(Lx)
    其中:

  • y ( x , t ) y(x, t) y(x,t): 弦在位置 x x x 和时间 t t t 的位移 (Displacement of the string at position x x x and time t t t),单位为米 (Meters, m)

  • A n A_n An: 第 n n n 个简谐波的余弦项振幅 (Amplitude of the cosine term of the n n n-th harmonic),单位为米 (Meters, m)

  • B n B_n Bn: 第 n n n 个简谐波的正弦项振幅 (Amplitude of the sine term of the n n n-th harmonic),单位为米 (Meters, m)

  • n π c L \frac{n\pi c}{L} Lc: 第 n n n 个简谐波的角频率 (Angular frequency of the n n n-th harmonic),单位为弧度每秒 (Radians per second, rad/s)

  • n π x L \frac{n\pi x}{L} Lx: 第 n n n 个简谐波的波数 (Wave number of the n n n-th harmonic),单位为弧度每米 (Radians per meter, rad/m)

  • c c c: 波速 (Wave Speed),单位为米每秒 (Meters per second, m/s)

  • L L L: 弦的长度 (Length of the string),单位为米 (Meters, m)

相关推荐

  1. matlab代码写函数

    2024-07-14 22:30:02       37 阅读

最近更新

  1. docker php8.1+nginx base 镜像 dockerfile 配置

    2024-07-14 22:30:02       67 阅读
  2. Could not load dynamic library ‘cudart64_100.dll‘

    2024-07-14 22:30:02       71 阅读
  3. 在Django里面运行非项目文件

    2024-07-14 22:30:02       58 阅读
  4. Python语言-面向对象

    2024-07-14 22:30:02       69 阅读

热门阅读

  1. C++数组

    2024-07-14 22:30:02       19 阅读
  2. 数仓工具—Hive语法之正则表达式函数

    2024-07-14 22:30:02       21 阅读
  3. 【Lua】lua实现C# continue效果

    2024-07-14 22:30:02       17 阅读
  4. 探索TASKCTL和 DataStage 的ETL任务调度协同

    2024-07-14 22:30:02       17 阅读
  5. TCPDump协议分析工具

    2024-07-14 22:30:02       20 阅读
  6. 神领物流项目第二天

    2024-07-14 22:30:02       17 阅读
  7. 日记

    2024-07-14 22:30:02       22 阅读