仅供参考
clc
clear
syms a11 a21 b11 b21 w
eq1=w^2*a11-1.9388*w*b11-852.64*a11+1.51256*w*b21+45403.08*(a21-a11)*((a21-a11)^2+(b21-b11)^2)+2.1316;
eq2=w^2*b11-1.9388*w*a11-852.64*b11+1.51256*w*a21+45403.08*(b21-b11)*((a21-a11)^2+(b21-b11)^2)-0.010658*w;
eq3=w^2*a21-15.1256*w*(b21-b11)-454030.8*(a21-a11)*((a21-a11)^2+(b21-b11)^2)-0.0025*w^2;
eq4=w^2*b21-15.1256*w*(a21-a11)-454030.8*(b21-b11)*((b21-b11)^2+(a21-a11)^2);
[a11s,a21s,b11s,b21s]=solve(eq1,eq2,eq3,eq4,a11,a21,b11,b21 );
Ax1=(a11s.^2+b11s.^2).^0.5;
Ax2=(a11s.^2+b11s.^2).^0.5;
w1=(0:0.1:60);
Ax1s=double(subs(Ax1,w,w1));
figure
plot(w1,Ax1s)
