24小时热门版块排行榜

返回列表

米欧米欧

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 48
散金: 50
红花: 1
帖子: 19
在线: 11.2小时
虫号: 3583297
注册: 2014-12-08
专业: 光学

[求助] 关于解25个线性方程组，急求大神们~~~ 已有2人参与

本人matlab小白，急求大神帮忙，
这个关于a11 a12 a13 a14 a15 a21 a22 a23 a24 a25 a31 a32 a33 a34 a35 a41 a42 a43 a44 a45 a51 a52 a53 a54 a55 的25个含有虚部的线性方程组，想求a31,
目前用这个代码算出来是空解。感激不尽~~~~

代码如下：

syms a11 a12 a13 a14 a15 a21 a22 a23 a24 a25 a31 a32 a33 a34 a35
syms a41 a42 a43 a44 a45 a51 a52 a53 a54 a55 deltap Omega2 Omega3
delta1=0;
delta2=0;
delta3=0;
r1=6*1000000;
r2=6*1000000;
r3=6*1000000;
r4=6*1000000;
r5=6*1000000;
Omegap=0.1*1000000;
r6=0.006*1000000;
r7=0.006*1000000;
Omega1=0.3*1000000;
P11=r6*a22+r1*a33+r7*a44+i*(-a13*Omegap+a31*Omegap);
P21=-((r6*a21)/2)+i*((-delta1+deltap)*a21+a31*Omega1+a51*Omega2-a23*Omegap);
P31=-(1/2)*(r1+r2+r3)*a31+i*(deltap*a31+a21*Omega1+a11*Omegap-a33*Omegap);
P41=-((r7*a41)/2)+i*((-delta1+delta2-delta3+deltap)*a41+a51*Omega3-a43*Omegap);
P51=-(1/2)*(r4+r5)*a51+i*((-delta1+delta2+deltap)*a51+a21*Omega2+a41*Omega3-a53*Omegap);
P12=-((r6*a12)/2)+i*((delta1-deltap)*a12-a13*Omega1-a15*Omega2+a32*Omegap);
P22=-r6*a22+r2*a33+r4*a55+i*(-a23*Omega1+a32*Omega1-a25*Omega2+a52*Omega2);
P32=-(1/2)*(r1+r2+r3+r6)*a32+i*(deltap*a32-(-delta1+deltap)*a32+a22*Omega1-a33*Omega1-a35*Omega2+a12*Omegap);
P42=-(1/2)*(r6+r7)*a42+i*((delta1-deltap)*a42+(-delta1+delta2-delta3+deltap)*a42-a43*Omega1-a45*Omega2+a52*Omega3);
P52=-(1/2)*(r4+r5+r6)*a52+i*((delta1-deltap)*a52+(-delta1+delta2+deltap)*a52-a53*Omega1+a22*Omega2-a55*Omega2+a42*Omega3);
P13=-(1/2)*(r1+r2+r3)*a13+i*(-deltap*a13-a12*Omega1-a11*Omegap+a33*Omegap);
P23=-(1/2)*(r1+r2+r3+r6)*a23+i*(-deltap*a23+(-delta1+deltap)*a23-a22*Omega1+a33*Omega1+a53*Omega2-a21*Omegap);
P33=(-r1-r2-r3)*a33+i*(a23*Omega1-a32*Omega1+a13*Omegap-a31*Omegap);
P43=-(1/2)*(r1+r2+r3+r7)*a43+i*(-deltap*a43+(-delta1+delta2-delta3+deltap)*a43-a42*Omega1+a53*Omega3-a41*Omegap);
P53=-(1/2)*(r1+r2+r3+r4+r5)*a53+i*(-deltap*a53+(-delta1+delta2+deltap)*a53-a52*Omega1+a23*Omega2+a43*Omega3-a51*Omegap);
P14=-((r7*a14)/2)+i*((delta1-delta2+delta3-deltap)*a14-a15*Omega3+a34*Omegap);
P24=-(1/2)*(r6+r7)*a24+i*((-delta1+deltap)*a24-(-delta1+delta2-delta3+deltap)*a24+a34*Omega1+a54*Omega2-a25*Omega3);
P34=-(1/2)*(r1+r2+r3+r7)*a34+i*(deltap*a34-(-delta1+delta2-delta3+deltap)*a34+a24*Omega1-a35*Omega3+a14*Omegap);
P44=r3*a33-r7*a44+r5*a55+i*(-a45*Omega3+a54*Omega3);
P54=-(1/2)*(r4+r5+r7)*a54+i*((-delta1+delta2+deltap)*a54-(-delta1+delta2-delta3+deltap)*a54+a24*Omega2+a44*Omega3-a55*Omega3);
P15=-(1/2)*(r4+r5)*a15+i*((delta1-delta2-deltap)*a15-a12*Omega2-a14*Omega3+a35*Omegap);
P25=-(1/2)*(r4+r5+r6)*a25+i*((-delta1+deltap)*a25-(-delta1+delta2+deltap)*a25+a35*Omega1-a22*Omega2+a55*Omega2-a24*Omega3);
P35=-(1/2)*(r1+r2+r3+r4+r5)*a35+i*(deltap*a35-(-delta1+delta2+deltap)*a35+a25*Omega1-a32*Omega2-a34*Omega3+a15*Omegap);
P45=-(1/2)*(r4+r5+r7)*a45+i*((delta1-delta2-deltap)*a45+(-delta1+delta2-delta3+deltap)*a45-a42*Omega2-a44*Omega3+a55*Omega3);
P55=a55+a11+a22+a33+a44-1;
x=solve(P11,P12,P13,P14,P15,P21,P22,P23,P24,P25,P31,P32,P33,P34,P35,P41,P42,P43,P44,P45,P51,P52,P53,P54,P55);
a31=x.a31

回复此楼