24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1542  |  回复: 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 的主题更新
信息提示
请填处理意见