 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
小木蟲論壇-學術科研互動平臺 » 計算模擬區(qū) » 仿真模擬 » MATLAB » 【求助】30金幣求一程序
| 查看: 1208 | 回復: 10 | |||
njhx505木蟲 (正式寫手)
|
[交流]
【求助】30金幣求一程序 已有4人參與
|
 請各位高手幫個忙，CA,CB是不同時間測得濃度值，求k1,k2  數據給出，希望高手可以幫忙，謝謝 [ Last edited by njhx505 on 2010-6-4 at 19:42 ] |
回復此樓» 猜你喜歡
一志愿中海洋材料工程專碩330分求調劑
已經有3人回復
-
328求調劑，英語六級551，有科研經歷
已經有9人回復
-
【同濟軟件】軟件（085405）考研求調劑
已經有3人回復
-
一志愿吉林大學材料學碩321求調劑
已經有4人回復
-
267一志愿南京工業(yè)大學0817化工求調劑
已經有7人回復
-
一志愿華中科技大學，080502，354分求調劑
已經有3人回復
-
一志愿中國海洋大學，生物學，301分，求調劑
已經有4人回復
-
一志愿武理材料305分求調劑
已經有6人回復
-
材料專碩英一數二306
已經有3人回復
-
085600材料與化工調劑 324分
已經有7人回復
» 本主題相關價值貼推薦，對您同樣有幫助:
- 30金幣求助：最長路徑問題，要求遍歷每個節(jié)點，總路徑最長
已經有7人回復
- 求高人指點用matlab求解非線性方程組，解決了追加100金幣；
已經有11人回復
- 【求助】50金幣求大家?guī)兔Ω囊欢纬绦?/a>
已經有13人回復
- 【求助】VS2005數值計算程序調試【急求�。。�！加金幣】
已經有9人回復
- 【求助】100金幣求助高人一個程序解方程
已經有15人回復
lijinfeng042
木蟲 (小有名氣)
Matlab
- 仿真EPI: 2
- 應助: 1 (幼兒園)
- 金幣: 2156.1
- 散金: 115
- 帖子: 291
- 在線: 31.5小時
- 蟲號: 1019062
- 注冊: 2010-05-15
- 性別: GG
- 專業(yè): 通信理論與系統
2樓2010-06-04 17:03:55
change0618
鐵桿木蟲 (著名寫手)
方丈大師
- 應助: 44 (小學生)
- 金幣: 17724.5
- 紅花: 17
- 帖子: 2413
- 在線: 546.7小時
- 蟲號: 496517
- 注冊: 2008-01-19
- 專業(yè): 化學反應工程
3樓2010-06-04 17:31:59
change0618
鐵桿木蟲 (著名寫手)
方丈大師
- 應助: 44 (小學生)
- 金幣: 17724.5
- 紅花: 17
- 帖子: 2413
- 在線: 546.7小時
- 蟲號: 496517
- 注冊: 2008-01-19
- 專業(yè): 化學反應工程
數據有的惡劣 4樓2010-06-04 20:28:04
njhx505
木蟲 (正式寫手)
- 應助: 0 (幼兒園)
- 金幣: 2163.9
- 散金: 6
- 帖子: 387
- 在線: 16.4小時
- 蟲號: 737675
- 注冊: 2009-04-01
- 性別: MM
- 專業(yè): 分離過程
5樓2010-06-04 20:39:35
★ ★
zzuwangshilei(金幣+2):辛苦了，多謝應助 2010-06-05 10:52:39
njhx505(金幣+3):謝謝 2010-06-06 21:57:51
zzuwangshilei(金幣+2):辛苦了，多謝應助 2010-06-05 10:52:39
njhx505(金幣+3):謝謝 2010-06-06 21:57:51
|
(2)式*2+(1)式得到 dCA(t)/dt+2*dCB(t)/dt=0; 在任意小區(qū)間[t0,t]積分得到：CA(t)+2*CB(t)=const ....(*)式 計算一下你表格中的數據，不滿足是常數。至于怎么處理，或者重新采集數據或者像你提到的滿足某種最優(yōu)。 另外： 將(*)式代入到上面其中一個方程比如(1)得到： dCA/dt=-k1*CA^2-k2*CA+const; 這個方程建議還是數值求解比較方便，比如向前歐拉，向后歐拉方法等等，要想得到高精度的就用Runge-Kutta。 |
6樓2010-06-05 03:24:37
njhx505
木蟲 (正式寫手)
- 應助: 0 (幼兒園)
- 金幣: 2163.9
- 散金: 6
- 帖子: 387
- 在線: 16.4小時
- 蟲號: 737675
- 注冊: 2009-04-01
- 性別: MM
- 專業(yè): 分離過程
7樓2010-06-05 12:50:54
change0618
鐵桿木蟲 (著名寫手)
方丈大師
- 應助: 44 (小學生)
- 金幣: 17724.5
- 紅花: 17
- 帖子: 2413
- 在線: 546.7小時
- 蟲號: 496517
- 注冊: 2008-01-19
- 專業(yè): 化學反應工程
8樓2010-06-05 14:58:06
njhx505
木蟲 (正式寫手)
- 應助: 0 (幼兒園)
- 金幣: 2163.9
- 散金: 6
- 帖子: 387
- 在線: 16.4小時
- 蟲號: 737675
- 注冊: 2009-04-01
- 性別: MM
- 專業(yè): 分離過程
9樓2010-06-05 16:38:43
change0618
鐵桿木蟲 (著名寫手)
方丈大師
- 應助: 44 (小學生)
- 金幣: 17724.5
- 紅花: 17
- 帖子: 2413
- 在線: 546.7小時
- 蟲號: 496517
- 注冊: 2008-01-19
- 專業(yè): 化學反應工程
★ ★
kuhailangyu(金幣+2):歡迎積極參與 2010-06-05 22:36:53
njhx505(金幣+24):謝謝你的程序，真是救了我的急 2010-06-06 21:59:25
kuhailangyu(金幣+2):歡迎積極參與 2010-06-05 22:36:53
njhx505(金幣+24):謝謝你的程序，真是救了我的急 2010-06-06 21:59:25
|
給你一個程序，但是數據不是你的，自己修改去吧 function Kinetics % 動力學ODE方程模型的參數估計 clear all clc ExpData = ... [ 0 0.1883 0.0100 0.2047 0.0200 0.2181 0.0300 0.2291 0.0400 0.2382 0.0500 0.2459 0.0600 0.2523 0.0700 0.2576 0.0800 0.2622 0.0900 0.2660 0.1000 0.2692 0.1100 0.2719 0.1200 0.2742 0.1300 0.2761 0.1400 0.2777 0.1500 0.2790 0.1600 0.2801 0.1700 0.2811 0.1800 0.2819 0.1900 0.2825 0.2000 0.2830 ]; t = ExpData(:,1); % ExpData第一列為時間 CB = ExpData(:,2); % ExpData第二列為組分Ｂ的濃度 CB0 = CB(1); % t=0時，組分B初始濃度 CA0 = 0.8; % t=0時，組分A初始濃度 CA = CA0-2*(ExpData(:,2)-CB0); % 由微分方程式導出來的CA與CB關系式 k0 = [20 50]; % 估值參數的猜想值 lb = [0 0]; % 設定的估值參數上限 ub = [+inf +inf]; % 設定的估值參數下限 [k,resnorm,residual,exitflag,output,lambda,jacobian] = ... lsqnonlin(@ObjFunc,k0,lb,ub,[],t,[CA0,CB0],[CA,CB]); ci = nlparci(k,residual,jacobian); % 計算置信區(qū)間 fprintf('\n\n估計參數值為:\n') fprintf('\tk1 = %.6f\t置信區(qū)間：[%.6f %.6f]\n',k(1),ci(1,: )) fprintf('\tk2 = %.6f\t置信區(qū)間：[%.6f %.6f]\n',k(2),ci(2,: )) tt = linspace(t(1),t(end),101); [tt C] = ode45(@KineticEqs,tt,[CA0,CB0],[],k); % 由估值k計算時間序列tt下的A，B濃度 figure(1) plot(t,CA,'o',tt,C(:,1),'r-') xlabel('t'); ylabel('C_A') figure(2) plot(t,CB,'o',tt,C(:,2),'r-') xlabel('t'); ylabel('C_B') % ------------------------------------------------------------------ function f = ObjFunc(k,tspan,x0,yexp) [t,y] = ode45(@KineticEqs,tspan,x0,[],k); f = y(: ) - yexp(: ); % ------------------------------------------------------------------ function dxdt = KineticEqs(t,x,k) dxdt =[-k(1)*x(1)^2+2*k(2)*x(2); 0.5*k(1)*x(1)^2-k(2)*x(2)]; [ Last edited by change0618 on 2010-6-5 at 18:23 ] |
10樓2010-06-05 18:15:24
相關版塊跳轉
我要訂閱樓主 njhx505 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085600材料與化工 +5 | 安全上岸！ 2026-03-16 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 化學工程321分求調劑 +15 | 大米飯！ 2026-03-15 | 18/900 |
|
|
|
[考研] 314求調劑 +8 | 無懈可擊的巨人 2026-03-12 | 8/400 |
|
|
|
[考研] 一志愿985，本科211，0817化學工程與技術319求調劑 +6 | Liwangman 2026-03-15 | 6/300 |
|
|
|
[考研] 302求調劑 +10 | 呼呼呼。。。。 2026-03-17 | 10/500 |
|
|
|
[考研] 303求調劑 +4 | 睿08 2026-03-17 | 6/300 |
|
|
|
[考研] 0703化學336分求調劑 +6 | zbzihdhd 2026-03-15 | 7/350 |
|
|
|
[考研] 環(huán)境工程調劑 +8 | 大可digkids 2026-03-16 | 8/400 |
|
|
|
[考研] 308求調劑 +4 | 是Lupa啊 2026-03-16 | 4/200 |
|
|
|
[考研] 【0856】化學工程（085602）313 分，本科學科評估A類院�；瘜W工程與工藝，誠求調劑 +7 | 小劉快快上岸 2026-03-11 | 8/400 |
|
|
|
[考研] 085600材料與化工求調劑 +5 | 緒幸與子 2026-03-17 | 5/250 |
|
|
|
[考研] 一志愿，福州大學材料專碩339分求調劑 +3 | 木子momo青爭 2026-03-15 | 3/150 |
|
|
|
[考研] 中科院材料273求調劑 +4 | yzydy 2026-03-15 | 4/200 |
|
|
|
[考研]
 機械專碩調劑
+3
|
笨笨兔子 2026-03-12 | 3/150 |
|
|
|
[考研] 26考研一志愿中國石油大學(華東)305分求調劑 +3 | 嘉年新程 2026-03-15 | 3/150 |
|
|
|
[考研] 330求調劑 +3 | ?醬給調劑跪了 2026-03-13 | 3/150 |
|
|
|
[考研] 招收0805（材料）調劑 +3 | 18595523086 2026-03-13 | 3/150 |
|
|
|
[考研] 311求調劑 +3 | 冬十三 2026-03-13 | 3/150 |
|
|
|
[考研] 求調劑 +3 | 程雨杭 2026-03-12 | 3/150 |
|
|
|
[考研] 081200-11408-276學碩求調劑 +3 | 崔wj 2026-03-12 | 4/200 |
|