 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 651 | 回復(fù): 3 | ||
| 當(dāng)前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | ||
[求助]
求大神幫忙看一下 我寫的積分法 求取動力學(xué)參數(shù)的 程式 已有2人參與
|
||
|
function KineticsEst1_int % 動力學(xué)參數(shù)辨識: 用積分法進(jìn)行反應(yīng)速率分析得到速率常數(shù)k和反應(yīng)級數(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); % 擬合效果圖(實驗與擬合的比較) figure(1) [t4plot X4plot] = ode45(@KineticsEqs,[tspan],X0,[],beta); plot(t,X(:,1),'bo',t4plot,X4plot,'k-') legend('Exp','Model') xlabel('時間 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ù)辨識結(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)) ]; 程序錯誤好多，但是不知道怎么改，尤其是 非線性擬合部分，關(guān)于矩陣beta，我是參照一本書上的一個例題寫的，但是他只有一個方程，我有5個，還請各位大神幫忙看一下。。萬分感謝�。� |
回復(fù)此樓» 猜你喜歡
285化工學(xué)碩求調(diào)劑（081700）
已經(jīng)有10人回復(fù)
-
求材料調(diào)劑
已經(jīng)有10人回復(fù)
-
本科鄭州大學(xué)物理學(xué)院，一志愿華科070200學(xué)碩，346求調(diào)劑
已經(jīng)有4人回復(fù)
-
國自科面上基金字體
已經(jīng)有8人回復(fù)
-
能源材料化學(xué)課題組招收碩士研究生8-10名
已經(jīng)有11人回復(fù)
-
材料專碩英一數(shù)二306
已經(jīng)有5人回復(fù)
-
085600材料與化工調(diào)劑 324分
已經(jīng)有9人回復(fù)
-
0703化學(xué)調(diào)劑
已經(jīng)有12人回復(fù)
-
0703化學(xué) 305求調(diào)劑
已經(jīng)有4人回復(fù)
-
0703化學(xué)調(diào)劑，求各位老師收留
已經(jīng)有10人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 航空動力學(xué)報審稿求大神幫助
已經(jīng)有3人回復(fù)
- 求ANSYS大神，熱力學(xué)問題，我都快崩潰了
已經(jīng)有8人回復(fù)
- 求助各位大神量子力學(xué)入門書籍，請知道的進(jìn)來說下
已經(jīng)有28人回復(fù)
4樓2016-11-09 11:58:17
dingd
鐵桿木蟲 (職業(yè)作家)
- 應(yīng)助: 1641 (講師)
- 金幣: 15037.3
- 散金: 101
- 紅花: 234
- 帖子: 3410
- 在線: 1223.7小時
- 蟲號: 291104
- 注冊: 2006-10-28
2樓2014-05-23 10:21:25
3樓2014-05-23 10:43:51
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿華中科技大學(xué)，080502，354分求調(diào)劑 +4 | 守候夕陽CF 2026-03-18 | 4/200 |
|
|---|---|---|---|---|
|
|
[考研] 311求調(diào)劑 +4 | 冬十三 2026-03-18 | 4/200 |
|
|
|
[考研] 267一志愿南京工業(yè)大學(xué)0817化工求調(diào)劑 +8 | SUICHILD 2026-03-12 | 8/400 |
|
|
[考研] 085600材料與化工 +5 | 安全上岸！ 2026-03-16 | 5/250 |
|
|
|
[考研] 08工科 320總分 求調(diào)劑 +5 | 梨花珞晚風(fēng) 2026-03-17 | 5/250 |
|
|
|
[考研] 297求調(diào)劑 +8 | 戲精丹丹丹 2026-03-17 | 8/400 |
|
|
|
[考研] 085601材料工程專碩求調(diào)劑 +6 | 慕寒mio 2026-03-16 | 6/300 |
|
|
|
[考研] 296求調(diào)劑 +5 | 大口吃飯 身體健 2026-03-13 | 5/250 |
|
|
|
[考研] 301求調(diào)劑 +4 | A_JiXing 2026-03-16 | 4/200 |
|
|
|
[考研] 26考研求調(diào)劑 +6 | 丶宏Sir 2026-03-13 | 6/300 |
|
|
|
[考研] 275求調(diào)劑 +4 | 太陽花天天開心 2026-03-16 | 4/200 |
|
|
|
[考研] 321求調(diào)劑 +5 | 大米飯！ 2026-03-15 | 5/250 |
|
|
|
[考研] 中科院材料273求調(diào)劑 +4 | yzydy 2026-03-15 | 4/200 |
|
|
|
[考研] 26考研一志愿中國石油大學(xué)(華東)305分求調(diào)劑 +3 | 嘉年新程 2026-03-15 | 3/150 |
|
|
|
[考研] 297求調(diào)劑 +4 | 學(xué)海漂泊 2026-03-13 | 4/200 |
|
|
|
[考研] 329求調(diào)劑 +3 | miaodesi 2026-03-12 | 4/200 |
|
|
|
[碩博家園] 085600 260分求調(diào)劑 +3 | 天空還下雨么 2026-03-13 | 5/250 |
|
|
|
[考研] 求調(diào)劑 +3 | 程雨杭 2026-03-12 | 3/150 |
|
|
|
[考研] 328化工專碩求調(diào)劑 +4 | 。，。，。，。i 2026-03-12 | 4/200 |
|
|
|
[考研] 材料301分求調(diào)劑 +5 | Liyouyumairs 2026-03-12 | 5/250 |
|