 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 534 | 回復(fù): 2 | ||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||
zidingxiang金蟲 (小有名氣)
|
[求助]
運(yùn)行結(jié)果與設(shè)置精度不符，結(jié)果應(yīng)與初值接近才對(duì) 已有1人參與
|
|
|
1.首先建立函數(shù)myfun,儲(chǔ)存方程組編程如下將myfun.m保存到工作路徑中: function F = myfun( x ) %this function designed for solving nonlinear equations a = x(1); b = x(2); c = x(3); d = x(4); e = x(5); f = x(6); F = [a+b*10+c*10^2+(d+e*10)/(0.6438^2-f^2)-2*1.620015/(1.620015^2-1)*(-0.87+2.6); a+b*10+c*10^2+(d+e*10)/(0.6328^2-f^2)-2*1.620507/(1.620507^2-1)*(-0.87+2.7); a+b*10+c*10^2+(d+e*10)/(0.5893^2-f^2)-2*1.622702/(1.622702^2-1)*(-0.87+2.8); a+b*10+c*10^2+(d+e*10)/(0.5461^2-f^2)-2*1.625401/(1.625401^2-1)*(-0.87+2.9); a+b*10+c*10^2+(d+e*10)/(0.4800^2-f^2)-2*1.631029/(1.631029^2-1)*(-0.88+3.3); a+b*10+c*10^2+(d+e*10)/(0.4358^2-f^2)-2*1.636369/(1.636369^2-1)*(-0.88+3.6)]; end 2.運(yùn)行程序 x0=[-3e-005,-1.5e-008,-5.0e-011,3.0e-007,2.0e-010,0]; options=optimset('Display','iter','TolFun',1e-020,'TolX',1e-020,'MaxIter',1e100); format long; [x,fval]=fsolve(@myfun,x0,options) 運(yùn)行結(jié)果顯示如下： Norm of First-order Trust-region Iteration Func-count f(x) step optimality radius 0 7 106.837 2.5e+003 1 1 14 1.58981 0.0393664 28.2 1 2 21 0.0425606 0.0513838 36.8 1 3 28 0.0198142 0.000578732 0.414 1 4 35 0.0194798 0.000755402 0.541 1 5 42 0.0194749 8.50804e-006 0.00609 1 6 49 0.0194748 1.11053e-005 0.00795 1 7 56 0.0194748 1.25077e-007 8.95e-005 1 8 63 0.0194748 1.6325e-007 0.000117 1 9 70 0.0194748 1.84156e-009 1.32e-006 1 10 77 0.0194748 2.41651e-009 1.58e-006 1 11 84 0.0194748 2.49273e-011 1.22e-008 1 Optimizer appears to be converging to a point which is not a root. Relative function value changing by less than max(options.TolFun^2,eps) but sum-of-squares of function values is greater than or equal to sqrt(options.TolFun) Try again with a new starting guess. x = 0.000167225746880 0.001972242468804 0.019722574638040 0.006232656757062 0.062323567700574 0.000000000137224 fval = 0.060290300323624 -0.083434220551602 -0.030535103096213 0.083729176350573 -0.030583681232919 0.000533528265654 注：函數(shù)精度遠(yuǎn)遠(yuǎn)小于設(shè)定精度值；參數(shù)值應(yīng)接近于設(shè)定值，但相差很遠(yuǎn)。 |
回復(fù)此樓» 猜你喜歡
08工學(xué)調(diào)劑
已經(jīng)有5人回復(fù)
-
281求調(diào)劑（0805）
已經(jīng)有25人回復(fù)
-
085600材料與化工
已經(jīng)有6人回復(fù)
-
265求調(diào)劑
已經(jīng)有8人回復(fù)
-
0817 化學(xué)工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章
已經(jīng)有21人回復(fù)
-
295材料求調(diào)劑，一志愿武漢理工085601專碩
已經(jīng)有4人回復(fù)
-
0856調(diào)劑，是學(xué)校就去
已經(jīng)有9人回復(fù)
-
298-一志愿中國(guó)農(nóng)業(yè)大學(xué)-求調(diào)劑
已經(jīng)有9人回復(fù)
-
求調(diào)劑
已經(jīng)有4人回復(fù)
-
學(xué)校已經(jīng)提交到NSFC，還能修改嗎？
已經(jīng)有7人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- matlab解非線性方程組程序糾錯(cuò)與正確運(yùn)行
已經(jīng)有6人回復(fù)
dingd
鐵桿木蟲 (職業(yè)作家)
- 應(yīng)助: 1641 (講師)
- 金幣: 15037.3
- 散金: 101
- 紅花: 234
- 帖子: 3410
- 在線: 1223.7小時(shí)
- 蟲號(hào): 291104
- 注冊(cè): 2006-10-28
3樓2014-07-08 09:14:23
zidingxiang
金蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 774.8
- 散金: 69
- 紅花: 1
- 帖子: 114
- 在線: 31.7小時(shí)
- 蟲號(hào): 958621
- 注冊(cè): 2010-03-01
- 性別: MM
- 專業(yè): 玻璃材料
2樓2014-07-08 08:23:19
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 295材料求調(diào)劑，一志愿武漢理工085601專碩 +4 | Charlieyq 2026-03-19 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 招收調(diào)劑碩士 +4 | lidianxing 2026-03-19 | 12/600 |
|
|
|
[考研] 一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑 +11 | Ymlll 2026-03-18 | 14/700 |
|
|
|
[考研] 材料專碩英一數(shù)二306 +6 | z1z2z3879 2026-03-18 | 6/300 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 +10 | 妮妮ninicgb 2026-03-15 | 14/700 |
|
|
|
[考研] 288求調(diào)劑，一志愿華南理工大學(xué)071005 +5 | ioodiiij 2026-03-17 | 5/250 |
|
|
|
[考研] 一志愿北京化工大學(xué)0703化學(xué)318分，有科研經(jīng)歷，求調(diào)劑 +3 | 一瓶苯甲酸 2026-03-14 | 3/150 |
|
|
|
[考研] 一志愿西安交通大學(xué)材料工程專業(yè) 282分求調(diào)劑 +5 | 楓橋ZL 2026-03-18 | 7/350 |
|
|
|
[考研] 286求調(diào)劑 +6 | lemonzzn 2026-03-16 | 10/500 |
|
|
|
[考研] 287求調(diào)劑 +3 | 晨昏線與星海 2026-03-19 | 4/200 |
|
|
|
[考研] 0703化學(xué)調(diào)劑，求各位老師收留 +10 | 秋有木北 2026-03-14 | 10/500 |
|
|
|
[考研] 一志愿華中科技大學(xué)，080502，354分求調(diào)劑 +4 | 守候夕陽CF 2026-03-18 | 4/200 |
|
|
|
[考研] 297求調(diào)劑 +8 | 戲精丹丹丹 2026-03-17 | 8/400 |
|
|
|
[考研]
 0854，計(jì)算機(jī)類招收調(diào)劑
+3
|
胡辣湯放糖 2026-03-15 | 6/300 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 +3 | 妮妮ninicgb 2026-03-17 | 3/150 |
|
|
|
[考研] 0703化學(xué)336分求調(diào)劑 +6 | zbzihdhd 2026-03-15 | 7/350 |
|
|
|
[考研] 268求調(diào)劑 +6 | 簡(jiǎn)單點(diǎn)0 2026-03-17 | 6/300 |
|
|
|
[考研] 304求調(diào)劑 +5 | 素年祭語 2026-03-15 | 5/250 |
|
|
|
[考研] 0856專碩279求調(diào)劑 +5 | 加油加油！? 2026-03-15 | 5/250 |
|
|
|
[考研] 本科南京大學(xué)一志愿川大藥學(xué)327 +3 | 麥田耕者 2026-03-14 | 3/150 |
|