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 - l

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

非线性单自由度结构的Newmark-Beta方法，matlab程序，Chopra, A.K. (2014) Dynamics of Structures, Theory and Applications to Earthquake Engineering: Global Edition. Pearson Education. 4th ed.

建立了关于中间变量的辅助微分方程，将其引入以重写色散手性介质的本构关系。 然后使用Newmark-Beta方法获得了离散的时域递归公式来计算中间变量。 最后，开发了用于分散手性介质的Newmark-FDTD配方。 仿真了色散手性平板的共极化和交叉极化以及色散手性球的双基地雷达截面的反射和透射系数，数值结果证明了该方法的有效性和可行性。