Newmark-beta法解微分方程的MATLAB程序

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