 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 1365 | 回復(fù): 10 | ||||
1135725495鐵桿木蟲 (著名寫手)
|
[求助]
求大家?guī)臀铱纯聪旅嬷衜atlab中的程序問題出在哪里？ 已有1人參與
|
|
求動(dòng)力學(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ù)此樓» 猜你喜歡
哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑
已經(jīng)有7人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑
已經(jīng)有5人回復(fù)
-
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)?
已經(jīng)有190人回復(fù)
-
求調(diào)劑
已經(jīng)有0人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑有論文
已經(jīng)有0人回復(fù)
-
請(qǐng)問還有沒有用Latex寫文章的小伙伴們？
已經(jīng)有0人回復(fù)
-
光學(xué)工程學(xué)碩調(diào)劑信息
已經(jīng)有31人回復(fù)
-
歡迎加入課題組
已經(jīng)有0人回復(fù)
-
散金幣，求好運(yùn)，祝面上順利！
已經(jīng)有70人回復(fù)
-
華東師范大學(xué)芯片設(shè)計(jì)徐瓏真實(shí)水平咋樣
已經(jīng)有1人回復(fù)
-
調(diào)劑
已經(jīng)有0人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 求助MATLAB中程序看不懂的簡(jiǎn)單問題！謝謝！
已經(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程序里哪里出錯(cuò)了，請(qǐng)大家?guī)兔纯矗?/a>
已經(jīng)有10人回復(fù)
- 請(qǐng)教MATLAB中的LMI求解問題，急！�。�
已經(jīng)有7人回復(fù)
- 求高手解答Matlab求解微分方程的初始值和邊界值問題
已經(jīng)有5人回復(fù)
- 求大神幫我改一下下面的程序！matlab迭代法的 不知道錯(cuò)在哪里！
已經(jīng)有7人回復(fù)
- matlab中遇到的畫圖小問題，程序沒問題。標(biāo)題要長(zhǎng)@@@@@@@@@@@@@@@@@@@@@@
已經(jīng)有8人回復(fù)
- 大家愛幫我看看這段MATLAB程序?yàn)樯冻鲥e(cuò)了
已經(jīng)有5人回復(fù)
- matlab程序問題，總是出錯(cuò)，麻煩大家看看問題在那里？
已經(jīng)有6人回復(fù)
- 【求助】幫我編個(gè)程啊~C語言，C++，Matlab，VB...都行~
已經(jīng)有14人回復(fù)
月只藍(lán)
主管區(qū)長(zhǎng) (職業(yè)作家)
-
專家經(jīng)驗(yàn): +1059 - 計(jì)算強(qiáng)帖: 8
- 應(yīng)助: 1712 (講師)
- 貴賓: 8.888
- 金幣: 68127.7
- 散金: 1938
- 紅花: 443
- 沙發(fā): 4
- 帖子: 4373
- 在線: 3291.4小時(shí)
- 蟲號(hào): 1122189
- 注冊(cè): 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í)
- 蟲號(hào): 2756771
- 注冊(cè): 2013-10-27
- 性別: MM
- 專業(yè): 有機(jī)化工
您好，您的計(jì)算結(jié)果我看到了，我想說的是我上面編的Matlab程序正確嗎？為什么會(huì)出現(xiàn)那樣的計(jì)算結(jié)果？還有就是您計(jì)算的方法是最小二乘嗎？我現(xiàn)在需要用matlab計(jì)算，因?yàn)閷懻撐男枰Ｋ阅懿荒軒兔匆幌鲁绦蚴欠裾_？ 3樓2015-05-26 15:38:21
月只藍(lán)
主管區(qū)長(zhǎng) (職業(yè)作家)
-
專家經(jīng)驗(yàn): +1059 - 計(jì)算強(qiáng)帖: 8
- 應(yīng)助: 1712 (講師)
- 貴賓: 8.888
- 金幣: 68127.7
- 散金: 1938
- 紅花: 443
- 沙發(fā): 4
- 帖子: 4373
- 在線: 3291.4小時(shí)
- 蟲號(hào): 1122189
- 注冊(cè): 2010-10-14
- 專業(yè): 宇宙學(xué)
- 管轄: 計(jì)算模擬區(qū)
4樓2015-05-26 16:25:47
dbb627
榮譽(yù)版主 (著名寫手)
-
專家經(jīng)驗(yàn): +4 - 計(jì)算強(qiáng)帖: 12
- 應(yīng)助: 289 (大學(xué)生)
- 貴賓: 0.589
- 金幣: 24640.4
- 散金: 551
- 紅花: 61
- 沙發(fā): 1
- 帖子: 1246
- 在線: 1794.8小時(shí)
- 蟲號(hào): 149791
- 注冊(cè): 2005-12-29
- 性別: GG
- 專業(yè): 污染控制化學(xué)
- 管轄: 計(jì)算模擬
|
程序有錯(cuò)，最后面函數(shù)寫錯(cuò)了 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ù)估計(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,'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
1135725495
鐵桿木蟲 (著名寫手)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 6868.7
- 散金: 18
- 紅花: 1
- 帖子: 1492
- 在線: 209.6小時(shí)
- 蟲號(hào): 2756771
- 注冊(cè): 2013-10-27
- 性別: MM
- 專業(yè): 有機(jī)化工
6樓2015-05-26 17:03:59
1135725495
鐵桿木蟲 (著名寫手)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 6868.7
- 散金: 18
- 紅花: 1
- 帖子: 1492
- 在線: 209.6小時(shí)
- 蟲號(hào): 2756771
- 注冊(cè): 2013-10-27
- 性別: MM
- 專業(yè): 有機(jī)化工
7樓2015-05-26 17:05:15
dbb627
榮譽(yù)版主 (著名寫手)
-
專家經(jīng)驗(yàn): +4 - 計(jì)算強(qiáng)帖: 12
- 應(yīng)助: 289 (大學(xué)生)
- 貴賓: 0.589
- 金幣: 24640.4
- 散金: 551
- 紅花: 61
- 沙發(fā): 1
- 帖子: 1246
- 在線: 1794.8小時(shí)
- 蟲號(hào): 149791
- 注冊(cè): 2005-12-29
- 性別: GG
- 專業(yè): 污染控制化學(xué)
- 管轄: 計(jì)算模擬
11357254958樓2015-05-26 17:08:44
1135725495
鐵桿木蟲 (著名寫手)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 6868.7
- 散金: 18
- 紅花: 1
- 帖子: 1492
- 在線: 209.6小時(shí)
- 蟲號(hào): 2756771
- 注冊(cè): 2013-10-27
- 性別: MM
- 專業(yè): 有機(jī)化工
9樓2015-05-26 17:09:52
1135725495
鐵桿木蟲 (著名寫手)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 6868.7
- 散金: 18
- 紅花: 1
- 帖子: 1492
- 在線: 209.6小時(shí)
- 蟲號(hào): 2756771
- 注冊(cè): 2013-10-27
- 性別: MM
- 專業(yè): 有機(jī)化工
10樓2015-05-26 17:13:33
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 1135725495 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 311求調(diào)劑 +3 | 26研0 2026-03-20 | 3/150 |
|
|---|---|---|---|---|
|
|
[考研]
 一志愿070300浙大化學(xué)358分，求調(diào)劑！
+3
|
酥酥魚.. 2026-03-21 | 3/150 |
|
|
[考研] 初試 317 +7 | 半拉月丙 2026-03-20 | 7/350 |
|
|
|
[考研] 材料學(xué)碩333求調(diào)劑 +3 | 北道巷 2026-03-18 | 3/150 |
|
|
|
[考研] 278求調(diào)劑 +9 | 煙火先于春 2026-03-17 | 9/450 |
|
|
|
[考研] 296求調(diào)劑 +4 | www_q 2026-03-20 | 4/200 |
|
|
|
[考研] 299求調(diào)劑 +5 | shxchem 2026-03-20 | 7/350 |
|
|
|
[考研] 求調(diào)劑 +3 | 白QF 2026-03-21 | 3/150 |
|
|
|
[考研] 求調(diào)劑 +6 | Mqqqqqq 2026-03-19 | 6/300 |
|
|
|
[考研] 299求調(diào)劑 +6 | △小透明* 2026-03-17 | 6/300 |
|
|
|
[考研] 一志愿西南交大，求調(diào)劑 +5 | 材化逐夢(mèng)人 2026-03-18 | 5/250 |
|
|
|
[考研] 22408 344分 求調(diào)劑 一志愿 華電計(jì)算機(jī)技術(shù) +4 | solanXXX 2026-03-20 | 4/200 |
|
|
|
[考研] 308求調(diào)劑 +3 | 阿姐阿姐家啊 2026-03-18 | 3/150 |
|
|
|
[考研] 350求調(diào)劑 +5 | weudhdk 2026-03-19 | 5/250 |
|
|
|
[考研] 環(huán)境工程調(diào)劑 +9 | 大可digkids 2026-03-16 | 9/450 |
|
|
|
[考研] 08工學(xué)調(diào)劑 +5 | 用戶573181 2026-03-20 | 5/250 |
|
|
|
[考研] 生物學(xué)調(diào)劑招人�。�！ +3 | 山海天嵐 2026-03-17 | 4/200 |
|
|
|
[考研] 286求調(diào)劑 +6 | lemonzzn 2026-03-16 | 10/500 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 +5 | pupcoco 2026-03-17 | 8/400 |
|
|
|
[考研] 085601求調(diào)劑 +4 | Du.11 2026-03-16 | 4/200 |
|