 |
 |
 |
| 24小時(shí)熱門(mén)版塊排行榜 |
- >論壇更新日志 (725)
- >考研 (52)
- >蟲(chóng)友互識(shí) (41)
- >導(dǎo)師招生 (34)
- >考博 (17)
- >休閑灌水 (11)
- >基金申請(qǐng) (6)
- >論文投稿 (5)
- >招聘信息布告欄 (4)
- >碩博家園 (4)
- >博后之家 (2)
- >經(jīng)濟(jì)學(xué) (1)
- >化工設(shè)備 (1)
- >機(jī)械 (1)
- >外文書(shū)籍求助 (1)
- >文獻(xiàn)求助 (1)
 |
 |
 |
||
 |
 |
 |
 |
 |
| 4 | 1/1 | 返回列表 |
| 查看: 648 | 回復(fù): 3 | ||
zxqworld新蟲(chóng) (初入文壇)
|
[求助]
求大神幫忙看一下 我寫(xiě)的積分法 求取動(dòng)力學(xué)參數(shù)的 程式 已有2人參與
|
|
function KineticsEst1_int % 動(dòng)力學(xué)參數(shù)辨識(shí): 用積分法進(jìn)行反應(yīng)速率分析得到速率常數(shù)k和反應(yīng)級(jí)數(shù)n clear all clc global X t =[30:30:120]; X1=[67.37 58.28 45.32 32.1]; X2=[7.88 17.88 17.45 26.43]; X3=[2.48 9.09 11.84 13.65]; X4=[1.58 0 0.94 0]; X5=[20.69 14.75 24.45 27.82]; % 非線性擬合 k0 = [0.00594 0.00391 0.00222 0.00810 0.00020 0.01071 ]; n0 =[1 1 1 1.43 1.38 1.26]; beta0=[k0; n0]; tspan = [30:120]; X0 = [67.37 7.88 2.48 1.58 20.69]; [beta,resnorm,residual,exitflag,output,lambda,jacobian] = ... lsqnonlin(@OptObjFunc,beta0,[],[],[],tspan,X0,X); ci = nlparci(beta,resid,jacobian); % 擬合效果圖(實(shí)驗(yàn)與擬合的比較) figure(1) [t4plot X4plot] = ode45(@KineticsEqs,[tspan],X0,[],beta); plot(t,X(:,1),'bo',t4plot,X4plot,'k-') legend('Exp','Model') xlabel('時(shí)間 min') ylabel('收率C_C, %') % 殘差關(guān)于擬合值的殘差圖 [t Xc] = ode45(@KineticsEqs,tspan,X0,[],k,n); figure plot(Xc,resid,'*') xlabel('收率擬合值（%）') ylabel('殘差R (%)') refline(0,0) % 參數(shù)辨識(shí)結(jié)果 fprintf('Estimated Parameters:\n') fprintf('\tk = %.4f ± %.4f\n',k,ci(1,2)-k) fprintf('\tn = %.2f ± %.2f\n',n,ci(2,2)-n) % ------------------------------------------------------------------ function f = OptObjFunc(k,n,tspan,X0,X) tspan=[30:30:120] [t Xc] = ode45(@KineticsEqs,tspan,X0,[],k,n); f1 = Xc(:,1) - X(:,1); f2 = Xc(:,2) - X(:,2); f3 = Xc(:,3) - X(:,3); f4 = Xc(:,4) - X(:,4); f5 = Xc(:,5) - X(:,5); f = f1+ f2+ f3+ f4+ f5; % ------------------------------------------------------------------ function dXdt = KineticsEqs(t,X,k,n) dXdt = ... [ ( -(k(1)+k(2)+k(3)+k(4))*X(1)^n(1) ) ( k(2)*X(1)^n(2)) ( k(1)*X(1)^n(1)+k(5)*X(4)^n(5)+k(6)*X(5)^n(6) ) ( k(3)^X(1)^n(3)-k(5)*X(4)^n(5)) ( k(4)*X(1)^n(4)-k(6)*X(5)^n(6)) ]; 程序錯(cuò)誤好多，但是不知道怎么改，尤其是 非線性擬合部分，關(guān)于矩陣beta，我是參照一本書(shū)上的一個(gè)例題寫(xiě)的，但是他只有一個(gè)方程，我有5個(gè)，還請(qǐng)各位大神幫忙看一下。。萬(wàn)分感謝！！ |
回復(fù)此樓» 猜你喜歡
281求調(diào)劑（0805）
已經(jīng)有16人回復(fù)
-
304求調(diào)劑
已經(jīng)有6人回復(fù)
-
材料工程專(zhuān)碩調(diào)劑
已經(jīng)有6人回復(fù)
-
一志愿天大材料與化工（085600）總分338
已經(jīng)有4人回復(fù)
-
085700資源與環(huán)境308求調(diào)劑
已經(jīng)有3人回復(fù)
-
求材料調(diào)劑
已經(jīng)有8人回復(fù)
-
294求調(diào)劑材料與化工專(zhuān)碩
已經(jīng)有5人回復(fù)
-
一志愿華中科技大學(xué)，080502，354分求調(diào)劑
已經(jīng)有4人回復(fù)
-
一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑
已經(jīng)有6人回復(fù)
-
085410人工智能專(zhuān)碩317求調(diào)劑（0854都可以）
已經(jīng)有3人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 航空動(dòng)力學(xué)報(bào)審稿求大神幫助
已經(jīng)有3人回復(fù)
- 求ANSYS大神，熱力學(xué)問(wèn)題，我都快崩潰了
已經(jīng)有8人回復(fù)
- 求助各位大神量子力學(xué)入門(mén)書(shū)籍，請(qǐng)知道的進(jìn)來(lái)說(shuō)下
已經(jīng)有28人回復(fù)
dingd
鐵桿木蟲(chóng) (職業(yè)作家)
- 應(yīng)助: 1641 (講師)
- 金幣: 15037.3
- 散金: 101
- 紅花: 234
- 帖子: 3410
- 在線: 1223.7小時(shí)
- 蟲(chóng)號(hào): 291104
- 注冊(cè): 2006-10-28
2樓2014-05-23 10:21:25
zxqworld
新蟲(chóng) (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 帖子: 8
- 在線: 2.1小時(shí)
- 蟲(chóng)號(hào): 2053599
- 注冊(cè): 2012-10-10
- 專(zhuān)業(yè): 有機(jī)合成
3樓2014-05-23 10:43:51
云少爺
新蟲(chóng) (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 3
- 帖子: 4
- 在線: 4.5小時(shí)
- 蟲(chóng)號(hào): 2063921
- 注冊(cè): 2012-10-15
4樓2016-11-09 11:58:17
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 zxqworld 的主題更新
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料工程專(zhuān)碩調(diào)劑 +5 | 204818@lcx 2026-03-17 | 6/300 |
|
|---|---|---|---|---|
|
|
[考研] 328求調(diào)劑，英語(yǔ)六級(jí)551，有科研經(jīng)歷 +3 | 生物工程調(diào)劑 2026-03-16 | 10/500 |
|
|
[考研] 344求調(diào)劑 +6 | knight344 2026-03-16 | 7/350 |
|
|
|
[考研] 0703化學(xué) 305求調(diào)劑 +3 | FY_yy 2026-03-14 | 3/150 |
|
|
|
[考研] 一志愿武理材料305分求調(diào)劑 +5 | 想上岸的鯉魚(yú) 2026-03-18 | 6/300 |
|
|
|
[考研] 生物學(xué)071000 329分求調(diào)劑 +3 | 我愛(ài)生物生物愛(ài)?/a> 2026-03-17 | 3/150 |
|
|
|
[考研] 301求調(diào)劑 +9 | yy要上岸呀 2026-03-17 | 9/450 |
|
|
|
[考博] 26博士申請(qǐng) +3 | 1042136743 2026-03-17 | 3/150 |
|
|
|
[考研] 考研化學(xué)學(xué)碩調(diào)劑，一志愿985 +4 | 張vvvv 2026-03-15 | 6/300 |
|
|
|
[考研] 26考研求調(diào)劑 +6 | 丶宏Sir 2026-03-13 | 6/300 |
|
|
|
[考研] 302求調(diào)劑 +4 | 小賈同學(xué)123 2026-03-15 | 8/400 |
|
|
|
[考研] 考研調(diào)劑 +3 | 淇ya_~ 2026-03-17 | 5/250 |
|
|
|
[考研] 085600材料與化工 求調(diào)劑 +13 | enenenhui 2026-03-13 | 14/700 |
|
|
|
[考研] 285求調(diào)劑 +6 | ytter 2026-03-12 | 6/300 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 290分有科研經(jīng)歷，論文在投 +7 | 膩膩gk 2026-03-14 | 7/350 |
|
|
|
[考研] 327求調(diào)劑 +6 | 拾光任染 2026-03-15 | 11/550 |
|
|
|
[考研] 材料080500調(diào)劑求收留 +3 | 一顆meteor 2026-03-13 | 3/150 |
|
|
|
[碩博家園] 085600 260分求調(diào)劑 +3 | 天空還下雨么 2026-03-13 | 5/250 |
|
|
|
[考研] 求調(diào)劑 +5 | 一定有學(xué)上- 2026-03-12 | 5/250 |
|
|
|
[考研] 材料301分求調(diào)劑 +5 | Liyouyumairs 2026-03-12 | 5/250 |
|