 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 3 | 1/1 | 返回列表 |
| 查看: 1870 | 回復(fù): 2 | |||
[求助]
窮蟲50金！ MATLAB的fsolve初值問題 已有1人參與
|
|
大家好，請(qǐng)教一個(gè)問題 在下用fsolve解非線性方程組，可是初值稍微改變，解就跟著變，我取的初值的絕對(duì)誤差應(yīng)該已經(jīng)小于1。 是不是fsolve這個(gè)算法不好，有沒有什么辦法優(yōu)化呢？ 先謝謝各位啦～  算法見下function f = emuch(x) B = 0.003; C = 0.0027; H = x(1)+75*(B+C)+0.5*x(11); alpha1 = 0.5+6*x(10)-H; beta1 = H*(0.01-1)+(0.5+6*x(10))*(1+0.01); gamma1 = 0.01*H; alpha2 =0.15-x(10); beta2 = 0.3; gamma2 = x(10); alpha3 = 0.3+0.2*x(10)-x(8); beta3 = x(8)*(0.1-1)+(0.3+0.2*x(10))*(1+0.1); gamma3 = x(6)*0.1; Q = 0.01*(1-x(9))+2*x(9)+0.05*x(6); P = 0.1+0.25*x(6); R = 0.3*(20*B+2*C+x(1)+x(11)+0.067*x(10)); J = 0.64*(B+x(1)+0.5*x(11)+x(10)); f1 = 0.05*(2*gamma1/(beta1+sqrt(beta1^2-4*alpha1*gamma1)))-0.1*x(1); f2 = 0.002+0.05*(2*gamma2/(beta2+sqrt(beta2^2-4*alpha2*gamma2)))-Q*x(2); f3 = 0.006+0.02*(2*gamma1/(beta1+sqrt(beta1^2-4*alpha1*gamma1)))-P*x(3); f4 =0.02+0.1*(2*gamma3/(beta3+sqrt(beta3^2-4*alpha3*gamma3)))-(0.01+R/(0.05+x(4)))*x(4); f5 = 50*x(10)*x(12)-(0.06+Q+R/(0.05+x(4)))*x(5); f6 = 50*x(11)*x(12)-(0.06+P+R/(0.05+x(4)))*x(6); f7 = (0.005+0.06*x(8))-0.08*x(7); f8 = x(7)-2*x(8); f9 = (0.04+2*x(8))*(1-x(9))/(1.05-x(9))-J*x(9)/(0.05+x(9)); f10 = x(2)-x(5); f11 = x(3)-x(6); f12 = x(4)-x(5)-x(6); f = [f1 f2 f3 f4 f5 f6 f7 f8 f9 f10 f11 f12]; end |
回復(fù)此樓» 猜你喜歡
求調(diào)劑
已經(jīng)有7人回復(fù)
-
272求調(diào)劑
已經(jīng)有7人回復(fù)
-
352分 化工與材料
已經(jīng)有3人回復(fù)
-
0703化學(xué)
已經(jīng)有6人回復(fù)
-
一志愿南昌大學(xué)324求調(diào)劑
已經(jīng)有3人回復(fù)
-
求調(diào)劑
已經(jīng)有3人回復(fù)
-
266分求材料化工冶金礦業(yè)等專業(yè)的調(diào)劑
已經(jīng)有4人回復(fù)
-
一志愿南師大0703化學(xué) 275求調(diào)劑
已經(jīng)有4人回復(fù)
-
085701環(huán)境工程求調(diào)劑
已經(jīng)有3人回復(fù)
-
266求調(diào)劑
已經(jīng)有10人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- matlab fsolve得出的解全部為復(fù)數(shù)解
已經(jīng)有15人回復(fù)
- 求助 MATLAB-fsolve
已經(jīng)有3人回復(fù)
- 使用matlab 中的fsolve 老提示出錯(cuò)
已經(jīng)有4人回復(fù)
- 求助matlab---fsolve解非線性方程組
已經(jīng)有6人回復(fù)
- 求助！matlab用fsolve函數(shù)求解非線性方程組的問題！
已經(jīng)有19人回復(fù)
- 【求助】matlab中fsolve解的精度控制【已解決】
已經(jīng)有4人回復(fù)
george1990
金蟲 (小有名氣)
- 應(yīng)助: 11 (小學(xué)生)
- 金幣: 1155.3
- 紅花: 2
- 帖子: 130
- 在線: 75.3小時(shí)
- 蟲號(hào): 2142478
- 注冊(cè): 2012-11-22
- 性別: GG
- 專業(yè): 機(jī)械動(dòng)力學(xué)
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
馮且右11: 金幣+10, ★有幫助, 謝謝啦～，不過這個(gè)的精度貌似也不太好，因?yàn)閺膄8顯然能看出x7和x8 2014-06-29 18:20:53
感謝參與，應(yīng)助指數(shù) +1
馮且右11: 金幣+10, ★有幫助, 謝謝啦～，不過這個(gè)的精度貌似也不太好，因?yàn)閺膄8顯然能看出x7和x8 2014-06-29 18:20:53
|
樓主可以使用1stOpt這個(gè)軟件來求解，相對(duì)于matlab求解此類方程時(shí)的，1stOpt省去了fsolve需要的初值的問題，而且精度一般都很高。 下面是程序： Parameter x(1:12); //Constant Constant B = 0.003; Constant C = 0.0027; //ConstStr ConstStr H = x1+75*(B+C)+0.5*x11; ConstStr alpha1 = 0.5+6*x10-H; ConstStr beta1 = H*(0.01-1)+(0.5+6*x10)*(1+0.01); ConstStr gamma1 = 0.01*H; ConstStr alpha2 =0.15-x10; ConstStr beta2 = 0.3; ConstStr gamma2 = x10; ConstStr alpha3 = 0.3+0.2*x10-x8; ConstStr beta3 = x8*(0.1-1)+(0.3+0.2*x10)*(1+0.1); ConstStr gamma3 = x6*0.1; ConstStr Q = 0.01*(1-x9)+2*x9+0.05*x6; ConstStr P = 0.1+0.25*x6; ConstStr R = 0.3*(20*B+2*C+x1+x11+0.067*x10); ConstStr J = 0.64*(B+x1+0.5*x11+x10); ConstStr f1 = 0.05*(2*gamma1/(beta1+sqrt(beta1^2-4*alpha1*gamma1)))-0.1*x1; ConstStr f2 = 0.002+0.05*(2*gamma2/(beta2+sqrt(beta2^2-4*alpha2*gamma2)))-Q*x2; ConstStr f3 = 0.006+0.02*(2*gamma1/(beta1+sqrt(beta1^2-4*alpha1*gamma1)))-P*x3; ConstStr f4 =0.02+0.1*(2*gamma2/(beta3+sqrt(beta3^2-4*alpha3*gamma3)))-(0.01+R/(0.05+x4))*x4; ConstStr f5 = 50*x10*x12-(0.06+Q+R/(0.05+x4))*x5; ConstStr f6 = 50*x11*x12-(0.06+P+R/(0.05+x4))*x6; ConstStr f7 = (0.005+0.06*x8)-0.08*x7; ConstStr f8 = x7-2*x8; ConstStr f9 = (0.04+2*x8)*(1-x9)/(1.05-x9)-J*x9/(0.05+x9); ConstStr f10 = x2-x5; ConstStr f11 = x3-x6; ConstStr f12 = x4-x5-x6; Function f1=0; f2=0; f3=0; f4=0; f5=0; f6=0; f7=0; f8=0; f9=0; f10=0; f11=0; f12=0; 求出來的x系列的值如下： x1: -0.00241408292469573 x2: -1.19963721286631E-5 x3: -0.00163511248524538 x4: -0.00164712704345382 x5: -1.19975380217937E-5 x6: -0.00163512908617766 x7: -7.04705240745422E-20 x8: 1.73072251386771E-18 x9: -0.00331312406855758 x10: -0.0122596988050143 x11: -1.65447541443873 x12: -0.000192121359156606 但是樓主注意，這個(gè)方程組好像答案不唯一，對(duì)于x是否還有別的約束條件呢，樓主可以在看看，這系列值只作為樓主參考。 |
2樓2014-06-28 23:47:34
3樓2014-06-29 18:47:11
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 馮且右11 的主題更新
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085701環(huán)境工程求調(diào)劑 +3 | 多久上課 2026-03-27 | 3/150 |
|
|---|---|---|---|---|
|
|
[考研] 求調(diào)劑推薦 材料 304 +13 | 荷包蛋hyj 2026-03-26 | 13/650 |
|
|
|
[考博] 26申博 +3 | 加油沖�。� 2026-03-26 | 3/150 |
|
|
[考研] 287求調(diào)劑 +10 | land xuxu 2026-03-26 | 10/500 |
|
|
|
[考研] 269專碩求調(diào)劑 +10 | 金恩貝 2026-03-21 | 10/500 |
|
|
|
[考研] 考研調(diào)劑 +10 | 呼呼？~+123456 2026-03-24 | 10/500 |
|
|
|
[考研] 調(diào)劑推薦 +5 | 清酒714 2026-03-26 | 6/300 |
|
|
|
[考研] 考研調(diào)劑 +9 | 小蠟新筆 2026-03-26 | 9/450 |
|
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +3 | 崔wj 2026-03-26 | 3/150 |
|
|
|
[考研] 284求調(diào)劑 +11 | junqihahaha 2026-03-26 | 12/600 |
|
|
|
[考研] 321求調(diào)劑 +6 | Ymlll 2026-03-24 | 6/300 |
|
|
|
[考研] 086000生物與醫(yī)藥292求調(diào)劑 +6 | 小小陳小小 2026-03-22 | 9/450 |
|
|
|
[考研] 材料277求調(diào)劑 +5 | min3 2026-03-24 | 5/250 |
|
|
|
[考研] 材料專碩 335 分求調(diào)劑 +4 | 拒絕冷暴力 2026-03-25 | 4/200 |
|
|
|
[考研] 293求調(diào)劑 +7 | 加一一九 2026-03-24 | 7/350 |
|
|
|
[考研] B區(qū)考研調(diào)劑 +4 | yqdszhdap－ 2026-03-22 | 5/250 |
|
|
|
[考研] 求調(diào)劑一志愿武漢理工大學(xué)材料工程（085601） +5 | WW.' 2026-03-23 | 7/350 |
|
|
|
[考研] 333求調(diào)劑 +3 | ALULU4408 2026-03-23 | 3/150 |
|
|
|
[考研] 306求調(diào)劑 +5 | 來好運(yùn)來來來 2026-03-22 | 5/250 |
|
|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境總分308求調(diào)劑 +7 | 墨墨漠 2026-03-20 | 7/350 |
|