版塊導航: 正在加載中...

登錄注冊

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

24小時熱門版塊排行榜

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

返回列表

當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖

杜?Ramos

鐵蟲 (小有名氣)

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

[交流] 怎樣設(shè)置求解精度，一直提醒我方程解接近于零應該怎么改呢已有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

回復此樓

» 猜你喜歡

材料求調(diào)劑已經(jīng)有11人回復
【求調(diào)劑】新能源材料本科，一志愿211，初試321 已經(jīng)有4人回復
336求調(diào)劑已經(jīng)有5人回復
材料求調(diào)劑已經(jīng)有8人回復
材料工程322分已經(jīng)有7人回復
材料考研調(diào)劑已經(jīng)有5人回復
生物學308分求調(diào)劑（一志愿華東師大）已經(jīng)有4人回復
求生物學調(diào)劑已經(jīng)有8人回復
08生物與醫(yī)藥專碩初試346找調(diào)劑已經(jīng)有6人回復
生物學327，求調(diào)劑已經(jīng)有6人回復

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

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

杜?Ramos

鐵蟲 (小有名氣)

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

引用回帖:

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

謝謝我試一下

回復此樓

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

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

查看全部 5 個回答

moguowei00

銅蟲 (小有名氣)

應助: 0 (幼兒園)
金幣: 324.6
紅花: 1
帖子: 197
在線: 25.1小時
蟲號: 1835287
注冊: 2012-05-26
性別: GG
專業(yè): 機械測試理論與技術(shù)

好難

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

回復此樓

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

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

杜?Ramos

鐵蟲 (小有名氣)

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

引用回帖:

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

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

贊一下

回復此樓

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

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

pdl9527

專家顧問 (小有名氣)

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

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

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

贊一下

回復此樓

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

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

查看全部 5 個回答

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

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

[考研] 336求調(diào)劑 +5	kiyy 2026-04-01	5/250	2026-04-02 07:14 by chixmc
[考研] 295材料工程專碩求調(diào)劑 +19	1428151015 2026-03-27	19/950	2026-04-01 22:34 by peike
[考研] 307分求調(diào)劑 +14	(o~o) 2026-03-31	15/750	2026-04-01 20:43 by longlotian
[考研] 一志愿鄭大材料工程290求調(diào)劑 +18	Youth_ 2026-03-30	18/900	2026-04-01 19:23 by 無際的草原
[考研] 309求調(diào)劑 +19	誰不是少年 2026-03-29	19/950	2026-04-01 15:47 by jp9609
[考研] 085600，321分求調(diào)劑 +13	大饞小子 2026-03-31	13/650	2026-04-01 12:35 by chemdavid
[考研] 309分085801求調(diào)劑 +7	學員Gtwj7W 2026-03-31	7/350	2026-04-01 02:36 by BruceLiu320
[考研] 375求調(diào)劑 +7	雨夏整夜 2026-03-29	7/350	2026-03-31 18:52 by xhai2011
[考研] 一志愿中海洋材料357 +4	麥恩莉. 2026-03-30	4/200	2026-03-31 14:35 by 記事本2026
[考研] 262求調(diào)劑 +7	ZZ..000 2026-03-30	8/400	2026-03-31 10:05 by cal0306
[考研] 一志愿華中師范化學332分求調(diào)劑 +3	Lyy930824@ 2026-03-29	3/150	2026-03-30 20:15 by DHUSHUAI
[考研] 085600 286分材料求調(diào)劑 +11	麻辣魷魚 2026-03-27	12/600	2026-03-30 19:33 by Wang200018
[考研] 抱歉 +3	田洪有 2026-03-30	3/150	2026-03-30 19:11 by 迷糊CCPs
[考研] 一志愿南開大學0710生物學359求調(diào)劑 +5	兔兔兔111223314 2026-03-29	7/350	2026-03-30 18:29 by 兔兔兔111223314
[考研] 292求調(diào)劑 +13	是妍子也是研子 2026-03-30	13/650	2026-03-30 18:01 by 小徐0109
[考研] 279求調(diào)劑 +4	蝶舞輕繞 2026-03-29	4/200	2026-03-29 09:45 by laoshidan
[考研] 305求調(diào)劑 +8	RuiFairyrui 2026-03-28	8/400	2026-03-29 08:22 by fmesaito
[考研] 304求調(diào)劑 +6	曼殊2266 2026-03-27	6/300	2026-03-28 14:10 by 唐沐兒
[考研] 復試調(diào)劑，一志愿南農(nóng)083200食品科學與工程 +5	XQTJZ 2026-03-26	5/250	2026-03-27 14:49 by 狂炫麥當當
[考研] 298調(diào)劑 +3	jiyingjie123 2026-03-27	3/150	2026-03-27 11:57 by wxiongid

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

24小時熱門版塊排行榜

杜?Ramos

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

» 猜你喜歡

杜?Ramos

moguowei00

杜?Ramos

pdl9527

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

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