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

登錄注冊

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

24小時熱門版塊排行榜

返回列表

是最后

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 9
帖子: 9
在線: 12.8小時
蟲號: 14633322
注冊: 2019-03-17
專業(yè): 傳動機械學(xué)

[求助] 多元非線性方程組求解，急求已有1人參與

這個4元非線性方程組怎么解啊？急求各位大神
function f=fun(x0)
syms x1 x2 x3 x4
x0=[x1 x2 x3 x4];

fi=0.515;
fe=0.53;
B=fi+fe-1;
d=pi/180;
a0=45*d;
Db=7.144;
Dm=40;
n=1200;
L=Db/Dm;
Z=atan(10/(pi*40));
p=7.801e-6;
E1=2.07e5;
E2=1.9e5;
u1=0.29;
u2=0.305;
Ro=17.15;
Ri=24.294;
ro=3.872;
ri=3.872;

g=x2/sin(x4)+x1/sin(x3);
A1=B*Db*sin(a0)+g;
A2=B*Db*cos(a0);
w1=sin(x3)*(0.03*Db+x1);
w2=A2-cos(x4)*(0.015*Db+x2);

v=atan(sin(x3)/(L+cos(x3)));
Wm=2*pi*n/(1+(1+L*cos(x3))*(cos(x4)+tan(v)*sin(x4))/((1-L*cos(x4))*(cos(x3)+tan(v)*sin(x3))));
Wr=-Wm*(1+L*cos(x3))*cos(Z)/(L*(sin(v)*sin(x3)+cos(v)*cos(x3)));

Fc=(pi/12)*p*(Db^3)*Dm*(Wm^2);
Mg=(1/60)*p*pi*(Db^5)*Wm*Wr*sin(v);

Ps11=2/Db;
Ps12=2/Db;
Pn11=2/Db;
Pn12=2/Db;
Ps21=1/Ro;
Pn21=-1/Ri;

Ps22=-ro;
Pn22=-ri;
S1=Ps11+Ps12+Ps21+Ps22;
S2=Pn11+Pn12+Pn21+Pn22;

Rs1=Ps11*Ps21/(Ps11+Ps21);
Rs2=Ps12*Ps22/(Ps12+Ps22);
Rn1=Pn11*Pn21/(Pn11+Pn21);
Rn2=Pn12*Pn22/(Pn12+Pn22);

K1=1.0339*(Rs2/Rs1)^0.636;
K2=1.0339*(Rn2/Rn1)^0.636;
Ke1=1.5277+0.6023*log(Rs2/Rs1);
Ke2=1.5277+0.6023*log(Rn2/Rn1);
Ee1=1.0003+0.5968*(Rs1/Rs2);
Ee2=1.0003+0.5968*(Rn1/Rn2);

e1=(1-1/(K1^2))^0.5;
e2=(1-1/(K2^2))^0.5;

y1=(2*Ke1/pi)*((1-(e1)^2)*pi/(2*Ee1))^(1/3);
y2=(2*Ke2/pi)*((1-(e2)^2)*pi/(2*Ee2))^(1/3);

E=0.5*((1-u1^2)/E1+(1-u2^2)/E2);
Q1=(8*(x1/y2)^3/((1.5*E)^2*S2))^0.5;
Q2=(8*(x2/y1)^3/((1.5*E)^2*S1))^0.5;

f1=(A1-w1)^2+(A2-w2)^2-(0.015*Db+x4)^2;
f2=w1^2+w2^2-(0.03*Db+x3);
f3=Q2*sin(x4)-Q1*sin(x3)-Mg/Db*(cos(x4)-cos(x3));
f4=Q2*cos(x4)-Q1*cos(x3)+Mg/Db*(sin(x4)-sin(x3))+Fc;

f=[f1; f2; f3 ;f4];

% format long;
% f=vpa(f,2);
% class(f);
end

回復(fù)此樓

» 猜你喜歡

重慶交通大學(xué)光子學(xué)微結(jié)構(gòu)與器件課題組招收2026年碩士研究生信息已經(jīng)有2人回復(fù)
一志愿鄭大材料學(xué)碩298分，求調(diào)劑已經(jīng)有6人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費? 已經(jīng)有220人回復(fù)
尋合作：應(yīng)力腐蝕多尺度模擬已經(jīng)有6人回復(fù)
考研交流已經(jīng)有0人回復(fù)
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有0人回復(fù)
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有0人回復(fù)
幫我的英語口語老師找學(xué)生已經(jīng)有3人回復(fù)

1樓 2019-08-11 14:55:31

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

獨孤神宇

版主 (知名作家)

應(yīng)助: 490 (碩士)
貴賓: 0.008
金幣: 31016.3
散金: 802
紅花: 122
沙發(fā): 1
帖子: 5600
在線: 856.5小時
蟲號: 3522474
注冊: 2014-11-06
性別: GG
專業(yè): 機械動力學(xué)
管轄: 計算模擬

直接用fsolve函數(shù)

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

贊一下

回復(fù)此樓

數(shù)值計算

2樓2019-08-11 15:05:19

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

是最后

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 9
帖子: 9
在線: 12.8小時
蟲號: 14633322
注冊: 2019-03-17
專業(yè): 傳動機械學(xué)

引用回帖:

2樓: Originally posted by 獨孤神宇 at 2019-08-11 15:05:19
直接用fsolve函數(shù)

fsolve，求解的不是自己想要的，跟初始值設(shè)置數(shù)值差不多，請問有沒有具體的解法可以求

贊一下

回復(fù)此樓

3樓2019-08-11 15:42:26

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

獨孤神宇

版主 (知名作家)

應(yīng)助: 490 (碩士)
貴賓: 0.008
金幣: 31016.3
散金: 802
紅花: 122
沙發(fā): 1
帖子: 5600
在線: 856.5小時
蟲號: 3522474
注冊: 2014-11-06
性別: GG
專業(yè): 機械動力學(xué)
管轄: 計算模擬

引用回帖:

3樓: Originally posted by 是最后 at 2019-08-11 15:42:26
fsolve，求解的不是自己想要的，跟初始值設(shè)置數(shù)值差不多，請問有沒有具體的解法可以求...

換初始值計算。

回復(fù)此樓

數(shù)值計算

4樓2019-08-11 15:56:56

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

是最后

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 9
帖子: 9
在線: 12.8小時
蟲號: 14633322
注冊: 2019-03-17
專業(yè): 傳動機械學(xué)

引用回帖:

4樓: Originally posted by 獨孤神宇 at 2019-08-11 15:56:56
換初始值計算。...

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.

<stopping criteria details>

x1 =

0.264988629986718 - 0.227162298915769i 0.157276044576307 - 0.503604785301162i 53.759850087022613 - 0.866858320147245i 34.658568519774903 + 0.590132528775326i

得出的是這個，這不算求出解啊

贊一下

回復(fù)此樓

5樓2019-08-11 16:51:12

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

獨孤神宇

版主 (知名作家)

應(yīng)助: 490 (碩士)
貴賓: 0.008
金幣: 31016.3
散金: 802
紅花: 122
沙發(fā): 1
帖子: 5600
在線: 856.5小時
蟲號: 3522474
注冊: 2014-11-06
性別: GG
專業(yè): 機械動力學(xué)
管轄: 計算模擬

初始值設(shè)置不對，多試幾次其他的值

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

贊一下

回復(fù)此樓

數(shù)值計算

6樓2019-08-11 19:05:08

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

dingd

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

計算強帖: 4
應(yīng)助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時
蟲號: 291104
注冊: 2006-10-28

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

★ ★
感謝參與，應(yīng)助指數(shù) +1
獨孤神宇: 金幣+2, 鼓勵交流 2019-08-11 22:12:35

公式有問題吧：
rs1 = 0.0482578901650419
rs2 = 0.301774232969708
rn1 = -0.048257890165042
rn2 = 0.301774232969708

Ke2=1.5277+0.6023*log(Rn2/Rn1)

上式中Rn2/Rn1為負數(shù)，再取對數(shù)？能計算嗎？

贊一下

回復(fù)此樓

7樓2019-08-11 21:43:49

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

southpark007

禁蟲 (著名寫手)

本帖內(nèi)容被屏蔽

8樓2019-08-11 22:58:13

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

dingd

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

計算強帖: 4
應(yīng)助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時
蟲號: 291104
注冊: 2006-10-28

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

公式邏輯錯誤，能否求解和軟件無關(guān)了。

贊一下

回復(fù)此樓

9樓2019-08-12 09:58:12

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

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主是最后的主題更新

返回列表

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

[考研] 求調(diào)劑，學(xué)校研究所都可以，材料與化工267分 +6	wmx1 2026-03-05	6/300	2026-03-05 20:01 by houyaoxu
[考研] 歡迎采礦、地質(zhì)、巖土、計算機、人工智能等專業(yè)的同學(xué)報考 +10	pin8023 2026-02-28	13/650	2026-03-05 18:54 by 3367308522
[考研] 304求調(diào)劑 +3	曼殊2266 2026-03-04	3/150	2026-03-05 10:39 by Iveryant
[考研] 化學(xué) 0703求調(diào)劑總分293 一志愿211 +9	土土小蟲 2026-03-03	11/550	2026-03-05 09:26 by vfdbghn
[考研] 材料復(fù)試調(diào)劑 +7	學(xué)材料的點 2026-03-01	8/400	2026-03-04 21:08 by wutongshun
[考研] 322分 085600求調(diào)劑，有互聯(lián)網(wǎng)+國金及主持省級大創(chuàng)經(jīng)歷 +6	熊境喆 2026-03-04	6/300	2026-03-04 20:32 by kakakapanpan
[考研] 一志愿中石油，化學(xué)工程，301求調(diào)劑 +3	豚饅枳枳 2026-03-04	3/150	2026-03-04 18:32 by houyaoxu
[考研] 085600 材料與化工 298 +14	小西笑嘻嘻 2026-03-03	14/700	2026-03-04 16:23 by sslc1985
[考研] 295求調(diào)劑 +4	小賽不吃香菜 2026-03-04	4/200	2026-03-04 16:13 by sslc1985
[考研] 0857求調(diào)劑 +4	CDDD1013 2026-03-04	4/200	2026-03-04 16:08 by zhukairuo
[考研] 本科太原理工采礦工程，求調(diào)劑 +3	onlx 2026-03-01	3/150	2026-03-04 15:57 by Stephen_ym
[考研] 322,求調(diào)劑 +3	菜菜愛玩 2026-03-04	3/150	2026-03-04 12:15 by xiongyaxuan
[考研] 085600求調(diào)劑 +4	LRZZZZZZ 2026-03-02	6/300	2026-03-04 09:11 by 無際的草原
[考研] 289求調(diào)劑 +7	BrightLL 2026-03-02	9/450	2026-03-03 16:02 by tgxtgxtgx9
[考研] 計算機學(xué)碩分數(shù)285求調(diào)劑 +4	glwshine 2026-03-02	5/250	2026-03-03 14:27 by king呀
[考研] 0856材料調(diào)劑 +5	沿岸有貝殼OUC 2026-03-02	5/250	2026-03-02 20:31 by hypershenger
[考研] 材料調(diào)劑 +3	恒順自然 2026-03-02	3/150	2026-03-02 18:49 by L135790
[考研] 一志愿東北大學(xué)材料專碩328，求調(diào)劑 +3	shs1083 2026-03-02	3/150	2026-03-02 17:27 by houyaoxu
[考研] 291 求調(diào)劑 +3	化工2026屆畢業(yè)?/a> 2026-03-02	3/150	2026-03-02 12:55 by houyaoxu
[考研] 材料類求調(diào)劑 +11	wana_kiko 2026-02-28	14/700	2026-03-02 08:46 by 聰明的大松鼠

24小時熱門版塊排行榜

[求助] 多元非線性方程組求解，急求 已有1人參與

» 猜你喜歡

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

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

[求助] 多元非線性方程組求解，急求已有1人參與