此代码的工作方式与 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] 的归一化二体问题的积分。

从猜测点可以获得下一个点，沿着由 Armijo 的步长指定的方向和距离测量。 在下一步计算要最小化的目标函数，如果它满足Armijo条件，那么它的值必须小于起点，它返回步的最优值。

 风力发电系统的输出功率受外界因数和风速的影响。为了提高小型风机发电机组的转换效率，文中采用一种最大功率优化跟踪算法。以变步长来跟踪风速变化，当功率变化小于一个阈值时停止搜索，来实现最大功率收索的快速性和稳定性。以带齿轮箱6 kW的鼠笼异步式风力发电并网为基础，通过Matlab/Simulink软件仿真结果证实此种方法与定步长爬山法相比，能够达到快速跟踪最大功率点和避免达到最大功率点附近的时候频繁波动。