这个方程解出来后有未​知数w，怎么再把它当​成自变量画图啊，各位​大佬。

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)

s=solve(eq1,eq2,eq3,eq4,a11,a21,b11,b21 )

Ax1=(a11^2+b11^2)^0.5

Ax2=(a11^2+b11^2)^0.5

w=(0:0.1:60)

figure

plot(w,Ax1)

Answer 1

[seqxtyxr]

仅供参考

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)