| 5 | 1/1 | 返回列表 |
| 查看: 1136 | 回復(fù): 4 | ||
| 當(dāng)前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | ||
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)該是在初值附近的,可為什么解不出來呢 |
鐵桿木蟲 (職業(yè)作家)
至尊木蟲 (著名寫手)
|
利用最小二乘法(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 |

新蟲 (初入文壇)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 311求調(diào)劑 +5 | zchqwer 2026-03-10 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 081200計算機科學(xué)與技術(shù)數(shù)一英一306求調(diào)劑 +3 | intankt 2026-03-08 | 3/150 |
|
|
[考研] 一志愿江南大學(xué)085701環(huán)境工程專碩總分287求調(diào)劑 +4 | 18266118446 2026-03-09 | 4/200 |
|
|
[考研] 0857 資源與環(huán)境 285分 +5 | 未名考生 2026-03-09 | 5/250 |
|
|
[考研] 材料工程,326分,求調(diào)劑 +3 | KRSLSR 2026-03-10 | 3/150 |
|
|
[考研] 313分生物學(xué)求調(diào)劑 +5 | Yyt楊1 2026-03-09 | 7/350 |
|
|
[考研] 求調(diào)劑 +4 | yfihxh 2026-03-09 | 4/200 |
|
|
[考研] 材料工程專碩調(diào)劑自薦信,初試323 +6 | 李白26 2026-03-07 | 6/300 |
|
|
[考研] A區(qū)一本交叉課題組,低分調(diào)劑,招收機械電子信息通信等交叉方向 +56 | lisimayy 2026-03-04 | 73/3650 |
|
|
[考研] 0817化學(xué)工程319求調(diào)劑 +7 | lv945 2026-03-08 | 9/450 |
|
|
[基金申請] 面上項目還需要AI說明嗎? +3 | liyundong 2026-03-08 | 3/150 |
|
|
[考研] 070300化學(xué)求調(diào)劑 +5 | 撲風(fēng)鈴的貓 2026-03-08 | 10/500 |
|
|
[考研] 0703化學(xué)求調(diào)劑,總分309分,一志愿華南師范 +3 | 花與葉@ 2026-03-08 | 3/150 |
|
|
[教師之家] 交大前校長王樹國:現(xiàn)在最先進的科技并不在大學(xué)實驗室,而是在企業(yè)研究院 +4 | zju2000 2026-03-08 | 6/300 |
|
|
[考研] 083000,總分284,求調(diào)劑 +5 | 徐yr 2026-03-04 | 5/250 |
|
|
[考研] 一志愿211 085600 280數(shù)二英二求調(diào)劑 +3 | 月山斜 2026-03-06 | 3/150 |
|
|
[考研] 267化工調(diào)劑求助 +7 | 聰少OZ 2026-03-04 | 7/350 |
|
|
[考研] 求調(diào)劑 +4 | 呼呼?~+123456 2026-03-05 | 5/250 |
|
|
[考研] 085602 293分求調(diào)劑 +3 | SivanNano. 2026-03-05 | 3/150 |
|
|
[考研] 274環(huán)境工程求調(diào)劑 +6 | 扶柳盈江 2026-03-05 | 6/300 |
|