查看: 1959  |  回复: 4

gyt3133

铁虫 (小有名气)

[求助] 求助:matlab解17个未知数,17个方程。

方程组如下,为17个未知数、17个方程,按理能求出其解。但提示:“Warning: Explicit solution could not be found.”不知为何,求助于各位大侠如何用matlab求解这种复杂方程组,或者推荐一个比较适合计算这种方程组的数学软件。
[a,b,c,d,e,f,g,h,i,j,k,l,m,n,o,p,q]=solve('0.4*0.03*0.01-0.99*a-0.99*54/102*i=0','0.62*0.02*0.01-0.99*b-0.99*(96/144*j+192/280*k+48/62*l+48/60*m+96/128*n)=0','0.23*0.3*0.01-0.99*c-0.99*28/60*o=0','0.23*1.6*0.01-0.99*d-0.99*(55/71*p+55/87*q)=0','0.81*0.4*0.01-0.99*e-0.99*(24/280*k+12/60*m)=0','0.75*0.001*0.01-0.99*f-0.99*14/62*l=0','0.95*0.1*0.01-0.99*g-0.99*(32/87*q+64/280*k+32/128*n)=0','0.97*0.001*0.01-0.99*h-0.99*(48/102*i+48/144*j+32/60*o+16/71*p)=0','1.09^2*0.37^3*(0.03+a)^2*(0.001+h)^3-2.72*10^-16=0','0.88^2*0.37^3*(0.02+b)^2*(0.001+h)^3-9.85*10^-15=0','1.34*0.37^2*(0.3+c)*(0.001+h)^2-2.38*10^-6=0','0.92*0.37*(1.6+d)*(0.001+h)-0.0216=0','0.92*1.04*(1.6+d)*(0.1+g)-0.84=0','0.88^4*1.17^2*1.04^2*(0.02+b)^4*(0.4+e)^2*(0.1+g)^2-1.64*10^-5=0','0.88*1.04*(0.02+b)*(0.001+f)-4.32*10^-4=0','0.88*1.17*(0.02+b)*(0.4+e)-0.639=0','0.88^2*1.04*(0.02+b)^2*(0.1+g)-5.35*10^-8=0','a','b','c','d','e','f','g','h','i','j','k','l','m','n','o','p','q')
在此,小弟拜谢各位高手!!!!!
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

鹰击长空
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

1stopt软件。
MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
2楼2013-05-07 08:55:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
dbb627: 金币+2, 专家考核, 谢谢指导 2013-05-08 09:52:01
gyt3133: 金币+3, ★★★很有帮助 2013-05-15 12:52:06
1:解析解不一定有,可以试试数值解
2:公示里有多处“10^-8”,是否应该写成“10^(-8)”

结果似乎不唯一,下面两组数值解供参考:
1:
a: 0.150336062520483
i: -0.283739163369587
b: -0.0198581738427485
j: 160218.463398046
k: -674651.464563644
l: -14.7339175704366
m: 267258.426372111
n: 189347.365315984
c: -0.299998949885594
o: 0.644348400594776
d: -1.58194654621143
p: -237006.335460672
q: 290418.722106044
e: 4375.58638947457
f: 3.32702122069614
g: 48.529278391489
h: 3.51383049202942

2:
a: -0.0299922858939372
i: 0.0568810518037125
b: -0.0199740454762181
j: 769.634667070334
k: -279543.909652118
l: -80.5373223843463
m: 245.622378547792
n: 254720.048302251
c: -0.299999963589788
o: 0.644350571334418
d: -1.60336151007992
p: -1056.29114435214
q: 1296.87069161139
e: 23911.7853386278
f: 18.1858545657712
g: -261.270247897996
h: -18.8778594365364
3楼2013-05-07 09:03:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gyt3133

铁虫 (小有名气)

引用回帖:
3楼: Originally posted by dingd at 2013-05-07 09:03:20
1:解析解不一定有,可以试试数值解
2:公示里有多处“10^-8”,是否应该写成“10^(-8)”

结果似乎不唯一,下面两组数值解供参考:
1:
a: 0.150336062520483
i: -0.283739163369587
b: -0.019858173842748 ...

谢谢帮助,有个不专业的问题,还望海涵。请问你在计算的时候用的数值解,具体的方法是什么(插值法、有限元),或者用了什么软件?
鹰击长空
4楼2013-05-07 11:47:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gyt3133

铁虫 (小有名气)

引用回帖:
2楼: Originally posted by 月只蓝 at 2013-05-07 08:55:54
1stopt软件。

好,学习用用,暂时不知道能不能
鹰击长空
5楼2013-05-07 11:58:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 gyt3133 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 极限拉扯,如何抉择 +6 Isak0320 2026-10-12 9/450 2026-10-12 21:13 by Isak0320
[基金申请] 入职后第一次管课题组的专利缴费,对账才发现官费比别人多掏了一倍 +6 13108017953 2026-10-08 6/300 2026-10-12 20:07 by w3w3y3
[有机交流] OpenAI又攻克数学难题了,离AI做有机合成还有多远? +4 wangyl_123 2026-10-10 8/400 2026-10-12 18:09 by wangyl_123
[教师之家] 课题组会帮你通过考核么? +3 好运来~9527 2026-10-10 6/300 2026-10-12 14:55 by 简单化xn
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:34 by XQNyDCa1Tk5n
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:26 by XQNyDCa1Tk5n
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 22:26 by XQNyDCa1Tk5n
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:15 by XQNyDCa1Tk5n
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 22:11 by XQNyDCa1Tk5n
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:06 by XQNyDCa1Tk5n
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 14:25 by VOMexwnDzOFZ
[教师之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 (EPI+-1)(金币-50) +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 14:05 by VOMexwnDzOFZ
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:51 by VOMexwnDzOFZ
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:45 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:54 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:34 by VOMexwnDzOFZ
[考博] 招收博士生一名 +3 mandolin 2026-10-10 3/150 2026-10-11 10:21 by 不是山谷!
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 qR2AdCWTVtkw 2026-10-09 4/200 2026-10-11 03:14 by M0zrYVSAnxpE
[有机交流] 各位大佬,路线设计求助 +6 zapen 2026-10-08 8/400 2026-10-10 15:59 by zapen
[有机交流] 我应该选哪条合成路线让实验成功率高一些? +5 jing816y 2026-10-08 9/450 2026-10-10 09:05 by zapen
信息提示
请填处理意见