 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 1206 | 回復(fù): 10 | |||
njhx505木蟲 (正式寫手)
|
[交流]
【求助】30金幣求一程序 已有4人參與
|
 請(qǐng)各位高手幫個(gè)忙，CA,CB是不同時(shí)間測(cè)得濃度值，求k1,k2  數(shù)據(jù)給出，希望高手可以幫忙，謝謝 [ Last edited by njhx505 on 2010-6-4 at 19:42 ] |
回復(fù)此樓» 猜你喜歡
一志愿武理材料305分求調(diào)劑
已經(jīng)有5人回復(fù)
-
材料專碩英一數(shù)二306
已經(jīng)有3人回復(fù)
-
085600材料與化工調(diào)劑 324分
已經(jīng)有7人回復(fù)
-
311求調(diào)劑
已經(jīng)有3人回復(fù)
-
26調(diào)劑/材料/英一數(shù)二/總分289/已過A區(qū)線
已經(jīng)有7人回復(fù)
-
295求調(diào)劑
已經(jīng)有5人回復(fù)
-
354求調(diào)劑
已經(jīng)有5人回復(fù)
-
材料專業(yè)求調(diào)劑
已經(jīng)有3人回復(fù)
-
求材料調(diào)劑
已經(jīng)有7人回復(fù)
-
一志愿天大材料與化工（085600）總分338
已經(jīng)有3人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 30金幣求助：最長(zhǎng)路徑問題，要求遍歷每個(gè)節(jié)點(diǎn)，總路徑最長(zhǎng)
已經(jīng)有7人回復(fù)
- 求高人指點(diǎn)用matlab求解非線性方程組，解決了追加100金幣；
已經(jīng)有11人回復(fù)
- 【求助】50金幣求大家?guī)兔Ω囊欢纬绦?/a>
已經(jīng)有13人回復(fù)
- 【求助】VS2005數(shù)值計(jì)算程序調(diào)試【急求！�。�！加金幣】
已經(jīng)有9人回復(fù)
- 【求助】100金幣求助高人一個(gè)程序解方程
已經(jīng)有15人回復(fù)
lijinfeng042
木蟲 (小有名氣)
Matlab
- 仿真EPI: 2
- 應(yīng)助: 1 (幼兒園)
- 金幣: 2156.1
- 散金: 115
- 帖子: 291
- 在線: 31.5小時(shí)
- 蟲號(hào): 1019062
- 注冊(cè): 2010-05-15
- 性別: GG
- 專業(yè): 通信理論與系統(tǒng)
2樓2010-06-04 17:03:55
change0618
鐵桿木蟲 (著名寫手)
方丈大師
- 應(yīng)助: 44 (小學(xué)生)
- 金幣: 17724.5
- 紅花: 17
- 帖子: 2413
- 在線: 546.7小時(shí)
- 蟲號(hào): 496517
- 注冊(cè): 2008-01-19
- 專業(yè): 化學(xué)反應(yīng)工程
3樓2010-06-04 17:31:59
change0618
鐵桿木蟲 (著名寫手)
方丈大師
- 應(yīng)助: 44 (小學(xué)生)
- 金幣: 17724.5
- 紅花: 17
- 帖子: 2413
- 在線: 546.7小時(shí)
- 蟲號(hào): 496517
- 注冊(cè): 2008-01-19
- 專業(yè): 化學(xué)反應(yīng)工程
數(shù)據(jù)有的惡劣 4樓2010-06-04 20:28:04
njhx505
木蟲 (正式寫手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 2163.9
- 散金: 6
- 帖子: 387
- 在線: 16.4小時(shí)
- 蟲號(hào): 737675
- 注冊(cè): 2009-04-01
- 性別: MM
- 專業(yè): 分離過程
5樓2010-06-04 20:39:35
★ ★
zzuwangshilei(金幣+2):辛苦了，多謝應(yīng)助 2010-06-05 10:52:39
njhx505(金幣+3):謝謝 2010-06-06 21:57:51
zzuwangshilei(金幣+2):辛苦了，多謝應(yīng)助 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 ....(*)式 計(jì)算一下你表格中的數(shù)據(jù)，不滿足是常數(shù)。至于怎么處理，或者重新采集數(shù)據(jù)或者像你提到的滿足某種最優(yōu)。 另外： 將(*)式代入到上面其中一個(gè)方程比如(1)得到： dCA/dt=-k1*CA^2-k2*CA+const; 這個(gè)方程建議還是數(shù)值求解比較方便，比如向前歐拉，向后歐拉方法等等，要想得到高精度的就用Runge-Kutta。 |
6樓2010-06-05 03:24:37
njhx505
木蟲 (正式寫手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 2163.9
- 散金: 6
- 帖子: 387
- 在線: 16.4小時(shí)
- 蟲號(hào): 737675
- 注冊(cè): 2009-04-01
- 性別: MM
- 專業(yè): 分離過程
7樓2010-06-05 12:50:54
change0618
鐵桿木蟲 (著名寫手)
方丈大師
- 應(yīng)助: 44 (小學(xué)生)
- 金幣: 17724.5
- 紅花: 17
- 帖子: 2413
- 在線: 546.7小時(shí)
- 蟲號(hào): 496517
- 注冊(cè): 2008-01-19
- 專業(yè): 化學(xué)反應(yīng)工程
★
zzuwangshilei(金幣+1):多謝積極參與 2010-06-05 15:16:11
njhx505(金幣+3): 2010-06-06 21:58:24
zzuwangshilei(金幣+1):多謝積極參與 2010-06-05 15:16:11
njhx505(金幣+3): 2010-06-06 21:58:24
|
你給的CA數(shù)據(jù)有問題 按照你的微分方程式可以得到一下關(guān)系： CA-CA0=-2*(CB-CB0) 所以我認(rèn)為你的數(shù)據(jù)擬合只需要 t=0時(shí)的 CA0 CB0 以及一組CA值或者CB值就可以了。 當(dāng)然CA，CB數(shù)據(jù)都有的話，可以驗(yàn)證你的數(shù)據(jù)的準(zhǔn)確性 [ Last edited by change0618 on 2010-6-5 at 15:02 ] |
8樓2010-06-05 14:58:06
njhx505
木蟲 (正式寫手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 2163.9
- 散金: 6
- 帖子: 387
- 在線: 16.4小時(shí)
- 蟲號(hào): 737675
- 注冊(cè): 2009-04-01
- 性別: MM
- 專業(yè): 分離過程
9樓2010-06-05 16:38:43
change0618
鐵桿木蟲 (著名寫手)
方丈大師
- 應(yīng)助: 44 (小學(xué)生)
- 金幣: 17724.5
- 紅花: 17
- 帖子: 2413
- 在線: 546.7小時(shí)
- 蟲號(hào): 496517
- 注冊(cè): 2008-01-19
- 專業(yè): 化學(xué)反應(yīng)工程
★ ★
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
|
給你一個(gè)程序，但是數(shù)據(jù)不是你的，自己修改去吧 function Kinetics % 動(dòng)力學(xué)ODE方程模型的參數(shù)估計(jì) 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第一列為時(shí)間 CB = ExpData(:,2); % ExpData第二列為組分Ｂ的濃度 CB0 = CB(1); % t=0時(shí)，組分B初始濃度 CA0 = 0.8; % t=0時(shí)，組分A初始濃度 CA = CA0-2*(ExpData(:,2)-CB0); % 由微分方程式導(dǎo)出來的CA與CB關(guān)系式 k0 = [20 50]; % 估值參數(shù)的猜想值 lb = [0 0]; % 設(shè)定的估值參數(shù)上限 ub = [+inf +inf]; % 設(shè)定的估值參數(shù)下限 [k,resnorm,residual,exitflag,output,lambda,jacobian] = ... lsqnonlin(@ObjFunc,k0,lb,ub,[],t,[CA0,CB0],[CA,CB]); ci = nlparci(k,residual,jacobian); % 計(jì)算置信區(qū)間 fprintf('\n\n估計(jì)參數(shù)值為:\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計(jì)算時(shí)間序列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
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 njhx505 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 354求調(diào)劑 +3 | Tyoumou 2026-03-18 | 5/250 |
|
|---|---|---|---|---|
|
|
[考研] 314求調(diào)劑 +8 | 無懈可擊的巨人 2026-03-12 | 8/400 |
|
|
|
[考研] 311求調(diào)劑 +11 | 冬十三 2026-03-15 | 12/600 |
|
|
[考研] 材料專碩306英一數(shù)二 +10 | z1z2z3879 2026-03-16 | 13/650 |
|
|
|
[考研] 0703化學(xué)336分求調(diào)劑 +6 | zbzihdhd 2026-03-15 | 7/350 |
|
|
|
[考研] 328求調(diào)劑，英語六級(jí)551，有科研經(jīng)歷 +3 | 生物工程調(diào)劑 2026-03-16 | 8/400 |
|
|
|
[考研] 301求調(diào)劑 +4 | A_JiXing 2026-03-16 | 4/200 |
|
|
|
[考研] 考研化學(xué)學(xué)碩調(diào)劑，一志愿985 +4 | 張vvvv 2026-03-15 | 6/300 |
|
|
|
[考研] 302求調(diào)劑 +4 | 小賈同學(xué)123 | 8/400 |
|
|
|
[考研] 070303一志愿西北大學(xué)學(xué)碩310找調(diào)劑 +5 | d如愿上岸 2026-03-12 | 8/400 |
|
|
|
[考研] 070303 總分349求調(diào)劑 +3 | LJY9966 2026-03-15 | 5/250 |
|
|
|
[考研] 289求調(diào)劑 +4 | 這么名字咋樣 2026-03-14 | 6/300 |
|
|
|
[基金申請(qǐng)] 現(xiàn)在如何回避去年的某一個(gè)專家，不知道名字 +3 | zk200107 2026-03-12 | 6/300 |
|
|
|
[考研] 學(xué)碩285求調(diào)劑 +13 | Wisjxn 2026-03-12 | 46/2300 |
|
|
|
[考研] 266求調(diào)劑 +4 | 學(xué)員97LZgn 2026-03-13 | 4/200 |
|
|
|
[考研] 0856材料與化工301求調(diào)劑 +5 | 奕束光 2026-03-13 | 5/250 |
|
|
|
[考研] 304求調(diào)劑 +7 | 7712b 2026-03-13 | 7/350 |
|
|
|
[考研] 311求調(diào)劑 +3 | 冬十三 2026-03-13 | 3/150 |
|
|
|
[考研] 307求調(diào)劑 +5 | 超級(jí)伊昂大王 2026-03-12 | 5/250 |
|
|
|
[考研] 工科材料085601 279求調(diào)劑 +8 | 困于星晨 2026-03-12 | 10/500 |
|