kzong.m

matlab编写结构计算振型分解法

function [K]=kzong(k0,dn)
%求解总刚K---输入变量：层间刚度矩阵k0
%                      结构层数dn
%            返回变量：总刚K
K=zeros(dn);
for i=1:dn-1
    K(i,i)=k0(i)+k0(i+1);
    K(i,i+1)=-k0(i+1);
    K(i+1,i)=-k0(i+1);
end 
K(dn,dn)=k0(dn);