| 7 | 1/1 | 返回列表 |
| 查看: 811 | 回復(fù): 6 | ||||
[求助]
MATLAB求解復(fù)雜方程組,求助!。
|
|
MATLAB新手。想用MATLAB求解復(fù)雜二元多次方程組。方程式如下所述: d/x=-1.131*[log(a/x)]^3+13.635*[log(a/x)^2]-30.594[[log(a/x)]]+29.267; b/ac=(-1.40557y^3+0.77526y^2+0.15830y-0.06831)[ log(a/x)]^3+(17.93006y^3-9.22091y^2-2.37733n+0.86295[[log(a/x)] ^2+(-79.99715y^3+40.55620y^2+9.00157y-2.54543[[log(a/x)]+(122.65069y^3-63.88418y^2-9.58936y+6.20045) 其中,a,b,c,d為參量,他們的數(shù)值在Excel表格中,x,y為未知數(shù),我想利用MATLAB,求解出上述方程的解。在a,b,c,d參量變化的情況下,求解出x,y的值。 我看到網(wǎng)上說可以使用m文件自定義函數(shù),也可以使用solve語句,但是我用solve語句寫出來如下, syms x y a b c d a= xlsread('20190718.xlsx','sheet2','B2:B6'); b= xlsread('20190718.xlsx','sheet2','C2:C6'); c= xlsread('20190718.xlsx','sheet2','A2:A6'); d= xlsread('20190718.xlsx','sheet2','D2 6');[x,y]=solve('-1.131*[log(a/x)]^3+13.635*[log(a/x)^2]-30.594[[log(a/x)]]+29.267==d/x,b/ac==(-1.40557y^3+0.77526y^2+0.15830y-0.06831)[ log(a/x)]^3+(17.93006y^3-9.22091y^2-2.37733n+0.86295[[log(a/x)] ^2+(-79.99715y^3+40.55620y^2+9.00157y-2.54543[[log(a/x)]+(122.65069y^3-63.88418y^2-9.58936y+6.20045)',x,y) x=simple(x),y=simple(y) 系統(tǒng)提示'"對無效變量名或定義數(shù)字的字符串的支持將在將來的版本中刪除。創(chuàng)建符號表達(dá)式,首先創(chuàng)建符號變量,然后對其使用操作。"" 我不知道該怎么辦了,求助,非常感謝! |
主管區(qū)長 (職業(yè)作家)
![]() |
專家經(jīng)驗(yàn): +1059 |
|
你的a b c d每個(gè)都是包含5個(gè)數(shù)的數(shù)組 你是想對于每組a b c d都解一次方程組 一共是解5次 對吧? 提供幾點(diǎn)思路: 1 不要用solve函數(shù) 這個(gè)函數(shù)不適合求復(fù)雜方程的數(shù)值解 改用fsolve函數(shù) 2 由于a b c d有5組 如果想一次求出所有的解 需要寫個(gè)循環(huán) 比如用for i=1:5 上述兩點(diǎn) 我以前寫過類似的代碼回在別人的帖子里 具體翻一下我回復(fù)記錄 3 MATLAB的fsolve函數(shù)需要提供初值 如果求解困難 可以試試1stOpt軟件來解 發(fā)自小木蟲Android客戶端 |

| 7 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0856材料與化工353分求調(diào)劑 +4 | NIFFFfff 2026-03-09 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 接受26屆調(diào)劑生 +20 | 豬豬豬毛 2026-03-06 | 20/1000 |
|
|
[考研] 一志愿中科大材料與化工。353分能調(diào)劑到哪? +3 | MayUxw1 2026-03-08 | 3/150 |
|
|
[考研] 一志愿武理314求調(diào)劑 +4 | ( ̄~ ̄;) 2026-03-08 | 5/250 |
|
|
[考研] 醫(yī)學(xué)檢驗(yàn)專業(yè)26考研求調(diào)劑 +3 | 琪ya 2026-03-03 | 5/250 |
|
|
[考研] 324求調(diào)劑 +9 | wxz2 2026-03-03 | 13/650 |
|
|
[考研] 22408 275分求調(diào)劑 +3 | 宇智波比 2026-03-03 | 4/200 |
|
|
[考研] 0817化學(xué)工程與技術(shù)312分求調(diào)劑 +7 | T123 tt 2026-03-04 | 7/350 |
|
|
[考研] 化學(xué)工程求調(diào)劑 +12 | 化工人999 2026-03-04 | 12/600 |
|
|
[考研] 346分材料求調(diào)劑 +5 | snow_反季節(jié)版 2026-03-07 | 5/250 |
|
|
[考研] 2026調(diào)劑】考試A區(qū)0703化學(xué)類323分 誠求接收 +3 | 卷柏卷柏 2026-03-05 | 4/200 |
|
|
[考研] 物理學(xué)求調(diào)劑 +3 | sx宋 2026-03-05 | 3/150 |
|
|
[考研] 085600材料與化工 298 調(diào)劑 +11 | 小西笑嘻嘻 2026-03-03 | 11/550 |
|
|
[考研] 268求調(diào)劑 +4 | 劉合華 2026-03-05 | 4/200 |
|
|
[考研] 287求調(diào)劑 +3 | 看看我. 2026-03-05 | 6/300 |
|
|
[考研] 080500材料科學(xué)與工程 +13 | 202114020319 2026-03-03 | 13/650 |
|
|
[考研] 求調(diào)劑 +7 | 博斯特525 2026-03-04 | 7/350 |
|
|
[考研] 武理材料工程302調(diào)劑 +3 | Doleres 2026-03-03 | 6/300 |
|
|
[考研] 293求調(diào)劑 +4 | 是樂渝哇 2026-03-03 | 4/200 |
|
|
[考研] 化工專碩調(diào)劑 +4 | 利好利好. 2026-03-03 | 7/350 |
|