此函数为显式和隐式方法（以及可选的自适应步长控制）实现了固定步长 Runge-Kutta 求解器。 该函数支持显式和隐式方法，也支持嵌入式方法。 任何 Runge-Kutta 方法都可以通过指定它们的屠夫表来简单地添加。 算法本身是通用的并且相对紧凑。 目前实施了大约 34 种方法。 MATLAB 的 ODE 求解器都是可变步长的，甚至不提供以固定步长运行的选项。 这是因为与固定步长相比，自适应步长可以使求解器更快、更精确。 但是，有时有充分的理由选择固定步长求解器： - 参数研究（比较不同模型参数的仿真结果） - 计算模拟结果的有限差分雅可比（自适应步长控制会引入明显的噪声） - 执行逐点计算，其中求解器输出和测量数据必须参考相同的时间向量- 具有用于模拟结果和固定计算时间的预分配数组 界面和选项在注释中进行了解释。 有两个例子： 示例 1 使用不同的方法和步长求解阻尼和驱动的谐振

瞬态动力学时间积分步长ΔT设置原则： ΔT 即两个时间点间的时间增量 ，它决定求解的精确度，必须采用相应的值才能得到分析现象。通常在每个循环子步中， ΔT 至少有20个时刻点应是足够的 即 ΔT=1/20(f) f是所关心的最高响应频率，而施加阶梯载荷时，为紧紧跟随载荷的阶跃变化， ΔT也许要小到和1/180f相近 接触分析中ΔT =1/30(fc) fc为接触频率 fc=1/2π(k/m)1/2 m 为有效质量 k 为间隙刚度

此代码的工作方式与 ode45、ode23 等系列非常相似，只是它使用固定步长 RK4 算法。 输入是函数句柄、时间跨度、初始条件和时间步长。 extraparameters 变量可用于将额外信息传递给派生例程，而不是使用全局变量。 next 是一个介于 1 和 100 之间的数字，用于通知用户模拟进度。 如果变量 quat = 'Quat'，则模拟将假设状态为 (x,y,z,q0,q1,q2,q3,u,v,w,p,q,r) 对四元数进行归一化。

变步长simpson积分c语言程序，期中f函数就是是被积函数。如果要改被积函数请改f函数

盲源分离中自然梯度算法批处理方法，他的主要步骤是将以前的一个点一个点的计算转化为多点计算。所谓的批处理方法就是迭代的过程不是用一个点去更新权值了，而是利用许多个点来更新权值。

为了快速准确地计算 GLONASS卫星坐标,通过对不同的卫星状态方程进行比较分析,给出一个准确的 GLONASS卫星运动状态方程。以该状态方程为基础,提出了古典形式、龙格库塔、基尔公式3种不同形式的龙格库塔轨道积分方法,比较分析得出3种积分方法精度相当,但基尔公式的舍入误差较小。然后,以基尔公式为基础,给出了定步长和自动选择步长 2种不同的积分方法的详细步骤。讨论了定步长的选择以及自动选择步长收敛值的取值,得到其最佳值在20～3 0之间。最后,利用 GLONASS的广播星历对基于基尔公式的自动选择积分步长

变步长的LMS自适应滤波算法matlab程序

变步长的LMS自适应滤波算法matlab程序

一种FxLMS自适应滤波算法

自启动 Radau IIA 积分器至少在速度和精度方面可与其他最好的积分器相媲美，而且通常更胜一筹，尤其是在高精度方面。 在这里，对于 e = 0.1 的偏心率，实现了从 t0 = 0 [s] 到 t = 3600 [s] 的归一化二体问题的积分。