 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 2497 | 回復(fù): 12 | ||
| 【懸賞金幣】回答本帖問(wèn)題，作者ZXR319將贈(zèng)送您 50 個(gè)金幣 | ||
ZXR319新蟲(chóng) (小有名氣)
|
[求助]
Matlab 擬合動(dòng)力學(xué)參數(shù)時(shí)設(shè)置的初值不同，輸出的結(jié)果不同怎么辦？ 已有3人參與
|
|
|
我在做 LH 機(jī)理動(dòng)力學(xué)參數(shù)估計(jì)，發(fā)現(xiàn)設(shè)置的參數(shù)初值不同輸出的結(jié)果就不同，我通過(guò)調(diào)整參數(shù)初值能使 lnk 和1/T 呈線性相關(guān)，但 lnK 和1/T 就不成線性了怎么辦？我該怎么調(diào)整參數(shù)初值？以下是我的代碼: % 主腳本 % 步驟1: 讀取Excel文件中的數(shù)據(jù) filename = 'D:\桌面文件\新60.xlsx'; data = readtable(filename); time = data{:, 't'}; % 時(shí)間數(shù)據(jù) CA_exp_333K = data{:, 'CA'}; % 333K下反應(yīng)物A的濃度數(shù)據(jù) CB_exp_333K = data{:, 'CB'}; % 333K下生成物B的濃度數(shù)據(jù) % 步驟2: 初始參數(shù)估計(jì) params_initial = [0.5,0.5,0.5,0.5,0.5]; % 需要根據(jù)實(shí)際情況進(jìn)行調(diào)整 lb=[0 0 0 0 0]; ub=[100 100 100 100 100]; % 使用非線性最小二乘法擬合參數(shù) options = optimoptions('lsqnonlin', 'Algorithm', 'levenberg-marquardt', 'Display', 'iter'); params_est = lsqnonlin(@(params) ObjFunc(params, time, CA_exp_333K, CB_exp_333K), params_initial, [], [], options); % 使用擬合參數(shù)重新計(jì)算濃度 options = odeset('RelTol', 1e-6, 'AbsTol', 1e-9); [t, C] = ode45(@(t, C) KineticModel2(t, C, params_est), time, [CA_exp_333K(1); CB_exp_333K(1)], options); CA_pred = C(:, 1); CB_pred = C(:, 2); % 步驟3: 計(jì)算總方差和 (SST) 并檢查 residuals = ObjFunc(params_est, time, CA_exp_333K, CB_exp_333K); SST = sum(residuals.^2); % 計(jì)算總方差和 % 步驟4: 根據(jù)SST值繪制圖形或重新設(shè)定初值 if SST < 1 % 如果SST小于1，繪制圖形并輸出參數(shù) fprintf('擬合的參數(shù)值為:\n'); fprintf('k: %f\n', params_est(1)); fprintf('K_A: %f\n', params_est(2)); fprintf('K_B: %f\n', params_est(3)); fprintf('K_H: %f\n', params_est(4)); fprintf('K_bc: %f\n', params_est(5)); fprintf('總方差和 (SST): %f\n', SST); % 輸出總方差和 % 繪制實(shí)驗(yàn)濃度和計(jì)算濃度相關(guān)圖 figure; plot(time, CA_exp_333K, 'b-', time, CA_pred, 'r--'); hold on; plot(time, CB_exp_333K, 'g-', time, CB_pred, 'k--'); legend('實(shí)驗(yàn)CA', '擬合CA', '實(shí)驗(yàn)CB', '擬合CB'); xlabel('時(shí)間 (min)'); ylabel('濃度 (mol/L)'); title('實(shí)驗(yàn)濃度與計(jì)算濃度對(duì)比'); hold off; else % 如果SST大于或等于1，提示需要重新設(shè)定初值 fprintf('SST值過(guò)大，需要重新設(shè)定參數(shù)初值。當(dāng)前SST: %f\n', SST); % 這里可以添加代碼來(lái)調(diào)整初始參數(shù)或采取其他措施 end function dCdt = KineticModel2(t, C, params) k = params(1); K_A = params(2); K_B = params(3); K_H = params(4); K_bc = params(5); P_H2 = 3e3; % 氫氣壓力，單位為KPa c_A = C(1); c_B = C(2); numerator = k * K_A * c_A * sqrt((K_H * P_H2) / K_bc); denominator = (1 + K_A * c_A + K_B * c_B+ sqrt((K_H * P_H2) / K_bc))^2; r = numerator / denominator; dC_A_dt = -r; dC_B_dt = r; dCdt = [dC_A_dt; dC_B_dt]; end function res = ObjFunc(params, time, CA_exp, CB_exp) options = odeset('RelTol', 1e-6, 'AbsTol', 1e-9); [t, C] = ode45(@(t, C) KineticModel2(t, C, params), time, [CA_exp(1); CB_exp(1)], options); CA_pred = C(:, 1); CB_pred = C(:, 2); % 計(jì)算總方差和 res = (CA_exp - CA_pred).^2 + (CB_exp - CB_pred).^2; end 以下圖片中是我的數(shù)據(jù)，跪求大佬指導(dǎo)? 參數(shù)初值通過(guò)查文獻(xiàn)給出的方法不現(xiàn)實(shí)，因?yàn)槲易鲞@個(gè)沒(méi)有人做過(guò)，類(lèi)似的文獻(xiàn)初值差距也很大，我覺(jué)得沒(méi)有什么參考價(jià)值???而且我不太明白，比如我有五個(gè)參數(shù) k、KA、KB、KH、Kbc，那這五個(gè)參數(shù)初值應(yīng)該一樣嗎？ 發(fā)自小木蟲(chóng)Android客戶(hù)端 |
回復(fù)此樓» 猜你喜歡
考研交流
已經(jīng)有0人回復(fù)
-
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有0人回復(fù)
-
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)?
已經(jīng)有174人回復(fù)
-
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有0人回復(fù)
-
幫我的英語(yǔ)口語(yǔ)老師找學(xué)生
已經(jīng)有3人回復(fù)
-
計(jì)算三階力常數(shù)用thirdorder擴(kuò)胞時(shí)參數(shù)如何選擇？臨近原子數(shù)怎么取值
已經(jīng)有2人回復(fù)
-
2026年循環(huán)經(jīng)濟(jì)功能材料國(guó)際會(huì)議（ICFMCE 2026）
已經(jīng)有0人回復(fù)
-
2026年第五屆電氣、電子與信息工程國(guó)際會(huì)議（ISEEIE 2026）
已經(jīng)有0人回復(fù)
xrs333
木蟲(chóng) (正式寫(xiě)手)
- 應(yīng)助: 151 (高中生)
- 金幣: 5117.2
- 紅花: 23
- 帖子: 390
- 在線: 239.5小時(shí)
- 蟲(chóng)號(hào): 1867708
- 注冊(cè): 2012-06-23
- 性別: GG
- 專(zhuān)業(yè): 內(nèi)流流體力學(xué)
2樓2024-10-07 09:10:05
johnliu1983
至尊木蟲(chóng) (著名寫(xiě)手)
- 應(yīng)助: 43 (小學(xué)生)
- 金幣: 15165.3
- 紅花: 8
- 帖子: 2446
- 在線: 149.9小時(shí)
- 蟲(chóng)號(hào): 370817
- 注冊(cè): 2007-05-14
- 性別: GG
- 專(zhuān)業(yè): 物理學(xué)II
3樓2024-10-08 18:39:53
ZXR319
新蟲(chóng) (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 74.6
- 帖子: 72
- 在線: 7.5小時(shí)
- 蟲(chóng)號(hào): 33303510
- 注冊(cè): 2023-03-13
- 專(zhuān)業(yè): 化學(xué)工程及工業(yè)化學(xué)
|
謝謝您的回復(fù)！關(guān)鍵問(wèn)題是我不知道合適的參數(shù)初值應(yīng)該是多少?我做的東西別人沒(méi)做過(guò)，類(lèi)似的文獻(xiàn)中的參數(shù)帶進(jìn)去擬合效果又不好，所以我應(yīng)該考慮全局最優(yōu)化算法嗎？ 發(fā)自小木蟲(chóng)Android客戶(hù)端 |
4樓2024-10-09 11:47:21
ZXR319
新蟲(chóng) (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 74.6
- 帖子: 72
- 在線: 7.5小時(shí)
- 蟲(chóng)號(hào): 33303510
- 注冊(cè): 2023-03-13
- 專(zhuān)業(yè): 化學(xué)工程及工業(yè)化學(xué)
|
謝謝您的回復(fù)！你說(shuō)的意思我明白，可我參數(shù)初值代0.5或1的時(shí)候得到的結(jié)果都不一樣，這個(gè)參數(shù)范圍實(shí)在是不好確定，而且就算確定了范圍，給出一個(gè)初值，擬合效果也不好，不知道為什么，真的愁死了 發(fā)自小木蟲(chóng)Android客戶(hù)端 |
5樓2024-10-09 11:50:23
xrs333
木蟲(chóng) (正式寫(xiě)手)
- 應(yīng)助: 151 (高中生)
- 金幣: 5117.2
- 紅花: 23
- 帖子: 390
- 在線: 239.5小時(shí)
- 蟲(chóng)號(hào): 1867708
- 注冊(cè): 2012-06-23
- 性別: GG
- 專(zhuān)業(yè): 內(nèi)流流體力學(xué)
6樓2024-10-09 18:12:01
johnliu1983
至尊木蟲(chóng) (著名寫(xiě)手)
- 應(yīng)助: 43 (小學(xué)生)
- 金幣: 15165.3
- 紅花: 8
- 帖子: 2446
- 在線: 149.9小時(shí)
- 蟲(chóng)號(hào): 370817
- 注冊(cè): 2007-05-14
- 性別: GG
- 專(zhuān)業(yè): 物理學(xué)II
ZXR3197樓2024-10-09 20:17:28
hzlhm
至尊木蟲(chóng) (著名寫(xiě)手)
- 應(yīng)助: 387 (碩士)
- 金幣: 17807
- 紅花: 53
- 帖子: 2879
- 在線: 606.3小時(shí)
- 蟲(chóng)號(hào): 1517335
- 注冊(cè): 2011-11-30
- 性別: GG
- 專(zhuān)業(yè): 常微分方程與動(dòng)力系統(tǒng)
【答案】應(yīng)助回帖
|
1、先用rand函數(shù)作為初值 2、然后修正初值 3、最后得到結(jié)果 鏈接: https://pan.baidu.com/s/1X3gWCp4EEbFsWrwCKLyeEA?pwd=wur3 提取碼: wur3 復(fù)制這段內(nèi)容后打開(kāi)百度網(wǎng)盤(pán)手機(jī)App，操作更方便哦 |
» 本帖已獲得的紅花（最新10朵）
8樓2024-10-09 21:36:21
hzlhm
至尊木蟲(chóng) (著名寫(xiě)手)
- 應(yīng)助: 387 (碩士)
- 金幣: 17807
- 紅花: 53
- 帖子: 2879
- 在線: 606.3小時(shí)
- 蟲(chóng)號(hào): 1517335
- 注冊(cè): 2011-11-30
- 性別: GG
- 專(zhuān)業(yè): 常微分方程與動(dòng)力系統(tǒng)
9樓2024-10-09 21:37:30
ZXR319
新蟲(chóng) (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 74.6
- 帖子: 72
- 在線: 7.5小時(shí)
- 蟲(chóng)號(hào): 33303510
- 注冊(cè): 2023-03-13
- 專(zhuān)業(yè): 化學(xué)工程及工業(yè)化學(xué)
10樓2024-10-10 11:14:27
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 ZXR319 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085600材料與化工調(diào)劑 280分 +13 | yyqqhh 2026-03-03 | 13/650 |
|
|---|---|---|---|---|
|
|
[基金申請(qǐng)]
 成果系統(tǒng)訪問(wèn)量大，請(qǐng)15分鐘后再?lài)L試。由此給您造成的不便，敬請(qǐng)諒解。
+15
|
xhuama 2026-03-02 | 17/850 |
|
|
|
[考研] 紡織、生物、化學(xué)、材料等專(zhuān)業(yè) +3 | Eember. 2026-03-05 | 7/350 |
|
|
|
[考研] 313求調(diào)劑 +3 | 水流年lc 2026-02-28 | 4/200 |
|
|
|
[考研] 0856材料與化工，270求調(diào)劑 +17 | YXCT 2026-03-01 | 20/1000 |
|
|
|
[考研] 材料與化工，291分，求調(diào)劑 +9 | 咕嚕咕嚕123123 2026-03-04 | 11/550 |
|
|
[考研] 316求調(diào)劑 +3 | 林小星發(fā)大財(cái) 2026-03-04 | 3/150 |
|
|
|
[考研] 材料學(xué)碩080500復(fù)試調(diào)劑294 +3 | 四葉zjz 2026-03-04 | 3/150 |
|
|
|
[考研]
085600 英一數(shù)二272求調(diào)劑
 5+6
|
vida_a 2026-03-01 | 47/2350 |
|
|
|
[考研]
材料325求調(diào)劑
30+5
|
mariusuki 2026-03-02 | 10/500 |
|
|
|
[考研] 課題組接收材料類(lèi)調(diào)劑研究生 +6 | gaoxiaoniuma 2026-02-28 | 9/450 |
|
|
|
[考研] 1 +7 | 黑！在干嘛 2026-02-28 | 8/400 |
|
|
|
[考研] 【2026 碩士調(diào)劑】課題組 招收調(diào)劑生 +3 | 考研版棒棒 2026-03-02 | 5/250 |
|
|
|
[考研] 0856求調(diào)劑285 +11 | 呂仔龍 2026-02-28 | 11/550 |
|
|
|
[考研] 291 求調(diào)劑 +3 | 化工2026屆畢業(yè)?/a> 2026-03-02 | 3/150 |
|
|
|
[考研] 哈工大計(jì)算機(jī)劉劼團(tuán)隊(duì)招生 +4 | hit_aiot 2026-03-01 | 6/300 |
|
|
|
[考研] 295求調(diào)劑 +8 | 19171856320 2026-02-28 | 8/400 |
|
|
|
[考研] 322求調(diào)劑 +3 | 熊境喆 2026-03-01 | 3/150 |
|
|
|
[考研] 272求調(diào)劑 +6 | 田智友 2026-02-28 | 6/300 |
|
|
|
[考研] 317一志愿華南理工電氣工程求調(diào)劑 +6 | Soliloquy_Q 2026-02-28 | 11/550 |
|