 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 1 | 1/1 | 返回列表 |
| 查看: 533 | 回復: 0 | ||
2011207156金蟲 (小有名氣)
|
[求助]
matlab擬合動力學參數(shù)，非常謝謝~
|
|
function butongshike201506171 clear all;clc format long tspan=(0:10:390); yexp=[0.046125761496241 0.078139373146412 0.117158327134899 0.159846402230616 0.199412869433285 0.231111460229405 0.259814231746068 0.282697037233637 0.306493661349134 0.326396561089024 0.345295684999190 0.360444423309371 0.377418750129898 0.395250807577438 0.414965633800510 0.430017989495138 0.445923212579888 0.457995753035510 0.469186259425871 0.480284344291621 0.491524402631167 0.501658351934599 0.510935126510884 0.520202070112304 0.531309455870521 0.537782255717000 0.547944130180314 0.558094213007469 0.565468701311900 0.573757577102420 0.581007002009939 0.589392949772071 0.597541754780127 0.607751809867969 0.613129267171725 0.620580135887054 0.626072239599544 0.632475341850705 0.638742438783615]; beta0=[67971 0.000839]; y0=0; lb=[-1 -1]*1e6; ub=[1e9 1e9]; yy=[0 yexp']; k0=beta0; [k,resnorm,residual,exitflag,output,lambda,jacobian] = ... lsqnonlin(@ObjFunc,k0,lb,ub,[],tspan,y0,yexp); fprintf('\n\n使用函數(shù)lsqnonlin()估計得到的參數(shù)值為:\n') fprintf('\t待擬合參數(shù) k1 = %.6f\n',k(1)) fprintf('\t待擬合參數(shù) k2= %.6f\n',k(2)) fprintf(' 殘差平方和= %.6f\n\n',resnorm) ts=0:10:max(tspan); [ts ys]=ode45(@KineticsEqs,ts,y0,[],k) [ttt XXsim] = ode45(@KineticsEqs,tspan,y0,[],k); y=XXsim(2:end); xexp=yexp; R2=1-sum((xexp-y).^2)./sum((xexp-mean(y)).^2); fprintf('\n\t決定系數(shù)R-Square = %.6f',R2); figure(1) plot(ts,ys,'b',tspan,yy,'or'),legend('計算值','實驗值','Location','best'); xlabel('x'),ylabel('y') %--------------------------------------------------------- function f = ObjFunc(k,tspan,y0,yexp) % 目標函數(shù) [t Xsim] = ode45(@KineticsEqs,tspan,y0,[],k) ; ysim = Xsim(2:end); size(ysim); size(yexp); f=ysim-yexp; %---------------------------------------------------------- function dydt = KineticsEqs(t,y,k) beta(1)=k(1); beta(2)=k(2); tspan=(0:10:390); h=0.068031; p=[3.472000000000000 3.519000000000000 3.548000000000000 3.589000000000000 3.634000000000000 3.676000000000000 3.710000000000000 3.740000000000000 3.765000000000000 3.789000000000000 3.811000000000000 3.832000000000000 3.847000000000000 3.865000000000000 3.885000000000000 3.906000000000000 3.921000000000000 3.938000000000000 3.952000000000000 3.964000000000000 3.977000000000000 3.988000000000000 4.000000000000000 4.010000000000000 4.020000000000000 4.032000000000000 4.039000000000000 4.050000000000000 4.061000000000000 4.069000000000000 4.078000000000000 4.087000000000000 4.095000000000000 4.105000000000000 4.115000000000000 4.122000000000000 4.129000000000000 4.135000000000000 4.142000000000000 4.150000000000000]'; D=p*h; sp=csaps(tspan,D,1); pp=fnder(sp); dD=fnval(pp,tspan); dydt= (83.38-0.0014*exp(1600*(1/277.55-1/298.15))... *interp1(tspan,p,t,'pchip')*10000-1/beta(2)... *interp1(tspan,D,t,'pchip')*4.229*(0.0598... -(y*0.0150*(54.2+17.68)+(0.236204-interp1(tspan,D,t,'pchip'))... *0.0598*4.229)/(45.07-interp1(tspan,D,t,'pchip')*4.229))... *interp1(tspan,dD,t,'pchip')/(-45.07+interp1(tspan,D,t,'pchip')... *4.229)/0.0150-(0.0598-(y*0.0150*(54.2+17.68)+(0.236204-... interp1(tspan,D,t,'pchip'))*0.0598*4.229)/... (45.07-interp1(tspan,D,t,'pchip')*4.229))*... interp1(tspan,dD,t,'pchip')/0.0150/beta(1))/... (1/beta(1)+1/beta(2)-1/beta(2)*interp1(tspan,D,t,'pchip')... *71.88/(-45.07+4.229*interp1(tspan,D,t,'pchip'))); 以上是求動力學參數(shù)beta(1)和beta(2)的代碼，p,D,和dD是隨時間變化的已知量，不知以上程序運行時interp1(tspan,D,t,'pchip')能否對應上不同時刻的D。非常謝謝。 |
回復此樓» 猜你喜歡
哈爾濱理工大學物理系招收物理學考研調(diào)劑
已經(jīng)有7人回復
-
0702一志愿吉大B區(qū)求調(diào)劑
已經(jīng)有5人回復
-
物理學I論文潤色/翻譯怎么收費?
已經(jīng)有96人回復
-
求調(diào)劑
已經(jīng)有0人回復
-
0702一志愿吉大B區(qū)求調(diào)劑有論文
已經(jīng)有0人回復
-
請問還有沒有用Latex寫文章的小伙伴們？
已經(jīng)有0人回復
-
光學工程學碩調(diào)劑信息
已經(jīng)有31人回復
-
歡迎加入課題組
已經(jīng)有0人回復
-
散金幣，求好運，祝面上順利！
已經(jīng)有68人回復
-
華東師范大學芯片設計徐瓏真實水平咋樣
已經(jīng)有1人回復
-
調(diào)劑
已經(jīng)有0人回復
找到一些相關的精華帖子，希望有用哦~
對反應動力學參數(shù)的擬合及表觀反應速率常數(shù)的計算
已經(jīng)有5人回復
- 動力學參數(shù)擬合
已經(jīng)有3人回復
- 這樣的動力學方程能夠有1stopt擬合出參數(shù)來么
已經(jīng)有4人回復
- 使用MATLAB擬合動力學參數(shù)問題
已經(jīng)有18人回復
- Matlab同時擬合活化能，反應級數(shù)，速率常數(shù)
已經(jīng)有26人回復
- 使用matlab最優(yōu)化方法擬合獲得多個動力學參數(shù)中的問題
已經(jīng)有4人回復
- 求助用matlab擬合動力學方程
已經(jīng)有9人回復
- 求助用matlab擬合動力學方程
已經(jīng)有13人回復
- matlab 擬合求助
已經(jīng)有8人回復
- 關于matlab擬合
已經(jīng)有4人回復
- 怎樣用一組參數(shù)同時擬合兩個曲線--matlab
已經(jīng)有5人回復
- MATLAB微分方程參數(shù)擬合問題，求大神
已經(jīng)有7人回復
- 誰有成熟的求動力學參數(shù)的MATLAB的程序。
已經(jīng)有11人回復
- matlab 擬合反應動力學參數(shù)結(jié)果很差。大家?guī)兔匆幌?/a>
已經(jīng)有14人回復
- 請教matlab反應動力學參數(shù)估計遇到的問題，謝謝
已經(jīng)有15人回復
- matlab如何編寫共享參數(shù)擬合程序
已經(jīng)有8人回復
- 非線性擬合估計方程參數(shù)（matlab，spss，origin都可以）
已經(jīng)有4人回復
- 【求助】用matlab最優(yōu)化方法進行參數(shù)擬合
已經(jīng)有17人回復
- 【求助】使用MATLAB怎么實現(xiàn)擬合力場參數(shù)的程序？【已解決】
已經(jīng)有15人回復
點擊這里搜索更多相關資源科研從小木蟲開始，人人為我，我為人人
相關版塊跳轉(zhuǎn)
我要訂閱樓主 2011207156 的主題更新
| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿北京化工大學070300 學碩336求調(diào)劑 +4 | vv迷 2026-03-21 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 297求調(diào)劑 +3 | 喜歡還是不甘心 2026-03-20 | 3/150 |
|
|
|
[考研] 277材料科學與工程080500求調(diào)劑 +6 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|
|
[考研] 306求0703調(diào)劑一志愿華中師范 +5 | 紙魚ly 2026-03-21 | 5/250 |
|
|
|
[考研] 能源材料化學課題組招收碩士研究生8-10名 +5 | 脫穎而出 2026-03-16 | 15/750 |
|
|
|
[考研] 324分 085600材料化工求調(diào)劑 +4 | llllkkkhh 2026-03-18 | 4/200 |
|
|
|
[考研] 295求調(diào)劑 +4 | 一志愿京區(qū)211 2026-03-18 | 6/300 |
|
|
|
[考研] 294求調(diào)劑材料與化工專碩 +15 | 陌の森林 2026-03-18 | 15/750 |
|
|
|
[考研] 324求調(diào)劑 +5 | lucky呀呀呀鴨 2026-03-20 | 5/250 |
|
|
|
[考研] 287求調(diào)劑 +7 | 晨昏線與星海 2026-03-19 | 8/400 |
|
|
|
[考研] 一志愿 西北大學 ，070300化學學碩，總分287，雙非一本，求調(diào)劑。 +4 | 晨昏線與星海 2026-03-19 | 4/200 |
|
|
|
[考研] 260求調(diào)劑 +3 | 朱芷琳 2026-03-20 | 3/150 |
|
|
|
[考研] 材料學求調(diào)劑 +4 | Stella_Yao 2026-03-20 | 4/200 |
|
|
|
[考研] 廣西大學家禽遺傳育種課題組2026年碩士招生（接收計算機專業(yè)調(diào)劑） +3 | 123阿標 2026-03-17 | 3/150 |
|
|
|
[考研] 312求調(diào)劑 +8 | 陌宸希 2026-03-16 | 9/450 |
|
|
|
[考博] 26博士申請 +3 | 1042136743 2026-03-17 | 3/150 |
|
|
|
[考研] 301求調(diào)劑 +4 | A_JiXing 2026-03-16 | 4/200 |
|
|
|
[碩博家園] 湖北工業(yè)大學 生命科學與健康學院-課題組招收2026級食品/生物方向碩士 +3 | 1喜春8 2026-03-17 | 5/250 |
|
|
|
[論文投稿] 有沒有大佬發(fā)小論文能帶我個二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|
|
|
[考研] 321求調(diào)劑 +5 | 大米飯！ 2026-03-15 | 5/250 |
|