论文“Matrix Formulation for Minimum Response of Undamped Structures”参考代码
clc;
clear;
M=[3 0 0;0 2 0;0 0 1];
K=[200 -100 0;-100 200 -100;0 -100 100];
P=[0;0;1];
C=[0;1;0];
G=[-1 0;0 0;0 1];
T=[0 -1;1 0;0 0];
K1=T'*K*G;
M1=T'*M*G;
[V,D,W]=eig(K1,M1);%W左特征向量,V右特征向量
[V0,D0]=eig(K,M);
[row,col]=find(D>0);
lamta=D(row,col);
circle_freq=sqrt(lamta);
%开始灵敏度计算
dK_dx=[1 -1 0;-1 1 0;0 0 0];
dK1_dx=T'*dK_dx*G;
dM1_dx=0;
dlamta_dx=V(:,col)'*(dK1_dx'-lamta*dM1_dx')*W(:,col)/(V(:,col)'*M1'*W(:,col));
dcircle_freq_dx=1/(2*circle_freq)*dlamta_dx;