%%K=globalStiffness; %K为刚度矩阵 900900矩阵Mc=globalMass; %Mc为质量矩阵 900900矩阵C=0.889Mc+0.000278K; %C为阻尼矩阵i=(-1)^(1/2);%应该怎样将w=1：80利用for循环输出80个D 使w=1时D=D1 ； w=2时 D=D2...D=K+iwC-w²Mc; %D为动态刚度矩阵
问题和相关代码如上