版塊導(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)有253人回復(fù)
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有0人回復(fù)
幫我的英語口語老師找學(xué)生已經(jīng)有3人回復(fù)
計算三階力常數(shù)用thirdorder擴(kuò)胞時參數(shù)如何選擇？臨近原子數(shù)怎么取值已經(jīng)有5人回復(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ā)表

[考研] 085701環(huán)境工程295求調(diào)劑，一志愿A區(qū)211 +7	等春來， 2026-03-06	8/400	2026-03-06 22:04 by 星空星月
[考研] 266材料化工求調(diào)劑 +5	哇塞王帥 2026-03-03	5/250	2026-03-06 16:44 by 學(xué)員FnSWZj
[考研] 材料工程調(diào)劑 +6	zbcm_zbcm 2026-03-05	6/300	2026-03-06 15:36 by @颯颯颯颯
[考研] 332材料求調(diào)劑 +6	zjy101327 2026-03-05	7/350	2026-03-06 00:18 by wutongshun
[考研] 【求調(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
[考研] 復(fù)試調(diào)劑 +5	呼呼？~+123456 2026-03-05	5/250	2026-03-05 19:43 by caszguilin
[考研] 歡迎采礦、地質(zhì)、巖土、計算機(jī)、人工智能等專業(yè)的同學(xué)報考 +10	pin8023 2026-02-28	13/650	2026-03-05 18:54 by 3367308522
[考研] 301求調(diào)劑 +4	李LJR 2026-03-04	4/200	2026-03-05 13:31 by Iveryant
[考研] 283求調(diào)劑 +10	鹿沫笙 2026-03-02	11/550	2026-03-05 13:12 by ZHANGTDK
[考研] 376求調(diào)劑 +3	王浩然sam 2026-03-04	3/150	2026-03-05 07:48 by bxbo
[考研] 0703化學(xué)求調(diào)劑，總分320分，一志愿中南民族大學(xué) +3	19286712440 2026-03-04	3/150	2026-03-04 18:30 by houyaoxu
[考研] 295求調(diào)劑 +4	小賽不吃香菜 2026-03-04	4/200	2026-03-04 16:13 by sslc1985
[考研] 085601 材料305分求助 +4	泡泡郵件 2026-03-03	6/300	2026-03-04 13:00 by 無際的草原
[考研] 291求調(diào)劑 +3	MuoLuo1312 2026-03-02	6/300	2026-03-03 12:13 by 熱情沙漠
[考研] 321求調(diào)劑一志愿東北林業(yè)大學(xué)材料與化工英二數(shù)二 +5	蟲蟲蟲蟲蟲7 2026-03-01	9/450	2026-03-02 21:33 by sunny81
[考研] 285求調(diào)劑 +9	滿頭大汗的學(xué)生 2026-02-28	9/450	2026-03-02 20:29 by hypershenger
[考研] 材料調(diào)劑 +6	愛擦汗的可樂冰 2026-02-28	7/350	2026-03-02 10:42 by Jy?
[考研] 材料類求調(diào)劑 +11	wana_kiko 2026-02-28	14/700	2026-03-02 08:46 by 聰明的大松鼠
[考研] 317一志愿華南理工電氣工程求調(diào)劑 +6	Soliloquy_Q 2026-02-28	11/550	2026-03-01 11:14 by 歌liekkas
[考研] 307求調(diào)劑 +4	73372112 2026-02-28	6/300	2026-03-01 00:04 by ll247

24小時熱門版塊排行榜

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

» 猜你喜歡

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

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

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