 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1837 | 回復(fù): 6 | ||
| 【懸賞金幣】回答本帖問題，作者Claire3741將贈(zèng)送您 5 個(gè)金幣 | ||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||
Claire3741新蟲 (初入文壇)
|
[求助]
求解高次方程組（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ù)此樓» 猜你喜歡
重慶交通大學(xué)光子學(xué)微結(jié)構(gòu)與器件課題組招收2026年碩士研究生信息
已經(jīng)有2人回復(fù)
-
一志愿鄭大材料學(xué)碩298分，求調(diào)劑
已經(jīng)有6人回復(fù)
-
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)?
已經(jīng)有117人回復(fù)
-
尋合作：應(yīng)力腐蝕多尺度模擬
已經(jīng)有6人回復(fù)
-
考研交流
已經(jīng)有0人回復(fù)
-
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有0人回復(fù)
-
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有0人回復(fù)
-
幫我的英語口語老師找學(xué)生
已經(jīng)有3人回復(fù)
Claire3741
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 58.6
- 帖子: 28
- 在線: 4.2小時(shí)
- 蟲號(hào): 23386947
- 注冊(cè): 2020-08-18
- 專業(yè): 電力電子學(xué)
3樓2020-08-24 23:38:23
獨(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客戶端 |
2樓2020-08-23 08:32:37
獨(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ì)算模擬
4樓2020-08-25 08:28:46
獨(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)助回帖
|
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； |
5樓2020-08-25 08:29:27
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 304求調(diào)劑 +3 | 52hz~~ 2026-03-05 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 一志愿武理085601專碩347分求調(diào)劑 +4 | 啊歐歐歐 2026-03-04 | 5/250 |
|
|
|
[考研]
材料學(xué)學(xué)碩308分/本科東北大學(xué)/一志愿西北工業(yè)大學(xué)/
 5+3
|
蘇堯幺幺 2026-03-03 | 8/400 |
|
|
|
[考研] 求調(diào)劑 +3 | 泡了個(gè)椒 2026-03-04 | 4/200 |
|
|
|
[考研] 302材料工程求調(diào)劑 +7 | Doleres 2026-03-01 | 8/400 |
|
|
|
[考研] 一志愿南京大學(xué)資源與環(huán)境專業(yè)283分求調(diào)劑 +3 | 追云?? 2026-03-04 | 3/150 |
|
|
|
[考研] 0817化學(xué)工程與技術(shù)312分求調(diào)劑 +3 | T123 tt 2026-03-04 | 3/150 |
|
|
|
[考研] 一志愿西交化工專碩288專業(yè)課93求調(diào)劑 +6 | 好運(yùn)好運(yùn)接接 2026-03-04 | 6/300 |
|
|
|
[基金申請(qǐng)] 沒有青基直接申請(qǐng)面上，感覺自己瘋了 +6 | kevin63t 2026-03-02 | 12/600 |
|
|
|
[考研] 264求調(diào)劑 +8 | 26調(diào)劑 2026-03-03 | 8/400 |
|
|
|
[考研] 282求調(diào)劑 +5 | 2103240126 2026-03-02 | 8/400 |
|
|
|
[考研] 085600求調(diào)劑 +4 | LRZZZZZZ 2026-03-02 | 6/300 |
|
|
|
[考研] 一志愿中科大能動(dòng)297求調(diào)劑，本科川大 +3 | 邵11 2026-03-03 | 3/150 |
|
|
|
[考研] 求調(diào)劑 +4 | Guo_yuxuan 2026-03-02 | 5/250 |
|
|
|
[考研] 求調(diào)劑 +7 | repeatt?t 2026-02-28 | 7/350 |
|
|
|
[考研] 調(diào)劑 +3 | 13853210211 2026-03-02 | 4/200 |
|
|
|
[考研] 291分工科求調(diào)劑 +9 | science餓餓 2026-03-01 | 10/500 |
|
|
|
[考研] 328求調(diào)劑 +3 | aaadim 2026-03-01 | 5/250 |
|
|
|
[考研] 295復(fù)試調(diào)劑 +3 | 簡(jiǎn)木ChuFront 2026-03-01 | 3/150 |
|
|
|
[考研] 307求調(diào)劑 +4 | 73372112 2026-02-28 | 6/300 |
|