 |
 |
 |
| 24小時(shí)熱門(mén)版塊排行榜 |
- >論壇更新日志 (4510)
- >導(dǎo)師招生 (635)
- >考研 (599)
- >基金申請(qǐng) (352)
- >文獻(xiàn)求助 (243)
- >考博 (171)
- >蟲(chóng)友互識(shí) (170)
- >休閑灌水 (85)
- >碩博家園 (68)
- >論文投稿 (52)
- >招聘信息布告欄 (50)
- >博后之家 (49)
- >論文道賀祈福 (32)
- >SciFinder/Reaxys (21)
- >綠色求助(高懸賞) (21)
- >找工作 (21)
 |
 |
 |
||
 |
 |
 |
 |
 |
| 1 | 1/1 | 返回列表 |
| 查看: 800 | 回復(fù): 0 | ||
泡花alkali新蟲(chóng) (小有名氣)
|
[求助]
關(guān)于自由基聚合參數(shù)擬合的一個(gè)matlab代碼，求大神看看錯(cuò)哪了？
|
|
我在跑一篇文獻(xiàn)上給出的動(dòng)力學(xué)模型，但是最后結(jié)果不太對(duì)，懷疑是其中的參數(shù)有問(wèn)題，于是我打算用文獻(xiàn)中的實(shí)驗(yàn)數(shù)據(jù)去擬合參數(shù)，模型中的物質(zhì)有6種，文獻(xiàn)中數(shù)據(jù)有濃度值的有一種（即單體轉(zhuǎn)化率），參數(shù)一共有9個(gè)，所以我打算固定其中8個(gè)參數(shù)不變，去擬合1個(gè)參數(shù)kp，代碼如下，但是跑不出來(lái)，求大神改改。 function KineticsMX180 clear all;clc global t kp kd ki kt ka km a b f kp0 tspan t=[0 7 10 20 30]; CAm=[4.45 3.231 2.759 1.936 1.317]; CA0=[4.45 0 0 0 0 0]'; kd=2.37e-4; % 初級(jí)自由基生成速率常數(shù),1/min ki=41.22; % 鏈引發(fā)速率常數(shù),L/(mol·min) kt=2.83e11; % 鏈終止速率常數(shù),L/(mol·min) ka=18.03; % 初級(jí)自由基遷入聚合物相速率,1/min km=42.83; % 鏈轉(zhuǎn)移速率常數(shù),L/(mol·min) a=0.4; b=0.4; f=0.8; kp0=[4.25e8]; %beta0=[kd ki kt ka km a b f kp]; tspan=t; [kp,resnorm,resid,exitflag,output,lambda,jacobian]=lsqnonlin(@OptObjFunc,kp0,[],[],[],tspan,CA0,CAm); ci=nlparci(kp,resid,jacobian); [t4plot CA4plot]=ode23(@KineticsEqs,[tspan(1) tspan(end)],CA0,[],kp); [t CAc]=ode23(@KineticsEqs,tspan,CA0,[],kp); function f=OptObjFunc(kp,tspan,CA0,CAm) global t kp kd ki kt ka km a b f kp0 tspan [t CAc]=ode23(@KineticsEqs,tspan,CA0,[],kp); CAc f=(CAm'-CAc); function dCAdt=KineticsEqs(t,CA,kp) global t kp kd ki kt ka km a b f kp0 tspan kd=2.37e-4; % 初級(jí)自由基生成速率常數(shù),1/min ki=41.22; % 鏈引發(fā)速率常數(shù),L/(mol·min) kt=2.83e11; % 鏈終止速率常數(shù),L/(mol·min) ka=18.03; % 初級(jí)自由基遷入聚合物相速率,1/min km=42.83; % 鏈轉(zhuǎn)移速率常數(shù),L/(mol·min) a=0.4; b=0.4; f=0.8; CA c1=CA(1,: ); c2=CA(2,: ); c3=CA(3,: ); % C=k2*c1+k3*c2+k4*c3+k5*c4 Co=c6+c7+c8+c9 c4=CA(4,: ); c5=CA(5,: ); c6=CA(6,: ); dc1dt=-kd*c1^a*c6^b; dc2dt=f*kd*c1^a*c6^b-ka*c2; dc3dt=ka*c2-ki*c3*c5-kt*c3*c4; dc4dt=ki*c3*c5-kt*c3*c4; dc5dt=kp*c5*c4; dc6dt=-3*kd*c1^a*c6^b; dCAdt=[dc1dt;dc2dt;dc3dt;dc4dt;dc5dt;dc6dt]; 其中c1~c6的：后加了空格是因?yàn)樗鼤?huì)自動(dòng)變成微笑表情。。 |
回復(fù)此樓» 猜你喜歡
2026年循環(huán)經(jīng)濟(jì)功能材料國(guó)際會(huì)議（ICFMCE 2026）
已經(jīng)有0人回復(fù)
-
2026年第五屆電氣、電子與信息工程國(guó)際會(huì)議（ISEEIE 2026）
已經(jīng)有0人回復(fù)
-
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)?
已經(jīng)有110人回復(fù)
-
哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑
已經(jīng)有5人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑
已經(jīng)有5人回復(fù)
-
求調(diào)劑
已經(jīng)有0人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑有論文
已經(jīng)有0人回復(fù)
-
請(qǐng)問(wèn)還有沒(méi)有用Latex寫(xiě)文章的小伙伴們？
已經(jīng)有0人回復(fù)
-
光學(xué)工程學(xué)碩調(diào)劑信息
已經(jīng)有30人回復(fù)
-
歡迎加入課題組
已經(jīng)有0人回復(fù)
-
散金幣，求好運(yùn)，祝面上順利！
已經(jīng)有59人回復(fù)
找到一些相關(guān)的精華帖子，希望有用哦~
使用MATLAB擬合動(dòng)力學(xué)參數(shù)問(wèn)題
已經(jīng)有18人回復(fù)
- 如何用Matlab曲線擬合一個(gè)多輸入?yún)?shù)，而且其中一個(gè)輸入?yún)?shù)與另一個(gè)參數(shù)有關(guān)呢？
已經(jīng)有6人回復(fù)
- 拜求！matlab擬合動(dòng)力學(xué)參數(shù)
已經(jīng)有14人回復(fù)
- ]緊急求助復(fù)介電常數(shù)的模型擬合-Matlab or 1stOpt
已經(jīng)有9人回復(fù)
- 關(guān)于matlab兩個(gè)自變量、一個(gè)因變量曲線擬合問(wèn)題
已經(jīng)有3人回復(fù)
- 求助用matlab 多元非線性擬合
已經(jīng)有4人回復(fù)
- matlab 擬合求助
已經(jīng)有8人回復(fù)
- 求助matlab擬合三維曲面
已經(jīng)有10人回復(fù)
- 求助，matlab自定義微分方程擬合實(shí)驗(yàn)數(shù)據(jù)來(lái)求方程中的參數(shù)
已經(jīng)有4人回復(fù)
- MATLAB微分方程參數(shù)擬合問(wèn)題，求大神
已經(jīng)有7人回復(fù)
- Matlab 含有兩個(gè)參數(shù)的公式擬合
已經(jīng)有4人回復(fù)
- matlab數(shù)據(jù)擬合問(wèn)題求助
已經(jīng)有7人回復(fù)
- 求助 用matlab擬合三元線性回歸方程及分析
已經(jīng)有9人回復(fù)
- 1stopt或matlab如何求解以下的非線性方程并擬合出相應(yīng)曲線？
已經(jīng)有13人回復(fù)
- matlab 擬合反應(yīng)動(dòng)力學(xué)參數(shù)結(jié)果很差。大家?guī)兔匆幌?/a>
已經(jīng)有14人回復(fù)
- matlab微分方程組參數(shù)擬合，以周為單位求解，匯總后以年為單位進(jìn)行數(shù)值比較
已經(jīng)有4人回復(fù)
- matlab擬合求值
已經(jīng)有17人回復(fù)
- matlab 求指點(diǎn) 動(dòng)力學(xué)方程擬合過(guò)程中導(dǎo)數(shù)的獲取
已經(jīng)有10人回復(fù)
- matlab非線性參數(shù)擬合問(wèn)題
已經(jīng)有7人回復(fù)
- 【求助】用matlab最優(yōu)化方法進(jìn)行參數(shù)擬合
已經(jīng)有17人回復(fù)
- 【求助】使用MATLAB擬合出按照公式的系數(shù)
已經(jīng)有19人回復(fù)
- 【求助】matlab曲面擬合表達(dá)式
已經(jīng)有10人回復(fù)
- 【求助】如何使得matlab擬合的曲線強(qiáng)制經(jīng)過(guò)一個(gè)點(diǎn)
已經(jīng)有8人回復(fù)
點(diǎn)擊這里搜索更多相關(guān)資源科研從小木蟲(chóng)開(kāi)始，人人為我，我為人人
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 泡花alkali 的主題更新
| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085600材料與化工調(diào)劑 324分 +10 | llllkkkhh 2026-03-18 | 10/500 |
|
|---|---|---|---|---|
|
|
[考研] 一志愿武漢理工材料工程專(zhuān)碩調(diào)劑 +3 | Doleres 2026-03-19 | 3/150 |
|
|
[考研] 274求調(diào)劑 +6 | S.H1 2026-03-18 | 6/300 |
|
|
|
[考研] 085410人工智能專(zhuān)碩317求調(diào)劑（0854都可以） +3 | xbxudjdn 2026-03-18 | 3/150 |
|
|
|
[考研] 一志愿天津大學(xué)化學(xué)工藝專(zhuān)業(yè)（081702）315分求調(diào)劑 +10 | yangfz 2026-03-17 | 10/500 |
|
|
|
[考研] 【同濟(jì)軟件】軟件（085405）考研求調(diào)劑 +3 | 2026eternal 2026-03-18 | 3/150 |
|
|
|
[考研] 298-一志愿中國(guó)農(nóng)業(yè)大學(xué)-求調(diào)劑 +7 | 手機(jī)用戶 2026-03-17 | 7/350 |
|
|
|
[考研] 331求調(diào)劑（0703有機(jī)化學(xué) +7 | ZY-05 2026-03-13 | 8/400 |
|
|
|
[考研] 070300化學(xué)319求調(diào)劑 +6 | 錦鯉0909 2026-03-17 | 6/300 |
|
|
|
[考研] 0703化學(xué)336分求調(diào)劑 +6 | zbzihdhd 2026-03-15 | 7/350 |
|
|
|
[考研] 278求調(diào)劑 +5 | 煙火先于春 2026-03-17 | 5/250 |
|
|
|
[考研] 268求調(diào)劑 +8 | 一定有學(xué)上- 2026-03-14 | 9/450 |
|
|
|
[考研] 308求調(diào)劑 +4 | 是Lupa啊 2026-03-16 | 4/200 |
|
|
|
[考研] 梁成偉老師課題組歡迎你的加入 +8 | 一鴨鴨喲 2026-03-14 | 10/500 |
|
|
|
[考研] 東南大學(xué)364求調(diào)劑 +5 | JasonYuiui 2026-03-15 | 5/250 |
|
|
|
[考研] 326求調(diào)劑 +3 | mlpqaz03 2026-03-15 | 3/150 |
|
|
|
[考研] 327求調(diào)劑 +6 | 拾光任染 2026-03-15 | 11/550 |
|
|
|
[考研] 26考研一志愿中國(guó)石油大學(xué)(華東)305分求調(diào)劑 +3 | 嘉年新程 2026-03-15 | 3/150 |
|
|
|
[考研] 266求調(diào)劑 +4 | 學(xué)員97LZgn 2026-03-13 | 4/200 |
|
|
|
[考研] 311求調(diào)劑 +3 | 冬十三 2026-03-13 | 3/150 |
|