Euler 方法被用作 Heun 方法的基础。 Euler 的方法使用函数在区间开始处的切线作为函数在区间上的斜率的估计,假设如果步长小,误差就会小。 然而,即使使用极小的步长,在大量的步长上,误差开始累积,估计与实际功能值不同。 例子: 输入 x 的初始值,即 x0: 0 输入 y 的初始值,即 y0:0.5 输入 x 的最终值:2 输入步长h:0.2 xy 0.000 0.500 0.200 0.630 0.400 0.833 0.600 1.124 0.800 1.523 1.000 2.054 1.200 2.746 1.400 3.634 1.600 4.762 1.800 6.181 2.000 7.957 >>
2022-05-17 11:04:55 1KB matlab
1
matlab的欧拉方法代码时间积分方案 这是一个有限元研究项目,其目的是针对梁单元在瞬态分析中研究时域中不同数值积分方案的阻尼效应。 特别是对于纵向自由度的冲击载荷。 研究的主要集成方案是Backward-Euler,Newmark方法,居中方案,一些Runge-kutta方案和广义alpha方案。 将所进行的仿真与商业代码Abaqus进行比较,并研究了不同的属性,例如准确性,稳定性,鲁棒性。这些方法和主要的有限元方法都是在Matlab中实现和编程的。 还对图形用户界面进行了编程,以允许用户轻松选择数值积分方法并查看位移结果。 GUI显示了在2m光束的末端进行3500 N冲击模拟的结果,持续时间为1e-4秒。 这些参数可以在 Newmark.m 和 dynamics.m 请享用 !
2022-05-11 16:24:51 88KB 系统开源
1
一旋LLG解算器 使用RK4,RK2或Euler方法求解LLG。 要更改集成方法,请转到main.c然后在#include "method"更改method找到所需的集成方法。 动画是由cria_anim.py创建的。 要选择动画时长,请转到cria_anim.py然后将调用方法faz(duration)更改为所需的时长。 在动画中,红色箭头是外部磁场,黑色箭头是自旋矢量(超过磁化范数的磁化矢量)。 外部字段随时间变化,要设置其变化方式,请转到method.c ,其中method是您正在使用的集成方法,并更改函数void* calccampo(void *k) 。 时间参数是k_ 。 例: 假设外部字段为H =(0.5t,t,0),则该函数为: void* calccampo(void *k){ double k_ = *(double*) k; H->x = 0.5 * k_; H->y
2022-04-08 15:43:57 7KB C
1
matlab的欧拉方法代码运动学/动力学 在为不同机器人系统的运动学和动力学进行课程分配时开发的代码。 内容: 用于5自由度机械手的FK / IK 对于正向运动学,所有关节位置和变换矩阵均使用DH方法进行计算。 对于逆运动学,该代码使用运动学去耦,并检查了关节极限,可行性并考虑了其他边缘情况。 2D摆动力学 dynamics.m使用Euler-Lagrange方法推导摆模型的闭合运动方程。 该符号表达式已复制到solution.m脚本中。 直接运行solution.m求解ODE并绘制wrt时间的联合变量值(假设扭矩= 0)。 它解决的另一种动力学问题是,当已知联合变量wrt time的函数时,运行torque.m来计算沿该轨迹的所需扭矩。 该代码解决了系统DOF = 2的动力学问题,但是很容易进行修改以解决高阶摆动力学模型。 但是,由于符号变量数量的增加,MATLAB中的微分运算会花费更多时间,因此可能需要事先进行一些矩阵运算以简化计算。
2022-03-18 11:21:06 9KB 系统开源
1
matlab的欧拉方法代码欧拉法 包含用于使用Euler方法和Modified Euler方法求解一阶ODE的代码。 理查森的外推法用于达到更高的精度要求。 所有代码都是使用Matlab的文件扩展名.m为Matlab编写的
2022-03-12 04:37:19 198KB 系统开源
1
matlab的欧拉方法代码PDE_Matlab_Code_solution- 作为论文的一部分,我必须用数值​​方法求解PDE。 PDE代表与微生物相关的项目在微观尺度上的梁状结构的位移。 它不是特别适合于解析数学解决方案。 BC取自相关文献,并采用Euler前向有限差​​分法。 PDE读取u_t =(beta / 4)*(s ^ 2 -1) u_ss-(1 / eta) (u_ssss)。
2022-02-28 11:36:13 2KB 系统开源
1
常微分方程的数值解法(C语言版本):包括用改进欧拉方法求解10x(1-y),四阶龙格—库塔方法求解10x(1-y)
1
matlab的欧拉方法代码神经元网络模型 动态耦合激发大脑中神经元的模型以产生复杂的网络同步 该项目提供了Matlab代码来模拟以下情况: 一个发射神经元细胞,使用三种不同的模型。 x个激发神经元的网络,使用静态耦合矩阵耦合在一起。 x个激发神经元的网络,使用基于神经元细胞之间突触模型的动态耦合功能耦合在一起。 在single_neuron_models目录中,运行着一些程序来模拟单个激发神经元的行为。每个程序的顶部都有一些示例运行。 OneNeuronTau.m:基于Tau常数的简单模型 OneNeuronIzhInF.m:伊兹凯维奇着名的“整合并发射”神经元模型 OneNeuronExpInF.m:更复杂的指数神经元模型 在Neuron_network_models目录中,运行程序“ Neuron Simulations”以打开一个GUI,该GUI允许配置和显示神经元网络。 NeuronSimulations.m:包含所有不同模型和行为的GUI NeuronNetworkTau.m:对通过耦合矩阵连接在一起的“ Tau激发”神经元网络进行建模。 使用正向Euler方法计算每个神经元
2021-12-12 12:38:03 461KB 系统开源
1
matlab的欧拉方法代码二维热扩散有限差分法 该项目的目的是使用有限差分法求解二维热方程。 这是一个MATLAB代码,可以解决不同的材料,例如(铜,铝,银等),或者允许用户通过输入导热系数,比热和密度来添加自己的材料。 可以应用许多随时间固定的“狄利克雷条件”的边界条件。 它还解决了印版的稳态温度,并告诉用户该印版达到该稳态所需的时间,并且用户选择了误差容限。 解决方案的精度将主要取决于在网格划分和开始迭代求解之前也可以选择的x和y方向上的节点数。 该代码可以用两种方式(Euler和二阶Runge-Kutte)以及具有中心有限差分的空间导数求解方程的时间导数部分。最后,在求解之后,时间上的图形模拟似乎表明了热量如何扩散到整个过程中在所选时间间隔内的板块。
2021-12-12 09:16:59 2.97MB 系统开源
1
matlab的欧拉方法代码我的硕士论文代码工作 标题:基于系统的微分平坦度特性和输入整形,具有悬浮载荷的四旋翼的轨迹滑移控制 摘要:在货运应用中对自动驾驶飞机的需求很高的情况下,这项工作提出了一种运动控制和轨迹生成解决方案,该解决方案由四轴转子组成,该系统具有悬架负载,旨在控制飞机的位置并减轻负载摇摆。 首先,使用Newton-Euler和Euler-Lagrange方法推导系统的动态模型,并将其分为两个部分:与机器人高度和偏航角关联的完全驱动子系统,以及与其他状态变量关联的欠驱动子系统。四旋翼。 每个子系统都由一个滑模控制器控制,该滑模控制器在Lyapunov的意义上被证明是稳定的,可以将系统驱动到滑模表面并停留在滑模表面上。 Routh-Hurwitz稳定性判据证明,与欠驱动子系统相关的滑动表面是局部稳定的,这是因为获得了一些使控制过程更容易的控制参数约束规则。 最后,提出了一种新的轨迹生成结构来抑制负载平衡,该结构包括建立负载的多项式轨迹,对其施加输入整形并利用系统的微分平坦度特性来计算飞机的期望状态。 在仿真中针对点对点轨迹和不同持续时间,对控制器和设计的轨迹生成器进行了测试。
2021-12-06 11:31:12 105KB 系统开源
1