 |
 |
 |
| 24小時(shí)熱門版塊排行榜 |
- >論壇更新日志 (1380)
- >考研 (176)
- >導(dǎo)師招生 (118)
- >蟲友互識(shí) (62)
- >考博 (26)
- >文獻(xiàn)求助 (25)
- >論文投稿 (23)
- >碩博家園 (16)
- >博后之家 (15)
- >文獻(xiàn)求助完結(jié)子版 (14)
- >找工作 (14)
- >基金申請(qǐng) (10)
- >論文道賀祈福 (8)
- >休閑灌水 (8)
- >教師之家 (7)
- >職場(chǎng)人生 (5)
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1016 | 回復(fù): 5 | ||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||
uklaoli新蟲 (初入文壇)
|
[求助]
Matlab 動(dòng)力學(xué)方程擬合的很差，郁悶，大俠們幫助一下。先謝了
|
|
| 我把程序和擬合的結(jié)果放在文件了。新手求助 |
回復(fù)此樓» 本帖附件資源列表
-
歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺(tái)，不對(duì)該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識(shí)產(chǎn)權(quán)問(wèn)題，其責(zé)任在于用戶本人，如對(duì)版權(quán)有異議，請(qǐng)聯(lián)系郵箱：xiaomuchong@tal.com - 附件 1 : 求助.pdf
2013-06-13 20:20:42, 143.12 K
» 猜你喜歡
一志愿大連理工大學(xué)材料求調(diào)劑
已經(jīng)有5人回復(fù)
-
材料與化工調(diào)劑一志愿大連海事085600，349
已經(jīng)有6人回復(fù)
-
一志愿中海洋320化學(xué)工程與技術(shù)學(xué)碩求調(diào)劑
已經(jīng)有4人回復(fù)
-
材料科學(xué)與工程調(diào)劑
已經(jīng)有11人回復(fù)
-
311求調(diào)劑一志愿合肥工業(yè)大學(xué)
已經(jīng)有4人回復(fù)
-
一志愿食品科學(xué)與工程083200求調(diào)劑
已經(jīng)有4人回復(fù)
-
材料科學(xué)與工程求調(diào)劑
已經(jīng)有10人回復(fù)
-
一志愿北化085600材料專碩275|有文章專利｜求調(diào)劑
已經(jīng)有11人回復(fù)
-
一志愿 南京航空航天大學(xué) ，080500材料科學(xué)與工程學(xué)碩
已經(jīng)有4人回復(fù)
-
085600材料與化工調(diào)劑
已經(jīng)有15人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- matlab 擬合反應(yīng)動(dòng)力學(xué)參數(shù)結(jié)果很差。大家?guī)兔匆幌?/a>
已經(jīng)有14人回復(fù)
- 為啥我的一級(jí)動(dòng)力學(xué)方程擬合出來(lái)的數(shù)據(jù)R是負(fù)數(shù)呢？
已經(jīng)有16人回復(fù)
- matlab 求指點(diǎn) 動(dòng)力學(xué)方程擬合過(guò)程中導(dǎo)數(shù)的獲取
已經(jīng)有10人回復(fù)
- 【原創(chuàng)】一級(jí)動(dòng)力學(xué)方程擬合求助
已經(jīng)有8人回復(fù)
- 【求助】使用Matlab擬合反應(yīng)動(dòng)力學(xué)方程問(wèn)題
已經(jīng)有7人回復(fù)
uklaoli
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 156
- 帖子: 22
- 在線: 72.7小時(shí)
- 蟲號(hào): 2134172
- 注冊(cè): 2012-11-18
- 專業(yè): 塑料
4樓2013-06-14 18:52:06
uklaoli
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 156
- 帖子: 22
- 在線: 72.7小時(shí)
- 蟲號(hào): 2134172
- 注冊(cè): 2012-11-18
- 專業(yè): 塑料
|
程序是這樣的 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
uklaoli
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 156
- 帖子: 22
- 在線: 72.7小時(shí)
- 蟲號(hào): 2134172
- 注冊(cè): 2012-11-18
- 專業(yè): 塑料
3樓2013-06-14 16:11:20
5樓2013-06-14 19:28:10
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 311求調(diào)劑一志愿合肥工業(yè)大學(xué) +4 | 秋二十二 2026-03-30 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 一志愿 南京航空航天大學(xué) ，080500材料科學(xué)與工程學(xué)碩 +4 | @taotao 2026-03-30 | 4/200 |
|
|
|
[考研] 274求調(diào)劑 +4 | xiao愛同學(xué) 2026-03-30 | 4/200 |
|
|
|
[考研] 303求調(diào)劑 +7 | DLkz1314. 2026-03-30 | 7/350 |
|
|
|
[考研] 一志愿哈爾濱工業(yè)大學(xué)材料與化工方向336分 +12 | 辰沐5211314 2026-03-26 | 12/600 |
|
|
|
[考研] 南京大學(xué)化學(xué)調(diào)劑 +10 | 景隨風(fēng) 2026-03-29 | 15/750 |
|
|
|
[考研] 282求調(diào)劑 +4 | wcq131415 2026-03-24 | 4/200 |
|
|
|
[考研] 337求調(diào)劑 +6 | 《樹》 2026-03-29 | 6/300 |
|
|
|
[考研] 321求調(diào)劑 +7 | 璞玉~~ 2026-03-25 | 8/400 |
|
|
|
[考研] 11408軟件工程求調(diào)劑 +3 | Qiu學(xué)ing 2026-03-28 | 3/150 |
|
|
|
[考研] 081200-314 +3 | LILIQQ 2026-03-27 | 4/200 |
|
|
|
[考研] 308求調(diào)劑 +7 | 墨墨漠 2026-03-27 | 7/350 |
|
|
|
[考研] 求調(diào)劑 +4 | 零八# 2026-03-27 | 4/200 |
|
|
|
[考研] 266求調(diào)劑 +11 | 陽(yáng)陽(yáng)哇塞 2026-03-27 | 12/600 |
|
|
|
[考研] 化學(xué)調(diào)劑 +4 | 愛吃番茄的旭 2026-03-24 | 5/250 |
|
|
|
[考研] 314求調(diào)劑 +3 | 溪云珂 2026-03-26 | 3/150 |
|
|
|
[考研]
 085602 289分求調(diào)劑
+8
|
WWW西西弗斯 2026-03-24 | 8/400 |
|
|
|
[考研] 303求調(diào)劑 +6 | 藍(lán)山月 2026-03-25 | 6/300 |
|
|
|
[考研] 086003食品工程求調(diào)劑 +6 | 淼淼111 2026-03-24 | 6/300 |
|
|
|
[考研] 080500求調(diào)劑 +3 | zzzzfan 2026-03-24 | 3/150 |
|