24小时热门版块排行榜    

CyRhmU.jpeg
南方科技大学公共卫生及应急管理学院2026级博士研究生招生报考通知(长期有效)
查看: 1730  |  回复: 7

zhaoyuan0426

至尊木虫 (著名写手)

[求助] Matlab求解多元多次方程组问题

求解的方程组如下:

1-(1/2*w1)^2/(w1+w2+w3+w4+w5)^2-1.34 /w1=0

1-(w1+1/2*w2)^2/(w1+w2+w3+w4+w5)^2-1.34 /w2=0

1-(w1+w2+1/2*w3)^2/(w1+w2+w3+w4+w5)^2-1.34 /w3=0

1-(w1+w2+w3+1/2*w4)^2/(w1+w2+w3+w4+w5)^2-1.34 /w4=0

1-(w1+w2+w3+w4+1/2*w5)^2/(w1+w2+w3+w4+w5)^2-1.34 /w5=0

我编的代码是:

syms w1 w2 w3 w4 w5;

eq1=1-(1/2*w1)^2/(w1+w2+w3+w4+w5)^2-1.34 /w1;

eq2=1-(w1+1/2*w2)^2/(w1+w2+w3+w4+w5)^2-1.34 /w2;

eq3=1-(w1+w2+1/2*w3)^2/(w1+w2+w3+w4+w5)^2-1.34 /w3;

eq4=1-(w1+w2+w3+1/2*w4)^2/(w1+w2+w3+w4+w5)^2-1.34 /w4;

eq5=1-(w1+w2+w3+w4+1/2*w5)^2/(w1+w2+w3+w4+w5)^2-1.34 /w5;

[w1 w2 w3 w4 w5]=solve(eq1, eq2, eq3, eq4, eq5)



为什么一直都是busy……
回复此楼

» 猜你喜欢

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

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

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
解析解不一定都能求出来的,用数值解试试。不过你的方程不加约束的话有多解:
1:
w1: -13.4538798456478
w2: 14.6773440575791
w3: 9.43746865097371
w4: -7.37476184493982
w5: 3.12888243346332

2:
w1: 1.35118332229718
w2: -18.0945159516444
w3: 19.1633085135825
w4: 1.65711370077759
w5: 3.34892665707608

3:
w1: 13.3106935818788
w2: -0.546696806714209
w3: -10.6339439083153
w4: 1.62593467173962
w5: 3.26198033012798
2楼2013-05-28 15:39:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhaoyuan0426

至尊木虫 (著名写手)

引用回帖:
2楼: Originally posted by dingd at 2013-05-28 15:39:10
解析解不一定都能求出来的,用数值解试试。不过你的方程不加约束的话有多解:
1:
w1: -13.4538798456478
w2: 14.6773440575791
w3: 9.43746865097371
w4: -7.37476184493982
w5: 3.12888243346332

2:
w ...

你好谢谢解答,但是我想要的都是正数解,就是w1-w5都要为正
3楼2013-05-28 16:30:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhaoyuan0426

至尊木虫 (著名写手)

引用回帖:
2楼: Originally posted by dingd at 2013-05-28 15:39:10
解析解不一定都能求出来的,用数值解试试。不过你的方程不加约束的话有多解:
1:
w1: -13.4538798456478
w2: 14.6773440575791
w3: 9.43746865097371
w4: -7.37476184493982
w5: 3.12888243346332

2:
w ...

具体怎么修改我的程序代码???
4楼2013-05-28 16:31:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhaoyuan0426

至尊木虫 (著名写手)

引用回帖:
2楼: Originally posted by dingd at 2013-05-28 15:39:10
解析解不一定都能求出来的,用数值解试试。不过你的方程不加约束的话有多解:
1:
w1: -13.4538798456478
w2: 14.6773440575791
w3: 9.43746865097371
w4: -7.37476184493982
w5: 3.12888243346332

2:
w ...

加个约束条件:
w1>0
w2>0
w3>0
w4>0
w5>0
w5>w4>w3>w2>w1>0
5楼2013-05-28 16:34:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

内容已删除
6楼2013-05-28 16:35:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhaoyuan0426

至尊木虫 (著名写手)

引用回帖:
6楼: Originally posted by dingd at 2013-05-28 16:35:45
都为正的话:

w1: 1.34614806140827
w2: 1.39902153592043
w3: 1.53006056380616
w4: 1.841052868708
w5: 3.84327694392641

我想看看代码?可以吗
7楼2013-05-28 18:17:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★
zhaoyuan0426: 金币+5, ★★★很有帮助 2013-05-29 10:33:12
1stOpt求解代码:
CODE:
ParameterDomain = [0,];
Function 1-(1/2*w1)^2/(w1+w2+w3+w4+w5)^2-1.34 /w1=0;
1-(w1+1/2*w2)^2/(w1+w2+w3+w4+w5)^2-1.34 /w2=0;
1-(w1+w2+1/2*w3)^2/(w1+w2+w3+w4+w5)^2-1.34 /w3=0;
1-(w1+w2+w3+1/2*w4)^2/(w1+w2+w3+w4+w5)^2-1.34 /w4=0;
1-(w1+w2+w3+w4+1/2*w5)^2/(w1+w2+w3+w4+w5)^2-1.34 /w5=0;

8楼2013-05-29 08:41:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zhaoyuan0426 的主题更新
信息提示
请填处理意见