版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

返回列表

【懸賞金幣】回答本帖問題，作者Claire3741將贈(zèng)送您 5 個(gè)金幣

Claire3741

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 58.6
帖子: 28
在線: 4.2小時(shí)
蟲號(hào): 23386947
注冊(cè): 2020-08-18
專業(yè): 電力電子學(xué)

[求助] 求解高次方程組（5次，matlab&1stopt）已有1人參與

我需要解一個(gè)五元方程組，首先使用matlab中的fsolve，我所設(shè)置的迭代次數(shù)為10000次，但是每次跑到400次就停了，而且還會(huì)報(bào)“solver stopped prematurely.fsolve stopped because it exceeded the iteration limit, options.MaxIterations=400(the default value).”我不知道原因。后面考慮到fsolve初值選取比較難，考慮用1stopt去解，仍然解不出來，我想著沒有精確解，有個(gè)看得過去的近似解也行啊。但是始終沒解。有沒有人能幫忙看看為什么前面兩種方法解不出來，或者是提供一些解答的新思路啊，謝謝大家了！
1、matlab代碼
x = sym('x', [1,2,3,4,5],'positive');
% 生成符號(hào)變量向量
w1=sqrt((x(1)+x(3)*x(4))./(2*x(2)*x(4)*x(5)));
z=((2*x(2)*x(4)+x(3)*x(4)*x(5))./(2*(x(1)+x(3)*x(4))))*w1;
w2=sqrt(1-z^2)*w1;
a=sqrt((1-2*x(5)*z*w1+x(5)^2*w1^2)./(1-z^2));
f1=atan((w2*x(5))./(1-z*w1*x(5)))-atan((sqrt(1-z^2)./(-z)));
% 參數(shù)組成的表達(dá)式
eq1=-0.5*x(4)*a./(x(1)+x(3)*x(4))+0.7959;
eq2=z*w1-1.275;
eq3=w2-2.284;
eq4=f1+0.4204;
eq5=-0.5*x(4)./(x(1)+x(3)*x(4))+0.3504;
f = matlabFunction([eq1;eq2;eq3;eq4;eq5], 'vars', {[x(1) x(2) x(3) x(4) x(5)]});
x0 = [0.5 2.5 1 0.05 0.2];
options=optimset('Display','iter','MaxFunEvals',10000,'TolFun',1e-10);
[x,fval,exitflag,output] = fsolve(f,x0,options);
2、1stopt
const delta=-0.5;
parameter x1[0,5],x2[0,5],x3[0,5],x4[0,1],x5[0,5];
conststr
w1=sqrt((x1+x3*x4)/(2*x2*x4*x5));
z=((2*x2*x4+x3*x4*x5)/(2*(x1+x3*x4)))*w1;
w2=sqrt(1-z^2)*w1;
a=sqrt((1-2*x5*z*w1+x5^2*w1^2)/(1-z^2));
f1=atan((w2*x5)/(1-z*w1*x5))-atan(sqrt(1-z^2)./(-z));
//參數(shù)組成的表達(dá)式
A=delta*x4*a/(x1+x3*x4)+0.7959;
B=z*w1-1.275;
C=w2-2.284;
E=f1+0.4204;
F=delta*x4/(x1+x3*x4)+0.3504；

matlab出錯(cuò)圖

回復(fù)此樓

» 猜你喜歡

2026年循環(huán)經(jīng)濟(jì)功能材料國(guó)際會(huì)議（ICFMCE 2026）已經(jīng)有0人回復(fù)
2026年第五屆電氣、電子與信息工程國(guó)際會(huì)議（ISEEIE 2026）已經(jīng)有0人回復(fù)
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有233人回復(fù)
哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑已經(jīng)有5人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑已經(jīng)有5人回復(fù)
求調(diào)劑已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有0人回復(fù)
請(qǐng)問還有沒有用Latex寫文章的小伙伴們？已經(jīng)有0人回復(fù)
光學(xué)工程學(xué)碩調(diào)劑信息已經(jīng)有23人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)

1樓 2020-08-22 20:49:39

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

獨(dú)孤神宇

版主 (知名作家)

應(yīng)助: 490 (碩士)
貴賓: 0.008
金幣: 31016.3
散金: 802
紅花: 122
沙發(fā): 1
帖子: 5600
在線: 856.5小時(shí)
蟲號(hào): 3522474
注冊(cè): 2014-11-06
性別: GG
專業(yè): 機(jī)械動(dòng)力學(xué)
管轄: 計(jì)算模擬

【答案】應(yīng)助回帖

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

解方程組經(jīng)量不要用符號(hào)運(yùn)算，直接使用數(shù)值計(jì)算

發(fā)自小木蟲Android客戶端

贊一下

回復(fù)此樓

數(shù)值計(jì)算

2樓2020-08-23 08:32:37

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

Claire3741

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 58.6
帖子: 28
在線: 4.2小時(shí)
蟲號(hào): 23386947
注冊(cè): 2020-08-18
專業(yè): 電力電子學(xué)

引用回帖:

2樓: Originally posted by 獨(dú)孤神宇 at 2020-08-23 08:32:37
解方程組經(jīng)量不要用符號(hào)運(yùn)算，直接使用數(shù)值計(jì)算

fsolve求的也是數(shù)值解吧

發(fā)自小木蟲Android客戶端

贊一下

回復(fù)此樓

3樓2020-08-24 23:38:23

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

獨(dú)孤神宇

版主 (知名作家)

應(yīng)助: 490 (碩士)
貴賓: 0.008
金幣: 31016.3
散金: 802
紅花: 122
沙發(fā): 1
帖子: 5600
在線: 856.5小時(shí)
蟲號(hào): 3522474
注冊(cè): 2014-11-06
性別: GG
專業(yè): 機(jī)械動(dòng)力學(xué)
管轄: 計(jì)算模擬

【答案】應(yīng)助回帖

引用回帖:

3樓: Originally posted by Claire3741 at 2020-08-24 23:38:23
fsolve求的也是數(shù)值解吧
...

你1stOpt 代碼寫的不對(duì)。。

x1: 0.0304532132539924
x2: 0.119504502058893
x3: 0.347292096631735
x4: 0.0266306732746401
x5: 0.911618798575995

贊一下

回復(fù)此樓

數(shù)值計(jì)算

4樓2020-08-25 08:28:46

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

獨(dú)孤神宇

版主 (知名作家)

應(yīng)助: 490 (碩士)
貴賓: 0.008
金幣: 31016.3
散金: 802
紅花: 122
沙發(fā): 1
帖子: 5600
在線: 856.5小時(shí)
蟲號(hào): 3522474
注冊(cè): 2014-11-06
性別: GG
專業(yè): 機(jī)械動(dòng)力學(xué)
管轄: 計(jì)算模擬

【答案】應(yīng)助回帖

引用回帖:

3樓: Originally posted by Claire3741 at 2020-08-24 23:38:23
fsolve求的也是數(shù)值解吧
...

5個(gè)參數(shù)，需要正版軟件才能計(jì)算。

Constant delta=-0.5;
Parameter x1=[0,5],x2=[0,5],x3=[0,5],x4=[0,1],x5=[0,5];
ConstStr
w1=sqrt((x1+x3*x4)/(2*x2*x4*x5)),
z=((2*x2*x4+x3*x4*x5)/(2*(x1+x3*x4)))*w1,
w2=sqrt(1-z^2)*w1,
aa=sqrt((1-2*x5*z*w1+x5^2*w1^2)/(1-z^2)),
f1=atan((w2*x5)/(1-z*w1*x5))-atan(sqrt(1-z^2)/(-z));
//參數(shù)組成的表達(dá)式
Function
delta*x4*aa/(x1+x3*x4)+0.7959=0;
z*w1-1.275=0;
w2-2.284=0;
f1+0.4204=0;
delta*x4/(x1+x3*x4)+0.3504=0；

贊一下

回復(fù)此樓

數(shù)值計(jì)算

5樓2020-08-25 08:29:27

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

Claire3741

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 58.6
帖子: 28
在線: 4.2小時(shí)
蟲號(hào): 23386947
注冊(cè): 2020-08-18
專業(yè): 電力電子學(xué)

引用回帖:

5樓: Originally posted by 獨(dú)孤神宇 at 2020-08-25 08:29:27
5個(gè)參數(shù)，需要正版軟件才能計(jì)算。

Constant delta=-0.5;
Parameter x1=,x2=,x3=,x4=,x5=;
ConstStr
w1=sqrt((x1+x3*x4)/(2*x2*x4*x5)),
z=((2*x2*x4+x3*x4*x5)/(2*(x1+x3*x4)))*w1,
w2=sqrt(1-z^2)*w1,
a ...

謝謝你的回復(fù)，我又重做了實(shí)驗(yàn)，用fsolve還是沒辦法解出來，我沒1stopt正版軟件，能不能麻煩你幫忙跑一下啊。我這個(gè)其實(shí)是6個(gè)參數(shù)，5個(gè)方程，我想弄兩組試看看能不能接出來。
1、6個(gè)參數(shù)，5個(gè)方程
Constant delta=-0.225；
Parameter x1=[0,5],x2=[0,5],x3=[0,5],x4=[0,5],x5=[0,5],x6=[0,5];
//生成變量K,F,T,H,D,R
ConstStr
w1=sqrt((x1+x5*x6)/(2*x3*x4*x6)),
z=((2*x6*x4+(x5*x6+x1*x2)*x3)/(2*(x1+x5*x6)))*w1,
w2=sqrt(1-z^2)*w1,
aa=sqrt((1-2*x3*z*w1+x3^2*w1^2)./(1-z^2)),
f1=atan((w2*x3)./(1-z*w1*x3))-atan((sqrt(1-z^2)./(-z)));
//參數(shù)組成的表達(dá)式
Function
delta*x6*aa/(x1+x5*x6)+0.02436=0;
z*w1-1.135=0;
w2-2.3=0;
f1+0.5595=0;
delta*x6./(x1+x5*x6)+0.01119=0;
2、5個(gè)參數(shù)，5個(gè)方程
Constant delta=-0.225；
Parameter x1=[0,5],x2=[0,5],x3=[0,5],x4=[0,5],x5=[0,5];
//生成變量K,F,T,H,D
ConstStr
w1=sqrt((x1+x5*0.05)/(2*x3*x4*0.05)),
z=((2*0.05*x4+(x5*0.05+x1*x2)*x3)/(2*(x1+x5*0.05)))*w1,
w2=sqrt(1-z^2)*w1,
aa=sqrt((1-2*x3*z*w1+x3^2*w1^2)./(1-z^2)),
f1=atan((w2*x3)./(1-z*w1*x3))-atan((sqrt(1-z^2)./(-z)));
//參數(shù)組成的表達(dá)式
Function
delta*0.05*aa/(x1+x5*0.05)+0.02436=0;
z*w1-1.135=0;
w2-2.3=0;
f1+0.5595=0;
delta*0.05./(x1+x5*0.05)+0.01119=0;
十分感謝！

贊一下

回復(fù)此樓

6樓2020-08-28 23:16:08

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

獨(dú)孤神宇

版主 (知名作家)

應(yīng)助: 490 (碩士)
貴賓: 0.008
金幣: 31016.3
散金: 802
紅花: 122
沙發(fā): 1
帖子: 5600
在線: 856.5小時(shí)
蟲號(hào): 3522474
注冊(cè): 2014-11-06
性別: GG
專業(yè): 機(jī)械動(dòng)力學(xué)
管轄: 計(jì)算模擬

【答案】應(yīng)助回帖

引用回帖:

6樓: Originally posted by Claire3741 at 2020-08-28 23:16:08
謝謝你的回復(fù)，我又重做了實(shí)驗(yàn)，用fsolve還是沒辦法解出來，我沒1stopt正版軟件，能不能麻煩你幫忙跑一下啊。我這個(gè)其實(shí)是6個(gè)參數(shù)，5個(gè)方程，我想弄兩組試看看能不能接出來。
1、6個(gè)參數(shù)，5個(gè)方程
Constant delt ...

目標(biāo)函數(shù)值(最小): 0.00886947526736046
x1: 0.963945105231679
x2: 0.134939569392818
x3: 0.858017777832085
x4: 1.79769482276588
x5: 1.67904225670153
x6: 0.0517955828557939

********************************************

目標(biāo)函數(shù)值(最小): 0.0088694752543157
x1: 0.829371852036702
x2: 0.0295612334196613
x3: 0.858017701890971
x4: 1.79745238455023
x5: 3.69943415831409

贊一下

回復(fù)此樓

數(shù)值計(jì)算

7樓2020-08-29 08:57:34

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

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 Claire3741 的主題更新

返回列表

不應(yīng)助 確定回帖應(yīng)助 (注意：應(yīng)助才可能被獎(jiǎng)勵(lì)，但不允許灌水，必須填寫15個(gè)字符以上)

普通表情龍兔虎貓

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

[考研] 268求調(diào)劑 +5	好運(yùn)連綿不絕 2026-03-12	5/250	2026-03-17 10:12 by Sammy2
[論文投稿] 有沒有大佬發(fā)小論文能帶我個(gè)二作 +3	增銳漏人 2026-03-17	4/200	2026-03-17 09:26 by xs74101122
[考研] 304求調(diào)劑 +7	小熊joy 2026-03-14	7/350	2026-03-17 08:53 by 霧散后相遇lc
[考研] 0703化學(xué)336分求調(diào)劑 +3	zbzihdhd 2026-03-15	4/200	2026-03-16 23:52 by zbzihdhd
[考研] 藥學(xué)383 求調(diào)劑 +3	藥學(xué)chy 2026-03-15	4/200	2026-03-16 20:51 by 元子^0^
[考研] 279求調(diào)劑 +3	抓著星星的女孩 2026-03-10	3/150	2026-03-13 23:47 by userper
[考研] 341求調(diào)劑 +4	番茄頭--- 2026-03-10	4/200	2026-03-13 23:12 by JourneyLucky
[考研] 279求調(diào)劑 +3	Dizzy123@ 2026-03-10	3/150	2026-03-13 23:02 by JourneyLucky
[考研] 材料專碩288分求調(diào)劑一志愿211 +4	在家想你 2026-03-11	4/200	2026-03-13 22:49 by JourneyLucky
[考研] 一志愿西南交大，材料專碩317求調(diào)劑 +5	lx8568 2026-03-11	5/250	2026-03-13 21:43 by peike
[考研] 四川大學(xué)085601材料工程專碩初試294求調(diào)劑 +4	祝我們好在冬天 2026-03-11	4/200	2026-03-13 21:39 by peike
[考研] 311求調(diào)劑 +3	冬十三 2026-03-13	3/150	2026-03-13 20:41 by JourneyLucky
[考研] 0703化學(xué)求調(diào)劑 +7	綠豆芹菜湯 2026-03-12	7/350	2026-03-13 17:25 by njzyff
[考研] 求b區(qū)學(xué)校調(diào)劑 +3	周56 2026-03-11	3/150	2026-03-13 16:20 by JourneyLucky
[考研] 274求調(diào)劑 +3	S.H1 2026-03-12	3/150	2026-03-13 15:15 by JourneyLucky
[考研] 0856化學(xué)工程280分求調(diào)劑 +4	shenzxsn 2026-03-11	4/200	2026-03-13 11:55 by ymwdoctor
[考研] 08食品或輕工求調(diào)劑，本科發(fā)表3篇sci一區(qū)top論文，一志愿南師大食品科學(xué)與工程 +3	我是一個(gè)兵， 2026-03-10	3/150	2026-03-13 10:21 by Yuyi.
[考研] 290求調(diào)劑 +3	ADT 2026-03-13	3/150	2026-03-13 10:19 by peike
[考博] 26讀博 +4	Rui135246 2026-03-12	10/500	2026-03-13 07:15 by gaobiao
[考研] 大連大學(xué)化學(xué)專業(yè)研究生調(diào)劑 +3	琪久. 2026-03-10	8/400	2026-03-11 10:02 by 琪久.

24小時(shí)熱門版塊排行榜

[求助] 求解高次方程組（5次，matlab&1stopt） 已有1人參與

» 猜你喜歡

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

[求助] 求解高次方程組（5次，matlab&1stopt）已有1人參與