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

登錄注冊

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

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

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

返回列表

當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖

杜?Ramos

鐵蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 6.1
散金: 387
紅花: 2
帖子: 53
在線: 22.2小時(shí)
蟲號: 4407039
注冊: 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ù)此樓

» 猜你喜歡

352分-085602-一志愿985 已經(jīng)有3人回復(fù)
291求調(diào)劑已經(jīng)有4人回復(fù)
330分求調(diào)劑已經(jīng)有4人回復(fù)
309求調(diào)劑已經(jīng)有4人回復(fù)
0856求調(diào)劑已經(jīng)有13人回復(fù)
0703化學(xué) 已經(jīng)有11人回復(fù)
329求調(diào)劑已經(jīng)有7人回復(fù)
321求調(diào)劑已經(jīng)有8人回復(fù)
0856材料化工調(diào)劑總分330 已經(jīng)有12人回復(fù)
070300化學(xué)354求調(diào)劑已經(jīng)有4人回復(fù)

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

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

pdl9527

專家顧問 (小有名氣)

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

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

你的方程的形式是 ax-b=0，所以這里提示你最后的結(jié)果能夠保證這個(gè)差值接近0了，如果你把求解的那一行改成 [n,fval,exitflag,output]=fsolve(@myfun,n0)，你會發(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的回帖

查看全部 5 個(gè)回答

moguowei00

銅蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 324.6
紅花: 1
帖子: 197
在線: 25.1小時(shí)
蟲號: 1835287
注冊: 2012-05-26
性別: GG
專業(yè): 機(jī)械測試?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í)
蟲號: 4407039
注冊: 2016-02-16
專業(yè): 聲學(xué)

引用回帖:

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

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

贊一下

回復(fù)此樓

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

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

杜?Ramos

鐵蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 6.1
散金: 387
紅花: 2
帖子: 53
在線: 22.2小時(shí)
蟲號: 4407039
注冊: 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)，你會發(fā)現(xiàn)exitflag=3，一般來說exitflag>0的情況下求得的解都是可靠的。exitf ...

謝謝我試一下

回復(fù)此樓

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

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

查看全部 5 個(gè)回答

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

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

[考研] 309求調(diào)劑 +4	誰不是少年 2026-03-29	4/200	2026-03-29 07:07 by wxiongid
[考研] 0703化學(xué) +11	妮妮ninicgb 2026-03-27	11/550	2026-03-29 06:45 by 544594351
[考研] 調(diào)劑310 +7	溫柔的晚安 2026-03-25	8/400	2026-03-29 01:09 by 我是小康
[考研] 298求調(diào)劑 +3	種圣賜 2026-03-28	3/150	2026-03-28 23:58 by wxiongid
[考研] 322求調(diào)劑 +7	宋明欣 2026-03-27	7/350	2026-03-28 21:27 by sanrepian
[考研] 壓國家一區(qū)線，求導(dǎo)師收留，有恩必謝！ +7	迷人的哈哈 2026-03-28	7/350	2026-03-28 16:47 by 催化大白
[考研] 081200-314 +3	LILIQQ 2026-03-27	4/200	2026-03-28 09:41 by 保護(hù)地球你我做?/a>
[考研] 286求調(diào)劑 +4	丟掉懶惰 2026-03-27	7/350	2026-03-28 08:07 by baoball
[考研] 求調(diào)劑推薦材料 304 +15	荷包蛋hyj 2026-03-26	15/750	2026-03-28 04:13 by fmesaito
[考研] 330一志愿中國海洋大學(xué) 化學(xué)工程 085602 有讀博意愿求調(diào)劑 +3	wywy.. 2026-03-27	4/200	2026-03-28 03:32 by fmesaito
[考研] 一志愿上海理工能源動(dòng)力（085800）310分求調(diào)劑 +3	zhangmingc 2026-03-27	4/200	2026-03-27 19:01 by 給你你注意休息
[考研] 08開頭275求調(diào)劑 +4	拉誰不重要 2026-03-26	4/200	2026-03-27 14:12 by Delta2012
[考研] 322求調(diào)劑 +4	我真的很想學(xué)習(xí) 2026-03-23	4/200	2026-03-27 13:51 by 楊楊楊紫
[考研] 317求調(diào)劑 +7	蛋黃咸肉粽 2026-03-26	7/350	2026-03-27 02:29 by fmesaito
[考研] 生物學(xué) 296 求調(diào)劑 +4	朵朵- 2026-03-26	6/300	2026-03-26 19:01 by 不吃魚的貓
[考研] 07化學(xué)303求調(diào)劑 +5	睿08 2026-03-25	5/250	2026-03-25 22:46 by 418490947
[考研] 0854AI CV方向招收調(diào)劑 +4	章小魚567 2026-03-23	4/200	2026-03-25 17:04 by CoderLoser
[考研] 【2026考研調(diào)劑】制藥工程 284分求相關(guān)專業(yè)調(diào)劑名額 +4	袁奐奐 2026-03-25	8/400	2026-03-25 14:32 by lbsjt
[考研] 生物學(xué)學(xué)碩求調(diào)劑 +7	小羊睡著了? 2026-03-23	10/500	2026-03-25 02:24 by 清風(fēng)拂揚(yáng)。 m
[考研] 一志愿國科過程所081700，274求調(diào)劑 +3	三水研0水立方 2026-03-23	3/150	2026-03-23 23:11 by MajorWen

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

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

杜?Ramos

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

» 猜你喜歡

pdl9527

moguowei00

杜?Ramos

杜?Ramos

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

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