 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 3582 | 回復(fù): 15 | |||
| 當(dāng)前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | |||
[求助]
matlab擬合方程參數(shù)時初值的選擇
|
|||
|
A+B=C A+W=R R+B=C+W 先考慮的三個反應(yīng)都是可逆反應(yīng)，有以下幾個式子：編寫程序如下 %30℃ clear all clc global keq1 keq2 t0 c0 keq1 = 3.94; % 平衡常數(shù)k1 keq3 = 1.08; % 平衡常數(shù)k3 t0 = [5,10, 15, 30, 45, 60,90,120,160,200,240,300,360]; c0 = [5.939601,5.203973,0,0.022484,0.078956] % A B C R W 組成的初始濃度 ci=[3.742508811 4.213566482 0.990406938 0.252498405 0.090415512; 3.045716838 3.387013322 1.816960097 0.215546741 0.081072022; 2.504754523 2.901991595 2.301981824 0.193427102 0.099462604; 1.65098702 2.036686653 3.167286766 0.170236582 0.095144044; 1.271633833 1.655460403 3.548513016 0.124730781 0.10334349; 1.129887499 1.630815394 3.573158025 0.129798695 0.097764543; 1.01239208 1.460748979 3.74322444 0.144649391 0.073614958; 1.007417603 1.429180292 3.774793127 0.123477124 0.094249307; 0.990005723 1.499216939 3.704756481 0.135524067 0.076072022; 0.970073726 1.493966669 3.71000675 0.130791851 0.076387812; 0.96903024 1.429705429 3.77426799 0.119510017 0.101202216; 0.980873385 1.457988124 3.745985296 0.124189274 0.097022161; 0.98477262 1.453402965 3.750570455 0.1200936 0.0762608]; % cA cC cR 對應(yīng)時間數(shù)值 k0=[1 1 1]; % k1+, k2, k3+ 初始向值 lb = [0 0 0]; ub = [1000 1000 500]; % 上下限 % 使用函數(shù)lsqnonlin()進(jìn)行參數(shù)估計 [k,resnorm,residual,exitflag] = lsqnonlin(@ObjFunc,k0,lb,ub,[],ci); k1plus=k(1); k1minus= k1plus/keq1; k2=k(2); k3plus=k(3); k3minus= k3plus/keq3; % ------------------------------------------------------------------ function f = ObjFunc(k,ci) % 目標(biāo)函數(shù) global t0 c0 [t,c_cal] = ode45(@Euqations,t0,c0,[],k); f= sum((c_cal-ci)^2); % ------------------------------------------------------------------ function dcdt = Euqations(t,c,k) global keq1 keq2 k1plus=k(1); k1minus= k1plus/keq1; k2=k(2); k3plus=k(3); k3minus= k3plus/keq3; cA=c(:,1); cB=c(:,2); cC=c(:,3); cR=c(:,4); cW=c(:,5); dcAdt = k1plus*cA*cB-k1minus*cC+k2*cA*cW; dcCdt = -k1plus*cA*cB+k1minus*cC+k3plus*cB*cR-k3minus*cC*cW; dcRdt = k2*cA*cW-k3plus*cB*cR+k3minus*cC*cW; dcdt = [dcAdt; dcCdt;dcRdt]; 但是k的初值一直影響結(jié)果，程序要怎么改一下初值對結(jié)果影響不大（這個程序有問題，） |
回復(fù)此樓» 猜你喜歡
2026年循環(huán)經(jīng)濟(jì)功能材料國際會議（ICFMCE 2026）
已經(jīng)有0人回復(fù)
-
2026年第五屆電氣、電子與信息工程國際會議（ISEEIE 2026）
已經(jīng)有0人回復(fù)
-
物理學(xué)I論文潤色/翻譯怎么收費(fèi)?
已經(jīng)有196人回復(fù)
-
哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑
已經(jīng)有5人回復(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ù)
-
請問還有沒有用Latex寫文章的小伙伴們？
已經(jīng)有0人回復(fù)
-
光學(xué)工程學(xué)碩調(diào)劑信息
已經(jīng)有24人回復(fù)
-
歡迎加入課題組
已經(jīng)有0人回復(fù)
-
散金幣，求好運(yùn)，祝面上順利！
已經(jīng)有17人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 在matlab中擬合曲面的表達(dá)式
已經(jīng)有29人回復(fù)
- matlab 求指點 動力學(xué)方程擬合過程中導(dǎo)數(shù)的獲取
已經(jīng)有10人回復(fù)
- matlab擬合擬均相動力學(xué)參數(shù)
已經(jīng)有6人回復(fù)
- matlab-常微分方程參數(shù)估計
已經(jīng)有12人回復(fù)
- 求教origin擬合簡單的滲流方程（又稱逾滲方程），做復(fù)合材料的應(yīng)該會遇到吧
已經(jīng)有20人回復(fù)
- 非線性擬合估計方程參數(shù)（matlab，spss，origin都可以）
已經(jīng)有4人回復(fù)
- matlab非線性參數(shù)擬合問題
已經(jīng)有7人回復(fù)
- MATLAB求解一個比較復(fù)雜的方程
已經(jīng)有14人回復(fù)
- 【求助】用matlab最優(yōu)化方法進(jìn)行參數(shù)擬合
已經(jīng)有17人回復(fù)
- 【求助】使用MATLAB擬合出按照公式的系數(shù)
已經(jīng)有19人回復(fù)
- 如何用MATLAB擬合指數(shù)函數(shù)
已經(jīng)有3人回復(fù)
- 【求助】最小二乘法編寫指數(shù)擬合方程
已經(jīng)有19人回復(fù)
- 【求助】使用Matlab擬合反應(yīng)動力學(xué)方程問題
已經(jīng)有7人回復(fù)
|
有的，我寫在附件里了，真是太謝謝了 |
» 本帖附件資源列表
-
歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺，不對該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識產(chǎn)權(quán)問題，其責(zé)任在于用戶本人，如對版權(quán)有異議，請聯(lián)系郵箱：xiaomuchong@tal.com - 附件 1 : 太殘忍doc (1)1.doc
2012-03-28 15:05:19, 32.5 K
7樓2012-03-28 15:05:21
dingd
鐵桿木蟲 (職業(yè)作家)
- 計算強(qiáng)帖: 4
- 應(yīng)助: 1641 (講師)
- 金幣: 15037.3
- 散金: 101
- 紅花: 234
- 帖子: 3410
- 在線: 1223.7小時
- 蟲號: 291104
- 注冊: 2006-10-28
2樓2012-03-28 10:37:30
|
我寫在附件里面了 |
» 本帖附件資源列表
-
歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺，不對該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識產(chǎn)權(quán)問題，其責(zé)任在于用戶本人，如對版權(quán)有異議，請聯(lián)系郵箱：xiaomuchong@tal.com - 附件 1 : 太殘忍doc.doc
2012-03-28 10:52:17, 30.5 K
3樓2012-03-28 10:52:39
dingd
鐵桿木蟲 (職業(yè)作家)
- 計算強(qiáng)帖: 4
- 應(yīng)助: 1641 (講師)
- 金幣: 15037.3
- 散金: 101
- 紅花: 234
- 帖子: 3410
- 在線: 1223.7小時
- 蟲號: 291104
- 注冊: 2006-10-28
4樓2012-03-28 11:07:18
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料專碩326求調(diào)劑 +6 | 墨煜姒莘 2026-03-15 | 7/350 |
|
|---|---|---|---|---|
|
|
[考研] 本人考085602 化學(xué)工程 專碩 +16 | 不知道叫什么！ 2026-03-15 | 18/900 |
|
|
|
[考研] 26考研求調(diào)劑 +6 | 丶宏Sir 2026-03-13 | 6/300 |
|
|
|
[考研] 268求調(diào)劑 +6 | 好運(yùn)連綿不絕 2026-03-12 | 7/350 |
|
|
|
[考博] 26申博 +4 | 八旬速覽 2026-03-16 | 4/200 |
|
|
|
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6 | 薛云鵬 2026-03-15 | 6/300 |
|
|
[考研] 302求調(diào)劑 +4 | 小賈同學(xué)123 2026-03-15 | 8/400 |
|
|
|
[考研] 267一志愿南京工業(yè)大學(xué)0817化工求調(diào)劑 +6 | SUICHILD 2026-03-12 | 6/300 |
|
|
|
[考研] 304求調(diào)劑 +7 | 小熊joy 2026-03-14 | 7/350 |
|
|
|
[考研] 283求調(diào)劑 +3 | 聽風(fēng)就是雨； 2026-03-16 | 3/150 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 290分有科研經(jīng)歷，論文在投 +7 | 膩膩gk 2026-03-14 | 7/350 |
|
|
|
[考研] 26考研一志愿中國石油大學(xué)(華東)305分求調(diào)劑 +3 | 嘉年新程 2026-03-15 | 3/150 |
|
|
|
[考研] 085601材料工程315分求調(diào)劑 +3 | yang_0104 2026-03-15 | 3/150 |
|
|
|
[考研] 中科大材料專碩319求調(diào)劑 +3 | 孟鑫材料 2026-03-13 | 3/150 |
|
|
|
[考研] 復(fù)試調(diào)劑 +9 | Copy267 2026-03-10 | 9/450 |
|
|
|
[考研] 求材料調(diào)劑 085600英一數(shù)二總分302 前三科235 精通機(jī)器學(xué)習(xí) 一志愿哈工大 +4 | 林yaxin 2026-03-12 | 4/200 |
|
|
|
[考研] 0703化學(xué)一志愿211 總分320求調(diào)劑 +5 | 瑪卡巴卡啊哈 2026-03-11 | 5/250 |
|
|
|
[考研] 281求調(diào)劑 +9 | Koxui 2026-03-12 | 11/550 |
|
|
|
[考研] 土木第一志愿276求調(diào)劑，科研和技能十分豐富，求新興方向的導(dǎo)師收留 +3 | 土木小天才 2026-03-12 | 3/150 |
|
|
|
[論文投稿]
投稿問題
 5+4
|
星光燦爛xt 2026-03-12 | 6/300 |
|