| 查看: 853 | 回复: 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 |
回复此楼» 猜你喜欢
材料与化工（0856）304求 B区 调剂
已经有3人回复
-
268求调剂
已经有8人回复
-
265求调剂
已经有14人回复
-
350 本科985求调剂，求老登收留
已经有3人回复
-
【考研调剂】化学专业 281分，一志愿四川大学，诚心求调剂
已经有9人回复
-
初试 317
已经有3人回复
-
279分求调剂 一志愿211
已经有15人回复
-
286求调剂
已经有4人回复
-
279求调剂
已经有4人回复
-
求调剂
已经有3人回复
FMStation
至尊木虫 (知名作家)
- 应助: 591 (博士)
- 贵宾: 0.03
- 金币: 18346.2
- 红花: 97
- 帖子: 8891
- 在线: 1489.9小时
- 虫号: 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
