 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 1117 | 回復(fù): 4 | ||
XiaoGuads新蟲 (初入文壇)
|
[求助]
數(shù)值解非線性方程組時(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)該是在初值附近的，可為什么解不出來呢 |
回復(fù)此樓» 猜你喜歡
A區(qū)一本交叉課題組，低分調(diào)劑，招收機(jī)械電子信息通信等交叉方向
已經(jīng)有27人回復(fù)
-
一志愿211，0860總分286食品工程專業(yè)求調(diào)劑
已經(jīng)有4人回復(fù)
-
282求調(diào)劑
已經(jīng)有8人回復(fù)
-
085600，一志愿鄭州大學(xué)，280分求調(diào)劑
已經(jīng)有5人回復(fù)
-
材料277分求調(diào)劑
已經(jīng)有11人回復(fù)
-
新疆大學(xué)地質(zhì)與礦業(yè)工程學(xué)院招生
已經(jīng)有5人回復(fù)
-
材料專碩290求調(diào)劑
已經(jīng)有4人回復(fù)
-
不限學(xué)校專業(yè)的調(diào)劑同學(xué)看過來
已經(jīng)有5人回復(fù)
-
招調(diào)劑
已經(jīng)有8人回復(fù)
-
一志愿蘇大材料工程專碩293求調(diào)劑
已經(jīng)有3人回復(fù)
2樓2019-01-30 17:49:49
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
月只藍(lán): 金幣+10, 感謝幫助！ 2019-01-31 18:04:04
感謝參與，應(yīng)助指數(shù) +1
月只藍(lán): 金幣+10, 感謝幫助！ 2019-01-31 18:04:04
|
利用最小二乘法（lsqnonlin），可以求得 x1=0.0109083348225495 ;x2=4.99999947195261;x3=0.0384214997768871;4=25.0000321470904 驗(yàn)證： 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 |
» 本帖已獲得的紅花（最新10朵）
XiaoGuads
3樓2019-01-31 14:39:33
XiaoGuads
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 51.2
- 帖子: 17
- 在線: 3小時(shí)
- 蟲號: 8840508
- 注冊: 2018-05-23
- 專業(yè): 物理學(xué)I
4樓2019-02-04 16:56:15
dingd
鐵桿木蟲 (職業(yè)作家)
- 應(yīng)助: 1641 (講師)
- 金幣: 15037.3
- 散金: 101
- 紅花: 234
- 帖子: 3410
- 在線: 1223.7小時(shí)
- 蟲號: 291104
- 注冊: 2006-10-28
5樓2019-02-28 18:14:16
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 XiaoGuads 的主題更新
| 5 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿211，0860總分286食品工程專業(yè)求調(diào)劑 +4 | 林林Winnie 2026-03-05 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 招調(diào)劑 +6 | 帆船哥 2026-03-04 | 8/400 |
|
|
|
[考研] 289求調(diào)劑 +3 | 馳光絕景 2026-03-04 | 3/150 |
|
|
|
[考研] 材料調(diào)劑 +4 | ounce. 2026-03-04 | 5/250 |
|
|
|
[考研] 材料調(diào)劑 +5 | ws 上岸鴨 2026-03-05 | 5/250 |
|
|
|
[考研] 【求調(diào)劑】293分環(huán)境工程求調(diào)劑材料/化工，服從調(diào)劑，抗壓能力強(qiáng)！ +8 | xiiiia 2026-03-04 | 8/400 |
|
|
|
[考研] 化學(xué)專業(yè)調(diào)劑 +3 | 好好好1233 2026-03-04 | 3/150 |
|
|
|
[考研] 一志愿天津大學(xué)085600 319分 材料與化工 金屬方向 求調(diào)劑 +6 | 青科11 2026-03-02 | 6/300 |
|
|
|
[考研] 材料328求調(diào)劑 +10 | 一個(gè)蘿卜02 2026-03-03 | 10/500 |
|
|
|
[考研] 求調(diào)劑院校 +6 | 云朵452 2026-03-02 | 12/600 |
|
|
|
[考研] 307求調(diào)劑 +4 | xuyuan111 2026-03-04 | 5/250 |
|
|
|
[考研]
085600 英一數(shù)二272求調(diào)劑
 5+6
|
vida_a 2026-03-01 | 47/2350 |
|
|
|
[基金申請] 請問大家，研究風(fēng)險(xiǎn)與應(yīng)對措施那里， 大家都怎么寫呢 ？ +3 | cauasen 2026-03-02 | 3/150 |
|
|
|
[考研] 276求調(diào)劑 +8 | 路lyh123 2026-02-28 | 10/500 |
|
|
|
[考研] 298求調(diào)劑 +10 | 人間唯你是清歡 2026-02-28 | 14/700 |
|
|
|
[考研] 0856材料調(diào)劑 +5 | 沿岸有貝殼OUC 2026-03-02 | 5/250 |
|
|
|
[考研] 0856化工專碩求調(diào)劑 +15 | 董boxing 2026-03-01 | 15/750 |
|
|
|
[考研] 264求調(diào)劑 +4 | 巴拉巴拉根556 2026-02-28 | 4/200 |
|
|
|
[考研] 295復(fù)試調(diào)劑 +3 | 簡木ChuFront 2026-03-01 | 3/150 |
|
|
|
[論文投稿]
求助coordination chemistry reviews 的寫作模板
10+3
|
ljplijiapeng 2026-02-27 | 4/200 |
|