Question

求问这个matlab解方程组哪里错了。。。运行之后不提示错误，但是就是没有办法运行处结果啊。

RT
求问啊！太感谢了！没什么米额。。能帮我看看真的是太感谢了。

clear;
%clc;
syms rho_33 rho_22 rho_13 rho_32 rho_12 rho_31 rho_23 rho_21
eqn1='gamma_32*rho_33-gamma_21*rho_22-(i/2)*(Omegac*rho_23-(conj(Omegac))*rho_32)=0';
eqn2='-(gamma_31+gamma_32)*rho_33+(i/2)*(Omegap*rho_13+Omegac*rho_23)-(i/2)*((conj(Omegap))*rho_31+(conj(Omegac))*rho_32)=0'; eqn3='-((gamma_31+gamma_32)/2-i*Deltap)*rho_13-(i/2)*((conj(Omegap))*(1-rho_22-rho_33)+(conj(Omegac))*rho_12)+(i/2)*(conj(Omegap))*rho_33=0'; eqn4='-((gamma_31+gamma_32)/2+i*Deltap)*rho_31+(i/2)*(Omegap*(1-rho_22-rho_33)+Omegac*rho_21)-(i/2)*Omegap*rho_33=0';
eqn5='-((gamma_31+gamma_32+gamma_21)/2+i*Deltac)*rho_32+(i/2)*(Omegap*rho_12+Omegac*rho_22)-(i/2)*Omegac*rho_33=0';
eqn6='-((gamma_31+gamma_32+gamma_21)/2-i*Deltac)*rho_23-(i/2)*((conj(Omegap))*rho_21+(conj(Omegac))*rho_22)+(i/2)*(conj(Omegac))*rho_33=0';
eqn7='-((gamma_21+gamma_ph)/2-i*(Deltap-Deltac))*rho_12-(i/2)*Omegac*rho_13+(i/2)*(conj(Omegap))*rho_32=0';
eqn8='-((gamma_21+gamma_ph)/2+i*(Deltap-Deltac))*rho_21+(i/2)*(conj(Omegac))*rho_31-(i/2)*Omegap*rho_23=0';

s=solve(eqn1,eqn2,eqn3,eqn4,eqn5,eqn6,eqn7,eqn8);

rho_33=s.rho_33;
rho_22=s.rho_22;
rho_13=s.rho_13;
rho_32=s.rho_32;
rho_12=s.rho_12;
rho_31=s.rho_31;
rho_21=s.rho_21;
rho_23=s.rho_23;

gamma_31=1e7; Omegap=10*gamma_31; Omegac=20*gamma_31; gamma_32=0.5*gamma_31;
gamma_21=0.01*gamma_31; gamma_ph=0.01*gamma_31; Deltac=0.3*gamma_31;

rho_22=subs(rho_22);
rho_33=subs(rho_33);
rho_13=subs(rho_13);
rho_32=subs(rho_32);
rho_12=subs(rho_12);

rho_31=subs(rho_31);
rho_23=subs(rho_23);
rho_21=subs(rho_21);

DeltaRatioGamma=-10: 0.02: 10;
Deltap=DeltaRatioGamma*gamma_31;

rho_22New=subs(rho_22);
rho_33New=subs(rho_33);
rho_13New=subs(rho_13);
rho_32New=subs(rho_32);
rho_12New=subs(rho_12);

rho_31New=subs(rho_31);
rho_23New=subs(rho_23);
rho_21New=subs(rho_21);

rho_11=1-rho_22-rho_33;

hold on;
plot(DeltaRatioGamma, real(rho_11), 'b');
hold on;
plot(DeltaRatioGamma, imag(rho_11), 'b:');
hold on;
plot(DeltaRatioGamma, real(rho_22New), 'y');
hold on;
plot(DeltaRatioGamma, imag(rho_22New), 'y:');
hold on;
plot(DeltaRatioGamma, real(rho_33New), 'r');
hold on;
plot(DeltaRatioGamma, imag(rho_33New), 'r:');
hold on;
plot(DeltaRatioGamma, real(rho_13New), 'b+');
hold on;
plot(DeltaRatioGamma, imag(rho_13New), 'b:+');
hold on;
plot(DeltaRatioGamma, real(rho_32New), 'g');
hold on;
plot(DeltaRatioGamma, imag(rho_32New), 'g:');
hold on;
plot(DeltaRatioGamma, real(rho_12New), 'g+');
hold on;
plot(DeltaRatioGamma, imag(rho_12New), 'g:+');

grid;

xlabel('\Delta_{p}/\gamma_{31}');
ylabel('\rho');
legend('Re\rho_{11}', 'Im\rho_{11}','Re\rho_{22}', 'Im\rho_{22}','Re\rho_{33}', 'Im\rho_{33}', 'Re\rho_{13}', 'Im\rho_{13}', 'Re\rho_{32}', 'Im\rho_{32}','Re\rho_{12}', 'Im\rho_{12}');
每次运行之后都是强制把进程关闭，每次运行都是卡在了solve那句。后面应该没有问题。

Answer 1

我个晕 这么长。。。检查一下，那个地方是不是用了汉字（非英文）的空格之类的，

追问

额，可不可以帮我就看看方程的部分，我已经弄了几天了，还是没有弄完成啊。太感谢了啊。
提示的是Warning: Explicit solution could not be found.

Answer 2

运行的时间是长些，提示里说，“没有显函数解。”

追问

也就是这个方程没有解？

追答

是的，至少在Matlab里没有解。