 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1015 | 回復(fù): 5 | ||
| 當(dāng)前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | ||
[求助]
Matlab 動力學(xué)方程擬合的很差，郁悶，大俠們幫助一下。先謝了
|
||
| 我把程序和擬合的結(jié)果放在文件了。新手求助 |
回復(fù)此樓» 本帖附件資源列表
-
歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺，不對該內(nèi)容負(fù)責(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
» 猜你喜歡
英一數(shù)一408，總分284，二戰(zhàn)真誠求調(diào)劑
已經(jīng)有5人回復(fù)
-
083000環(huán)境科學(xué)與工程調(diào)劑，總分281
已經(jīng)有4人回復(fù)
-
求調(diào)劑
已經(jīng)有3人回復(fù)
-
食品工程專碩一志愿中海洋309求調(diào)劑
已經(jīng)有10人回復(fù)
-
085600材料與化工調(diào)劑
已經(jīng)有14人回復(fù)
-
一志愿北化085600材料專碩275|有文章專利｜求調(diào)劑
已經(jīng)有10人回復(fù)
-
材料科學(xué)與工程求調(diào)劑
已經(jīng)有9人回復(fù)
-
一志愿 南京航空航天大學(xué) ，080500材料科學(xué)與工程學(xué)碩
已經(jīng)有3人回復(fù)
-
材料科學(xué)與工程調(diào)劑
已經(jīng)有10人回復(fù)
-
274求調(diào)劑
已經(jīng)有4人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- matlab 擬合反應(yīng)動力學(xué)參數(shù)結(jié)果很差。大家?guī)兔匆幌?/a>
已經(jīng)有14人回復(fù)
- 為啥我的一級動力學(xué)方程擬合出來的數(shù)據(jù)R是負(fù)數(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ù)
3樓2013-06-14 16:11:20
|
程序是這樣的 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
4樓2013-06-14 18:52:06
5樓2013-06-14 19:28:10
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿中海洋材料357 +3 | 麥恩莉. 2026-03-30 | 3/150 |
|
|---|---|---|---|---|
|
|
[考研] 一志愿哈爾濱工業(yè)大學(xué)材料與化工方向336分 +12 | 辰沐5211314 2026-03-26 | 12/600 |
|
|
[考研] 材料專碩 085600求調(diào)劑 +7 | BBQ233 2026-03-30 | 7/350 |
|
|
|
[考研] 310求調(diào)劑 +10 | 爭取九點睡 2026-03-30 | 10/500 |
|
|
|
[考研] 材料專碩調(diào)劑 +11 | 椰椰。 2026-03-29 | 11/550 |
|
|
|
[考研] 求調(diào)劑 +7 | 青春裁為三截 2026-03-29 | 7/350 |
|
|
|
[考研] 317求調(diào)劑 +10 | 蛋黃咸肉粽 2026-03-26 | 10/500 |
|
|
|
[考研] 318一志愿吉林大學(xué)生物與醫(yī)藥 求調(diào)劑 +5 | 篤行致遠(yuǎn). 2026-03-28 | 5/250 |
|
|
|
[考研] 070305高分子化學(xué)與物理 304分求調(diào)劑 +12 | c297914 2026-03-28 | 12/600 |
|
|
|
[考研] 356求調(diào)劑 +4 | gysy?s?a 2026-03-28 | 4/200 |
|
|
|
[考研] 11408軟件工程求調(diào)劑 +3 | Qiu學(xué)ing 2026-03-28 | 3/150 |
|
|
|
[考研] 299求調(diào)劑 +7 | 嗯嗯嗯嗯2 2026-03-27 | 7/350 |
|
|
|
[考研] 085701環(huán)境工程，267求調(diào)劑 +16 | minht 2026-03-26 | 16/800 |
|
|
|
[考研] 0856，材料與化工321分求調(diào)劑 +12 | 大饞小子 2026-03-27 | 13/650 |
|
|
|
[考研] 308求調(diào)劑 +7 | 墨墨漠 2026-03-27 | 7/350 |
|
|
|
[考研] 285求調(diào)劑 +4 | AZMK 2026-03-27 | 7/350 |
|
|
|
[考研] 266求調(diào)劑 +11 | 陽陽哇塞 2026-03-27 | 12/600 |
|
|
|
[考研] 085600，材料與化工321分調(diào)劑 +4 | 大饞小子 2026-03-27 | 6/300 |
|
|
|
[考研] 348求調(diào)劑 +4 | 小懶蟲不懶了 2026-03-27 | 5/250 |
|
|
|
[考研] 0854人工智能方向招收調(diào)劑 +4 | 章小魚567 2026-03-24 | 4/200 |
|