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

登錄注冊

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

24小時熱門版塊排行榜

返回列表

是最后

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 9
帖子: 9
在線: 12.8小時
蟲號: 14633322
注冊: 2019-03-17
專業(yè): 傳動機(jī)械學(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ù)此樓

» 猜你喜歡

考研交流已經(jīng)有0人回復(fù)
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有0人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費? 已經(jīng)有128人回復(fù)
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有0人回復(fù)
幫我的英語口語老師找學(xué)生已經(jīng)有3人回復(fù)
計算三階力常數(shù)用thirdorder擴(kuò)胞時參數(shù)如何選擇？臨近原子數(shù)怎么取值已經(jīng)有3人回復(fù)
2026年循環(huán)經(jīng)濟(jì)功能材料國際會議（ICFMCE 2026）已經(jīng)有0人回復(fù)
2026年第五屆電氣、電子與信息工程國際會議（ISEEIE 2026）已經(jīng)有0人回復(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è): 機(jī)械動力學(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è): 傳動機(jī)械學(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è): 機(jī)械動力學(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è): 傳動機(jī)械學(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è): 機(jī)械動力學(xué)
管轄: 計算模擬

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

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

贊一下

回復(fù)此樓

數(shù)值計算

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

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

dingd

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

計算強(qiáng)帖: 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為負(fù)數(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è)作家)

計算強(qiáng)帖: 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ā)表

[考研] 歡迎211本科同學(xué)，過A區(qū)國家線，A區(qū)非偏遠(yuǎn)一本，交叉學(xué)科課題組 +18	lisimayy 2026-03-04	30/1500	2026-03-06 15:42 by lisimayy
[考研] 0703 東華大學(xué) 理學(xué)化學(xué)均過 a 線，總分 281 求調(diào)劑 +6	香香 Lu 2026-03-03	11/550	2026-03-06 15:37 by @颯颯颯颯
[考研] 08工科求調(diào)劑 +3	隆LLL 2026-03-06	4/200	2026-03-06 15:26 by 隆LLL
[考研] 085602高分子方向求調(diào)劑 +7	tlgudy 2026-03-04	7/350	2026-03-06 14:17 by guosr9609
[考研] 332材料求調(diào)劑 +6	zjy101327 2026-03-05	7/350	2026-03-06 00:18 by wutongshun
[考博] 26申博-目前4篇SCI一作 +4	chen_2024 2026-03-02	4/200	2026-03-05 22:59 by sicilyl0001
[考研] 【求調(diào)劑】293分環(huán)境工程求調(diào)劑材料/化工，服從調(diào)劑，抗壓能力強(qiáng)！ +8	xiiiia 2026-03-04	8/400	2026-03-05 20:01 by wangjihu
[考研] 085602 293分求調(diào)劑 +3	SivanNano. 2026-03-05	3/150	2026-03-05 19:41 by caszguilin
[考研] 304求調(diào)劑 +4	曼殊2266 2026-03-05	4/200	2026-03-05 17:10 by zhukairuo
[考研] 283求調(diào)劑 +10	鹿沫笙 2026-03-02	11/550	2026-03-05 13:12 by ZHANGTDK
[考研] 材料學(xué)碩080500復(fù)試調(diào)劑294 +3	四葉zjz 2026-03-04	3/150	2026-03-05 07:09 by kunm555
[考研] 求調(diào)劑院校 +6	云朵452 2026-03-02	12/600	2026-03-04 23:17 by 云朵452
[考研] 成績276，專業(yè)代碼0856求調(diào)劑 +10	小陳朵 2026-03-03	10/500	2026-03-04 20:16 by 1124861332
[考研] 298求調(diào)劑一志愿中海洋 +3	lour. 2026-03-03	3/150	2026-03-03 20:41 by gxg2025
[考研] 江蘇省農(nóng)科院招調(diào)劑1名 +5	Qwertyuop 2026-03-01	5/250	2026-03-03 17:50 by 梅子鴨鴨
[考研] 材料270求調(diào)劑 6+6	Eiiiio 2026-03-01	11/550	2026-03-03 15:04 by Phoebelala
[考研] 材料工程求調(diào)劑 +3	1431251 2026-03-03	3/150	2026-03-03 11:58 by EBSD
[考研] 288求調(diào)劑 +3	少71.8 2026-03-02	5/250	2026-03-03 06:01 by tgxtgxtgx9
[考研] 295求調(diào)劑 +8	19171856320 2026-02-28	8/400	2026-03-02 11:19 by yuchj
[考研] 調(diào)劑 +3	簡木ChuFront 2026-02-28	3/150	2026-03-01 11:46 by 王偉要上岸啊

24小時熱門版塊排行榜

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

» 猜你喜歡

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

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

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