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

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

24小時(shí)熱門版塊排行榜

北京石油化工學(xué)院2026年研究生招生接收調(diào)劑公告

返回列表

杜?Ramos

鐵蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 6.1
散金: 387
紅花: 2
帖子: 53
在線: 22.2小時(shí)
蟲號(hào): 4407039
注冊(cè): 2016-02-16
專業(yè): 聲學(xué)

[交流] 怎樣設(shè)置求解精度，一直提醒我方程解接近于零應(yīng)該怎么改呢已有2人參與

function equation()
global Rg M T Veff P
Rg=8.31;
M=5.8*10^-3;
T=273.15;
Veff=1;
P=10000;
n0=[53.156*10^-3 0 0 5.98*10^-9 0 0 0 0];
n=fsolve(@myfun,n0);
end

function F=myfun(n)
global Rg M T Veff P
n1=n(1);
n2=n(2);
n3=n(3);
n4=n(4);
n5=n(5);
n6=n(6);
n7=n(7);
n8=n(8);
kf1=1.912*10^(23)*(T^(-1.83))*exp(-1.185*10^5/(8.31*T));
kr1=2.2*10^(22)*(T^(-2))*exp(0/(8.31*T));
kf2=4.515*10^(17)*(T^(-0.64))*exp(-1.189*10^5/(8.31*T));
kr2=6.165*10^(15)*(T^(-0.5))*exp(0/(8.31*T));
kf3=2.97*10^(6)*(T^(2.02))*exp(-1.34*10^4/(8.31*T));
kr3=1.465*10^(5)*(T^(2.11))*exp(2.904*10^3/(8.31*T));
kf4=1.475*10^(12)*(T^(0.6))*exp(0/(8.31*T));
kr4=3.09*10^(12)*(T^(0.53))*exp(-4.887*10^4/(8.31*T));
kf5=2.667*10^(4)*(T^(2.65))*exp(-4.88*10^3/(8.31*T));
kr5=3.82*10^(12)*(T^(0))*exp(-7.948*10^3/(8.31*T));
kf6=2.298*10^(9)*(T^(1.4))*exp(-1.832*10^4/(8.31*T));
kr6=2.16*10^(8)*(T^(1.52))*exp(-3.45*10^3/(8.31*T));
kf7=3.25*10^(13)*(T^(0))*exp(0/(8.31*T));
kr7=3.252*10^(12)*(T^(0.33))*exp(-5.328*10^4/(8.31*T));
kf8=7.079*10^(13)*(T^(0))*exp(-2.95*10^2/(8.31*T));
kr8=2.027*10^(10)*(T^(-0.72))*exp(-3.684*10^4/(8.31*T));
kf9=1.838*10^(10)*(T^(0.59))*exp(-3.089*10^4/(8.31*T));
kr9=1*10^(12)*(T^(0))*exp(0/(8.31*T));
kf10=8.66*10^(3)*(T^(2.68))*exp(-1.856*10^4/(8.31*T));
kr10=9.550*10^(6)*(T^(2))*exp(-3.97*10^3/(8.31*T));
kf11=1.66*10^(13)*(T^(0))*exp(-8.23*10^2/(8.31*T));
kr11=3.164*10^(12)*(T^(0.35))*exp(-5.551*10^4/(8.31*T));
kf12=1*10^(14)*(T^(-0.37))*exp(0/(8.31*T));
kr12=2.951*10^(14)*(T^(0))*exp(-4.843*10^4/(8.31*T));
kf13=4.2*10^(14)*(T^(0))*exp(-1.198*10^4/(8.31*T));
kr13=4.624*10^(16)*(T^(-0.35))*exp(-5.067*10^4/(8.31*T));
nt=n1+n2+n3+n4+n5+n6+n7+n8;
eq1=-(nt*Rg*T)/P*((kf1*(n1/Veff)*(M/Veff)-kr1*(n2/Veff)*(n3/Veff)*(M/Veff))+(kf3*(n1/Veff)*(n5/Veff)-kr3*(n3/Veff)*(n3/Veff))+(kf6*(n1/Veff)*(n2/Veff)-kr6*(n7/Veff)*(n3/Veff))+(kf9*(n1/Veff)*(n6/Veff)-kr9*(n8/Veff)*(n3/Veff)));
eq2=-(nt*Rg*T)/P*((kr1*(n2/Veff)*(n3/Veff)*(M/Veff)-kf1*(n1/Veff)*(M/Veff))+(kf4*(n2/Veff)*(n4/Veff)*(M/Veff)-kr4*(n6/Veff)*(M/Veff))+(kf5*(n2/Veff)*(n3/Veff)-kr5*(n5/Veff)*(n7/Veff))+(kf6*(n2/Veff)*(n1/Veff)-kr6*(n3/Veff)*(n7/Veff))+(kf8*(n2/Veff)*(n6/Veff)-kr8*(n3/Veff)*(n3/Veff))+(kf11*(n6/Veff)*(n2/Veff)-kr11*(n7/Veff)*(n4/Veff)));
eq3=-(nt*Rg*T)/P*((kr1*(n2/Veff)*(n3/Veff)*(M/Veff)-kf1*(n1/Veff)*(M/Veff))+(2*kr3*(n3/Veff)*(n3/Veff)-kr3*(n1/Veff)*(n5/Veff))+(kf5*(n2/Veff)*(n3/Veff)-kr5*(n5/Veff)*(n7/Veff))+(kr6*(n3/Veff)*(n7/Veff)-kf6*(n2/Veff)*(n1/Veff))+(kr7*(n3/Veff)*(n4/Veff)-kf7*(n6/Veff)*(n5/Veff))+(2*kr8*(n3/Veff)*(n3/Veff)-kf8*(n6/Veff)*(n2/Veff))+(kr9*(n8/Veff)*(n3/Veff)-kf9*(n1/Veff)*(n6/Veff))+(kf10*(n3/Veff)*(n6/Veff)-kr10*(n8/Veff)*(n5/Veff))+(2*kf12*(n3/Veff)*(n3/Veff)*(M/Veff)-kr11*(n8/Veff)*(M/Veff)));
eq4=-(nt*Rg*T)/P*((kf2*(n4/Veff)*(M/Veff)-kr2*(n5/Veff)*(n5/Veff)*(M/Veff))+(kf4*(n2/Veff)*(n4/Veff)*(M/Veff)-kr4*(n6/Veff)*(M/Veff))+(kr7*(n3/Veff)*(n4/Veff)-kf7*(n6/Veff)*(n5/Veff))+(kr11*(n7/Veff)*(n4/Veff)-kf11*(n6/Veff)*(n2/Veff))+(kr13*(n8/Veff)*(n4/Veff)-kf13*(n6/Veff)*(n6/Veff)));
eq5=-(nt*Rg*T)/P*((2*kr2*(n5/Veff)*(n5/Veff)*(M/Veff)-kf2*(n4/Veff)*(M/Veff))+(kf3*(n1/Veff)*(n5/Veff)-kr3*(n3/Veff)*(n3/Veff))+(kr5*(n5/Veff)*(n7/Veff)-kf5*(n2/Veff)*(n3/Veff))+(kf7*(n6/Veff)*(n5/Veff)-kr7*(n3/Veff)*(n4/Veff))+(kr10*(n8/Veff)*(n5/Veff)-kf10*(n3/Veff)*(n6/Veff)));
eq6=-(nt*Rg*T)/P*((kr4*(n6/Veff)*(M/Veff)-kf4*(n2/Veff)*(n4/Veff)*(M/Veff))+(kf7*(n6/Veff)*(n5/Veff)-kr7*(n3/Veff)*(n4/Veff))+(kf8*(n6/Veff)*(n2/Veff)-kr8*(n3/Veff)*(n3/Veff))+(kf9*(n1/Veff)*(n6/Veff)-kr9*(n3/Veff)*(n8/Veff))+(kf10*(n3/Veff)*(n6/Veff)-kr10*(n8/Veff)*(n5/Veff))+(kf11*(n6/Veff)*(n2/Veff)-kr11*(n7/Veff)*(n4/Veff))+(2*kf13*(n6/Veff)*(n6/Veff)-kr13*(n8/Veff)*(n4/Veff)));
eq7=-(nt*Rg*T)/P*((kr5*(n5/Veff)*(n7/Veff)-kf5*(n2/Veff)*(n3/Veff))+(kr11*(n4/Veff)*(n7/Veff)-kf11*(n6/Veff)*(n2/Veff)));
eq8=-(nt*Rg*T)/P*((kr9*(n3/Veff)*(n8/Veff)-kf9*(n1/Veff)*(n6/Veff))+(kr10*(n8/Veff)*(n5/Veff)-kf10*(n3/Veff)*(n6/Veff))+(kr12*(n8/Veff)*(M/Veff)-kf12*(n3/Veff)*(n3/Veff)*(M/Veff))+(kr13*(n8/Veff)*(n4/Veff)-kf13*(n6/Veff)*(n6/Veff)));
F=[eq1;eq2;eq3;eq4;eq5;eq6;eq7;eq8];
end

回復(fù)此樓

» 猜你喜歡

生物學(xué)學(xué)碩，一志愿湖南大學(xué)，初試成績(jī)338 已經(jīng)有6人回復(fù)
321求調(diào)劑已經(jīng)有7人回復(fù)
343求調(diào)劑已經(jīng)有4人回復(fù)
305求調(diào)劑已經(jīng)有3人回復(fù)
本科新能源科學(xué)與工程，一志愿華理能動(dòng)285求調(diào)劑已經(jīng)有6人回復(fù)
332求調(diào)劑已經(jīng)有3人回復(fù)
學(xué)碩274求調(diào)劑已經(jīng)有8人回復(fù)
壓國(guó)家一區(qū)線，求導(dǎo)師收留，有恩必謝！已經(jīng)有7人回復(fù)
081200-11408-276學(xué)碩求調(diào)劑已經(jīng)有4人回復(fù)
22408 359分調(diào)劑已經(jīng)有4人回復(fù)

1樓 2017-10-23 14:21:31

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

moguowei00

銅蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 324.6
紅花: 1
帖子: 197
在線: 25.1小時(shí)
蟲號(hào): 1835287
注冊(cè): 2012-05-26
性別: GG
專業(yè): 機(jī)械測(cè)試?yán)碚撆c技術(shù)

好難

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

回復(fù)此樓

2樓2017-10-23 14:34:08

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

杜?Ramos

鐵蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 6.1
散金: 387
紅花: 2
帖子: 53
在線: 22.2小時(shí)
蟲號(hào): 4407039
注冊(cè): 2016-02-16
專業(yè): 聲學(xué)

引用回帖:

2樓: Originally posted by moguowei00 at 2017-10-23 14:34:08
好難

函數(shù)表達(dá)式有點(diǎn)長(zhǎng)，就是求解屬于很簡(jiǎn)單。我不知道怎樣設(shè)置求解的取值精度

贊一下

回復(fù)此樓

3樓2017-10-23 14:47:34

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

pdl9527

專家顧問 (小有名氣)

專家經(jīng)驗(yàn): +8
應(yīng)助: 100 (初中生)
金幣: 2110.3
紅花: 23
帖子: 282
在線: 112小時(shí)
蟲號(hào): 1227333
注冊(cè): 2011-03-09
性別: GG
專業(yè): 分離過程
管轄: 計(jì)算模擬

★
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖

你的方程的形式是 ax-b=0，所以這里提示你最后的結(jié)果能夠保證這個(gè)差值接近0了，如果你把求解的那一行改成 [n,fval,exitflag,output]=fsolve(@myfun,n0)，你會(huì)發(fā)現(xiàn)exitflag=3，一般來說exitflag>0的情況下求得的解都是可靠的。exitflag具體代表的意思你可以在命令框輸入doc fsolve自行了解。

贊一下

回復(fù)此樓

4樓2017-10-23 16:52:40

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

杜?Ramos

鐵蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 6.1
散金: 387
紅花: 2
帖子: 53
在線: 22.2小時(shí)
蟲號(hào): 4407039
注冊(cè): 2016-02-16
專業(yè): 聲學(xué)

引用回帖:

4樓: Originally posted by pdl9527 at 2017-10-23 16:52:40
你的方程的形式是 ax-b=0，所以這里提示你最后的結(jié)果能夠保證這個(gè)差值接近0了，如果你把求解的那一行改成 =fsolve(@myfun,n0)，你會(huì)發(fā)現(xiàn)exitflag=3，一般來說exitflag>0的情況下求得的解都是可靠的。exitf ...

謝謝我試一下

回復(fù)此樓

5樓2017-10-23 17:11:55

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

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主杜?Ramos 的主題更新

返回列表

普通表情龍兔虎貓高級(jí)回復(fù) (可上傳附件)

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

[考研] 調(diào)劑310 +6	溫柔的晚安 2026-03-25	7/350	2026-03-28 16:04 by qq2000
[考研] 食品工程專碩一志愿中海洋309求調(diào)劑 +3	小張zxy張 2026-03-26	6/300	2026-03-28 15:55 by 小張zxy張
[考研] 化學(xué)0703 調(diào)劑 306分一志愿211 +4	26要上岸 2026-03-28	4/200	2026-03-28 15:30 by 1018329917
[考研] 347求調(diào)劑 +3	山頂見α 2026-03-25	3/150	2026-03-28 14:13 by 唐沐兒
[考研] 308求調(diào)劑 +7	墨墨漠 2026-03-27	7/350	2026-03-28 07:43 by 熱情沙漠
[有機(jī)交流] 高溫高壓反應(yīng)求助 10+4	chibby 2026-03-25	4/200	2026-03-27 21:08 by BT20230424
[考研] 348求調(diào)劑 +4	小懶蟲不懶了 2026-03-27	5/250	2026-03-27 12:47 by 果果媽咪
[考研] 315調(diào)劑 +4	0860求調(diào)劑 2026-03-26	5/250	2026-03-27 11:23 by wangjy2002
[考研] 調(diào)劑 +3	李嘉圖·S·路 2026-03-27	3/150	2026-03-27 11:19 by wangjy2002
[考研] 352求調(diào)劑 +4	大米飯！ 2026-03-22	4/200	2026-03-26 16:40 by 不吃魚的貓
[考研] 332求調(diào)劑 +6	032500 2026-03-25	6/300	2026-03-25 22:45 by 418490947
[考研] 材料與化工328分調(diào)劑 +6	。，。，。，。i 2026-03-23	6/300	2026-03-25 22:30 by 418490947
[考研] 材料專碩331求調(diào)劑 +4	鮮當(dāng)牛 2026-03-24	4/200	2026-03-24 15:58 by JourneyLucky
[考研] 求調(diào)劑一志愿武漢理工大學(xué)材料工程（085601） +5	WW.' 2026-03-23	7/350	2026-03-24 14:50 by sprinining
[考研] 一志愿山東大學(xué)藥學(xué)學(xué)碩求調(diào)劑 +3	開開心心沒煩惱 2026-03-23	4/200	2026-03-24 00:06 by 開開心心沒煩惱
[考研] 材料/農(nóng)業(yè)專業(yè)，07/08開頭均可，過線就行 +3	呵唔哦豁 2026-03-23	4/200	2026-03-23 22:30 by 汪��？！
[考研] 361求調(diào)劑 +3	Glack 2026-03-22	3/150	2026-03-23 22:03 by fuyu_
[考研] 一志愿070300浙大化學(xué)358分，求調(diào)劑！ +4	酥酥魚.. 2026-03-21	4/200	2026-03-23 08:12 by Iveryant
[考研] 293求調(diào)劑 +3	濤濤Wjt 2026-03-22	5/250	2026-03-22 22:21 by jiangpengfei
[考研] 308求調(diào)劑 +3	墨墨漠 2026-03-21	3/150	2026-03-22 16:54 by i_cooler

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产 高清 中文字幕,99re热久久亚洲综合精品成人,熟妇 一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

24小時(shí)熱門版塊排行榜

杜?Ramos

[交流] 怎樣設(shè)置求解精度，一直提醒我方程解接近于零 應(yīng)該怎么改呢 已有2人參與

» 猜你喜歡

moguowei00

杜?Ramos

pdl9527

杜?Ramos

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产高清中文字幕,99re热久久亚洲综合精品成人,熟妇一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

[交流] 怎樣設(shè)置求解精度，一直提醒我方程解接近于零應(yīng)該怎么改呢已有2人參與