版塊導(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í)熱門(mén)版塊排行榜

北京石油化工學(xué)院2026年研究生招生接收調(diào)劑公告

返回列表

qs001

新蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 6.5
帖子: 2
在線: 1.9小時(shí)
蟲(chóng)號(hào): 2219462
注冊(cè): 2013-01-03
專業(yè): 防災(zāi)工程

[求助] fsolve計(jì)算結(jié)果問(wèn)題

計(jì)算結(jié)果感覺(jué)誤差好大，是由于初值取值造成的影響嗎，還是要設(shè)置option來(lái)控制精度，第一次使用matlab一竅不通，請(qǐng)大家?guī)兔聪隆?br />                                     First-Order                   Norm of
Iteration  Func-count Residual    optimality    Lambda          step
   0          8       39.2026    1.13e+003       0.01
   1       16    0.00103045          0.705       0.001    0.294456
   2       24 3.55079e-005       0.0727    0.0001    0.237055
   3       32 7.28141e-006          0.19    1e-005    0.144512
   4       40 2.91808e-006       0.0123    1e-006    0.0248491
   5       48 2.90763e-006    3.58e-007    1e-007 0.000338087
   6       56 2.90763e-006    1.02e-009    1e-008 1.90914e-006

Equation solved, fsolve stalled.

fsolve stopped because the relative size of the current step is less than the
default value of the step size tolerance and the vector of function values
is near zero as measured by the selected value of the function tolerance.

y =

-2.2413 2.2083 -0.1310 0.1318 1.0596 0.0211 -0.0272

fval =

0.0000
0.0008
-0.0008
0.0003
0.0000
0.0000
-0.0011
0.0006

下面是我要解的方程組
function F=myfun(x)
b=1;a=1;
F=[-0.0575*x(3)/a-x(7)*a*0.5-x(6);
x(3)/a^2+x(4)+17.39*x(7)*log(a)+0.8695652174*b*log(a);
-x(1)/(x(5)^2)-x(2)+x(3)/(x(5)^2)+x(4)-92.49880554*x(7)*log(x(5))+0.8695652174*b*log(x(5));
0.0065*x(1)/x(5)-0.0026*x(2)*x(5)-0.0575*x(3)/x(5)+x(7)*x(5)*((0.7-0.4*log(x(5)))/1.4-0.9347826087);
-2*x(1)/(x(5)^2)+109.8901099*x(7)-b;
2*x(1)/(1.5^2)-109.8901099*x(7)-b;
0.444444444444444*x(1)/a^2+x(2)+109.8901099*x(7)*log(1.5*a);
-x(1)*log(1.5*a/x(5))+0.5*(1.5^2*a^2-x(5)^2)*x(2)-x(3)*log(x(5)/a)+0.5*(x(5)^2-a^2)*x(4)+17.39130435*x(7)*(0.5*x(5)^2*(log(x(5))+0.5)-0.5*a^2*(log(a)+0.5))+109.8901099*x(7)*(0.5*1.5*1.5*a^2*(0.5+log(1.5*a))-0.5*x(5)^2*(log(x(5))+0.5))+0.0209742*1.5^2*a^2*b+0.8695652174*b*(0.5*x(5)^2*(log(x(5))+0.5)-0.5*a^2*(log(a)+0.5))];

選取的初值是x0=[-2,2,-0.5,0.5,1.05,0.01,-0.01]

回復(fù)此樓

» 猜你喜歡

本科新能源科學(xué)與工程，一志愿華理能動(dòng)285求調(diào)劑已經(jīng)有5人回復(fù)
一志愿華北電力大學(xué)能動(dòng)專碩，293，求調(diào)劑已經(jīng)有5人回復(fù)
調(diào)劑310 已經(jīng)有7人回復(fù)
食品工程專碩一志愿中海洋309求調(diào)劑已經(jīng)有6人回復(fù)
275求調(diào)劑已經(jīng)有14人回復(fù)
283求調(diào)劑已經(jīng)有3人回復(fù)
085602 化工專碩 338分求調(diào)劑已經(jīng)有12人回復(fù)
312，生物學(xué)求調(diào)劑已經(jīng)有3人回復(fù)
化學(xué)0703 調(diào)劑 306分一志愿211 已經(jīng)有4人回復(fù)
311（085601）求調(diào)劑已經(jīng)有3人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

求助 MATLAB-fsolve 已經(jīng)有3人回復(fù)
使用fsolve求解非線性方程問(wèn)題已經(jīng)有9人回復(fù)
MATLAB求解代數(shù)方程組問(wèn)題，請(qǐng)幫忙已經(jīng)有9人回復(fù)
各位朋友，如何用matlab編寫(xiě)這個(gè)方程，并求解？已經(jīng)有15人回復(fù)
求助 MATLAB解方程組-fslove 已經(jīng)有7人回復(fù)
vasp計(jì)算態(tài)密度結(jié)果出問(wèn)題已經(jīng)有11人回復(fù)
NBO計(jì)算結(jié)果問(wèn)題已經(jīng)有8人回復(fù)
求助！matlab用fsolve函數(shù)求解非線性方程組的問(wèn)題！已經(jīng)有19人回復(fù)
Ansys載荷步計(jì)算結(jié)果的保存問(wèn)題已經(jīng)有8人回復(fù)
各位大俠，請(qǐng)問(wèn)下我在linux計(jì)算，結(jié)果老是出現(xiàn)下面的問(wèn)題，請(qǐng)問(wèn)怎么辦呀。已經(jīng)有23人回復(fù)
matlab 問(wèn)題求助，是否計(jì)算結(jié)果錯(cuò)誤，如何改正？已經(jīng)有3人回復(fù)
非線性方程組matlab擬合已經(jīng)有9人回復(fù)
mathematica請(qǐng)教計(jì)算結(jié)果問(wèn)題已經(jīng)有4人回復(fù)
為什么計(jì)算區(qū)域的大小對(duì)計(jì)算結(jié)果的精度影響這么大？顆粒擾流問(wèn)題已經(jīng)有6人回復(fù)
氣相色譜問(wèn)題，出峰但是不計(jì)算結(jié)果已經(jīng)有5人回復(fù)
Abinit 新舊版計(jì)算結(jié)果差別多大可以認(rèn)為正常，沒(méi)有安裝錯(cuò)誤之類的問(wèn)題。已經(jīng)有4人回復(fù)
【求助】請(qǐng)高手幫我解答個(gè)比較疑惑的問(wèn)題（UNIFAC估算）已經(jīng)有4人回復(fù)
【求助】vasp5 光學(xué)計(jì)算結(jié)果問(wèn)題已經(jīng)有4人回復(fù)
【求助】請(qǐng)教一個(gè)關(guān)于解方程的問(wèn)題已經(jīng)有8人回復(fù)
【求助】fluent的傳熱問(wèn)題，計(jì)算結(jié)果有沒(méi)有問(wèn)題，請(qǐng)大家?guī)兔纯�。【已解決】已經(jīng)有10人回復(fù)

1樓 2013-01-07 14:10:08

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

nono2009

超級(jí)版主 (文學(xué)泰斗)

No gains, no pains.

專家經(jīng)驗(yàn): +21105
仿真EPI: 1
應(yīng)助: 28684 (院士)
貴賓: 513.911
金幣: 2555230
散金: 27828
紅花: 2148
沙發(fā): 66666
帖子: 1602255
在線: 65200.9小時(shí)
蟲(chóng)號(hào): 827383
注冊(cè): 2009-08-13
性別: GG
專業(yè): 工程熱物理與能源利用
管轄: 科研家籌備委員會(huì)

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

★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
qs001: 金幣+5, ★有幫助 2013-01-08 20:44:22
xiegangmai: 金幣+1, 謝謝建議 2013-01-08 21:06:39

建議樓主采用 format long（雙精度）再試試

贊一下

回復(fù)此樓

2樓2013-01-07 14:23:41

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

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