版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

返回列表

lxyy

銀蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 237.1
散金: 480
紅花: 9
帖子: 182
在線: 210.5小時(shí)
蟲號(hào): 1131417
注冊(cè): 2010-10-25
專業(yè): 化學(xué)反應(yīng)工程

[求助] matlab擬合擬均相動(dòng)力學(xué)參數(shù)

模型如附件所示，要求里面的k+,k-,m,n,p，我自己編了個(gè)程序，但是一直有問(wèn)題
M文件
function f=fortyfour_D(c,t,rD)
t=[5,10,15,20,30,45,60,80,100,120,150,180,210,240,300,360];
CD=[749.16 ,994.49 ,1382.00,1589.18,1909.81,2222.86,2525.68,2721.90,2896.67,2951.42,2960.67,2999.60,3017.97,3048.91,3072.27,3063.85];
knots = 3;K = 3;    %三次B樣條
sp=spap2(knots,K,t,CD);
sp=spap2(newknt(sp),K,t,CD);
pp=fnder(sp);       % 計(jì)算B樣條函數(shù)的導(dǎo)函數(shù)
dCDdt=fnval(pp,t);    % 計(jì)算t處的導(dǎo)函數(shù)值
Wcat=85.6312 %催化劑濃度
rD=dCDdt./Wcat;
f=c(1)*CA.^c(2)*CB.^c(3)-c(4)*CD.^c(5)))-rD;
命令如下
c0=[1,1,1,1,1];
for i=1:100
c=lsqnonlin('fortyfour_D ',c0);
c0=c;
end
c
% 繪制圖形
ti = linspace(t(1),t(end),200); %y=linspace（a，b，n）生成一個(gè)行向量，該向量將a與b之間平分為n個(gè)點(diǎn)，包含端點(diǎn)a和b。
CDi = fnval(sp,ti); %計(jì)算樣條函數(shù)區(qū)間內(nèi)任意一點(diǎn)的值
plot(t,CD,'ro',ti,CDi,'b-')，xlabel('t'),ylabel('C_D')
legend('實(shí)驗(yàn)數(shù)據(jù)','擬合曲線')
c=[1,1,1,1,1];for i=1:100;c=lsqnonlin('fortyfour_D ',c);c;end
CD=[749.16 ,994.49 ,1382.00,1589.18,1909.81,2222.86,2525.68,2721.90,2896.67,2951.42,2960.67,2999.60,3017.97,3048.91,3072.27,3063.85];
CA = [2732.43 2468.85 2223.97 2023.64 1696.41 1343.42 1054.34 902.20 756.91 717.86 703.53 631.00 652.41 649.53 624.92 744.95];
CB= [3127.04 2722.48 2414.51 2267.55 1927.70 1500.17 1285.07 1111.91 924.27 886.32 830.86 821.82 801.97 780.39 794.41 783.73];
rD=[0.5631,0.5282,0.4934,0.4585,0.3888,0.2842,0.1796,0.0482,0.0404,0.0326,0.0209,0.0092,0.0046,0.0040,0.0028,0.0015];
>> c0=[1,1,1,1,1];
for i=1:1000
c=lsqcurvefit('ffD ',c0,CD,rD);
c0=c;
end
c

請(qǐng)大家?guī)兔纯次页绦蚰睦镉袉?wèn)題[ Last edited by lxyy on 2011-12-30 at 16:05 ]

回復(fù)此樓

» 本帖附件資源列表

歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺(tái)，不對(duì)該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識(shí)產(chǎn)權(quán)問(wèn)題，其責(zé)任在于用戶本人，如對(duì)版權(quán)有異議，請(qǐng)聯(lián)系郵箱：xiaomuchong@tal.com
附件 1 : clip_image002.gif

2011-12-30 16:05:36, 503 bytes

» 猜你喜歡

2026年第五屆電氣、電子與信息工程國(guó)際會(huì)議（ISEEIE 2026）已經(jīng)有0人回復(fù)
哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑已經(jīng)有6人回復(fù)
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有272人回復(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)還有沒有用Latex寫文章的小伙伴們？已經(jīng)有0人回復(fù)
光學(xué)工程學(xué)碩調(diào)劑信息已經(jīng)有30人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
散金幣，求好運(yùn)，祝面上順利！已經(jīng)有65人回復(fù)
華東師范大學(xué)芯片設(shè)計(jì)徐瓏真實(shí)水平咋樣已經(jīng)有1人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

matlab擬合動(dòng)力學(xué)參數(shù) 已經(jīng)有4人回復(fù)
非線性擬合估計(jì)方程參數(shù)（matlab，spss，origin都可以）已經(jīng)有4人回復(fù)
matlab非線性參數(shù)擬合問(wèn)題已經(jīng)有7人回復(fù)
請(qǐng)教matlab擬合已經(jīng)有9人回復(fù)
求助：化學(xué)反應(yīng)動(dòng)力學(xué)數(shù)據(jù)擬合得到Arrenius方程參數(shù)（A，belta，Ea）已經(jīng)有4人回復(fù)
動(dòng)力學(xué)數(shù)據(jù)擬合已經(jīng)有5人回復(fù)
動(dòng)力學(xué)方程擬合參數(shù)求解已經(jīng)有4人回復(fù)
【求助】用matlab最優(yōu)化方法進(jìn)行參數(shù)擬合已經(jīng)有17人回復(fù)
求助matlab擬合參數(shù)的問(wèn)題已經(jīng)有2人回復(fù)
【求助】使用MATLAB擬合出按照公式的系數(shù) 已經(jīng)有19人回復(fù)
【求助】使用Matlab擬合反應(yīng)動(dòng)力學(xué)方程問(wèn)題已經(jīng)有7人回復(fù)

1樓 2011-12-30 16:02:07

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

dingd

鐵桿木蟲 (職業(yè)作家)

計(jì)算強(qiáng)帖: 4
應(yīng)助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時(shí)
蟲號(hào): 291104
注冊(cè): 2006-10-28

★
dbb627(金幣+1): 歡迎交流 2011-12-31 08:41:48

不知理解的對(duì)否，用1stOpt試下：

CODE:

Variable CD,CA,CB,rD;

Function rD=k1*Ca^m*Cb^n-k2*Cd^p;

Data;

CD=[749.16,994.49,1382.00,1589.18,1909.81,2222.86,2525.68,2721.90,2896.67,2951.42,2960.67,2999.60,3017.97,3048.91,3072.27,3063.85];

CA=[2732.43,2468.85,2223.97,2023.64,1696.41,1343.42,1054.34,902.20,756.91,717.86,703.53,631.00,652.41,649.53,624.92,744.95];

CB=[3127.04,2722.48,2414.51,2267.55,1927.70,1500.17,1285.07,1111.91,924.27,886.32,830.86,821.82,801.97,780.39,794.41,783.73];

rD=[0.5631,0.5282,0.4934,0.4585,0.3888,0.2842,0.1796,0.0482,0.0404,0.0326,0.0209,0.0092,0.0046,0.0040,0.0028,0.0015];

均方差(RMSE): 0.0187338627500751
殘差平方和(SSE): 0.00561532181661841
相關(guān)系數(shù)(R): 0.996169068786095
相關(guān)系數(shù)之平方(R^2): 0.992352813606156
決定系數(shù)(DC): 0.99235277918993

參數(shù) 最佳估算
---------- -------------
k1 14.3465921423275
m 0.0595653946350665
n -0.452413851540345
k2 1.79201130280021E-8
p 2.22793339858939

贊一下(1人)

回復(fù)此樓

2樓2011-12-30 22:02:01

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

lxyy

銀蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 237.1
散金: 480
紅花: 9
帖子: 182
在線: 210.5小時(shí)
蟲號(hào): 1131417
注冊(cè): 2010-10-25
專業(yè): 化學(xué)反應(yīng)工程

引用回帖:

樓: Originally posted by dingd at 2011-12-30 22:02:01:
不知理解的對(duì)否，用1stOpt試下：
[code]
Variable CD,CA,CB,rD;
Function rD=k1*Ca^m*Cb^n-k2*Cd^p;
Data;
CD=[749.16,994.49,1382.00,1589.18,1909.81,2222.86,2525.68,2721.90,2896.67,2951.42,2960.67,29 ...

作為一個(gè)菜鳥來(lái)說(shuō)我不知道你這樣算對(duì)不對(duì)，但是以我的經(jīng)驗(yàn)來(lái)看mn應(yīng)該接近1，p算出來(lái)應(yīng)該接近2才對(duì)

贊一下

回復(fù)此樓

3樓2011-12-30 22:39:49

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

dingd

鐵桿木蟲 (職業(yè)作家)

計(jì)算強(qiáng)帖: 4
應(yīng)助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時(shí)
蟲號(hào): 291104
注冊(cè): 2006-10-28

【答案】應(yīng)助回帖

感謝參與，應(yīng)助指數(shù) +1
lxyy(金幣+1): ★有幫助 2012-01-05 16:34:35

如果公式和數(shù)據(jù)沒錯(cuò)，參數(shù)也沒有范圍限制的話，2樓給出的結(jié)果從數(shù)學(xué)角度應(yīng)該是最好的了，可以將該結(jié)果和你知道的結(jié)果代入驗(yàn)證一下就知道了。

贊一下

回復(fù)此樓

4樓2011-12-31 08:18:50

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

lxyy

銀蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 237.1
散金: 480
紅花: 9
帖子: 182
在線: 210.5小時(shí)
蟲號(hào): 1131417
注冊(cè): 2010-10-25
專業(yè): 化學(xué)反應(yīng)工程

M文件
function f=q(c,t,CD,rD)
t=[0 0.08333 0.166667 0.25 0.5 0.75 1 1.5 2 2.666667 3.333333 4 5];
CD=[0 0.53259 1.154323 1.621033 2.644088 3.163488 3.461446 3.710468 3.784974 3.746928 3.761759 3.713404 3.643825];
rD=[0.0710291014584411 0.0700187830917774 0.0655437729479848 0.0576582625377526 0.0329826401194250 0.0166083513736405 0.00912074378901541 0.00278932296049053 0.000202497370868135 0 0 0 0]；
f= c(1)*((4.6026-CD).^c(2)*(5.064335-CD).^c(3)-c(4)*CD.^c(5))-rD;
命令如下
c0=[0.1,0.1,0.1,0.1,0.1];
for i=1:50
c=lsqnonlin('q ',c0);
c0=c;
end
c
我把程序簡(jiǎn)化了下，還是一直運(yùn)行不了

贊一下

回復(fù)此樓

5樓2011-12-31 11:56:24

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

bao2007pin

銅蟲 (正式寫手)

應(yīng)助: 2 (幼兒園)
金幣: 2047.6
散金: 207
帖子: 307
在線: 170.5小時(shí)
蟲號(hào): 1119550
注冊(cè): 2010-10-11
性別: GG
專業(yè): 有機(jī)合成

表示我不懂！坐等高人吧�。。。。�！

贊一下

回復(fù)此樓

6樓2011-12-31 12:43:36

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

lxyy

銀蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 237.1
散金: 480
紅花: 9
帖子: 182
在線: 210.5小時(shí)
蟲號(hào): 1131417
注冊(cè): 2010-10-25
專業(yè): 化學(xué)反應(yīng)工程

引用回帖:

樓: Originally posted by bao2007pin at 2011-12-31 12:43:36:
表示我不懂！坐等高人吧�。。。。�！

謝謝

回復(fù)此樓

7樓2011-12-31 14:00:47

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 lxyy 的主題更新

返回列表

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[基金申請(qǐng)] 學(xué)校已經(jīng)提交到NSFC，還能修改嗎？ 40+3	babangida 2026-03-19	7/350	2026-03-20 14:13 by hejicker
[考研] 296求調(diào)劑 +4	www_q 2026-03-18	8/400	2026-03-20 12:12 by 學(xué)員8dgXkO
[考研] 26調(diào)劑/材料/英一數(shù)二/總分289/已過(guò)A區(qū)線 +8	步川酷紫123 2026-03-13	8/400	2026-03-20 08:47 by xingguangj
[考研] 307求調(diào)劑 +9	冷笙123 2026-03-17	9/450	2026-03-19 22:44 by 學(xué)員8dgXkO
[考研] 梁成偉老師課題組歡迎你的加入 +9	一鴨鴨喲 2026-03-14	11/550	2026-03-19 17:22 by ！本暗一次！
[考研] 085600材料與化工求調(diào)劑 +6	緒幸與子 2026-03-17	6/300	2026-03-19 13:27 by houyaoxu
[考研] 287求調(diào)劑 +3	晨昏線與星海 2026-03-19	4/200	2026-03-19 12:32 by peike
[考研] 材料考研調(diào)劑 +3	xwt。 2026-03-19	3/150	2026-03-19 11:22 by w沐陽(yáng)w
[考研] 085600材料與化工 +5	安全上岸！ 2026-03-16	5/250	2026-03-18 15:33 by cmz0325
[考研] 302求調(diào)劑 +10	呼呼呼。。。。 2026-03-17	10/500	2026-03-18 12:45 by Linda Hu
[考研] 301求調(diào)劑 +9	yy要上岸呀 2026-03-17	9/450	2026-03-18 08:58 by 無(wú)際的草原
[考研] 308求調(diào)劑 +4	是Lupa啊 2026-03-16	4/200	2026-03-17 17:12 by ruiyingmiao
[考研] 一志愿蘇州大學(xué)材料工程（085601）專碩有科研經(jīng)歷三項(xiàng)國(guó)獎(jiǎng)兩個(gè)實(shí)用型專利一項(xiàng)省級(jí)立項(xiàng) +6	大火山小火山 2026-03-16	8/400	2026-03-17 15:05 by 無(wú)懈可擊111
[考研] 一志愿，福州大學(xué)材料專碩339分求調(diào)劑 +3	木子momo青爭(zhēng) 2026-03-15	3/150	2026-03-17 07:52 by laoshidan
[考研] 東南大學(xué)364求調(diào)劑 +5	JasonYuiui 2026-03-15	5/250	2026-03-16 21:28 by 木瓜膏
[基金申請(qǐng)] 今年的國(guó)基金是打分制嗎？ 50+3	zhanghaozhu 2026-03-14	3/150	2026-03-16 17:07 by 北京萊茵潤(rùn)色
[考研] 304求調(diào)劑 +3	曼殊2266 2026-03-14	3/150	2026-03-16 16:39 by houyaoxu
[考研] 326求調(diào)劑 +3	mlpqaz03 2026-03-15	3/150	2026-03-16 07:33 by Iveryant
[考研] 085601材料工程315分求調(diào)劑 +3	yang_0104 2026-03-15	3/150	2026-03-15 10:58 by peike
[考研] 288求調(diào)劑 +4	奇點(diǎn)0314 2026-03-14	4/200	2026-03-14 23:04 by JourneyLucky

24小時(shí)熱門版塊排行榜

[求助] matlab擬合擬均相動(dòng)力學(xué)參數(shù)

» 本帖附件資源列表

» 猜你喜歡

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

【答案】應(yīng)助回帖

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助: