 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 6 | 1/1 | 返回列表 |
| 查看: 1004 | 回復(fù): 5 | ||
[求助]
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
» 猜你喜歡
271材料工程求調(diào)劑
已經(jīng)有5人回復(fù)
-
281求調(diào)劑（0805）
已經(jīng)有16人回復(fù)
-
304求調(diào)劑
已經(jīng)有6人回復(fù)
-
材料工程專碩調(diào)劑
已經(jīng)有6人回復(fù)
-
一志愿天大材料與化工（085600）總分338
已經(jīng)有4人回復(fù)
-
085700資源與環(huán)境308求調(diào)劑
已經(jīng)有3人回復(fù)
-
求材料調(diào)劑
已經(jīng)有8人回復(fù)
-
294求調(diào)劑材料與化工專碩
已經(jīng)有5人回復(fù)
-
一志愿華中科技大學(xué)，080502，354分求調(diào)劑
已經(jīng)有4人回復(fù)
-
一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑
已經(jīng)有6人回復(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ù)
|
程序是這樣的 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ā)表 | |
|---|---|---|---|---|
|
[考研] 328求調(diào)劑，英語六級551，有科研經(jīng)歷 +3 | 生物工程調(diào)劑 2026-03-17 | 7/350 |
|
|---|---|---|---|---|
|
|
[考研] 070303一志愿西北大學(xué)學(xué)碩310找調(diào)劑 +6 | d如愿上岸 2026-03-12 | 9/450 |
|
|
|
[考研] 0703化學(xué) 305求調(diào)劑 +3 | FY_yy 2026-03-14 | 3/150 |
|
|
|
[考研] 【同濟軟件】軟件（085405）考研求調(diào)劑 +3 | 2026eternal 2026-03-18 | 3/150 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 +3 | 妮妮ninicgb 2026-03-17 | 3/150 |
|
|
[考研] 307求調(diào)劑 +3 | 冷笙123 2026-03-17 | 3/150 |
|
|
|
[考研] 268求調(diào)劑 +7 | 好運連綿不絕 2026-03-12 | 8/400 |
|
|
|
[考研] 本人考085602 化學(xué)工程 專碩 +16 | 不知道叫什么！ 2026-03-15 | 18/900 |
|
|
|
[考研] 332求調(diào)劑 +6 | Zz版 2026-03-13 | 6/300 |
|
|
|
[論文投稿] 有沒有大佬發(fā)小論文能帶我個二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|
|
|
[考研] 一志愿，福州大學(xué)材料專碩339分求調(diào)劑 +3 | 木子momo青爭 2026-03-15 | 3/150 |
|
|
|
[考研] 283求調(diào)劑 +3 | 聽風(fēng)就是雨； 2026-03-16 | 3/150 |
|
|
|
[考研] 藥學(xué)383 求調(diào)劑 +3 | 藥學(xué)chy 2026-03-15 | 4/200 |
|
|
|
[考研] 304求調(diào)劑 +3 | 曼殊2266 2026-03-14 | 3/150 |
|
|
|
[考研] 中科大材料與化工319求調(diào)劑 +3 | 孟鑫材料 2026-03-14 | 3/150 |
|
|
|
[考研] 本科南京大學(xué)一志愿川大藥學(xué)327 +3 | 麥田耕者 2026-03-14 | 3/150 |
|
|
|
[考研] 289求調(diào)劑 +4 | 這么名字咋樣 2026-03-14 | 6/300 |
|
|
|
[碩博家園] 085600 260分求調(diào)劑 +3 | 天空還下雨么 2026-03-13 | 5/250 |
|
|
|
[考研] 328化工專碩求調(diào)劑 +4 | 。，。，。，。i 2026-03-12 | 4/200 |
|
|
|
[考研] 321求調(diào)劑（食品/專碩） +3 | xc321 2026-03-12 | 6/300 |
|