| 2 | 1/1 | 返回列表 |
| 查看: 2552 | 回復(fù): 1 | |||
qs001新蟲(chóng) (初入文壇)
|
[求助]
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] |
超級(jí)版主 (文學(xué)泰斗)
No gains, no pains.
![]() |
專家經(jīng)驗(yàn): +21105 |
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料求調(diào)劑一志愿哈工大324 +5 | 閆旭東 2026-03-28 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 0703化學(xué) +6 | 妮妮ninicgb 2026-03-27 | 6/300 |
|
|
[考研] 材料與化工(0856)304求B區(qū)調(diào)劑 +7 | 邱gl 2026-03-27 | 7/350 |
|
|
[考研] 085600材料與化工調(diào)劑 +10 | A-哆啦Z夢(mèng) 2026-03-23 | 16/800 |
|
|
[考研] 一志愿太原理工安全工程300分,求調(diào)劑 +4 | 0857求調(diào)劑. 2026-03-24 | 4/200 |
|
|
[考研] 化學(xué)308分求調(diào)劑 +8 | 你好明天你好 2026-03-23 | 9/450 |
|
|
[考研] 304求調(diào)劑 +3 | 曼殊2266 2026-03-27 | 3/150 |
|
|
[考研]
|
nnnnnnn5 2026-03-25 | 5/250 |
|
|
[考研] 351求調(diào)劑 +4 | 麥克阿磊 2026-03-24 | 4/200 |
|
|
[考研] 342求調(diào)劑 +3 | 加油a李zs 2026-03-26 | 3/150 |
|
|
[考研] 求調(diào)劑 一志愿 本科 北科大 化學(xué) 343 +6 | 13831862839 2026-03-24 | 7/350 |
|
|
[考研] 一志愿華理,數(shù)一英一285求A區(qū)調(diào)劑 +8 | AZMK 2026-03-25 | 10/500 |
|
|
[考研] 機(jī)械學(xué)碩310分,數(shù)一英一,一志愿211本科雙非找調(diào)劑信息 +3 | @357 2026-03-25 | 3/150 |
|
|
[考研] 297求調(diào)劑 +6 | 田洪有 2026-03-26 | 6/300 |
|
|
[考研] 材料專碩 335 分求調(diào)劑 +4 | 拒絕冷暴力 2026-03-25 | 4/200 |
|
|
[考研] 340求調(diào)劑 +5 | 話梅糖111 2026-03-24 | 5/250 |
|
|
[考研] 一志愿南航材料專317分求調(diào)劑 +5 | 炸呀炸呀炸薯?xiàng)l 2026-03-23 | 5/250 |
|
|
[考研] 336化工調(diào)劑 +4 | 王大坦1 2026-03-23 | 5/250 |
|
|
[考研]
|
酥酥魚(yú).. 2026-03-21 | 4/200 |
|
|
[考研] 308求調(diào)劑 +3 | 墨墨漠 2026-03-21 | 3/150 |
|