对随机微分方程与乘法噪声的直接数值解决方案
1Department of Chemistry, University of Colorado Boulder, Boulder, Colorado 80309, USA.
Physical review letters
|July 12, 2024
概括
我们开发了一种新的数值方法,用于乘法噪声的随机微分方程,避免轨迹平均. 这种方法准确地预测了振荡器的分叉,并且在计算上比传统的模拟更有效.
科学领域:
- 计算物理学的计算物理.
- 数字分析 数字分析
- 非线性动力学是一种非线性动力学.
背景情况:
- 经典的随机微分方程 (SDEs) 与乘法噪声带来了重要的计算挑战.
- 传统方法通常依赖于轨迹平均,这可能是计算上昂贵的,可能无法准确地捕捉像分叉这样的关键现象.
- 量子力学的途径积分解决方案提供了替代的理论框架.
研究的目的:
- 开发一种新的数值方法,用于用乘数噪声解决经典的SDEs.
- 为了避免轨迹平均方法的计算负担和潜在的不准确性.
- 为了准确地模拟接近分叉制度的系统.
主要方法:
- 灵感来自量子放松路径的整体解决方案.
- 为经典的SDEs开发一个无轨道平均数值的数值方法.
- 模拟与非马科夫噪声相结合的经典振荡器.
- 使用张量因子化技术加速方法.
主要成果:
- 开发的方法准确地估计了经典振荡器的过渡到双叉模式.
- 该方法在准确度上明显优于轨迹平均模拟.
- 与传统方法相比,计算成本是数量级较低的.
结论:
- 这种新的数值方法提供了一种高效和准确的方法来模拟用乘数噪声模拟经典的SDEs.
- 这种方法特别有利于研究临界点附近的系统,例如分叉.
- 张量因子化技术为这些模拟提供了大量的计算加速度.
相关概念视频
Mechanistic Models: Compartment Models in Algorithms for Numerical Problem Solving
48
Mechanistic models play a crucial role in algorithms for numerical problem-solving, particularly in nonlinear mixed effects modeling (NMEM). These models aim to minimize specific objective functions by evaluating various parameter estimates, leading to the development of systematic algorithms. In some cases, linearization techniques approximate the model using linear equations.
In individual population analyses, different algorithms are employed, such as Cauchy's method, which uses a...
In individual population analyses, different algorithms are employed, such as Cauchy's method, which uses a...
48
Second Order systems II
96
In an underdamped second-order system, where the damping ratio ζ is between 0 and 1, a unit-step input results in a transfer function that, when transformed using the inverse Laplace method, reveals the output response. The output exhibits a damped sinusoidal oscillation, and the difference between the input and output is termed the error signal. This error signal also demonstrates damped oscillatory behavior. Eventually, as the system reaches a steady state, the error diminishes to zero.
96
Propagation of Uncertainty from Random Error
664
An experiment often consists of more than a single step. In this case, measurements at each step give rise to uncertainty. Because the measurements occur in successive steps, the uncertainty in one step necessarily contributes to that in the subsequent step. As we perform statistical analysis on these types of experiments, we must learn to account for the propagation of uncertainty from one step to the next. The propagation of uncertainty depends on the type of arithmetic operation performed on...
664
Difference Equation Solution using z-Transform
282
The z-transform is a powerful tool for analyzing practical discrete-time systems, often represented by linear difference equations. Solving a higher-order difference equation requires knowledge of the input signal and the initial conditions up to one term less than the order of the equation.
The z-transform facilitates handling delayed signals by shifting the signal in the z-domain, which corresponds to delaying the signal in the time domain, and advancing signals by similarly shifting in the...
The z-transform facilitates handling delayed signals by shifting the signal in the z-domain, which corresponds to delaying the signal in the time domain, and advancing signals by similarly shifting in the...
282
Linear Approximation in Time Domain
81
Nonlinear systems often require sophisticated approaches for accurate modeling and analysis, with state-space representation being particularly effective. This method is especially useful for systems where variables and parameters vary with time or operating conditions, such as in a simple pendulum or a translational mechanical system with nonlinear springs.
For a simple pendulum with a mass evenly distributed along its length and the center of mass located at half the pendulum's length,...
For a simple pendulum with a mass evenly distributed along its length and the center of mass located at half the pendulum's length,...
81
Multimachine Stability
150
Multimachine stability analysis is crucial for understanding the dynamics and stability of power systems with multiple synchronous machines. The objective is to solve the swing equations for a network of M machines connected to an N-bus power system.
In analyzing the system, the nodal equations represent the relationship between bus voltages, machine voltages, and machine currents. The nodal equation is given by:
In analyzing the system, the nodal equations represent the relationship between bus voltages, machine voltages, and machine currents. The nodal equation is given by:
150


