| 查看: 986 | 回复: 4 | |||
[求助]
关于解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 |
回复此楼» 猜你喜欢
关于filecode，很负责任的告诉大家
已经有7人回复
-
【2027博士申请】纳米药物递送方向
已经有5人回复
-
关于豆爷回答的JTJC与%2F数量
已经有7人回复
-
2026国自然放榜时间
已经有3人回复
-
国基金的申报应该改成非等额制，评价高的钱多评价低的钱少，但是增加资助率
已经有7人回复
-
好奇怪的filecode
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有5人回复
FMStation
至尊木虫 (知名作家)
- 应助: 591 (博士)
- 贵宾: 0.03
- 金币: 18701.7
- 红花: 97
- 帖子: 8891
- 在线: 1493.2小时
- 虫号: 2400059
- 注册: 2013-04-04
- 专业: 计算机应用技术
2楼2016-08-23 07:41:11
3楼2016-08-24 15:38:21
Mr__Right
专家顾问 (著名写手)
-
专家经验: +31 - 应助: 317 (大学生)
- 金币: 14456.3
- 散金: 500
- 红花: 54
- 帖子: 2716
- 在线: 950.6小时
- 虫号: 1972612
- 注册: 2012-09-04
- 性别: GG
- 专业: 应用数学方法
- 管辖: 程序语言
4楼2016-08-25 08:27:27
Mr__Right
专家顾问 (著名写手)
-
专家经验: +31 - 应助: 317 (大学生)
- 金币: 14456.3
- 散金: 500
- 红花: 54
- 帖子: 2716
- 在线: 950.6小时
- 虫号: 1972612
- 注册: 2012-09-04
- 性别: GG
- 专业: 应用数学方法
- 管辖: 程序语言
5楼2016-08-25 08:45:21
100