 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 3 | 1/1 | 返回列表 |
| 查看: 349 | 回復(fù): 2 | ||
清輝12365新蟲 (初入文壇)
|
[求助]
弄動力學(xué)參數(shù)擬合參數(shù)，不知道怎么辦了，哪個大神幫忙弄下代碼
|
| 下面是我上傳的問題和數(shù)據(jù)，看了幾個月的實用化工計算機模擬沒有進展，請高手幫忙一下， |
回復(fù)此樓» 本帖附件資源列表
-
歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺，不對該內(nèi)容負責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識產(chǎn)權(quán)問題，其責(zé)任在于用戶本人，如對版權(quán)有異議，請聯(lián)系郵箱：xiaomuchong@tal.com - 附件 1 : 問題.docx
- 附件 2 : 數(shù)據(jù).xlsx
2015-03-02 08:46:05, 19.53 K
2015-03-02 08:46:11, 10.7 K
» 猜你喜歡
求調(diào)劑
已經(jīng)有3人回復(fù)
-
265求調(diào)劑
已經(jīng)有4人回復(fù)
-
085700資源與環(huán)境308求調(diào)劑
已經(jīng)有6人回復(fù)
-
一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑
已經(jīng)有12人回復(fù)
-
286分人工智能專業(yè)請求調(diào)劑愿意跨考！
已經(jīng)有3人回復(fù)
-
329求調(diào)劑
已經(jīng)有5人回復(fù)
-
申請回稿延期一個月，編輯同意了。但系統(tǒng)上的時間沒變，給編輯又寫郵件了，沒回復(fù)
已經(jīng)有4人回復(fù)
-
材料學(xué)碩318求調(diào)劑
已經(jīng)有5人回復(fù)
-
一志愿中國海洋大學(xué)，生物學(xué)，301分，求調(diào)劑
已經(jīng)有6人回復(fù)
-
081700化工學(xué)碩調(diào)劑
已經(jīng)有3人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 這樣的動力學(xué)方程能夠有1stopt擬合出參數(shù)來么
已經(jīng)有4人回復(fù)
- 使用MATLAB擬合動力學(xué)參數(shù)問題
已經(jīng)有18人回復(fù)
- 1stopt數(shù)據(jù)擬合求助
已經(jīng)有8人回復(fù)
- 跪求高手幫忙看看這組數(shù)據(jù)如何進行擬合
已經(jīng)有7人回復(fù)
- 使用matlab最優(yōu)化方法擬合獲得多個動力學(xué)參數(shù)中的問題
已經(jīng)有4人回復(fù)
- 求助1stopt擬合動力學(xué)參數(shù)
已經(jīng)有4人回復(fù)
- 求大神幫忙擬合一個非線性方程，求出模型參數(shù)
已經(jīng)有15人回復(fù)
- 拜求！matlab擬合動力學(xué)參數(shù)
已經(jīng)有14人回復(fù)
- 關(guān)于兩個參數(shù)的曲線擬合問題，向大神們求助
已經(jīng)有3人回復(fù)
- 求助用matlab擬合動力學(xué)方程
已經(jīng)有9人回復(fù)
- 使用lsqnonlin函數(shù)優(yōu)化動力學(xué)參數(shù)，總是得不到合理的結(jié)果
已經(jīng)有6人回復(fù)
- 求大家?guī)蛶兔纯磎atlab擬合動力學(xué)方程的問題，已經(jīng)走投無路了，求大家?guī)蛶兔Π。?/a>
已經(jīng)有19人回復(fù)
- 怎樣用一組參數(shù)同時擬合兩個曲線--matlab
已經(jīng)有5人回復(fù)
- MATLAB微分方程參數(shù)擬合問題，求大神
已經(jīng)有7人回復(fù)
- 急求高手幫寫個matlab程序擬合參數(shù)求動力學(xué)方程。
已經(jīng)有4人回復(fù)
- 動力學(xué)參數(shù)擬合
已經(jīng)有26人回復(fù)
- matlab 擬合反應(yīng)動力學(xué)參數(shù)結(jié)果很差。大家?guī)兔匆幌?/a>
已經(jīng)有14人回復(fù)
- 【求助】用matlab最優(yōu)化方法進行參數(shù)擬合
已經(jīng)有17人回復(fù)
清輝12365
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 23
- 帖子: 21
- 在線: 53.2小時
- 蟲號: 2804052
- 注冊: 2013-11-15
- 專業(yè): 化學(xué)反應(yīng)工程
2樓2015-03-02 09:14:22
清輝12365
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 23
- 帖子: 21
- 在線: 53.2小時
- 蟲號: 2804052
- 注冊: 2013-11-15
- 專業(yè): 化學(xué)反應(yīng)工程
|
1-初始量，2-溫度，3-壓力，4-甲醇摩爾分數(shù)x1，5-CO摩爾分數(shù)x2，6-O2摩爾分數(shù)X3，8-產(chǎn)物DMC摩爾數(shù)y1 9-產(chǎn)物CO2摩爾分數(shù)y2 上面的1、2、3、4、5、6、8、9在Donglixue.mat中 我的編程是：function yyh clear all;clc load Donglixue; %輸入數(shù)據(jù) NO=Donglixue(:,1); %輸入初始量 ym0=Donglixue(:,4); %輸入甲醇摩爾比 yco0=Donglixue(:,5); %co的摩爾比 yo2=Donglixue(:,6); %o2的摩爾比 ydmc=Donglixue(:,8); %實驗的DMC的摩爾比 yco2=Donglixue(:,9); %實驗的CO2的摩爾比 t=Donglixue(:,2); %輸入溫度 T=t+273.15; %溫度換成K P=Donglixue(:,2); %輸入壓力 %lsqnonlin函數(shù)擬合 beta0=[10000 500 0 0 0 10000 10000 0.2 -0.3]; lb=[0 0 -2 -2 -2 0 0 -2 -2]; %參數(shù)的下限 Y0=[0 0]; tspan=[0,0.375]; %積分限 [beta,resonrm,residual,exitflag,output,lambda,jacobian]=... lsqnonlin(@fun1,beta0,lb,[],[],tspan,Y0,NO,ym0,yco0,yo2,ydmc,yco2,T,P) function f=fun1(beta,tspan,Y0,NO,ym0,yco0,yo2,ydmc,yco2,T,P) [t,Y]=ode45(@fun2,tspan,Y0,[],beta); f=[Y(:,1)-ydmc;Y(:,2)-yco2]; function ff=fun2(t,Y,NO,ym0,yco0,yo2,ydmc,yco2,T,P) ff(1)=((1+1.5*Y(:,1)+0.5*Y(:,2)).*(1+1.5*Y(:,1))*... (ym0+Y(:,1).*(1.5*ym0-2)+0.5*ym0.*Y(:,2)).^beta(3)... .*P.^(beta(3)+beta(4)+beta(5)).*(yco0+Y(:,1).*(1.5*yco0-1)+Y(:,2).*(0.5*yco0-1)).^beta(4).*... (yo2+0.5.*Y(:,1)*(3*yco0-1)+0.5*Y(:,2)*(yo2-1)).^beta(5)... .*beta(1)*exp(-beta(2)./T)) +0.5*Y(:,1).*beta(6)*exp(-beta(7)./T).*P.^(beta(8)+beta(9)).*... (yco0+Y(:,1).*(1.5*yco0-1)+Y(:,2).*(0.5*yco0-1)).^beta(8).*... (yo2+0.5.*Y(:,1)*(3*yco0-1)+0.5*Y(:,2)*(yo2-1)).^beta(9)))./NO ff(2)=(1+1.5*Y(:,1)+0.5*Y(:,2)).* ((0.5*Y(:,2)+1).*beta(6)*exp(-beta(7)./T).*P.^(beta(8)+beta(9)).*... (yco0+Y(:,1).*(1.5*yco0-1)+Y(:,2).*(0.5*yco0-1)).^beta(8).*... (yo2+0.5.*Y(:,1)*(3*yco0-1)+0.5*Y(:,2)*(yo2-1)).^beta(9)+... 1.5*Y(:,2)*... (ym0+Y(:,1).*(1.5*ym0-2)+0.5*ym0.*Y(:,2)).^beta(3).*... (yco0+Y(:,1).*(1.5*yco0-1)+Y(:,2).*(0.5*yco0-1)).^beta(4).*... (yo2+0.5.*Y(:,1)*(3*yco0-1)+0.5*Y(:,2)*(yo2-1)).^beta(5).*P.^(beta(3)+beta(4)+beta(5).*... beta(1)*exp(-beta(2)./T)./NO 因為r1、r2式子比較麻煩，所以看的式子比較亂，請大神們指導(dǎo)一下， 1.jpg 2.jpg |
» 本帖附件資源列表
-
歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺，不對該內(nèi)容負責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識產(chǎn)權(quán)問題，其責(zé)任在于用戶本人，如對版權(quán)有異議，請聯(lián)系郵箱：xiaomuchong@tal.com - 附件 1 : Donglixue.mat
2015-03-02 16:23:20, 1 K
3樓2015-03-02 16:23:24
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 清輝12365 的主題更新
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料學(xué)碩318求調(diào)劑 +5 | February_Feb 2026-03-19 | 5/250 |
|
|---|---|---|---|---|
|
[考研]
 320求調(diào)劑0856
+3
|
不想起名字112 2026-03-19 | 3/150 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 ，六級已過，有科研經(jīng)歷 +12 | 曦熙兮 2026-03-15 | 12/600 |
|
|
|
[考研] 266求調(diào)劑 +5 | 陽陽哇塞 2026-03-14 | 10/500 |
|
|
|
[考研] 304求調(diào)劑 +3 | 曼殊2266 2026-03-18 | 3/150 |
|
|
|
[考研] 求調(diào)劑 +3 | Mqqqqqq 2026-03-19 | 3/150 |
|
|
|
[考研] 一志愿天大材料與化工（085600）總分338 +5 | 蔡大美女 2026-03-13 | 5/250 |
|
|
|
[考研] 材料080500調(diào)劑求收留 +4 | 一顆meteor 2026-03-13 | 4/200 |
|
|
|
[考研] 一志愿武理材料305分求調(diào)劑 +5 | 想上岸的鯉魚 2026-03-18 | 6/300 |
|
|
|
[考研] 085601專碩，總分342求調(diào)劑，地區(qū)不限 +5 | share_joy 2026-03-16 | 5/250 |
|
|
|
[考研] 299求調(diào)劑 +5 | △小透明* 2026-03-17 | 5/250 |
|
|
|
[考研] 085601求調(diào)劑 +4 | Du.11 2026-03-16 | 4/200 |
|
|
|
[考研] 一志愿蘇州大學(xué)材料工程（085601）專碩有科研經(jīng)歷三項國獎兩個實用型專利一項省級立項 +6 | 大火山小火山 2026-03-16 | 8/400 |
|
|
|
[考研] 302求調(diào)劑 +4 | 小賈同學(xué)123 2026-03-15 | 8/400 |
|
|
|
[論文投稿] 有沒有大佬發(fā)小論文能帶我個二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|
|
|
[考研] 藥學(xué)383 求調(diào)劑 +3 | 藥學(xué)chy 2026-03-15 | 4/200 |
|
|
|
[考研] 304求調(diào)劑 +4 | ahbd 2026-03-14 | 4/200 |
|
|
|
[考研] 326求調(diào)劑 +3 | mlpqaz03 2026-03-15 | 3/150 |
|
|
|
[考研] 26考研一志愿中國石油大學(xué)(華東)305分求調(diào)劑 +3 | 嘉年新程 2026-03-15 | 3/150 |
|
|
|
[考研] 288求調(diào)劑 +4 | 奇點0314 2026-03-14 | 4/200 |
|