kweewu

銅蟲 (小有名氣)

應助: 0 (幼兒園)
金幣: 2905.9
帖子: 295
在線: 84.2小時
蟲號: 1568403
注冊: 2012-01-07
性別: GG
專業(yè): 物理電子學

[求助] 求助大神解一個非線性方程組已有1人參與

方程組形式如下：

x^2+y^2-2xycosω1=a^2
y^2+z^2-2yzcosω2=b^2
z^2+x^2-2zxcosω3=c^2

除去x，y，z外，其他均為常數(shù)
cosω1=0.37461
cosω2=0.24192
cosω3=0.17365

a=7.2
b=14.2
c=18.1

小弟菜鳥一個，使用excel規(guī)劃求解得出x=9.759132，y=3.746134，z=15.8171，
目標函數(shù)選擇為f=（x^2+y^2-2xycosω1-a^2）^2+(y^2+z^2-2yzcosω2-b^2)^2+(z^2+x^2-2zxcosω3=c^2)^2
最終的函數(shù)值f應該為0，但當取以上解時f=3333.691；
所以應該還有更精確的解。
因為沒用過MATLAB，所以希望求得能解該方程組的MATLAB代碼。
謝謝！

eq.png

回復此樓

PEMFC

1樓 2015-08-11 11:38:36

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

月只藍

主管區(qū)長 (職業(yè)作家)

專家經(jīng)驗: +1059
應助: 1712 (講師)
貴賓: 8.888
金幣: 68128.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線: 3291.4小時
蟲號: 1122189
注冊: 2010-10-14
專業(yè): 宇宙學
管轄: 計算模擬區(qū)

【答案】應助回帖

感謝參與，應助指數(shù) +1

不存在滿足f=0的解。

使得f盡量小的解：
8.063894135668461
0.049339384803544
-14.560697883414663

或者：
  -8.063894144363893
  -0.049339340193223
  14.560697881727988

MATLAB代碼：

CODE:

function solve_eq

clear all;clc

format long

global a b c cos1 cos2 cos3

cos1=0.37461;

cos2=0.24192;

cos3=0.17365;

a=7.2;

b=14.2;

c=18.1;

x0=[-1 0 0];

options=optimset('TolFun',1e-9);

[x fval]=fsolve(@fun,x0,options);

[x;fval]'

x'

function f=fun(u)

global a b c cos1 cos2 cos3

x=u(1);

y=u(2);

z=u(3);

f(1)=x^2+y^2-2*x*y*cos1-a^2;

f(3)=y^2+z^2-2*y*z*cos2-b^2;

f(2)=z^2+x^2-2*z*x*cos3-c^2;

贊一下

回復此樓

» 本帖已獲得的紅花（最新10朵）

kweewu

MATLAB、MS小問題、普通問題請發(fā)帖求助！時間精力有限，恕不接受無償私信求助。

2樓2015-08-11 13:09:22

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

kweewu

銅蟲 (小有名氣)

應助: 0 (幼兒園)
金幣: 2905.9
帖子: 295
在線: 84.2小時
蟲號: 1568403
注冊: 2012-01-07
性別: GG
專業(yè): 物理電子學

送紅花一朵

引用回帖:

2樓: Originally posted by 月只藍 at 2015-08-11 13:09:22
不存在滿足f=0的解。

使得f盡量小的解：
8.063894135668461
0.049339384803544
-14.560697883414663

或者：
  -8.063894144363893
  -0.049339340193223
  14.560697881727988

MATLAB代碼：
...

您好！忘了說了，x，y，z均是正數(shù)，這樣怎樣相應的修改一下代碼呢？

贊一下

回復此樓

PEMFC

3樓2015-08-11 15:46:12

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

月只藍

主管區(qū)長 (職業(yè)作家)

專家經(jīng)驗: +1059
應助: 1712 (講師)
貴賓: 8.888
金幣: 68128.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線: 3291.4小時
蟲號: 1122189
注冊: 2010-10-14
專業(yè): 宇宙學
管轄: 計算模擬區(qū)

【答案】應助回帖

引用回帖:

3樓: Originally posted by kweewu at 2015-08-11 15:46:12
您好！忘了說了，x，y，z均是正數(shù)，這樣怎樣相應的修改一下代碼呢？...

如果要求x y z大于0，我想你的計算結(jié)果已經(jīng)是最優(yōu)了：
x= 9.75912727167921
y= 3.74632253182401
z=15.8171084695522
此時
(x^2+y^2-2*x*y*cos1-a^2)^2+(y^2+z^2-2*y*z*cos2-b^2)^2+(z^2+x^2-2*z*x*cos3-c^2)^2 =3333.69114238958；

贊一下

回復此樓

MATLAB、MS小問題、普通問題請發(fā)帖求助！時間精力有限，恕不接受無償私信求助。

4樓2015-08-11 16:00:32

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

kweewu

銅蟲 (小有名氣)

應助: 0 (幼兒園)
金幣: 2905.9
帖子: 295
在線: 84.2小時
蟲號: 1568403
注冊: 2012-01-07
性別: GG
專業(yè): 物理電子學

引用回帖:

4樓: Originally posted by 月只藍 at 2015-08-11 16:00:32
如果要求x y z大于0，我想你的計算結(jié)果已經(jīng)是最優(yōu)了：
x= 9.75912727167921
y= 3.74632253182401
z=15.8171084695522
此時
(x^2+y^2-2*x*y*cos1-a^2)^2+(y^2+z^2-2*y*z*cos2-b^2)^2+(z^2+x^2-2*z*x*cos3-c^2) ...

謝謝！那如此在您的代碼上要怎樣使x，y，z都是正數(shù)呢？

贊一下

回復此樓

PEMFC

5樓2015-08-11 16:03:34

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

月只藍

主管區(qū)長 (職業(yè)作家)

專家經(jīng)驗: +1059
應助: 1712 (講師)
貴賓: 8.888
金幣: 68128.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線: 3291.4小時
蟲號: 1122189
注冊: 2010-10-14
專業(yè): 宇宙學
管轄: 計算模擬區(qū)

【答案】應助回帖

引用回帖:

5樓: Originally posted by kweewu at 2015-08-11 16:03:34
謝謝！那如此在您的代碼上要怎樣使x，y，z都是正數(shù)呢？...

不再用fsolve函數(shù)了，改用有約束的非線性最優(yōu)化函數(shù)，求目標函數(shù);(x^2+y^2-2*x*y*cos1-a^2)^2+(y^2+z^2-2*y*z*cos2-b^2)^2+(z^2+x^2-2*z*x*cos3-c^2)^2的最小值，然后附加x y z大于0的約束。

贊一下

回復此樓

MATLAB、MS小問題、普通問題請發(fā)帖求助！時間精力有限，恕不接受無償私信求助。

6樓2015-08-11 16:09:25

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

kweewu

銅蟲 (小有名氣)

應助: 0 (幼兒園)
金幣: 2905.9
帖子: 295
在線: 84.2小時
蟲號: 1568403
注冊: 2012-01-07
性別: GG
專業(yè): 物理電子學

引用回帖:

6樓: Originally posted by 月只藍 at 2015-08-11 16:09:25
不再用fsolve函數(shù)了，改用有約束的非線性最優(yōu)化函數(shù)，求目標函數(shù);(x^2+y^2-2*x*y*cos1-a^2)^2+(y^2+z^2-2*y*z*cos2-b^2)^2+(z^2+x^2-2*z*x*cos3-c^2)^2的最小值，然后附加x y z大于0的約束。...

謝謝！我應該如何操作才能采用您的應助？找了一會兒沒找到這個選項o(╯□╰)o

贊一下

回復此樓

PEMFC

7樓2015-08-11 16:26:51

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

kweewu

銅蟲 (小有名氣)

應助: 0 (幼兒園)
金幣: 2905.9
帖子: 295
在線: 84.2小時
蟲號: 1568403
注冊: 2012-01-07
性別: GG
專業(yè): 物理電子學

引用回帖:

改用有約束的非線性最優(yōu)化函數(shù)，附加x y z大于0的約束，這個代碼該怎么寫呢？

贊一下

回復此樓

PEMFC

8樓2015-08-11 16:30:11

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

月只藍

主管區(qū)長 (職業(yè)作家)

專家經(jīng)驗: +1059
應助: 1712 (講師)
貴賓: 8.888
金幣: 68128.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線: 3291.4小時
蟲號: 1122189
注冊: 2010-10-14
專業(yè): 宇宙學
管轄: 計算模擬區(qū)

【答案】應助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
kweewu: 金幣+30, ★★★★★最佳答案, 非常感謝您的耐心幫助！ 2015-08-11 17:07:24

引用回帖:

8樓: Originally posted by kweewu at 2015-08-11 16:30:11
改用有約束的非線性最優(yōu)化函數(shù)，附加x y z大于0的約束，這個代碼該怎么寫呢？...

CODE:

function youhua

clear all;clc

format long

A=[-1 0 0

    0 -1 0

    0 0 -1];

b=[0 0 0];

x0=[1 1 1];

[x fval]=fmincon(@fun,x0,A,b)

function f=fun(u)

cos1=0.37461;

cos2=0.24192;

cos3=0.17365;

a=7.2;

b=14.2;

c=18.1;

x=u(1);

y=u(2);

z=u(3);

f=(x^2+y^2-2*x*y*cos1-a^2)^2+(y^2+z^2-2*y*z*cos2-b^2)^2+(z^2+x^2-2*z*x*cos3-c^2)^2;

計算結(jié)果：
x y z分別=

9.759125101118036 3.746323040238880 15.817108773681586

目標函數(shù)值fval =

3.333691142391446e+003

上述代碼須完整復制進一個新建的m文件中，運行即可，不能在主程序窗口運行。

贊一下

回復此樓

MATLAB、MS小問題、普通問題請發(fā)帖求助！時間精力有限，恕不接受無償私信求助。

9樓2015-08-11 16:47:55

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 一志愿南開大學0710生物學359求調(diào)劑 +3	兔兔兔111223314 2026-03-29	3/150	2026-03-29 23:14 by Evan_Liu
[考研] 085600 286分材料求調(diào)劑 +10	麻辣魷魚 2026-03-27	11/550	2026-03-29 21:52 by nanaliuyun
[考研] 總分293求調(diào)劑 +8	加一一九 2026-03-25	11/550	2026-03-29 19:53 by 無際的草原
[考研] 327求調(diào)劑 +6	汲亦昊 2026-03-29	6/300	2026-03-29 13:40 by peike
[考研] 279求調(diào)劑 +4	蝶舞輕繞 2026-03-29	4/200	2026-03-29 09:45 by laoshidan
[考研] 材料學碩333求調(diào)劑 +11	北道巷 2026-03-24	11/550	2026-03-29 08:31 by fmesaito
[考研] 0856求調(diào)劑 +7	楒桉 2026-03-28	7/350	2026-03-29 08:28 by fmesaito
[考研] 312，生物學求調(diào)劑 +3	小譯同學abc 2026-03-28	3/150	2026-03-28 15:32 by 落�？伤�
[考研] 394求調(diào)劑 +3	好事多磨靜候佳?/a> 2026-03-26	5/250	2026-03-28 14:24 by 唐沐兒
[考研] 352分化工與材料 +5	海納百川Ly 2026-03-27	5/250	2026-03-28 03:39 by fmesaito
[考研] 考研化學308分求調(diào)劑 +10	你好明天你好 2026-03-23	12/600	2026-03-27 14:43 by shangxh
[考研] 一志愿華東理工大學081700，初試分數(shù)271 +6	kotoko_ik 2026-03-23	7/350	2026-03-27 12:29 by 惠州彭于晏
[考研] 316求調(diào)劑 +5	Pigcasso 2026-03-24	5/250	2026-03-27 12:10 by zhshch
[考研] 0703化學一志愿南京師范大學303求調(diào)劑 +3	zzffylgg 2026-03-24	3/150	2026-03-27 10:42 by shangxh
[考研] 324求調(diào)劑 +5	hanamiko 2026-03-26	5/250	2026-03-27 10:33 by wangjy2002
[考研] 324求調(diào)劑 +8	hanamiko 2026-03-26	10/500	2026-03-27 08:06 by hypershenger
[考研] 321求調(diào)劑 +6	Ymlll 2026-03-24	6/300	2026-03-26 20:50 by 不吃魚的貓
[考研] 081200-11408-276學碩求調(diào)劑 +3	崔wj 2026-03-26	3/150	2026-03-26 19:57 by nihaoar
[考研] 求b區(qū)院校調(diào)劑 +4	周56 2026-03-24	5/250	2026-03-25 17:12 by yishunmin
[考研] 材料/農(nóng)業(yè)專業(yè)，07/08開頭均可，過線就行 +3	呵唔哦豁 2026-03-23	4/200	2026-03-23 22:30 by 汪�。�！