海神之光上传的视频是由对应的完整代码运行得来的，完整代码皆可运行，亲测可用，适合小白； 1、从视频里可见完整代码的内容 主函数：main.m； 调用函数：其他m文件；无需运行 运行结果效果图； 2、代码运行版本 Matlab 2019b；若运行有误，根据提示修改；若不会，私信博主； 3、运行操作步骤 步骤一：将所有文件放到Matlab的当前文件夹中； 步骤二：双击打开main.m文件； 步骤三：点击运行，等程序运行完得到结果； 4、仿真咨询 如需其他服务，可私信博主； 4.1 博客或资源的完整代码提供 4.2 期刊或参考文献复现 4.3 Matlab程序定制 4.4 科研合作

三维空间车轨耦合动力学程序：基于Newmark-Beta法的车辆轨道耦合动力学MATLAB代码实现，已嵌入轨道不平顺激励。,根据翟书编写的三维空间车轨耦合动力学程序 通过newmark-beta法求解的车辆-轨道空间耦合动力学matlab代码 已在代码里面加入轨道不平顺激励使用即可，无需动脑 ,翟书编写;三维空间车轨耦合动力学程序;Newmark-beta法;车辆-轨道空间耦合动力学Matlab代码;轨道不平顺激励。,翟书编写的三维空间车轨耦合动力学程序——Newmark-beta法求解车辆轨道耦合动力学MATLAB代码

求解粘滞阻尼质量弹簧系统的常微分运动方程，该系统在受到任意外力时会表现出非线性，周期性，力-位移行为。 导数通过使用[1]中介绍的Newmark-beta方法在时间上隐式积分。 然后，将控制方程的完整数值形式表示为残差，并使用从[2]中获得的Newton-Raphson算法找到解决方案。 弹簧的材料特性可以是线性的也可以是非线性的，因为牛顿法应该以任何一种方式收敛。 通过使用主干曲线定义恢复力与位移之间的非线性关系。 数据作为一组必须严格为正的横坐标对输入。 然后，用户提供的骨架通过水平轴和垂直轴反射，从而使关系变为各向同性，即，无论弹簧是拉伸还是压缩，该关系都是相同的。 假定第一个数据点为初始弹性屈服点。 如果系统在屈服时开始恢复，则主干将从其原始屈服点移至当前位移。 这种移动产生了应变硬化的非常基本的形式。 但是，它无法捕捉到循环应变硬化和降解的更实际效果。 改善弹簧滞后性能的

利用newmark-β的方法求解振动方程，设置步长，积分常数

Newmark-beta法解微分方程的MATLAB程序。 function [q,v,a]=newmark(M,C,K,F,q0,v0,dt,nt) % newmark-beta method % [q,v,a]=newmark(M,C,K,F,q0,v0,dt,nt) % obtain the response of the dynamic system % M - mass matrix % K - stiffness matrix % C - damping matrix % F - loads matrix(nt columns) % q0 - initial displacement % v0 - initial velocity % dt - interval(time step) % nt - number of sampling points % [q,v,a] - disp,velocity,acceleration