MATLAB中solve的解为什么全是零，求大神指点

omega = (50:20:2*wnmax/2/pi)'*2*pi; %扫描频率
syms Gamma1

%% 由公式联立计算出gamma和周期
for n=1:length(omega);
Gamma(n,1)=solve(tan(30*omega(n)*R*Gamma1/(Kpd*Nt)-1.5*pi)==(-omega(n)*Gamma1*Gc(n)*Gpd(n)+Hc(n)*(1-Hpd(n)*omega(n)*Gamma1))/(Gc(n)*(1-Hpd(n)*omega(n)*Gamma1)+Hc(n)*Gpd(n)*omega(n)*Gamma1));
Gamma(n)=Gamma1;
end
这里的变量是Gamma1，其余的全是已知的量（矩阵）
pi和0的值带进去做差便可
FRFc=subs(FRFc1,phi,pi)-subs(FRFc1,phi,0);
FRFpd=subs(FRFpd1,phi,pi)-subs(FRFpd1,phi,0);
Hpd=imag(FRFpd);%过程阻尼传递函数的虚部
Gpd=real(FRFpd);%过程阻尼传递函数的实部
Gc=real(FRFc);%结构传递函数的实部
Hc=imag(FRFc);%结构传递函数的虚部
omega = (50:20:2*wnmax/2/pi)'*2*pi; %扫描频率

相关了解……