版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門(mén)版塊排行榜

返回列表

當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖

XiaoGuads

新蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 51.2
帖子: 17
在線: 3小時(shí)
蟲(chóng)號(hào): 8840508
注冊(cè): 2018-05-23
專(zhuān)業(yè): 物理學(xué)I

[求助] 數(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)該是在初值附近的，可為什么解不出來(lái)呢

回復(fù)此樓

» 猜你喜歡

一志愿安徽大學(xué)材料工程專(zhuān)碩313分，求調(diào)劑的學(xué)校已經(jīng)有3人回復(fù)
考研一志愿長(zhǎng)安大學(xué)材料與化工309分請(qǐng)求調(diào)劑已經(jīng)有8人回復(fù)
282分材料專(zhuān)業(yè)求調(diào)劑院校已經(jīng)有14人回復(fù)
085600材料與化工 326 求調(diào)劑已經(jīng)有4人回復(fù)
085600材料與化工，一志愿廣州985，求調(diào)劑已經(jīng)有15人回復(fù)
環(huán)境調(diào)劑已經(jīng)有4人回復(fù)
材料調(diào)劑，307分已經(jīng)有9人回復(fù)
0817一志愿蘇州大學(xué)280復(fù)試調(diào)劑已經(jīng)有5人回復(fù)
化工284求調(diào)劑已經(jīng)有11人回復(fù)
2026級(jí)碩士研究生招生/調(diào)劑已經(jīng)有7人回復(fù)

1樓 2019-01-30 16:45:55

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

XiaoGuads

新蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 51.2
帖子: 17
在線: 3小時(shí)
蟲(chóng)號(hào): 8840508
注冊(cè): 2018-05-23
專(zhuān)業(yè): 物理學(xué)I

送紅花一朵

引用回帖:

3樓: Originally posted by hzlhm at 2019-01-31 14:39:33
利用最小二乘法（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 ...

我在這方面是小白，請(qǐng)問(wèn)大神能否貼出代碼呢，謝謝！

贊一下

回復(fù)此樓

4樓2019-02-04 16:56:15

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

查看全部 5 個(gè)回答

竹一拿下

銅蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 0 (幼兒園)
金幣: 15.9
散金: 1367
紅花: 10
帖子: 942
在線: 260.9小時(shí)
蟲(chóng)號(hào): 7380745
注冊(cè): 2017-10-27
性別: GG
專(zhuān)業(yè): 力學(xué)

★
jjdg: 金幣+1, 感謝參與 2019-01-31 00:28:29

試試牛頓

發(fā)自小木蟲(chóng)Android客戶端

回復(fù)此樓

2樓2019-01-30 17:49:49

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

hzlhm

至尊木蟲(chóng) (著名寫(xiě)手)

應(yīng)助: 387 (碩士)
金幣: 17808.5
紅花: 53
帖子: 2879
在線: 606.3小時(shí)
蟲(chóng)號(hào): 1517335
注冊(cè): 2011-11-30
性別: GG
專(zhuān)業(yè): 常微分方程與動(dòng)力系統(tǒng)

【答案】應(yīng)助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(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

贊一下(1人)

回復(fù)此樓

» 本帖已獲得的紅花（最新10朵）

XiaoGuads

QQ：2120156492

3樓2019-01-31 14:39:33

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

dingd

鐵桿木蟲(chóng) (職業(yè)作家)

應(yīng)助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時(shí)
蟲(chóng)號(hào): 291104
注冊(cè): 2006-10-28

【答案】應(yīng)助回帖

1stOpt試試，不需初值，精度更高：

CODE:

Parameter x(4);

Constant alpha=0.125*pi/180,xx1=0,z1=5,xx2=0,z2=25,L1=5,L2=5,f=30,h0=20,n2=1.49386;

Function

tan(alpha) - x(1)/x(2);

n2*sin(atan((xx1 - x(1))/(z1 - x(2))) - atan((x(2) - x(4))/(x(1) - x(3)))) - sin(pi/2 - alpha + atan((xx1 - x(1))/(z1 - x(2)))) ;

n2*sin(atan((xx2 - x(3))/(z2 - x(4))) - atan((x(2) - x(4))/(x(1) - x(3)))) - sin(pi/2 + atan((xx2 - x(3))/(z2 - x(4))));

L1 + L2 - f + h0 - h0*n2 + n2*((x(1) - x(3))^2 + (x(2) - x(4))^2)^(1/2);

有周期性三角函數(shù)，多解：
1：
x1: 0.0108298192893512
x2: 4.96401563027804
x3: 11.024819043045
x4: -11.7300564127835

2：
x1: 0.0109088181961899
x2: 5.0000209101071
x3: -0.00562667845017508
x4: 25.00001407242

贊一下

回復(fù)此樓

5樓2019-02-28 18:14:16

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

查看全部 5 個(gè)回答

普通表情龍兔虎貓

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 一志愿安徽大學(xué)材料工程專(zhuān)碩313分，求調(diào)劑的學(xué)校 +3	Yu先生 2026-03-10	3/150	2026-03-10 08:34 by Hongeo93
[考研] 材料調(diào)劑，307分 +9	張泳銘1 2026-03-09	9/450	2026-03-10 08:21 by Linda Hu
[考研] 320求調(diào)劑 +4	魏zy 2026-03-08	4/200	2026-03-09 16:14 by ruiyingmiao
[考研] 0703化學(xué) 學(xué)碩理工科均可不區(qū)分研究方向總分279求調(diào)劑 +9	1一11 2026-03-03	9/450	2026-03-09 13:51 by macy2011
[考研] 313分生物學(xué)求調(diào)劑 +4	Yyt楊1 2026-03-09	5/250	2026-03-09 13:42 by macy2011
[考研] 一志愿211 材料與化工 280求調(diào)劑 +3	Sanity蔣 2026-03-08	3/150	2026-03-09 06:35 by houyaoxu
[考研] 招調(diào)劑 +20	帆船哥 2026-03-04	23/1150	2026-03-08 16:25 by taycat
[考研] 085700資環(huán)求調(diào)劑，初始279，六級(jí)已過(guò)，英語(yǔ)能力強(qiáng) +4	085700資環(huán)調(diào)劑 2026-03-03	5/250	2026-03-08 13:35 by 30660438
[考研] 化學(xué)工程求調(diào)劑 +12	化工人999 2026-03-04	12/600	2026-03-08 11:17 by 我的船我的海
[考研] 材料工程專(zhuān)碩調(diào)劑自薦信，初試323 +4	李白26 2026-03-07	4/200	2026-03-07 21:59 by 星空星月
[考研] 0307化學(xué)求調(diào)劑 +6	0ok0k 2026-03-07	6/300	2026-03-07 20:10 by pies112
[考研] 278求調(diào)劑 +5	Gale1314 2026-03-06	5/250	2026-03-07 14:41 by 2735147993
[考研] 一志愿南京大學(xué)資源與環(huán)境專(zhuān)業(yè)283分求調(diào)劑 +8	追云?? 2026-03-04	8/400	2026-03-07 09:23 by 斬魂滴兔子！
[考研] 求調(diào)劑 +4	呼呼？~+123456 2026-03-05	5/250	2026-03-06 23:15 by L135790
[考研] 求調(diào)劑 +3	泡了個(gè)椒 2026-03-04	4/200	2026-03-05 14:37 by 泡了個(gè)椒
[考研] 化學(xué) 0703求調(diào)劑總分293 一志愿211 +9	土土小蟲(chóng) 2026-03-03	11/550	2026-03-05 09:26 by vfdbghn
[考研] 376求調(diào)劑 +3	王浩然sam 2026-03-04	3/150	2026-03-05 07:48 by bxbo
[考研] 322分 085600求調(diào)劑，有互聯(lián)網(wǎng)+國(guó)金及主持省級(jí)大創(chuàng)經(jīng)歷 +6	熊境喆 2026-03-04	6/300	2026-03-04 20:32 by kakakapanpan
[考研] 學(xué)碩材料275調(diào)劑 +9	路三三 2026-03-03	9/450	2026-03-04 17:02 by 夢(mèng)天888
[考研] 復(fù)試調(diào)劑 +4	杰尼龜aaa 2026-03-03	4/200	2026-03-03 23:06 by zhukairuo