| 4 | 1/1 | 返回列表 |
| 查看: 460 | 回復(fù): 3 | ||
liaoliao8545新蟲 (小有名氣)
|
[求助]
matlab
|
|
我用ezplot畫函數(shù)圖找零點運(yùn)行后出現(xiàn)下面問題 Current plot held Warning: Function failed to evaluate on array inputs; vectorizing the function may speed up its evaluation and avoid the need to loop over array elements. 什么原因呢 還有我用fsolve解非線性方程組的解 x30 x40 x50 x60 為初值 運(yùn)行 出現(xiàn)下面的情況 x30 = 0.0158 x40 = 0.6910 x50 = 1.3824e-035 x60 = 2.8181e-032 Optimization terminated: no further progress can be made. Trust-region radius less than 2*eps. Problem may be ill-conditioned or Jacobian may be inaccurate. Try using exact Jacobian or check Jacobian for errors. z = 0.0158 0.6910 -0.0000 0.0000 什么原因呢 這是非線性方程組m文件 function f=fun_fun(x) global C1 C2 C7 C8 C9 C10 d1 d2 d3 n m;%%%%這些是由主程序求得的 f=[2*C2*x(4)*sqrt(x(1))+2*x(4)+C7*sqrt(x(4))*sqrt(x(4))+C9*sqrt(x(4))-(C1*x(3)*sqrt(x(1))+x(3))*d1; 2*C1*x(3)*sqrt(x(1))+C2*x(4)*sqrt(x(3))+2*x(1)+x(3)+C7*sqrt(x(4))*sqrt(x(1))+C8*sqrt(x(1))+C10*sqrt(x(1))*sqrt(x(2))-(C1*x(3)*sqrt(x(1))+x(3))*d2; 2*x(2)+C10*sqrt(x(1))*sqrt(x(2))-(C1*x(3)*sqrt(x(1))+x(3))*d3; C1*x(3)*sqrt(x(1))+C2*x(4)*sqrt(x(1))+x(3)+x(2)+x(3)+x(4)+C7*sqrt(x(4))*sqrt(x(3))+C8*sqrt(x(3))+C9*sqrt(x(4))+C10*sqrt(x(1))*sqrt(x(2))+1/(n*(C1*x(3)*sqrt(x(1))+x(3)))-1;]; |
新蟲 (小有名氣)
|
補(bǔ)充一下非線性方程組如下 :比較容易看清 2C_2 x_4 √(x_1 )+2x_4+C_7 √(x_4 ) √(x_1 )+C_9 √(x_4 )-(C_1 x_3 √(x_1 )+x_3 ) d_1=0 2C_1 x_3 √(x_1 )+C_2 x_4 √(x_1 )+2x_1+x_3+C_7 √(x_4 ) √(x_1 )+C_8 √(x_1 )+C_10 √(x_1 ) √(x_2 )-(C_1 x_3 √(x_1 )+x_3 ) d_2=0 2x_2+C_10 √(x_1 ) √(x_2 )-(C_1 x_3 √(x_1 )+x_3 ) d_3=0 C_1 x_3 √(x_1 )+C_2 x_4 √(x_1 )+x_1+x_2+x_3+x_4+C_7 √(x_4 ) √(x_1 )+C_8 √(x_1 )+C_9 √(x_4 )+C_10 √(x_1 ) √(x_2 )+1/n (C_1 x_3 √(x_1 )+x_3 )=1 |
新蟲 (小有名氣)
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑 +7 | Ymlll 2026-03-18 | 10/500 |
|
|---|---|---|---|---|
|
[考研] 294求調(diào)劑材料與化工專碩 +13 | 陌の森林 2026-03-18 | 13/650 |
|
|
[考研] 本人考085602 化學(xué)工程 專碩 +17 | 不知道叫什么! 2026-03-15 | 19/950 |
|
|
[考研] 一志愿天津大學(xué)化學(xué)工藝專業(yè)(081702)315分求調(diào)劑 +11 | yangfz 2026-03-17 | 11/550 |
|
|
[考研] 304求調(diào)劑 +3 | 曼殊2266 2026-03-18 | 3/150 |
|
|
[考研] 287求調(diào)劑 +3 | 晨昏線與星海 2026-03-19 | 4/200 |
|
|
[考研] 材料080500調(diào)劑求收留 +4 | 一顆meteor 2026-03-13 | 4/200 |
|
|
[考研] 0703化學(xué) 305求調(diào)劑 +4 | FY_yy 2026-03-14 | 4/200 |
|
|
[考研] 工科材料085601 279求調(diào)劑 +6 | 困于星晨 2026-03-17 | 6/300 |
|
|
[考研] 301求調(diào)劑 +9 | yy要上岸呀 2026-03-17 | 9/450 |
|
|
[考研] 268求調(diào)劑 +8 | 一定有學(xué)上- 2026-03-14 | 9/450 |
|
|
[考研] 308求調(diào)劑 +4 | 是Lupa啊 2026-03-16 | 4/200 |
|
|
[考研]
|
zhouzhen654 2026-03-16 | 3/150 |
|
|
[考研] [導(dǎo)師推薦]西南科技大學(xué)國防/材料導(dǎo)師推薦 +3 | 尖角小荷 2026-03-16 | 6/300 |
|
|
[考研] 機(jī)械專碩325,尋找調(diào)劑院校 +3 | y9999 2026-03-15 | 5/250 |
|
|
[考研] 0854控制工程 359求調(diào)劑 可跨專業(yè) +3 | 626776879 2026-03-14 | 9/450 |
|
|
[考研] 326求調(diào)劑 +3 | mlpqaz03 2026-03-15 | 3/150 |
|
|
[考研] 中科大材料專碩319求調(diào)劑 +3 | 孟鑫材料 2026-03-13 | 3/150 |
|
|
[基金申請] 現(xiàn)在如何回避去年的某一個專家,不知道名字 +3 | zk200107 2026-03-12 | 6/300 |
|
|
[考研] 070303一志愿西北大學(xué)學(xué)碩310找調(diào)劑 +3 | d如愿上岸 2026-03-13 | 3/150 |
|