 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 7 | 1/1 | 返回列表 |
| 查看: 1829 | 回復: 6 | ||
| 【懸賞金幣】回答本帖問題，作者Claire3741將贈送您 5 個金幣 | ||
[求助]
求解高次方程組（5次，matlab&1stopt） 已有1人參與
|
||
|
我需要解一個五元方程組，首先使用matlab中的fsolve，我所設置的迭代次數(shù)為10000次，但是每次跑到400次就停了，而且還會報“solver stopped prematurely.fsolve stopped because it exceeded the iteration limit, options.MaxIterations=400(the default value).”我不知道原因。后面考慮到fsolve初值選取比較難，考慮用1stopt去解，仍然解不出來，我想著沒有精確解，有個看得過去的近似解也行啊。但是始終沒解。有沒有人能幫忙看看為什么前面兩種方法解不出來，或者是提供一些解答的新思路啊，謝謝大家了！ 1、matlab代碼 x = sym('x', [1,2,3,4,5],'positive'); % 生成符號變量向量 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ù)組成的表達式 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ù)組成的表達式 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出錯圖 |
回復此樓» 猜你喜歡
德國亥姆霍茲Hereon中心招收兩位醫(yī)用鎂合金腐蝕與LPSO相變方向2026公派博士生
已經(jīng)有0人回復
-
推薦一款可以AI輔助寫作的Latex編輯器SmartLatexEditor，超級好用，推薦試試
已經(jīng)有11人回復
-
物理學I論文潤色/翻譯怎么收費?
已經(jīng)有289人回復
-
2026-CJ開始申報了
已經(jīng)有1人回復
-
西安電子科學大學杭州研究院劉麗香教授招收智能多模態(tài)傳感器和微型儲能器件方向博士
已經(jīng)有13人回復
-
重慶交通大學光子學微結(jié)構(gòu)與器件課題組招收2026年碩士研究生信息
已經(jīng)有1人回復
-
一志愿鄭大材料學碩298分，求調(diào)劑
已經(jīng)有6人回復
-
尋合作：應力腐蝕多尺度模擬
已經(jīng)有3人回復
-
考研交流
已經(jīng)有0人回復
獨孤神宇
版主 (知名作家)
- 應助: 490 (碩士)
- 貴賓: 0.008
- 金幣: 31016.3
- 散金: 802
- 紅花: 122
- 沙發(fā): 1
- 帖子: 5600
- 在線: 856.5小時
- 蟲號: 3522474
- 注冊: 2014-11-06
- 性別: GG
- 專業(yè): 機械動力學
- 管轄: 計算模擬
2樓2020-08-23 08:32:37
3樓2020-08-24 23:38:23
獨孤神宇
版主 (知名作家)
- 應助: 490 (碩士)
- 貴賓: 0.008
- 金幣: 31016.3
- 散金: 802
- 紅花: 122
- 沙發(fā): 1
- 帖子: 5600
- 在線: 856.5小時
- 蟲號: 3522474
- 注冊: 2014-11-06
- 性別: GG
- 專業(yè): 機械動力學
- 管轄: 計算模擬
4樓2020-08-25 08:28:46
獨孤神宇
版主 (知名作家)
- 應助: 490 (碩士)
- 貴賓: 0.008
- 金幣: 31016.3
- 散金: 802
- 紅花: 122
- 沙發(fā): 1
- 帖子: 5600
- 在線: 856.5小時
- 蟲號: 3522474
- 注冊: 2014-11-06
- 性別: GG
- 專業(yè): 機械動力學
- 管轄: 計算模擬
【答案】應助回帖
|
5個參數(shù)，需要正版軟件才能計算。 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ù)組成的表達式 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
|
謝謝你的回復，我又重做了實驗，用fsolve還是沒辦法解出來，我沒1stopt正版軟件，能不能麻煩你幫忙跑一下啊。我這個其實是6個參數(shù)，5個方程，我想弄兩組試看看能不能接出來。 1、6個參數(shù)，5個方程 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ù)組成的表達式 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個參數(shù)，5個方程 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ù)組成的表達式 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; 十分感謝！ |
6樓2020-08-28 23:16:08
獨孤神宇
版主 (知名作家)
- 應助: 490 (碩士)
- 貴賓: 0.008
- 金幣: 31016.3
- 散金: 802
- 紅花: 122
- 沙發(fā): 1
- 帖子: 5600
- 在線: 856.5小時
- 蟲號: 3522474
- 注冊: 2014-11-06
- 性別: GG
- 專業(yè): 機械動力學
- 管轄: 計算模擬
【答案】應助回帖
|
目標函數(shù)值(最小): 0.00886947526736046 x1: 0.963945105231679 x2: 0.134939569392818 x3: 0.858017777832085 x4: 1.79769482276588 x5: 1.67904225670153 x6: 0.0517955828557939 ******************************************** 目標函數(shù)值(最小): 0.0088694752543157 x1: 0.829371852036702 x2: 0.0295612334196613 x3: 0.858017701890971 x4: 1.79745238455023 x5: 3.69943415831409 |
7樓2020-08-29 08:57:34
相關版塊跳轉(zhuǎn)
我要訂閱樓主 Claire3741 的主題更新
| 7 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 312求調(diào)劑 +5 | 醋精華了一下發(fā)?/a> 2026-03-03 | 7/350 |
|
|---|---|---|---|---|
|
|
[考研] 一志愿西交大材料學碩 346 求調(diào)劑 +3 | zju51 2026-03-04 | 3/150 |
|
|
|
[考研] 331求調(diào)劑 +3 | zzZ&zZ 2026-03-03 | 3/150 |
|
|
[考博] 26超級電容器申博 +4 | dhdjdjend 2026-02-25 | 4/200 |
|
|
|
[考博] 26申博 求博導 +3 | 愛讀書的小帥 2026-02-28 | 5/250 |
|
|
|
[考研] 中國林科院林化所（南京）2026年招收化學/材料/環(huán)境工程等背景碩士研究生3名 +3 | realstar2006 2026-02-27 | 3/150 |
|
|
|
[考研] 290求調(diào)劑 +9 | ErMiao1020 2026-03-02 | 9/450 |
|
|
|
[考研] 085600材料與化工 298 調(diào)劑 +3 | 小西笑嘻嘻 2026-03-03 | 3/150 |
|
|
|
[考研] 283求調(diào)劑 +9 | 鹿沫笙 2026-03-02 | 9/450 |
|
|
|
[考研] 299求調(diào)劑 +5 | kkcoco25 2026-03-02 | 9/450 |
|
|
|
[考研]
材料工程專碩283求調(diào)劑
 5+8
|
，��？ 2026-03-02 | 10/500 |
|
|
|
[考研] 324求調(diào)劑 +4 | wxz2 2026-03-03 | 5/250 |
|
|
|
[考研] 11408，學碩276求調(diào)劑 +3 | 崔wj 2026-03-02 | 5/250 |
|
|
|
[考研] 285求調(diào)劑 +9 | 滿頭大汗的學生 2026-02-28 | 9/450 |
|
|
|
[考研] 265分求調(diào)劑不調(diào)專業(yè)和學校有行學上就 +6 | 禮堂丁真258 2026-02-28 | 9/450 |
|
|
|
[考研] 求調(diào)劑 +3 | 熬夜的貓頭鷹 2026-03-02 | 3/150 |
|
|
|
[考研] 281求調(diào)劑 +5 | 2026計算機_誠心 2026-03-01 | 8/400 |
|
|
|
[考研] 材料類求調(diào)劑 +11 | wana_kiko 2026-02-28 | 14/700 |
|
|
|
[基金申請] 成果系統(tǒng)訪問量大，請一小時后再嘗試。---NSFC啥時候好哦，已經(jīng)兩天這樣了 +4 | NSFC2026我來了 2026-02-28 | 4/200 |
|
|
|
[基金申請]
剛錄用，沒有期刊號，但是在線可看的論文可以放為代表作嗎
10+3
|
arang1 2026-03-01 | 3/150 |
|