 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1380 | 回復(fù): 10 | ||||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||||
1135725495鐵桿木蟲 (著名寫手)
|
[求助]
求大家?guī)臀铱纯聪旅嬷衜atlab中的程序問題出在哪里？ 已有1人參與
|
|||
|
求動力學(xué)參數(shù) 方程為：r1 =(k(1)*k(2)*c^2-0.25*k(3)*k(4)*(9.404-c)^2)/(1+k(2)*c+0.5*k(4)*(9.404-c)); t c 0 8.769713879 3 5.909461709 6 5.736047766 9 4.756936544 12 4.505799773 15 4.354160734 18 4.124025932 21 3.994725485 31 3.732949397 41 3.680613046 51 3.537445782 61 3.283577643 求k1,k2,k3,,k4 下面是寫的程序，大家?guī)兔纯�，為什么結(jié)果是直線呢？ function parafit % % r1 =(k(1)*k(2)*c^2-0.25*k(3)*k(4)*(9.404-c)^2)/(1+k(2)*c+0.5*k(4)*(9.404-c)); % % dCAdt = - r1; clear all clc % t/min CA / mol/L Kinetics=[0 7.029559717 1 3.960706476 2 3.163862085 3 2.963057756 6 2.820434159 9 2.698109429 12 2.556212898 15 2.447100967 18 2.207246045 21 2.13921916 31 2.000403189 41 1.917172016 51 1.851299168 61 1.636287388]; c=Kinetics( :,1)'; k0 = [0.671 0.02253 0.641 0.04028]; % 參數(shù)初值 lb = [0 0 0 0]; % 參數(shù)下限 ub = [100 100 100 100]; % 參數(shù)上限 x0 = [7.029559717]; yexp = Kinetics; warning off % 使用函數(shù)fmincon()進(jìn)行參數(shù)估計(jì) [k,fval,flag] = fmincon(@ObjFunc7Fmincon,k0,[],[],[],[],lb,ub,[],[],x0,yexp); fprintf('\n使用函數(shù)fmincon()估計(jì)得到的參數(shù)值為:\n') fprintf('\tk1 = %.11f\n',k(1)) fprintf('\tk2 = %.11f\n',k(2)) fprintf('\tk3 = %.11f\n',k(3)) fprintf('\tk4 = %.11f\n',k(4)) fprintf(' The sum of the squares is: %.1e\n\n',fval) k_fm= k; warning off % 使用函數(shù)lsqnonlin()進(jìn)行參數(shù)估計(jì) [k,resnorm,residual,exitflag,output,lambda,jacobian] = ... lsqnonlin(@ObjFunc7LNL,k0,lb,ub,[],x0,yexp); ci = nlparci(k,residual,jacobian); fprintf('\n\n使用函數(shù)lsqnonlin()估計(jì)得到的參數(shù)值為:\n') fprintf('\tk1 = %.11f\n',k(1)) fprintf('\tk2 = %.11f\n',k(2)) fprintf('\tk3 = %.11f\n',k(3)) fprintf('\tk4 = %.11f\n',k(4)) fprintf(' The sum of the squares is: %.1e\n\n',resnorm) k_ls = k; output warning off % 以函數(shù)fmincon()估計(jì)得到的結(jié)果為初值，使用函數(shù)lsqnonlin()進(jìn)行參數(shù)估計(jì) k0 = k_fm; [k,resnorm,residual,exitflag,output,lambda,jacobian] = ... lsqnonlin(@ObjFunc7LNL,k0,lb,ub,[],x0,yexp); ci = nlparci(k,residual,jacobian); fprintf('\n\n以fmincon()的結(jié)果為初值，使用函數(shù)lsqnonlin()估計(jì)得到的參數(shù)值為:\n') fprintf('\tk1 = %.11f\n',k(1)) fprintf('\tk2 = %.11f\n',k(2)) fprintf('\tk3 = %.11f\n',k(3)) fprintf('\tk4 = %.11f\n',k(4)) fprintf(' The sum of the squares is: %.1e\n\n',resnorm) k_fmls = k; output tspan = [0 1 2 3 6 9 12 15 18 21 31 41 51 61]; [t x] = ode45(@KineticEqs,tspan,x0,[],k_fmls); figure; plot(t,x(:,1),t,yexp(:,2),'*');legend('ca-pr','ca-real') x(:,1) figure;plot(t,x); hold on % ------------------------------------------------------------------ function f = ObjFunc7Fmincon(k,x0,yexp) tspan = [0 1 2 3 6 9 12 15 18 21 31 41 51 61]; [t x] = ode45(@KineticEqs,tspan,x0,[],k); y(:,2) = x(:,1); f = sum((y(:,2)-yexp(:,2)).^2) ; % ------------------------------------------------------------------ function f = ObjFunc7LNL(k,x0,yexp) tspan = [0 1 2 3 6 9 12 15 18 21 31 41 51 61]; [t x] = ode45(@KineticEqs,tspan,x0,[],k); y(:,2) = x(:,1); f1 = y(:,2) - yexp(:,2); f = [f1]; % ------------------------------------------------------------------ function dxdt = KineticEqs(t,x,k) Kinetics=[0 7.029559717 1 3.960706476 2 3.163862085 3 2.963057756 6 2.820434159 9 2.698109429 12 2.556212898 15 2.447100967 18 2.207246045 21 2.13921916 31 2.000403189 41 1.917172016 51 1.851299168 61 1.636287388]; c=Kinetics(:,1)'; for i=1:14 c(i)=c(:,i); r1 =(k(1)*k(2)*c(i)^2-0.25*k(3)*k(4)*(9.404-c(i))^2)/(1+k(2)*c(i)+0.5*k(4)*(9.404-c(i))); dCAdt = - r1; dxdt = [dCAdt]; end 結(jié)果如下： H8_DJFZ~6NZC}07J0EXL1]L.png |
回復(fù)此樓» 猜你喜歡
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有0人回復(fù)
-
中德博士后2026項(xiàng)目5月30日截止- 國內(nèi)博士后入站、德國兩年全職
已經(jīng)有0人回復(fù)
-
物理學(xué)I論文潤色/翻譯怎么收費(fèi)?
已經(jīng)有284人回復(fù)
-
面向全球申請人的中德博士后2026項(xiàng)目5月30日截止- 國內(nèi)博后入站德國兩年全職
已經(jīng)有0人回復(fù)
-
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有2人回復(fù)
-
重慶交大材料學(xué)院26年碩士招生復(fù)試實(shí)施細(xì)則已出。歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。
已經(jīng)有0人回復(fù)
-
北工大宋曉艷教授團(tuán)隊(duì)招收2026年9月入學(xué)博士【學(xué)博和專博名額均有】！
已經(jīng)有1人回復(fù)
-
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收10人
已經(jīng)有0人回復(fù)
-
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收調(diào)劑，10人
已經(jīng)有0人回復(fù)
-
計(jì)算與實(shí)驗(yàn)科學(xué)工程國際會議(8月7-11日,香港) 二維材料分會S5與聲學(xué)換能器分會S8
已經(jīng)有0人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對您同樣有幫助:
- 求助MATLAB中程序看不懂的簡單問題！謝謝！
已經(jīng)有3人回復(fù)
- 如何求這個(gè)matlab程序畫出的legend曲線的法矢量？
已經(jīng)有6人回復(fù)
- 高分求將我的MATLAB程序改寫為FORTRAN程序（1000BB）
已經(jīng)有13人回復(fù)
- 求在matlab中進(jìn)行分類的程序
已經(jīng)有5人回復(fù)
- 用matlab畫一個(gè)函數(shù)，能幫我看看程序哪里有問題嗎
已經(jīng)有8人回復(fù)
- 求教如何用matlab繪制下面兩圖
已經(jīng)有5人回復(fù)
- 求大神用matlab給我處理一組數(shù)據(jù)，有人沒有啊？
已經(jīng)有15人回復(fù)
- matlab一運(yùn)行就出現(xiàn)的問題
已經(jīng)有7人回復(fù)
- 有個(gè)matlab程序，但是運(yùn)行不下來，求高手
已經(jīng)有23人回復(fù)
- matlab程序中調(diào)用m文件的問題
已經(jīng)有13人回復(fù)
- 幫忙詳解示范一下MATLAB中g(shù)input函數(shù)的應(yīng)用吧
已經(jīng)有3人回復(fù)
- 實(shí)在是看不出這段matlab程序里哪里出錯了，請大家?guī)兔纯矗?/a>
已經(jīng)有10人回復(fù)
- 請教MATLAB中的LMI求解問題，急�。。�
已經(jīng)有7人回復(fù)
- 求高手解答Matlab求解微分方程的初始值和邊界值問題
已經(jīng)有5人回復(fù)
- 求大神幫我改一下下面的程序！matlab迭代法的 不知道錯在哪里！
已經(jīng)有7人回復(fù)
- matlab中遇到的畫圖小問題，程序沒問題。標(biāo)題要長@@@@@@@@@@@@@@@@@@@@@@
已經(jīng)有8人回復(fù)
- 大家愛幫我看看這段MATLAB程序?yàn)樯冻鲥e了
已經(jīng)有5人回復(fù)
- matlab程序問題，總是出錯，麻煩大家看看問題在那里？
已經(jīng)有6人回復(fù)
- 【求助】幫我編個(gè)程啊~C語言，C++，Matlab，VB...都行~
已經(jīng)有14人回復(fù)
1135725495
鐵桿木蟲 (著名寫手)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 6868.7
- 散金: 18
- 紅花: 1
- 帖子: 1492
- 在線: 209.6小時(shí)
- 蟲號: 2756771
- 注冊: 2013-10-27
- 性別: MM
- 專業(yè): 有機(jī)化工
9樓2015-05-26 17:09:52
月只藍(lán)
主管區(qū)長 (職業(yè)作家)
-
專家經(jīng)驗(yàn): +1059 - 計(jì)算強(qiáng)帖: 8
- 應(yīng)助: 1712 (講師)
- 貴賓: 8.888
- 金幣: 68128.7
- 散金: 1938
- 紅花: 443
- 沙發(fā): 4
- 帖子: 4373
- 在線: 3291.4小時(shí)
- 蟲號: 1122189
- 注冊: 2010-10-14
- 專業(yè): 宇宙學(xué)
- 管轄: 計(jì)算模擬區(qū)
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
1135725495: 金幣+10, ★有幫助 2015-05-26 15:35:53
感謝參與，應(yīng)助指數(shù) +1
1135725495: 金幣+10, ★有幫助 2015-05-26 15:35:53
|
兩組參考結(jié)果： Correlation Coef. (R): 0.986609345793065 R-Square: 0.97339800120622 Determination Coef. (DC): 0.973395608126228 F-Statistic: 73.3952829093829 Parameters Best Estimate -------------------- ------------- k1 0.0633244071641639 k2 -0.123667552854961 k3 0.0140433904054541 k4 -0.681989563148558 ====== Output Results ===== File: Data file - 1 No Obs. c Cal. c 1 5.909461709 6.03887691479301 2 5.736047766 5.36051185683508 3 4.756936544 4.91507500618295 4 4.505799773 4.59634816701972 5 4.354160734 4.35722735701414 6 4.124025932 4.17236164824689 7 3.994725485 4.02646466838814 8 3.732949397 3.71546906582311 9 3.680613046 3.55315135651946 10 3.537445782 3.46435427377623 11 3.283577643 3.41460522937653 Correlation Coef. (R): 0.98660934579373 R-Square: 0.973398001207533 Determination Coef. (DC): 0.973395608126228 F-Statistic: 73.3952829093831 Parameters Best Estimate -------------------- ------------- k1 0.481614078106072 k2 1.61663483596038 k3 0.0205007668163867 k4 46.4476555663092 ====== Output Results ===== File: Data file - 1 No Obs. c Cal. c 1 5.909461709 6.03887691237705 2 5.736047766 5.36051185622609 3 4.756936544 4.91507500636805 4 4.505799773 4.5963481675709 5 4.354160734 4.35722735771124 6 4.124025932 4.17236164896553 7 3.994725485 4.02646466905497 8 3.732949397 3.71546906608907 9 3.680613046 3.55315135632873 10 3.537445782 3.46435427319561 11 3.283577643 3.41460522849613 |
2樓2015-05-26 13:56:42
1135725495
鐵桿木蟲 (著名寫手)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 6868.7
- 散金: 18
- 紅花: 1
- 帖子: 1492
- 在線: 209.6小時(shí)
- 蟲號: 2756771
- 注冊: 2013-10-27
- 性別: MM
- 專業(yè): 有機(jī)化工
您好，您的計(jì)算結(jié)果我看到了，我想說的是我上面編的Matlab程序正確嗎？為什么會出現(xiàn)那樣的計(jì)算結(jié)果？還有就是您計(jì)算的方法是最小二乘嗎？我現(xiàn)在需要用matlab計(jì)算，因?yàn)閷懻撐男枰�。所以您能不能幫忙看一下程序是否正確？ 3樓2015-05-26 15:38:21
月只藍(lán)
主管區(qū)長 (職業(yè)作家)
-
專家經(jīng)驗(yàn): +1059 - 計(jì)算強(qiáng)帖: 8
- 應(yīng)助: 1712 (講師)
- 貴賓: 8.888
- 金幣: 68128.7
- 散金: 1938
- 紅花: 443
- 沙發(fā): 4
- 帖子: 4373
- 在線: 3291.4小時(shí)
- 蟲號: 1122189
- 注冊: 2010-10-14
- 專業(yè): 宇宙學(xué)
- 管轄: 計(jì)算模擬區(qū)
4樓2015-05-26 16:25:47
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿：西北大學(xué)，英一數(shù)一408-284分求調(diào)劑 +4 | 12.27 2026-03-27 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 一志愿北京化工大學(xué)材料與化工（085600）296求調(diào)劑 +14 | 稻妻小編 2026-03-26 | 14/700 |
|
|
|
[考研]
 一志愿武漢理工，總分321，英一數(shù)二，求老師收留。
+10
|
nnnnnnn5 2026-03-25 | 10/500 |
|
|
|
[考研] 本科新能源科學(xué)與工程，一志愿華理能動285求調(diào)劑 +7 | AZMK 2026-03-28 | 11/550 |
|
|
|
[考研] 本科新能源科學(xué)與工程，一志愿華理能動285求調(diào)劑 +3 | AZMK 2026-03-27 | 5/250 |
|
|
|
[考研]
一志愿廈門大學(xué)化學(xué)學(xué)碩307求調(diào)劑
+10
|
y7czhao 2026-03-26 | 10/500 |
|
|
|
[考研] 299求調(diào)劑 +7 | 嗯嗯嗯嗯2 2026-03-27 | 7/350 |
|
|
|
[考研] 085600 286分 材料求調(diào)劑 +7 | 麻辣魷魚 2026-03-27 | 8/400 |
|
|
|
[考研] 材料求調(diào)劑 一志愿哈工大總分298分，前三科223分 +5 | dongfang59 2026-03-27 | 5/250 |
|
|
|
[考研] 265求調(diào)劑11408 +3 | 劉小鹿lu 2026-03-27 | 3/150 |
|
|
|
[考研] 一志愿211院校 344分 東北農(nóng)業(yè)大學(xué)生物學(xué)學(xué)碩，求調(diào)劑 +5 | 丶風(fēng)雪夜歸人丶 2026-03-26 | 8/400 |
|
|
|
[考研] 0856調(diào)劑 +5 | 求求讓我有書讀?/a> 2026-03-26 | 6/300 |
|
|
|
[考研] 085600，材料與化工321分調(diào)劑 +4 | 大饞小子 2026-03-27 | 6/300 |
|
|
|
[考研] 305求調(diào)劑 +5 | 哇盧卡庫 2026-03-26 | 5/250 |
|
|
|
[考研] 調(diào)劑推薦 +5 | 清酒714 2026-03-26 | 6/300 |
|
|
|
[碩博家園] 北京林業(yè)大學(xué)碩導(dǎo)招生廣告 +6 | kongweilin 2026-03-26 | 8/400 |
|
|
|
[考研] 341求調(diào)劑 +7 | 青檸檬1 2026-03-26 | 7/350 |
|
|
|
[考研] 0703化學(xué)求調(diào)劑 +3 | 丹青奶蓋 2026-03-26 | 5/250 |
|
|
|
[考研] 一志愿南航 335分 | 0856材料化工 | GPA 4.07 | 有科研經(jīng)歷 +6 | cccchenso 2026-03-23 | 6/300 |
|
|
|
[考研] 335求調(diào)劑 +4 | yuyu宇 2026-03-23 | 5/250 |
|