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
2019-12-21 19:40:44
1KB
MATLAB
1