 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1398 | 回復(fù): 10 | ||||
| 當(dāng)前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | ||||
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ù)估計 [k,fval,flag] = fmincon(@ObjFunc7Fmincon,k0,[],[],[],[],lb,ub,[],[],x0,yexp); fprintf('\n使用函數(shù)fmincon()估計得到的參數(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ù)估計 [k,resnorm,residual,exitflag,output,lambda,jacobian] = ... lsqnonlin(@ObjFunc7LNL,k0,lb,ub,[],x0,yexp); ci = nlparci(k,residual,jacobian); fprintf('\n\n使用函數(shù)lsqnonlin()估計得到的參數(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()估計得到的結(jié)果為初值，使用函數(shù)lsqnonlin()進(jìn)行參數(shù)估計 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()估計得到的參數(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ù)此樓» 猜你喜歡
計算與實驗科學(xué)工程國際會議(8月7-11日,香港) 二維材料分會S5與聲學(xué)換能器分會S8
已經(jīng)有0人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑有論文
已經(jīng)有4人回復(fù)
-
物理學(xué)I論文潤色/翻譯怎么收費?
已經(jīng)有102人回復(fù)
-
歡迎調(diào)劑到贛南師范大學(xué)智能制造與未來能源學(xué)院
已經(jīng)有28人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑 本科期間發(fā)表一篇Sci
已經(jīng)有4人回復(fù)
-
歡迎調(diào)劑到贛南師范大學(xué)智能制造與未來能源學(xué)院
已經(jīng)有26人回復(fù)
-
重慶交大26年碩士生招生擬調(diào)劑通知已出！歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。
已經(jīng)有0人回復(fù)
-
UJN物理學(xué)專業(yè)調(diào)劑
已經(jīng)有3人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑 本科期間發(fā)表一篇Sci
已經(jīng)有0人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 求助MATLAB中程序看不懂的簡單問題！謝謝！
已經(jīng)有3人回復(fù)
- 如何求這個matlab程序畫出的legend曲線的法矢量？
已經(jīng)有6人回復(fù)
- 高分求將我的MATLAB程序改寫為FORTRAN程序（1000BB）
已經(jīng)有13人回復(fù)
- 求在matlab中進(jìn)行分類的程序
已經(jīng)有5人回復(fù)
- 用matlab畫一個函數(shù)，能幫我看看程序哪里有問題嗎
已經(jīng)有8人回復(fù)
- 求教如何用matlab繪制下面兩圖
已經(jīng)有5人回復(fù)
- 求大神用matlab給我處理一組數(shù)據(jù)，有人沒有啊？
已經(jīng)有15人回復(fù)
- matlab一運行就出現(xiàn)的問題
已經(jīng)有7人回復(fù)
- 有個matlab程序，但是運行不下來，求高手
已經(jīng)有23人回復(fù)
- matlab程序中調(diào)用m文件的問題
已經(jīng)有13人回復(fù)
- 幫忙詳解示范一下MATLAB中g(shù)input函數(shù)的應(yīng)用吧
已經(jīng)有3人回復(fù)
- 實在是看不出這段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程序為啥出錯了
已經(jīng)有5人回復(fù)
- matlab程序問題，總是出錯，麻煩大家看看問題在那里？
已經(jīng)有6人回復(fù)
- 【求助】幫我編個程啊~C語言，C++，Matlab，VB...都行~
已經(jīng)有14人回復(fù)
月只藍(lán)
無蟲
-
專家經(jīng)驗: +1059 - 計算強帖: 8
- 應(yīng)助: 1712 (講師)
- 貴賓: 8.888
- 金幣: 68128.7
- 散金: 1938
- 紅花: 443
- 沙發(fā): 4
- 帖子: 4373
- 在線: 3291.4小時
- 蟲號: 1122189
- 注冊: 2010-10-14
- 專業(yè): 宇宙學(xué)
- 管轄: 計算模擬區(qū)
4樓2015-05-26 16:25:47
月只藍(lán)
禁蟲
-
專家經(jīng)驗: +1059 - 計算強帖: 8
- 應(yīng)助: 1712 (講師)
- 貴賓: 8.888
- 金幣: 68128.7
- 散金: 1938
- 紅花: 443
- 沙發(fā): 4
- 帖子: 4373
- 在線: 3291.4小時
- 蟲號: 1122189
- 注冊: 2010-10-14
- 專業(yè): 宇宙學(xué)
- 管轄: 計算模擬區(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
主管區(qū)長
- 應(yīng)助: 1 (幼兒園)
- 金幣: 6868.7
- 散金: 18
- 紅花: 1
- 帖子: 1492
- 在線: 209.6小時
- 蟲號: 2756771
- 注冊: 2013-10-27
- 性別: MM
- 專業(yè): 有機化工
您好，您的計算結(jié)果我看到了，我想說的是我上面編的Matlab程序正確嗎？為什么會出現(xiàn)那樣的計算結(jié)果？還有就是您計算的方法是最小二乘嗎？我現(xiàn)在需要用matlab計算，因為寫論文需要。所以您能不能幫忙看一下程序是否正確？ 3樓2015-05-26 15:38:21
dbb627
新蟲
-
專家經(jīng)驗: +4 - 計算強帖: 12
- 應(yīng)助: 289 (大學(xué)生)
- 貴賓: 0.589
- 金幣: 24640.4
- 散金: 551
- 紅花: 61
- 沙發(fā): 1
- 帖子: 1246
- 在線: 1794.8小時
- 蟲號: 149791
- 注冊: 2005-12-29
- 性別: GG
- 專業(yè): 污染控制化學(xué)
- 管轄: 計算模擬
|
程序有錯，最后面函數(shù)寫錯了 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]; k0 = [0.0671 1.253 0.02 48.028]; % 參數(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ù)估計 [k,fval,flag] = fmincon(@ObjFunc7Fmincon,k0,[],[],[],[],lb,ub,[],[],x0,yexp); fprintf('\n使用函數(shù)fmincon()估計得到的參數(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ù)估計 [k,resnorm,residual,exitflag,output,lambda,jacobian] = ... lsqnonlin(@ObjFunc7LNL,k0,lb,ub,[],x0,yexp); ci = nlparci(k,residual,jacobian); fprintf('\n\n使用函數(shù)lsqnonlin()估計得到的參數(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()估計得到的結(jié)果為初值，使用函數(shù)lsqnonlin()進(jìn)行參數(shù)估計 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()估計得到的參數(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,'o-',t,yexp(:,2),'*');legend('ca-pr','ca-real') % ------------------------------------------------------------------ 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); f = sum((x-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); f1 = x - yexp(:,2); f = [f1]; % ------------------------------------------------------------------ function dxdt = KineticEqs(t,x,k) dxdt=-(k(1)*k(2)*x^2-0.25*k(3)*k(4)*(9.404-x)^2)/(1+k(2)*x+0.5*k(4)*(9.404-x)); |
5樓2015-05-26 16:56:17
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 265求調(diào)劑 +4 | 林深溫瀾 2026-04-01 | 6/300 |
|
|---|---|---|---|---|
|
|
[考研] 319求調(diào)劑 +9 | 太容易1018 2026-04-01 | 9/450 |
|
|
[考研] 環(huán)境工程297分求調(diào)劑一志愿杭高院 +10 | GENJIOW 2026-03-31 | 11/550 |
|
|
|
[考研] 307分求調(diào)劑 +14 | (o~o) 2026-03-31 | 15/750 |
|
|
|
[考研] 286求調(diào)劑 +5 | lim0922 2026-03-26 | 5/250 |
|
|
|
[考研] 材料科學(xué)與工程339求調(diào)劑 +11 | hyz0119 2026-03-31 | 12/600 |
|
|
|
[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +8 | 邱gl 2026-03-30 | 16/800 |
|
|
|
[考研] 311求調(diào)劑一志愿合肥工業(yè)大學(xué) +14 | 秋二十二 2026-03-30 | 14/700 |
|
|
|
[考研] 材料調(diào)劑 +11 | 一樣YWY 2026-03-31 | 11/550 |
|
|
|
[考研] 314求調(diào)劑 +6 | 1xiaojun23 2026-03-31 | 6/300 |
|
|
|
[考研] 346求調(diào)劑 一志愿070303有機化學(xué) +11 | 蘿卜燉青菜 2026-03-28 | 12/600 |
|
|
|
[考研] 343求調(diào)劑 +8 | 愛羈絆 2026-03-28 | 8/400 |
|
|
|
[考研] 085600材料與化工329分求調(diào)劑 +5 | 這是有太陽哇 2026-03-25 | 5/250 |
|
|
|
[考研] 085601一志愿西北工業(yè)大學(xué)初試346 +4 | 085601初試346 2026-03-30 | 4/200 |
|
|
|
[考研] 一志愿雙一流機械285分求調(diào)劑 +4 | 幸運的三木 2026-03-29 | 5/250 |
|
|
|
[考研] 一志愿北京理工大學(xué)本科211材料工程294求調(diào)劑 +8 | mikasa的圍巾 2026-03-28 | 8/400 |
|
|
|
[考研] 復(fù)試調(diào)劑 +3 | raojunqi0129 2026-03-28 | 3/150 |
|
|
|
[考研] 394求調(diào)劑 +3 | 好事多磨靜候佳?/a> 2026-03-26 | 5/250 |
|
|
|
[考研] 藥學(xué)105500求調(diào)劑 +3 | Ssun。。 2026-03-28 | 3/150 |
|
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +4 | 崔wj 2026-03-26 | 4/200 |
|