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 %

function [qva]=newmark（MCKFq0v0dtnt）
% newmark-beta method
%           [qva]=newmark（MCKFq0v0dtnt）
%           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
%           [qva] - dispvelocityacceleration
    if nt~= length（F）
        disp（‘The size of force matrix not match the number of sampling points ‘）
        return
    end
    dof=length（M）;
    alpha=0