matlab求解模态,[转载]MATLAB与模态叠加法求解高等动力学问题，两者误差比较

图示表示四层楼的抗剪模型，其剪切刚度系数及楼板质量均表示在图中，在顶层受一水平力的简谐激振力pcos(Ωt)，仅考虑其稳态响应。求顶层的响应函数。

%%%%%%%%%% 运动方程求解 %%%%%% 运动方程求解 %%%%

T=input('Time is ');

tb=input('Buchang is ');

%---------------中心差分法------------------

t=0:tb:T;

M=[1 0 0 0; 0 2 0 0; 0 0 2 0; 0 0 0 3];

K=[800 -800 0 0; -800 2400 -1600 0; 0 -1600 4000

-2400; 0 0 -2400 5600];

a0=(tb^2)1;

a1=(2*tb)1;

a2=2*a0;

u0=[0.00 0.00 0.00 0.00]';v0=[0.00 0.00 0.00

0.00]';ac0=[100.00 0.00 0.00 0.00]';

u00=u0-tb*v0+a2ac0;

%---------------(-tb时刻的位移)------------------

DM=a0*M;

%---------------等效质量阵------------------

n=fix(tbT);

u=zeros(4,n);

%---------------定义４x N的位移矩阵------------------

u(:,1)=u00;u(:,2)=u0;