[jianbaoliao]

程序如下：
[x,d]=eig(ik,m); %ik和im分别为结构的刚度矩阵和质量矩阵
d=diag(sqrt(d));
for i=1:cn %cn为结构的层数，即质点数
[d1(i),j]=min(d); %请问这句的j怎么取值，是不是循环啊xgd(:,i)=x(:,j);
d(j)=max(d)+1; 这句又怎么跳出循环呢？
end