版塊導航: 正在加載中...

登錄注冊

應《網絡安全法》要求，自2017年10月1日起，未進行實名認證將不得使用互聯(lián)網跟帖服務。為保障您的帳號能夠正常使用，請盡快對帳號進行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

北京石油化工學院2026年研究生招生接收調劑公告

返回列表

當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖

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

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

月只藍

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

專家經驗: +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;

計算結果：
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

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

查看全部 9 個回答

月只藍

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

專家經驗: +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

已閱回復此樓關注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

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

月只藍

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

專家經驗: +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，我想你的計算結果已經是最優(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

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

查看全部 9 個回答

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

[考研] 張芳銘-中國農業(yè)大學-環(huán)境工程專碩-298 +6	手機用戶 2026-03-26	6/300	2026-03-31 09:43 by jp9609
[考研] 085600材料與化工329分求調劑 +4	這是有太陽哇 2026-03-25	4/200	2026-03-31 09:39 by Huaxue_Wang
[考研] 物理學調劑 +3	小羊36 2026-03-30	3/150	2026-03-31 09:20 by jp9609
[考研] 求調劑 +3	圖鑒212 2026-03-30	3/150	2026-03-31 00:31 by jp9609
[考研] 求調劑，一志愿南京航空航天大學，080500材料科學與工程學碩，總分289分 +9	@taotao 2026-03-29	9/450	2026-03-30 22:29 by 我是小康
[考研] 085701環(huán)境工程求調劑 +11	多久上課 2026-03-27	12/600	2026-03-30 21:21 by 研究僧導導
[考研] 085600，專業(yè)課化工原理，320分求調劑 +5	大饞小子 2026-03-29	5/250	2026-03-30 19:02 by 源_2020
[考研] 070300化學354求調劑 +7	101次希望 2026-03-28	7/350	2026-03-30 13:22 by Iveryant
[考研] 南京大學化學調劑 +10	景隨風 2026-03-29	15/750	2026-03-30 11:21 by limeifeng
[考研] 349求調劑 +6	李木子啊哈哈 2026-03-25	6/300	2026-03-29 12:47 by 無際的草原
[考研] 283求調劑 +3	A child 2026-03-28	3/150	2026-03-28 15:41 by ms629
[考研] 320分，材料與化工專業(yè)，求調劑 +9	一定上岸aaa 2026-03-27	13/650	2026-03-28 15:00 by 神馬都不懂
[考研] 0703化學求調劑，各位老師看看我�。�！ +5	祁祺祺 2026-03-25	5/250	2026-03-27 21:44 by 東方豬豬
[考研] 274求調劑 +17	顧九笙要謙虛 2026-03-24	23/1150	2026-03-27 15:16 by caszguilin
[考研] 305求調劑 +5	哇盧卡庫 2026-03-26	5/250	2026-03-27 14:01 by laoshidan
[考研] 314求調劑 +3	溪云珂 2026-03-26	3/150	2026-03-27 11:20 by sanrepian
[考研] 一志愿鄭大085600，310分求調劑 +5	李瀟可 2026-03-26	5/250	2026-03-27 11:14 by 不吃魚的貓
[考研] 292求調劑 +4	求求了收下我吧?/a> 2026-03-26	4/200	2026-03-27 10:37 by zhshch
[考研] 303求調劑 +6	藍山月 2026-03-25	6/300	2026-03-25 22:47 by 418490947
[考研] 材料與化工304求B區(qū)調劑 +3	邱gl 2026-03-25	3/150	2026-03-25 19:03 by Ainin_

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产 高清 中文字幕,99re热久久亚洲综合精品成人,熟妇 一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

24小時熱門版塊排行榜

kweewu

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

» 收錄本帖的淘帖專輯推薦

» 猜你喜歡

» 本主題相關價值貼推薦，對您同樣有幫助:

月只藍

【答案】應助回帖

月只藍

【答案】應助回帖

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

kweewu

月只藍

【答案】應助回帖

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产高清中文字幕,99re热久久亚洲综合精品成人,熟妇一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

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