 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 902 | 回復(fù): 5 | ||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||
805592836_01鐵蟲 (初入文壇)
|
[求助]
未完成的matlab 已有2人參與
|
|
|
有實(shí)驗(yàn)測(cè)得的一系列時(shí)間點(diǎn) t=（0,180,300,420,600,900,1200,1800,2400） CA=（1.845，1.414,1.237,1.065,0.873,0.629,0.440,0.226,0.124） 微分方程：k(1)*k(2)*k(3)*CA0*(CB.^3)./(1+k(2)*CA+k(3)^(1/6)*CB.^(1/2)).^7 CB=0.4為定值 希望由此得到三個(gè)參數(shù)k(1),k(2),k(3)，以及實(shí)驗(yàn)值與擬合值的比較圖�，F(xiàn)已經(jīng)編了部分代碼，但不是很確定對(duì)否，希望高手可以修改一下 clear all; clc k0=[1,1,1]; lb=[0,0,0]; ub=[+inf,+inf,+inf]; tspan=[0,180,300,420,600,900,1200,1800,2400]; yexp=[1.414,1.237,1.065,0.873,0.629,0.440,0.226,0.124]'; y0=1.845; % 使用fmincon()進(jìn)行參數(shù)估計(jì) [k,fval,flag]=fmincon(@ObjFunc4Fmincon,k0,lb,ub,option,tspan,y0,yexp); fprintf('\tk1=%.4f\n',k(1)), %本征動(dòng)力學(xué)參數(shù) fprintf('\tk2=%.4f\n',k(2)), % A的吸附平衡常數(shù)KA fprintf('\tk3=%.4f\n',k(3)), % H2的吸附平衡常數(shù)KB fprintf('The sum of the squares is:%.1e\n\n',fval), k_fmincon=k; % 使用函數(shù)Isqnonlin()進(jìn)行參數(shù)估計(jì) yy=[y0 yexp']; options = optimset('MaxFunEvals',100000) [k,resmorm,residual,exitflag,output,lambda,jacobian]=... Isqnonlin(@ObjFunc4LNL,k0,lb,ub,options,tspan,y0,yexp); ci=nlparci(k,residual,jacobian); fprintf('n\n使用Iaqunonlin()估計(jì)得到的參數(shù)值:\n'),Output % 以函數(shù)fmincon()估計(jì)得到的結(jié)果為初值，使用Isqnonlin()進(jìn)行參數(shù)估計(jì) k0=k_fmincon; [k,resmorm,residual,exitflag,output,lambda,jacobian]=... Isqnonlin(@ObjFunc4LNL,k0,lb,ub,options,tspan,y0,yexp); ci=nlparci(k,residual,jacobian); fprintf('n\n以fmincon()結(jié)果為初值，使用函數(shù)Isqnonlin()估計(jì)得到的參數(shù)值為:\n') Output %------------------------------------------------------------ function f=ObjFunc4Fmincon(k,tspan,y0,yexp) [t y]=ode45(@KineticsEqs,tspan,y0,yexp,k); %---------------------------------------------------------- function f=Objfunc4LNL(k,tspan,yexp) [t y]=ode45(@KineticEqs,tspan,tspan,[],k); %---------------------------------------------------------- function dydt = KineticsEqs(t,y,k) y0=1.845; CB=y; CB0=0.4; CA0=y0; k1=k(1); KA=k(2); KB=k(3); dxdt=k(1)*k(2)*k(3)*CA0*(CB.^3)./(1+k(2)*CA0+k(3)^(1/6)*CB.^(1/2)).^7 |
回復(fù)此樓» 猜你喜歡
歡迎調(diào)劑到贛南師范大學(xué)智能制造與未來(lái)能源學(xué)院
已經(jīng)有27人回復(fù)
-
重慶交大26年碩士生招生擬調(diào)劑通知已出！歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。
已經(jīng)有0人回復(fù)
-
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)?
已經(jīng)有104人回復(fù)
-
UJN物理學(xué)專業(yè)調(diào)劑
已經(jīng)有12人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑 本科期間發(fā)表一篇Sci
已經(jīng)有0人回復(fù)
-
北京納米能源與系統(tǒng)研究所王中林院士/曹南穎研究員課題組招收2026碩士研究生1名
已經(jīng)有0人回復(fù)
-
山西大同大學(xué)物理學(xué)碩士研究生招收調(diào)劑生
已經(jīng)有0人回復(fù)
-
[調(diào)劑信息]211智能人工感知方向國(guó)家青年特聘專家課題組招收調(diào)劑研究生
已經(jīng)有0人回復(fù)
-
［調(diào)劑信息］211智能人工感知方向國(guó)家青年特聘專家課題組招收調(diào)劑研究生
已經(jīng)有0人回復(fù)
-
[調(diào)劑信息]211智能人工感知方向國(guó)家青年特聘專家課題組招收調(diào)劑研究生
已經(jīng)有0人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- matlab中的計(jì)算問(wèn)題？
已經(jīng)有8人回復(fù)
- MATLAB求助，關(guān)于importdata的問(wèn)題
已經(jīng)有6人回復(fù)
- 用matlab制作的M文件保存在哪里才能直接調(diào)用？
已經(jīng)有9人回復(fù)
- matlab一運(yùn)行就出現(xiàn)的問(wèn)題
已經(jīng)有7人回復(fù)
- matlab的m文件關(guān)聯(lián)問(wèn)題
已經(jīng)有17人回復(fù)
- 求助一個(gè)matlab關(guān)于分類討論的問(wèn)題
已經(jīng)有9人回復(fù)
- matlab編程中遇到的問(wèn)題
已經(jīng)有4人回復(fù)
- 幫忙修改Matlab中出現(xiàn)的錯(cuò)誤
已經(jīng)有4人回復(fù)
- matlab讀取圖像的問(wèn)題
已經(jīng)有12人回復(fù)
- 分享matlab的一些代碼（入門貼）
已經(jīng)有19人回復(fù)
- 求MATLAB清華版的教學(xué)視頻
已經(jīng)有14人回復(fù)
- 關(guān)于Matlab作圖的問(wèn)題
已經(jīng)有14人回復(fù)
- MATLAB初學(xué)者，這樣的式子怎么就是畫不了圖，總是出錯(cuò)!求指導(dǎo)!!!!!!!!!!!!!
已經(jīng)有8人回復(fù)
- matlab程序?qū)D中的數(shù)據(jù)記錄下來(lái)
已經(jīng)有4人回復(fù)
- Matlab在地學(xué)中的應(yīng)用
已經(jīng)有24人回復(fù)
- matlab中dydt的編寫要求
已經(jīng)有12人回復(fù)
- 求助svm在matlab中的應(yīng)用
已經(jīng)有7人回復(fù)
- 關(guān)于matlab 2010b 的一個(gè)小問(wèn)題
已經(jīng)有7人回復(fù)
- MATLAB安裝后出現(xiàn)一堆java的問(wèn)題
已經(jīng)有12人回復(fù)
- 【求助】關(guān)于matlab的高維數(shù)組
已經(jīng)有6人回復(fù)
805592836_01
鐵蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 823
- 帖子: 10
- 在線: 23.4小時(shí)
- 蟲號(hào): 2393210
- 注冊(cè): 2013-04-01
- 性別: MM
- 專業(yè): 能源化工
3樓2014-03-12 10:02:43
dingd
鐵桿木蟲 (職業(yè)作家)
- 計(jì)算強(qiáng)帖: 4
- 應(yīng)助: 1641 (講師)
- 金幣: 15037.3
- 散金: 101
- 紅花: 234
- 帖子: 3410
- 在線: 1223.7小時(shí)
- 蟲號(hào): 291104
- 注冊(cè): 2006-10-28
2樓2014-03-12 09:56:13
dingd
鐵桿木蟲 (職業(yè)作家)
- 計(jì)算強(qiáng)帖: 4
- 應(yīng)助: 1641 (講師)
- 金幣: 15037.3
- 散金: 101
- 紅花: 234
- 帖子: 3410
- 在線: 1223.7小時(shí)
- 蟲號(hào): 291104
- 注冊(cè): 2006-10-28
4樓2014-03-12 10:35:53
805592836_01
鐵蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 823
- 帖子: 10
- 在線: 23.4小時(shí)
- 蟲號(hào): 2393210
- 注冊(cè): 2013-04-01
- 性別: MM
- 專業(yè): 能源化工
5樓2014-03-12 11:31:29
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085801 總分275 本科新能源 求調(diào)劑 +16 | bradoner 2026-04-01 | 20/1000 |
|
|---|---|---|---|---|
|
|
[考研] 366求調(diào)劑一志愿東北大學(xué) +7 | 運(yùn)氣來(lái)得若有似?/a> 2026-04-02 | 7/350 |
|
|
[考研] 材料求調(diào)劑一志愿哈工大324 +10 | 閆旭東 2026-03-28 | 12/600 |
|
|
|
[考研] 266分，一志愿電氣工程，本科材料，求材料專業(yè)調(diào)劑 +4 | 哇呼哼呼哼 2026-04-02 | 4/200 |
|
|
|
[考研] 327求調(diào)劑 +9 | 汲亦昊 2026-03-29 | 9/450 |
|
|
|
[考研] 材料專碩322分 +11 | 哈哈哈吼吼吼哈 2026-04-01 | 11/550 |
|
|
|
[考研] 085600，材料與化工321分調(diào)劑 +8 | 大饞小子 2026-03-27 | 10/500 |
|
|
|
[考研] 考研調(diào)劑 +12 | Amber00 2026-03-31 | 12/600 |
|
|
|
[考研] 343求調(diào)劑085601 +4 | 要努力學(xué)習(xí)x 2026-03-29 | 5/250 |
|
|
|
[考研] 化學(xué)工程專碩324分，一志愿中國(guó)礦業(yè)大學(xué)求調(diào)劑 +7 | 耿耿1314 2026-04-01 | 7/350 |
|
|
|
[考研] 349求調(diào)劑 +6 | 吃的不少 2026-04-01 | 6/300 |
|
|
|
[考研] 309求調(diào)劑 +19 | 誰(shuí)不是少年 2026-03-29 | 19/950 |
|
|
|
[考研] 288資源與環(huán)境專碩求調(diào)劑，不限專業(yè)，有學(xué)上就行 +25 | lllllos 2026-03-30 | 26/1300 |
|
|
|
[考研] 【調(diào)劑】一志愿廈大生物與醫(yī)藥調(diào)劑 +3 | Echo蝦米 2026-03-31 | 3/150 |
|
|
|
[考研] 土木304求調(diào)劑 +5 | 頂級(jí)擦擦 2026-03-31 | 5/250 |
|
|
|
[考研] 材料與化工272求調(diào)劑 +25 | 阿斯蒂芬2004 2026-03-28 | 25/1250 |
|
|
|
[考研] 279求調(diào)劑 +12 | j的立方 2026-03-29 | 12/600 |
|
|
|
[考研] 一志愿南開大學(xué)0710生物學(xué)359求調(diào)劑 +5 | 兔兔兔111223314 2026-03-29 | 7/350 |
|
|
|
[考研] 332求調(diào)劑 +6 | @MZB382400 2026-03-28 | 6/300 |
|
|
|
[考研] 一志愿北京工業(yè)大學(xué)，324分求調(diào)劑 +6 | 零八# 2026-03-28 | 6/300 |
|