应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » MATLAB/Mathematica » 关于解25个线性方程组，急求大神们~~~
51/1返回列表
查看: 968  |  回复: 4
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

米欧米欧

新虫 (初入文坛)

[求助] 关于解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
回复此楼

» 猜你喜欢
1楼 2016-08-22 11:57:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

FMStation

至尊木虫 (知名作家)

【答案】应助回帖

★ ★
感谢参与，应助指数 +1
米欧米欧(月只蓝代发): 金币+2, 感谢应助 2017-08-14 09:44:53
Warning: Explicit solution could not be found.
Please check your P11.. P55
赞一下(1人) 回复此楼
2楼2016-08-23 07:41:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

这里有颗瓜

铁虫 (小有名气)
如果方程没问题的话，建议尝试数值近似解

发自小木虫Android客户端
赞一下(2人) 回复此楼
3楼2016-08-24 15:38:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Mr__Right

专家顾问 (著名写手)

【答案】应助回帖

感谢参与，应助指数 +1
齐次线性方程组系数矩阵满秩，只有平凡解。
赞一下(1人) 回复此楼
文章乃身外之物,要多考虑编辑、审稿人和读者的感受。
4楼2016-08-25 08:27:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Mr__Right

专家顾问 (著名写手)
★ ★
米欧米欧(月只蓝代发): 金币+2, 感谢回复。 2017-08-14 09:45:04
引用回帖:
4楼: Originally posted by Mr__Right at 2016-08-25 08:27:27
齐次线性方程组系数矩阵满秩，只有平凡解。

我看错了。如果系数矩阵中必须包含变量，则结果非常复杂。最好把系数矩阵中的变量都改为常量再求解
赞一下(1人) 回复此楼
文章乃身外之物,要多考虑编辑、审稿人和读者的感受。
5楼2016-08-25 08:45:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 米欧米欧 的主题更新
51/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 小木虫没落了，除了祈祷帖子，几乎看不到有价值的帖子 +14 botar 2026-07-14 15/750 2026-07-19 00:15 by beefly
[基金申请] 面上的水很深 +13 布布和一二 2026-07-18 18/900 2026-07-18 20:37 by 日月天地
[基金申请] 国社科通讯评审 +4 潜水冠军123 2026-07-18 5/250 2026-07-18 19:58 by chenmodeyu
[基金申请] 5个%2F，非固定段22个字母 +8 jaasxulu 2026-07-15 9/450 2026-07-18 17:40 by dxcgbm
[基金申请] 生命口C13面上会评时间 +7 luckinging 2026-07-14 11/550 2026-07-17 20:43 by 虫卜1988
[教师之家] 内心匮乏 +16 水冰月月野兔 2026-07-12 16/800 2026-07-17 19:25 by otani
[基金申请] 生命口面上会评结束了吗大家有知道结果吗？现在还没有结果正常不？ +6 碧水00 2026-07-16 9/450 2026-07-17 15:21 by jaasxulu
[论文投稿] 投英文期刊如何查重 10+5 不知道取什么名� 2026-07-15 6/300 2026-07-17 09:07 by bobvan
[基金申请] 不要再数国自然申请书的 filecode 的分隔符个数了 +14 sadpencil 2026-07-12 21/1050 2026-07-17 00:15 by wangze12014
[有机交流] 伯胺上甲基求助 100+4 002tmumu 2026-07-14 7/350 2026-07-16 19:47 by 002tmumu
[论文投稿] 跨出版社商投稿 20+5 倪好520 2026-07-13 5/250 2026-07-16 17:24 by 无视军团
[基金申请] 生生命学部会评 +4 wuke100666 2026-07-15 4/200 2026-07-16 00:12 by 啦啦小妮子
[基金申请] H口，面上，时间戳7月7日。 +3 plumage5 2026-07-15 3/150 2026-07-15 21:06 by 6543yes
[基金申请] 会评 +5 无名者登山 2026-07-15 5/250 2026-07-15 18:15 by 杠就是你对
[基金申请] E口会评 +11 布布和一二 2026-07-13 13/650 2026-07-15 16:34 by Kittylucky
[基金申请] E口青C会评结束了吗？ 40+3 小小书虫c 2026-07-14 9/450 2026-07-15 16:04 by 安静的知了
[论文投稿] 农业机械学报投稿周期 +3 G0DLIN 2026-07-14 7/350 2026-07-15 09:26 by xs74101122
[基金申请] 国自然面上D口祈祷 +4 monkeynana 2026-07-14 4/200 2026-07-15 09:14 by 白水煮
[基金申请] 没收到消息，再次陪跑 +3 天空星辰fly 2026-07-14 3/150 2026-07-14 19:14 by 6543yes
[基金申请] 明天E口面上会评 +8 布布和一二 2026-07-12 9/450 2026-07-13 22:59 by QTSorange
信息提示
关闭
请填处理意见
关闭