people891108

木蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 4137.1
散金: 38
紅花: 7
帖子: 466
在線: 192.6小時(shí)
蟲號(hào): 2122748
注冊(cè): 2012-11-12
專業(yè): 化學(xué)反應(yīng)工程

[求助] 用1stopt解非線性方程組問題已有1人參與

本人剛開始接觸1stopt,最開始用matlab求解結(jié)果不好，打算用1stopt求解非線性方程組，下面是我的程序，解的方法選用包維爾法，每次解一次方程組需要迭代很長時(shí)間，求大神幫我優(yōu)化這個(gè)程序和算法，第二個(gè)方程0.7491=x3+x4+x5中0.7491是個(gè)變量，每次我都手動(dòng)修改數(shù)值，希望大神能幫我將其他數(shù)值[0.5549; 0.6009; 0.7181; 0.7491; 1.0403]也能編進(jìn)去，謝謝

程序
Title "concentration";
Parameters x(1:8)[0,];
//Variable ;
Function 1.678=x1+x2;
0.749151=x3+x4+x5;
50=x6+x7+x8+x4+x5;
0=x2+x7-x4-2*x5-x8;
3.82024*10^(-13)= x7*x8;
0.43996*10^(-6)=(x7*x4/(x3));
0.74128*10^(-10)=(x7*x5/(x4));
0.039833*10^(-6)=(x7*x1/(x2));

回復(fù)此樓

» 猜你喜歡

2026年第五屆電氣、電子與信息工程國際會(huì)議（ISEEIE 2026）已經(jīng)有0人回復(fù)
哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑已經(jīng)有6人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費(fèi)? 已經(jīng)有58人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑已經(jīng)有5人回復(fù)
求調(diào)劑已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有0人回復(fù)
請(qǐng)問還有沒有用Latex寫文章的小伙伴們？已經(jīng)有0人回復(fù)
光學(xué)工程學(xué)碩調(diào)劑信息已經(jīng)有30人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
散金幣，求好運(yùn)，祝面上順利！已經(jīng)有65人回復(fù)
華東師范大學(xué)芯片設(shè)計(jì)徐瓏真實(shí)水平咋樣已經(jīng)有1人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

用1stopt自定義函數(shù)來擬合非線性曲線，函數(shù)中有復(fù)數(shù)，為什么得不到結(jié)果啊？？？已經(jīng)有5人回復(fù)
用1stOpt解出非線性方程組的解作為初值用Matlab求解，解不出來已經(jīng)有4人回復(fù)
看了下1stOpt的介紹，感覺略顯浮夸已經(jīng)有14人回復(fù)
MATLAB解非線性方程組已經(jīng)有9人回復(fù)
請(qǐng)教高手，如何解非線性方程組�。�！已經(jīng)有13人回復(fù)
非線性方程的線性化問題，大神們進(jìn)來指教下�。。�！感謝了�。。。。。。。�！已經(jīng)有9人回復(fù)
1stopt解方程組已經(jīng)有12人回復(fù)
牛頓迭代法解非線性矩陣方程已經(jīng)有4人回復(fù)
用matlab求解一個(gè)非線性方程組的解已經(jīng)有4人回復(fù)
用matlab求解非線性方程組說無解，一定是方程組本身無解，還是有可能程序有問題呢？已經(jīng)有11人回復(fù)
sos!!1stopt解非線性方程組，一直出現(xiàn)常數(shù)定義錯(cuò)誤的問題！求大神幫助！已經(jīng)有12人回復(fù)
matlab解非線性方程組解不開已經(jīng)有3人回復(fù)
用1stopt擬合非線性方程結(jié)果與其他軟件擬合結(jié)果差異大已經(jīng)有5人回復(fù)
1stopt運(yùn)行程序，總是出現(xiàn)字符函數(shù)錯(cuò)誤的提示，請(qǐng)高手幫助已經(jīng)有5人回復(fù)
用1stopt或者origin非線性擬合磁導(dǎo)率已經(jīng)有5人回復(fù)
1stopt4階非線性常微分方程，幫小弟看一下，跪求。已經(jīng)有7人回復(fù)
1stopt或matlab如何求解以下的非線性方程并擬合出相應(yīng)曲線？已經(jīng)有13人回復(fù)
【求助】求教matlab解非線性方程組已經(jīng)有9人回復(fù)

1樓 2015-05-12 17:21:57

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

dingd

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

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

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

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
people891108: 金幣+20, ★★★很有幫助, 非常感謝 2015-05-13 09:04:34

循環(huán)常數(shù) a x1 x2 x3 x4 x5 x6 x7 x8
0.5549 1.14173198863636 0.536268011363644 0.0224757584546091 0.528580867483308 0.00384337406208307 49.4675753252834 1.87075303955625E-8 4.14463700340161E-7
0.6009 1.10180642808357 0.576193571916431 0.027075916130922 0.571866314491743 0.00195776937733538 49.4257641557996 2.08306045808516E-8 0.000411739500621744
0.7181 0.999546948870126 0.678453051129874 0.0414689554283631 0.674809322011355 0.00182172256028892 49.3233686175059 2.70367977046525E-8 3.11034738591439E-7
0.7491 0.973269657335869 0.704730342664131 0.0459880433975497 0.701493763743082 0.00161819285935604 49.2968877925099 2.88425929463931E-8 2.22044929482337E-7
1.0403 0.743218579157401 0.934781420842599 0.106317357324403 0.933709132758929 0.000273509916853229 49.0654919888832 5.00963406693268E-8 0.000525318346303248

贊一下

回復(fù)此樓

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

people891108

2樓2015-05-12 20:14:30

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

people891108

木蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 4137.1
散金: 38
紅花: 7
帖子: 466
在線: 192.6小時(shí)
蟲號(hào): 2122748
注冊(cè): 2012-11-12
專業(yè): 化學(xué)反應(yīng)工程

送紅花一朵

引用回帖:

2樓: Originally posted by dingd at 2015-05-12 20:14:30
循環(huán)常數(shù) a x1 x2 x3 x4 x5 x6 x7 x8
0.5549 1.14173198863636 0.536268011363644 0.0224757584546091 0.528580867483308 0.00384337406208307 49.4675753252834 1.87075303955625E-8 4.14463700340161E-7
0.6009 ...

請(qǐng)問能把程序代碼發(fā)給我么，我用循環(huán)常數(shù)試了一下，結(jié)果只出一個(gè)解

贊一下

回復(fù)此樓

3樓2015-05-13 09:42:26

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

dingd

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

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

引用回帖:

3樓: Originally posted by people891108 at 2015-05-13 09:42:26
請(qǐng)問能把程序代碼發(fā)給我么，我用循環(huán)常數(shù)試了一下，結(jié)果只出一個(gè)解...

如果你用的是老版本就別試了，無法進(jìn)行這種循環(huán)計(jì)算。

贊一下

回復(fù)此樓

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

情弦風(fēng)

4樓2015-05-13 09:51:11

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

情弦風(fēng)

新蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 907.2
散金: 145
帖子: 95
在線: 59.6小時(shí)
蟲號(hào): 3754702
注冊(cè): 2015-03-21
專業(yè): 磁共振結(jié)構(gòu)成像與疾病診斷

送紅花一朵

引用回帖:

4樓: Originally posted by dingd at 2015-05-13 09:51:11
如果你用的是老版本就別試了，無法進(jìn)行這種循環(huán)計(jì)算。...

大神我需要求解5600個(gè)循環(huán)的非線性方程組，可以把代碼貼出來嗎

贊一下

回復(fù)此樓

5樓2016-02-25 20:48:55

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

仰望星空7

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 1237.7
帖子: 42
在線: 21.3小時(shí)
蟲號(hào): 4127268
注冊(cè): 2015-10-08
專業(yè): 工程熱物理相關(guān)交叉領(lǐng)域

function F=myfun(X)
T1=X(1);T2=X(2);T3=X(3);T4=X(4);Tf=X(5);q=X(6);
c1=0.003;
c2=0.003;
h=0.002;
w1=0.002;
w2=0.002;
v=5;
T0=3287;
%燃?xì)鈱?duì)流
syms sigma At A Pr hc1
dt=0.013;d=0.03;Pc=3.5*10^6;c=2000;C0=9.18*10^5;Cp=1500;u=1.56*10^(-5);
m=0.022;k=1.26;Ma=3;
sigma=(T1*(1+(k-1)*Ma*Ma/2)/(2*T0)+0.5)^(-0.68)*(1+(k-1)*Ma*Ma/2)^(-0.12);
At=pi*dt*dt/4;A=pi*d*d/4;
Pr=u*Cp/m;
hc1=(0.026/dt^0.2)*(u^0.2*Cp/Pr^0.6)*(Pc/c)^0.8*(At/A)^0.9*sigma;
F(1)=hc1*(T0-T1)+C0-q;
%內(nèi)壁導(dǎo)熱
syms k2
m1=387;
k2=m1/c1;
F(2)=k2*(T1-T2)-q;
%冷卻劑對(duì)流
syms Cpl ml lf lt pf pt uf ut Ref Prf hc20 hc2 yita
Cpl=4.1*Tf+1892.2;
ml=-0.0001*Tf+0.1383;lf=10^(-8)*(-5.3*Tf+683.7);lt=10^(-8)*(-5.3*T2+683.7);
pf=-0.8102*Tf+831.4661;pt=-0.8102*T2+831.4661;uf=pf*lf;ut=pt*lt;
de=2*h*w1/(h+w1);Ref=v*w1/lf;Prf=uf*Cpl/ml;
hc20=0.023*Ref^0.8*Prf^0.34*(uf/ut)*ml/de;
m2=202.4;m=(2*hc20/(m2*w2))^0.5;yita=tanh(m*h)/(m*h);
hc2=hc20*yita;
F(3)=hc2*(T2-Tf)-q;
F(4)=T2+T3-2*Tf;
%外壁導(dǎo)熱
m3=387;
F(5)=(T3-T4)/(c2/m3)-q;
%空氣對(duì)流
syms av Gra hc3
g=9.8;Ta=293;av=1/((T4+Ta)/2);d=0.2;pa=1.06;Pra=0.7;ma=0.04;la=4*10^(-5);
Gra=g*av*(T4-Ta)*d^3/la^2;hc3=0.48*(Gra*Pra)^0.25*ma/d;
F(6)=hc3*(T4-Ta)-q;
end

贊一下

回復(fù)此樓

6樓2016-03-10 15:22:13

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

小雙子藍(lán)胖子

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 9.5
帖子: 1
在線: 1.8小時(shí)
蟲號(hào): 2635402
注冊(cè): 2013-09-04
專業(yè): 控制理論與方法

引用回帖:

4樓: Originally posted by dingd at 2015-05-13 09:51:11
如果你用的是老版本就別試了，無法進(jìn)行這種循環(huán)計(jì)算。...

你好有沒有新的版本下載地址呀求助謝謝需要用到循環(huán)

贊一下

回復(fù)此樓

7樓2017-04-13 23:00:23

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

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

[考研] 【考研調(diào)劑】化學(xué)專業(yè) 281分，一志愿四川大學(xué)，誠心求調(diào)劑 +6	吃吃吃才有意義 2026-03-19	6/300	2026-03-20 10:47 by 盡舜堯1
[考研] 081700化工學(xué)碩調(diào)劑 +3	【1】 2026-03-16	3/150	2026-03-19 23:40 by edmund7
[考研] 一志愿福大288有機(jī)化學(xué)，求調(diào)劑 +3	小木蟲200408204 2026-03-18	3/150	2026-03-19 13:31 by houyaoxu
[考研] 085600材料與化工求調(diào)劑 +6	緒幸與子 2026-03-17	6/300	2026-03-19 13:27 by houyaoxu
[考研] 一志愿武理材料305分求調(diào)劑 +5	想上岸的鯉魚 2026-03-18	6/300	2026-03-18 17:53 by 無際的草原
[考研] 收復(fù)試調(diào)劑生 +4	雨后秋荷 2026-03-18	4/200	2026-03-18 14:16 by elevennnne
[考研] 312求調(diào)劑 +8	陌宸希 2026-03-16	9/450	2026-03-18 12:39 by Linda Hu
[考研] 303求調(diào)劑 +4	睿08 2026-03-17	6/300	2026-03-18 11:01 by Iveryant
[考研] 0703化學(xué)調(diào)劑 +3	妮妮ninicgb 2026-03-17	3/150	2026-03-18 10:29 by macy2011
[考研] 268求調(diào)劑 +6	簡單點(diǎn)0 2026-03-17	6/300	2026-03-18 09:04 by 無際的草原
[考研] 334求調(diào)劑 +3	志存高遠(yuǎn)意在機(jī)?/a> 2026-03-16	3/150	2026-03-18 08:34 by lm4875102
[考研] 一志愿蘇州大學(xué)材料工程（085601）專碩有科研經(jīng)歷三項(xiàng)國獎(jiǎng)兩個(gè)實(shí)用型專利一項(xiàng)省級(jí)立項(xiàng) +6	大火山小火山 2026-03-16	8/400	2026-03-17 15:05 by 無懈可擊111
[考研] 一志愿南京大學(xué)，080500材料科學(xué)與工程，調(diào)劑 +4	Jy? 2026-03-16	4/200	2026-03-17 11:02 by gaoqiong
[考研] 考研調(diào)劑 +3	淇ya_~ 2026-03-17	5/250	2026-03-17 09:25 by Winj1e
[考研] 283求調(diào)劑 +3	聽風(fēng)就是雨； 2026-03-16	3/150	2026-03-17 07:41 by 熱情沙漠
[考研] 304求調(diào)劑 +3	曼殊2266 2026-03-14	3/150	2026-03-16 16:39 by houyaoxu
[考研] 327求調(diào)劑 +6	拾光任染 2026-03-15	11/550	2026-03-15 22:47 by 拾光任染
[考研] 求老師收留調(diào)劑 +4	jiang姜66 2026-03-14	5/250	2026-03-15 20:11 by Winj1e
[考研] 一志愿哈工大材料324分求調(diào)劑 +5	閆旭東 2026-03-14	5/250	2026-03-14 14:53 by 木瓜膏
[考研] 311求調(diào)劑 +3	冬十三 2026-03-13	3/150	2026-03-13 20:41 by JourneyLucky