 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1892 | 回復(fù): 8 | ||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||
是最后新蟲 (初入文壇)
|
[求助]
多元非線性方程組求解，急求 已有1人參與
|
|
|
這個(gè)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é)生
已經(jīng)有4人回復(fù)
-
計(jì)算三階力常數(shù)用thirdorder擴(kuò)胞時(shí)參數(shù)如何選擇？臨近原子數(shù)怎么取值
已經(jīng)有7人回復(fù)
-
物理學(xué)I論文潤色/翻譯怎么收費(fèi)?
已經(jīng)有58人回復(fù)
-
2026年循環(huán)經(jīng)濟(jì)功能材料國際會(huì)議（ICFMCE 2026）
已經(jīng)有0人回復(fù)
-
2026年第五屆電氣、電子與信息工程國際會(huì)議（ISEEIE 2026）
已經(jīng)有0人回復(fù)
-
哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑
已經(jīng)有0人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑
已經(jīng)有2人回復(fù)
-
求調(diào)劑
已經(jīng)有0人回復(fù)
獨(dú)孤神宇
版主 (知名作家)
- 應(yīng)助: 490 (碩士)
- 貴賓: 0.008
- 金幣: 31016.3
- 散金: 802
- 紅花: 122
- 沙發(fā): 1
- 帖子: 5600
- 在線: 856.5小時(shí)
- 蟲號: 3522474
- 注冊: 2014-11-06
- 性別: GG
- 專業(yè): 機(jī)械動(dòng)力學(xué)
- 管轄: 計(jì)算模擬
4樓2019-08-11 15:56:56
獨(dú)孤神宇
版主 (知名作家)
- 應(yīng)助: 490 (碩士)
- 貴賓: 0.008
- 金幣: 31016.3
- 散金: 802
- 紅花: 122
- 沙發(fā): 1
- 帖子: 5600
- 在線: 856.5小時(shí)
- 蟲號: 3522474
- 注冊: 2014-11-06
- 性別: GG
- 專業(yè): 機(jī)械動(dòng)力學(xué)
- 管轄: 計(jì)算模擬
2樓2019-08-11 15:05:19
是最后
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 9
- 帖子: 9
- 在線: 12.8小時(shí)
- 蟲號: 14633322
- 注冊: 2019-03-17
- 專業(yè): 傳動(dòng)機(jī)械學(xué)
3樓2019-08-11 15:42:26
是最后
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 9
- 帖子: 9
- 在線: 12.8小時(shí)
- 蟲號: 14633322
- 注冊: 2019-03-17
- 專業(yè): 傳動(dòng)機(jī)械學(xué)
|
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 得出的是這個(gè)，這不算求出解啊 |
5樓2019-08-11 16:51:12
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料調(diào)劑，307分 +9 | 張泳銘1 2026-03-09 | 9/450 |
|
|---|---|---|---|---|
|
|
[考研] 一志愿北京化工大學(xué)材料與化工296分求調(diào)劑 +8 | 稻妻小編 2026-03-09 | 8/400 |
|
|
|
[考研] 320材料與化工，求調(diào)劑 +11 | 鶴遨予卿 2026-03-04 | 13/650 |
|
|
[考研] 313求調(diào)劑 +4 | Yyt楊1 2026-03-07 | 5/250 |
|
|
|
[考研] 313分生物學(xué)求調(diào)劑 +4 | Yyt楊1 2026-03-09 | 5/250 |
|
|
|
[考研] 中科大材料299求調(diào)劑 +10 | DAIjiayo 2026-03-05 | 16/800 |
|
|
|
[考研] 材料調(diào)劑 +4 | xxxcm 2026-03-08 | 7/350 |
|
|
|
[考研] 0856求調(diào)劑 +3 | squirtle11 2026-03-07 | 3/150 |
|
|
|
[考研] 0703化學(xué)求調(diào)劑，總分309分，一志愿華南師范 +3 | 花與葉@ 2026-03-08 | 3/150 |
|
|
|
[考研] 理學(xué)287分求調(diào)劑 +4 | 看看我. 2026-03-05 | 4/200 |
|
|
|
[考研] 求調(diào)劑 +4 | 呼呼？~+123456 2026-03-05 | 5/250 |
|
|
|
[考研] 306求調(diào)劑 +7 | Bahati 2026-03-05 | 7/350 |
|
|
|
[考研] 不限學(xué)校專業(yè)的調(diào)劑同學(xué)看過來 +5 | 啊擺啊擺 2026-03-05 | 9/450 |
|
|
|
[考研] 材料考研339求調(diào)劑 +3 | Karry*^_^* 2026-03-04 | 3/150 |
|
|
|
[考研] 328求調(diào)劑 +6 | vuzhdkfjkx 2026-03-04 | 7/350 |
|
|
|
[考研] 347求調(diào)劑 +6 | 啊歐歐歐 2026-03-03 | 8/400 |
|
|
|
[考研] 材料085600 303求調(diào)劑 +7 | 1bygone 2026-03-04 | 7/350 |
|
|
|
[考研] 267調(diào)劑求助 +5 | 聰少OZ 2026-03-04 | 5/250 |
|
|
|
[考研] 322分 085600求調(diào)劑，有互聯(lián)網(wǎng)+國金及主持省級大創(chuàng)經(jīng)歷 +6 | 熊境喆 2026-03-04 | 6/300 |
|
|
|
[考研] 266求調(diào)劑 +7 | 哇塞王帥 2026-03-03 | 7/350 |
|