 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 6 | 1/1 | 返回列表 |
| 查看: 1006 | 回復(fù): 5 | |||
[求助]
Matlab 動力學(xué)方程擬合的很差，郁悶，大俠們幫助一下。先謝了
|
| 我把程序和擬合的結(jié)果放在文件了。新手求助 |
回復(fù)此樓» 本帖附件資源列表
-
歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺，不對該內(nèi)容負責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識產(chǎn)權(quán)問題，其責(zé)任在于用戶本人，如對版權(quán)有異議，請聯(lián)系郵箱：xiaomuchong@tal.com - 附件 1 : 求助.pdf
2013-06-13 20:20:42, 143.12 K
» 猜你喜歡
梁成偉老師課題組歡迎你的加入
已經(jīng)有11人回復(fù)
-
復(fù)試調(diào)劑
已經(jīng)有6人回復(fù)
-
0703化學(xué)調(diào)劑 ，六級已過，有科研經(jīng)歷
已經(jīng)有11人回復(fù)
-
0856調(diào)劑，是學(xué)校就去
已經(jīng)有4人回復(fù)
-
321求調(diào)劑
已經(jīng)有10人回復(fù)
-
東華理工大學(xué)化材專業(yè)26屆碩士博士申請
已經(jīng)有8人回復(fù)
-
288求調(diào)劑
已經(jīng)有11人回復(fù)
-
【考研調(diào)劑】化學(xué)專業(yè) 281分，一志愿四川大學(xué)，誠心求調(diào)劑
已經(jīng)有5人回復(fù)
-
能源材料化學(xué)課題組招收碩士研究生8-10名
已經(jīng)有12人回復(fù)
-
0703化學(xué)調(diào)劑
已經(jīng)有4人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- matlab 擬合反應(yīng)動力學(xué)參數(shù)結(jié)果很差。大家?guī)兔匆幌?/a>
已經(jīng)有14人回復(fù)
- 為啥我的一級動力學(xué)方程擬合出來的數(shù)據(jù)R是負數(shù)呢？
已經(jīng)有16人回復(fù)
- matlab 求指點 動力學(xué)方程擬合過程中導(dǎo)數(shù)的獲取
已經(jīng)有10人回復(fù)
- 【原創(chuàng)】一級動力學(xué)方程擬合求助
已經(jīng)有8人回復(fù)
- 【求助】使用Matlab擬合反應(yīng)動力學(xué)方程問題
已經(jīng)有7人回復(fù)
|
程序是這樣的 function myparafit clear all clc warning off close all k0 = [0.0000000005 0.0000000005 0.0000000005]; % initial value lb = [0 0 0 0]; % ???? ub = [1 1 1 1]; % ???? x0 = [100 0 0 0 25]; yexp=[0 100 0 0 0 20 74.1 5 3.1 17.8 40 41.4 5 12.3 41.3 60 7.2 7.5 20.9 64.4 90 2.9 15 20.4 62]; % [k,fval,flag] = fmincon(@ObjFunc7Fmincon,k0,[],[],[],[],lb,ub,[],[],x0,yexp); k_fm= k; [k,resnorm,residual,exitflag,output,lambda,jacobian] = ... lsqnonlin(@ObjFunc7LNL,k0,lb,ub,[],x0,yexp); ci = nlparci(k,residual,jacobian); k_ls = k; k0 = k_fm; [k,resnorm,residual,exitflag,output,lambda,jacobian] = ... lsqnonlin(@ObjFunc7LNL,k0,lb,ub,[],x0,yexp); ci = nlparci(k,residual,jacobian); k_fmls = k; tspan=[0 20 40 60 90]; [t x] = ode45(@KineticEqs,tspan,x0,[],k_fmls); figure; plot(t,x(:,1),'*-',t,yexp(:,2),'o-'); figure; plot(t,x(:,2:4),'*-'); hold on plot(t,yexp(:,3:5),'o:'); figure hold on k_fmls; % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% function f = ObjFunc7Fmincon(k,x0,yexp) tspan=[0 20 40 60 90]; [t x] = ode45(@KineticEqs,tspan,x0,[],k); y(:,2:5) = x(:,1:4); f =sum((y(:,2)-yexp(:,1)).^2)+sum((y(:,3)-yexp(:,2)).^2)+sum((y(:,4)-yexp(:,3)).^2)+sum((y(:,5)-yexp(:,4)).^2); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% function f = ObjFunc7LNL(k,x0,yexp) tspan=[0 20 40 60 90]; [t x] = ode45(@KineticEqs,tspan,x0,[],k); y(:,2) = x(:,1); y(:,3) = x(:,2); y(:,4) = x(:,3); y(:,5) = x(:,4); f1 = y(:,2) - yexp(:,2); f2 = y(:,3) - yexp(:,3); f3 = y(:,4) - yexp(:,4); f4 = y(:,5) - yexp(:,5); f = [f1; f2; f3; f4]; %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% function dxdt = KineticEqs(t,x,k) % global Alpha k_d W=25; m0=100; k_d=0.01; P=x(1); G=x(2); L=x(3); W=x(4); Alpha=x(5); dxdt1=-(k(1)+k(2)+k(3))*P.*Alpha; dxdt2=k(1)*P.*Alpha; dxdt3=k(2)*P.*Alpha; dxdt4=k(3)*P.*Alpha; dxdt5=-k_d*exp(-k_d*x(3)*(100/25)); dxdt=[dxdt1;dxdt2;dxdt3;dxdt4;dxdt5]; |
2樓2013-06-13 23:17:57
3樓2013-06-14 16:11:20
4樓2013-06-14 18:52:06
5樓2013-06-14 19:28:10
6樓2013-06-14 22:18:41
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 uklaoli 的主題更新
| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 【考研調(diào)劑】化學(xué)專業(yè) 281分，一志愿四川大學(xué)，誠心求調(diào)劑 +5 | 吃吃吃才有意義 2026-03-19 | 5/250 |
|
|---|---|---|---|---|
|
|
[考研] 招收調(diào)劑碩士 +4 | lidianxing 2026-03-19 | 10/500 |
|
|
|
[考研] 一志愿武漢理工材料工程專碩調(diào)劑 +4 | Doleres 2026-03-19 | 4/200 |
|
|
[考研] 307求調(diào)劑 +6 | 冷笙123 2026-03-17 | 6/300 |
|
|
|
[考研] 一志愿 西北大學(xué) ，070300化學(xué)學(xué)碩，總分287，雙非一本，求調(diào)劑。 +3 | 晨昏線與星海 2026-03-19 | 3/150 |
|
|
|
[考研] 一志愿華中科技大學(xué)，080502，354分求調(diào)劑 +4 | 守候夕陽CF 2026-03-18 | 4/200 |
|
|
|
[考研] 0854可跨調(diào)劑，一作一項核心論文五項專利，省、國級證書40+數(shù)一英一287 +8 | 小李0854 2026-03-16 | 8/400 |
|
|
|
[考研] 材料專碩306英一數(shù)二 +10 | z1z2z3879 2026-03-16 | 13/650 |
|
|
|
[考研] 312求調(diào)劑 +8 | 陌宸希 2026-03-16 | 9/450 |
|
|
|
[考博] 環(huán)境領(lǐng)域全國重點實驗室招收博士1-2名 +3 | QGZDSYS 2026-03-13 | 5/250 |
|
|
|
[考研] 265求調(diào)劑 +3 | 梁梁校校 2026-03-17 | 3/150 |
|
|
|
[考研] 268求調(diào)劑 +7 | 好運連綿不絕 2026-03-12 | 8/400 |
|
|
|
[考研] 一志愿蘇州大學(xué)材料工程（085601）專碩有科研經(jīng)歷三項國獎兩個實用型專利一項省級立項 +6 | 大火山小火山 2026-03-16 | 8/400 |
|
|
|
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6 | 薛云鵬 2026-03-15 | 6/300 |
|
|
|
[考研] 278求調(diào)劑 +3 | Yy7400 2026-03-13 | 3/150 |
|
|
|
[考研]
 11408 一志愿西電，277分求調(diào)劑
+3
|
zhouzhen654 2026-03-16 | 3/150 |
|
|
|
[考研] 304求調(diào)劑 +5 | 素年祭語 2026-03-15 | 5/250 |
|
|
|
[考研] 材料與化工 323 英一+數(shù)二+物化，一志愿：哈工大 本人本科雙一流 +4 | 自由的_飛翔 2026-03-13 | 5/250 |
|
|
|
[考研] 297求調(diào)劑 +4 | 學(xué)海漂泊 2026-03-13 | 4/200 |
|
|
|
[考研] 295求調(diào)劑 +3 | 小匕仔汁 2026-03-12 | 3/150 |
|