| 5 | 1/1 | 返回列表 |
| 查看: 1138 | 回復(fù): 4 | |||
XiaoGuads新蟲 (初入文壇)
|
[求助]
數(shù)值解非線性方程組時找不到解,希望各位能幫幫我 已有2人參與
|
|
function y=fun(x) alpha=0.125*pi/180;x1=0;z1=5;x2=0;z2=25;L1=5;L2=5;f=30;h0=20;n2=1.49386; y(1) = tan(alpha) - x(1)/x(2); y(2) = n2*sin(atan((x1 - x(1))/(z1 - x(2))) - atan((x(2) - x(4))/(x(1) - x(3)))) - sin(pi/2 - alpha + atan((x1 - x(1))/(z1 - x(2)))) ; y(3) = n2*sin(atan((x2 - x(3))/(z2 - x(4))) - atan((x(2) - x(4))/(x(1) - x(3)))) - sin(pi/2 + atan((x2 - x(3))/(z2 - x(4)))); y(4) = L1 + L2 - f + h0 - h0*n2 + n2*((x(1) - x(3))^2 + (x(2) - x(4))^2)^(1/2); end ---------------------------------------------------------------------- 命令行輸入 >> fsolve('fun',[0.1 5.1 0.1 25.1]) ---------------------------------------------------------------------- No solution found. fsolve stopped because the relative size of the current step is less than the default value of the step size tolerance squared, but the vector of function values is not near zero as measured by the default value of the function tolerance. 合理的解應(yīng)該是在初值附近的,可為什么解不出來呢 |
至尊木蟲 (著名寫手)
|
利用最小二乘法(lsqnonlin),可以求得 x1=0.0109083348225495 ;x2=4.99999947195261;x3=0.0384214997768871;4=25.0000321470904 驗證: y1 = tan(alpha) - x(1)/x(2)=-2.16860202711944e-09 y2 = n2*sin(atan((x1 - x(1))/(z1 - x(2))) - atan((x(2) - x(4))/(x(1) - x(3)))) - sin(pi/2 - alpha + atan((x1 - x(1))/(z1 - x(2)))) =5.90246506399546e-06 y3 = n2*sin(atan((x2 - x(3))/(z2 - x(4))) - atan((x(2) - x(4))/(x(1) - x(3)))) - sin(pi/2 + atan((x2 - x(3))/(z2 - x(4))))=-3.15643589759254e-05 |

新蟲 (初入文壇)
鐵桿木蟲 (職業(yè)作家)
| 5 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085600 材料與化工 295 求調(diào)劑 +9 | dream…… 2026-03-10 | 10/500 |
|
|---|---|---|---|---|
|
[考研] 304求調(diào)劑(085602一志愿985) +9 | 化工人999 2026-03-09 | 9/450 |
|
|
[考研] 材料工程085601調(diào)劑求老師收留 +9 | 強木木木 2026-03-07 | 11/550 |
|
|
[考研] 0817化學(xué)工程319求調(diào)劑 +8 | lv945 2026-03-08 | 10/500 |
|
|
[考研] 一志愿北京化工大學(xué)材料與化工296分求調(diào)劑 +15 | 稻妻小編 2026-03-09 | 17/850 |
|
|
[考研]
一志愿天津大學(xué)材料與化工275求調(diào)劑
10+5
|
穿只靴子 2026-03-07 | 22/1100 |
|
|
[考研] 材料工程,326分,求調(diào)劑 +5 | KRSLSR 2026-03-10 | 5/250 |
|
|
[考研] 298求調(diào)劑 +3 | Vv呀! 2026-03-10 | 3/150 |
|
|
[考研] 收調(diào)劑 +7 | 調(diào)劑的考研學(xué)生 2026-03-10 | 7/350 |
|
|
[考博] 求材料讀博院校 +9 | yanglei131 2026-03-08 | 9/450 |
|
|
[考研] 材料與化工 一志愿山大 321分 求調(diào)劑 +6 | 每天散步 2026-03-09 | 7/350 |
|
|
[考研] 化工284求調(diào)劑 +12 | 小黑想變白 2026-03-07 | 12/600 |
|
|
[考研] 一志愿山東大學(xué),總分327,英語二79,有論文,有競賽,已過四六級 +3 | 木木目目1 2026-03-09 | 3/150 |
|
|
[考研] 320求調(diào)劑 +4 | 魏zy 2026-03-08 | 4/200 |
|
|
[考研] 0856材料與化工290求調(diào)劑 +7 | Nebulala 2026-03-08 | 8/400 |
|
|
[考研] 理學(xué)287分求調(diào)劑 +4 | 看看我. 2026-03-05 | 4/200 |
|
|
[考研] 083000環(huán)境科學(xué)與工程調(diào)劑 +5 | 加油呀fxy 2026-03-07 | 6/300 |
|
|
[考研] 301求調(diào)劑 +5 | 一二LV 2026-03-07 | 5/250 |
|
|
[考研] 070300化學(xué)求調(diào)劑292分 +3 | 打烊eee 2026-03-07 | 3/150 |
|
|
[考研] 275求調(diào)劑 +4 | 大爆炸難民 2026-03-06 | 5/250 |
|