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

登錄注冊

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

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

返回列表

【懸賞金幣】回答本帖問題，作者ZXR319將贈(zèng)送您 40 個(gè)金幣

ZXR319

新蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 74.6
帖子: 72
在線: 7.5小時(shí)
蟲號(hào): 33303510
注冊: 2023-03-13
專業(yè): 化學(xué)工程及工業(yè)化學(xué)

[求助] 請會(huì)1stopt 的大佬幫幫我已有2人參與

我在做 lh 機(jī)理動(dòng)力學(xué)，已經(jīng)有一部分?jǐn)?shù)據(jù)了，目前已經(jīng)用 matlab 處理了，擬合初值是隨便給的，我現(xiàn)在想用1stopt 軟件擬合一下，擬合的結(jié)果作為 matlab 的參數(shù)初值，以便更準(zhǔn)確。求大佬幫我擬合一下，最好有源代碼，金幣我可以盡量給，有多少都可以給，好商量。
數(shù)據(jù)如下：t       ca       cb
0       0.1679       0
15    0.1511 0.0164
30 0.1322 0.0376
45 0.115    0.0512
60       0.1037      0.06
90 0.1003    0.069
120    0.09486    0.07368
還需要提供什么可隨時(shí)聯(lián)系我，拜謝各位！
圖2是我用 matlab 擬合的思路，目標(biāo)函數(shù)就是計(jì)算值和實(shí)驗(yàn)測量值的差（比如殘差平方和）最小。

@月只藍(lán) 發(fā)自小木蟲Android客戶端

回復(fù)此樓

» 猜你喜歡

重慶交通大學(xué)光子學(xué)微結(jié)構(gòu)與器件課題組招收2026年碩士研究生信息已經(jīng)有2人回復(fù)
一志愿鄭大材料學(xué)碩298分，求調(diào)劑已經(jīng)有6人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費(fèi)? 已經(jīng)有97人回復(fù)
尋合作：應(yīng)力腐蝕多尺度模擬已經(jīng)有5人回復(fù)
考研交流已經(jīng)有0人回復(fù)
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有0人回復(fù)
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有0人回復(fù)
幫我的英語口語老師找學(xué)生已經(jīng)有2人回復(fù)

1樓 2024-09-09 11:24:26

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

獨(dú)孤神宇

版主 (知名作家)

應(yīng)助: 490 (碩士)
貴賓: 0.008
金幣: 31016.3
散金: 802
紅花: 122
沙發(fā): 1
帖子: 5600
在線: 856.5小時(shí)
蟲號(hào): 3522474
注冊: 2014-11-06
性別: GG
專業(yè): 機(jī)械動(dòng)力學(xué)
管轄: 計(jì)算模擬

模型缺少 r的數(shù)據(jù)。

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

贊一下

回復(fù)此樓

數(shù)值計(jì)算

2樓2024-09-09 13:31:42

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

ZXR319

新蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 74.6
帖子: 72
在線: 7.5小時(shí)
蟲號(hào): 33303510
注冊: 2023-03-13
專業(yè): 化學(xué)工程及工業(yè)化學(xué)

引用回帖:

2樓: Originally posted by 獨(dú)孤神宇 at 2024-09-09 13:31:42
模型缺少 r的數(shù)據(jù)。

老師您的意思是數(shù)據(jù)里每一個(gè)時(shí)間下應(yīng)該寫出來對應(yīng)的反應(yīng)速率嗎？我通過實(shí)驗(yàn)只能得到不同時(shí)間的反應(yīng)物濃度和生成物濃度，那在1stopt 中可以通過代碼得到不同時(shí)刻的反應(yīng)速率嗎？我在 Matlab 中是通過 ODE45函數(shù)求反應(yīng)速率的

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

贊一下

回復(fù)此樓

3樓2024-09-09 14:57:27

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

ZXR319

新蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 74.6
帖子: 72
在線: 7.5小時(shí)
蟲號(hào): 33303510
注冊: 2023-03-13
專業(yè): 化學(xué)工程及工業(yè)化學(xué)

Function r = (k * KA * cA * sqrt(KH * PH2 / Kbc)) / ((1 + KA * cA + KB * cB + sqrt(KH * PH2 / Kbc))^2);
這是模型方程的代碼形式

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

贊一下

回復(fù)此樓

4樓2024-09-09 15:28:45

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

hzlhm

至尊木蟲 (著名寫手)

應(yīng)助: 387 (碩士)
金幣: 17807
紅花: 53
帖子: 2879
在線: 606.3小時(shí)
蟲號(hào): 1517335
注冊: 2011-11-30
性別: GG
專業(yè): 常微分方程與動(dòng)力系統(tǒng)

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

感謝參與，應(yīng)助指數(shù) +1

你的問題還得提供反應(yīng)速率的微分方程，才能做擬合系數(shù)的工作。

贊一下

回復(fù)此樓

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

ZXR319

QQ：2120156492

5樓2024-09-09 22:28:55

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

ZXR319

新蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 74.6
帖子: 72
在線: 7.5小時(shí)
蟲號(hào): 33303510
注冊: 2023-03-13
專業(yè): 化學(xué)工程及工業(yè)化學(xué)

送紅花一朵

引用回帖:

5樓: Originally posted by hzlhm at 2024-09-09 22:28:55
你的問題還得提供反應(yīng)速率的微分方程，才能做擬合系數(shù)的工作。

謝謝老師我做的是加氫反應(yīng)，1molA 和3mol 氫氣生成1molB，A 和 B 摩爾比是1:1，所以是不是-dcA/dt=r；dcB/dt=r。不知道這樣可以計(jì)算了嗎？

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

贊一下

回復(fù)此樓

6樓2024-09-10 08:18:34

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

獨(dú)孤神宇

版主 (知名作家)

應(yīng)助: 490 (碩士)
貴賓: 0.008
金幣: 31016.3
散金: 802
紅花: 122
沙發(fā): 1
帖子: 5600
在線: 856.5小時(shí)
蟲號(hào): 3522474
注冊: 2014-11-06
性別: GG
專業(yè): 機(jī)械動(dòng)力學(xué)
管轄: 計(jì)算模擬

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

感謝參與，應(yīng)助指數(shù) +1

引用回帖:

6樓: Originally posted by ZXR319 at 2024-09-10 08:18:34
謝謝老師我做的是加氫反應(yīng)，1molA 和3mol 氫氣生成1molB，A 和 B 摩爾比是1:1，所以是不是-dcA/dt=r；dcB/dt=r。不知道這樣可以計(jì)算了嗎？
...

假設(shè) dcB/dt=r

模型公式: cb' = (k*ka*ca*sqrt(kh*3/kbc))/((1+ka*ca+kb*cb+sqrt(kh*3/kbc))^2)

殘差平方和(SSE): 0.000100907681851596
均方差(RMSE): 0.00410096902881087
相關(guān)系數(shù)(R): 0.978424525403226
相關(guān)系數(shù)之平方(R^2): 0.957314551910527
修正R平方(Adj. R^2): 0.928857586517545
確定系數(shù)(DC): 0.956042109592344
F統(tǒng)計(jì)(F-Statistic): 5.07803520426272

參數(shù) 最佳估算
---- --------
k    3947916360160.45
ka    4.19389648022836E15
kh    0.406855254346946
kbc 8.49921244458393
kb    1.28426131584115E16

====== 結(jié)果輸出 ======

文件: 數(shù)據(jù)文件-1
No t 目標(biāo) cb 計(jì)算 cb
1 15 0.0164 0.023012428915109
2 30 0.0376 0.0370171975730106
3 45 0.0512 0.0473206201899215
4 60 0.06 0.05546914388147
5 90 0.069 0.068244316554898
6 120 0.07368 0.0782291213083908

贊一下

回復(fù)此樓

數(shù)值計(jì)算

7樓2024-09-11 21:54:51

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

獨(dú)孤神宇

版主 (知名作家)

應(yīng)助: 490 (碩士)
貴賓: 0.008
金幣: 31016.3
散金: 802
紅花: 122
沙發(fā): 1
帖子: 5600
在線: 856.5小時(shí)
蟲號(hào): 3522474
注冊: 2014-11-06
性別: GG
專業(yè): 機(jī)械動(dòng)力學(xué)
管轄: 計(jì)算模擬

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

引用回帖:

假設(shè) dcA/dt=-r

模型公式: ca' = -(k*ka*ca*sqrt(kh*3/kbc))/((1+ka*ca+kb*cb+sqrt(kh*3/kbc))^2)

殘差平方和(SSE): 0.000192586899315775
均方差(RMSE): 0.005665493495948
相關(guān)系數(shù)(R): 0.959232883604991
相關(guān)系數(shù)之平方(R^2): 0.920127724989146
修正R平方(Adj. R^2): 0.866879541648577
確定系數(shù)(DC): 0.91769457654023
F統(tǒng)計(jì)(F-Statistic): 2.54023664039648

參數(shù) 最佳估算
---- --------
k    2837450.21546848
ka    2.67743713648235E-7
kh    67.0303631498266
kbc 0.100651239217611
kb    687.010099643109

====== 結(jié)果輸出 ======

文件: 數(shù)據(jù)文件-1
No t 目標(biāo) ca 計(jì)算 ca
1 15 0.1511 0.143509362937322
2 30 0.1322 0.129868469396974
3 45 0.115 0.120093597136612
4 60 0.1037 0.112222162940757
5 90 0.1003 0.0997387665702967
6 120 0.09486 0.0893241064493534

贊一下

回復(fù)此樓

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

ZXR319

數(shù)值計(jì)算

8樓2024-09-11 21:58:39

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

ZXR319

新蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 74.6
帖子: 72
在線: 7.5小時(shí)
蟲號(hào): 33303510
注冊: 2023-03-13
專業(yè): 化學(xué)工程及工業(yè)化學(xué)

送紅花一朵

引用回帖:

8樓: Originally posted by 獨(dú)孤神宇 at 2024-09-11 21:58:39
假設(shè) dcA/dt=-r

模型公式: ca' = -(k*ka*ca*sqrt(kh*3/kbc))/((1+ka*ca+kb*cb+sqrt(kh*3/kbc))^2)

殘差平方和(SSE): 0.000192586899315775
均方差(RMSE): 0.005665493495948
相關(guān)系數(shù)(R): 0.9592328836049 ...

謝謝老師，您辛苦了，太感謝您了

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

贊一下

回復(fù)此樓

9樓2024-09-13 09:34:45

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

ZXR319

新蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 74.6
帖子: 72
在線: 7.5小時(shí)
蟲號(hào): 33303510
注冊: 2023-03-13
專業(yè): 化學(xué)工程及工業(yè)化學(xué)

送紅花一朵

引用回帖:

老師，能不能麻煩您給個(gè)源代碼？我自己弄的有問題，擬合出來反應(yīng)速率和濃度都是0。Title "動(dòng)力學(xué)模型2";
Parameters k,KA,KB,KH,Kbc;
Variable t,CA,CB;
Constant PH2=3;
Function r = (k * K_A * c_A * sqrt(K_H * P_H2 / K_bc)) / ((1 + K_A * c_A + K_B * c_B + sqrt(K_H * P_H2 / K_bc))^2);
Data;
t  CA  CB
0  0.1679  0
30  0.1314  0.0367
60  0.0964  0.0719
90  0.0759  0.096
120  0.0575  0.1113
STARTPROGRAM
var i: integer;
begin
  for i = 1 to length(t)-1 do
r = -(c_A - c_A[i-1]) / (t - t[i-1]);
end
ENDPROGRAM

FIT
這是Ai 給出的代碼，我不知道怎么改

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

贊一下

回復(fù)此樓

10樓2024-09-15 09:34:32

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

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

返回列表

不應(yīng)助 確定回帖應(yīng)助 (注意：應(yīng)助才可能被獎(jiǎng)勵(lì)，但不允許灌水，必須填寫15個(gè)字符以上)

普通表情龍兔虎貓

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

[考研] 復(fù)試調(diào)劑 +3	呼呼？~+123456 2026-03-05	5/250	2026-03-05 08:19 by ms629
[考研] 316求調(diào)劑 +3	林小星發(fā)大財(cái) 2026-03-04	3/150	2026-03-05 07:49 by bxbo
[考研] 376求調(diào)劑 +3	王浩然sam 2026-03-04	3/150	2026-03-05 07:48 by bxbo
[考研] 求調(diào)劑院校 +6	云朵452 2026-03-02	12/600	2026-03-04 23:17 by 云朵452
[考研] 085600材料與化工（高分子）290分求調(diào)劑 +7	wengyujian 2026-03-04	7/350	2026-03-04 20:50 by ms629
[考研] 320材料與化工，求調(diào)劑 +6	鶴遨予卿 2026-03-04	8/400	2026-03-04 20:47 by wutongshun
[考研] 一志愿鄭州大學(xué)，學(xué)碩，物理化學(xué)， 333求調(diào)劑 +5	李魔女斗篷 2026-03-04	5/250	2026-03-04 18:33 by houyaoxu
[考研] 085701環(huán)境工程求調(diào)劑 +6	xiiiia 2026-03-04	6/300	2026-03-04 16:45 by yuningshan
[考研] 292求調(diào)劑 +9	yhk_819 2026-02-28	9/450	2026-03-04 16:06 by sslc1985
[考研] 歡迎采礦、地質(zhì)、巖土、計(jì)算機(jī)、人工智能等專業(yè)的同學(xué)報(bào)考 +9	pin8023 2026-02-28	12/600	2026-03-04 15:27 by rorr
[考研] 325求調(diào)劑 +5	學(xué)家科 2026-03-04	5/250	2026-03-04 14:04 by kakakapanpan
[考研] 347分電子信息專碩求調(diào)劑 +3	咕嘟咕嘟1 2026-03-03	3/150	2026-03-04 12:32 by 熱情沙漠
[考研] 0703化學(xué)調(diào)劑 +4	G212 2026-03-03	5/250	2026-03-04 09:34 by 每天只擺一小會(huì)
[考研] 085602化學(xué)工程350，調(diào)劑，有沒有211的 +5	利好利好. 2026-03-02	9/450	2026-03-03 17:06 by 利好利好.
[考研] 0805總分292，求調(diào)劑 +12	幻想之殤 2026-03-01	12/600	2026-03-03 15:58 by tgxtgxtgx9
[考研] 課題組接收材料類調(diào)劑研究生 +6	gaoxiaoniuma 2026-02-28	9/450	2026-03-03 09:27 by gaoxiaoniuma
[考研] 338求調(diào)劑 +5	18162027187 2026-03-02	6/300	2026-03-03 07:52 by njzyff
[考博] 博士自薦 +4	kkluvs 2026-02-28	5/250	2026-03-02 19:19 by 輕松不少隨
[考研] 一志愿華南理工大學(xué)材料與化工326分，求調(diào)劑 +3	wujinrui1 2026-02-28	3/150	2026-03-02 16:36 by chuocheng
[考研] 299求調(diào)劑 +3	Y墨明棋妙Y 2026-02-28	5/250	2026-03-01 21:01 by tangxiaotian

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

[求助] 請會(huì)1stopt 的大佬幫幫我 已有2人參與

» 猜你喜歡

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

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

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

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

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

[求助] 請會(huì)1stopt 的大佬幫幫我已有2人參與