 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 4 | 1/1 | 返回列表 |
| 查看: 910 | 回復(fù): 3 | ||
| 【懸賞金幣】回答本帖問題，作者修竹好衰將贈(zèng)送您 20 個(gè)金幣 | ||
修竹好衰新蟲 (初入文壇)
|
[求助]
初學(xué)用MATLAB求解動(dòng)力學(xué)模型，但是在非線性二乘法處出現(xiàn)問題 ，麻煩各位大神幫忙看看 已有2人參與
|
|
|
出現(xiàn)的問題是： 錯(cuò)誤使用 lsqncommon (line 64) The Levenberg-Marquardt algorithm does not handle bound constraints and the trust-region-reflective algorithm requires at least as many equations as variables; aborting. 出錯(cuò) lsqnonlin (line 262) lsqncommon(funfcn,xCurrent,lb,ub,options,defaultopt,optimgetFlag,caller,... 出錯(cuò) kinetic (line 16) lsqnonlin(@fun,k0,lb,ub,[],yexp); 其中反應(yīng)數(shù)據(jù)是引用文件：‘reaction.xls’中每7行為一個(gè)反應(yīng)溫度下的各個(gè)時(shí)間點(diǎn)數(shù)據(jù)，總共3個(gè)溫度下的數(shù)據(jù)，共計(jì)21行，已在后面圖片里面顯示 程序： clc clear all; global a b tspan=[0 1 3 5 7 9 11]; %各個(gè)反應(yīng)時(shí)間點(diǎn)，0h，1h，3h…… yexp=xlsread('reaction.xls') %y的賦值從reaction這個(gè)excel文件中調(diào)取 k0=[2*10^8 60000 6.75*10^10 72000 3*10^11 90000]; %每個(gè)反應(yīng)速率常數(shù)k所對(duì)應(yīng)的A、E，總計(jì)6個(gè)參數(shù)初值設(shè)定（正反應(yīng)、逆反應(yīng)、副反應(yīng)） lb=[0 0 0 0 0 0]; ub=[+inf +inf +inf +inf +inf +inf]; [k,resnorm,residual,exitflag,output,lambda,jacobian]=lsqnonlin(@fun,k0,lb,ub,[],yexp); %===================================================================================== function f=fun(k,yexp) %這個(gè)就是擬合的目標(biāo)函數(shù),得到對(duì)應(yīng)最合適的k f=[]; for i=1:3 %對(duì)應(yīng)三個(gè)溫度下的數(shù)據(jù)， tspan=[0 1 3 5 7 9 11]; a=7*(i-1)+1; x0=yexp(a,  ;[t,x]=ode45(@func,tspan,x0,[],k,i); d=a+6; yc1=x(:,1); yc2=x(:,2); yc3=x(:,3); yc4=x(:,4); f11=yexp(a:d,1)-yc1; f12=yexp(a:d,2)-yc2; f13=yexp(a:d,3)-yc3; f14=yexp(a:d,4)-yc4; ff=[f11 f12 f13 f14]; f=[f:ff]; end %===================================================================================== function dxdt=func(t,x,k,i) switch i case 1 m=0.6; V=0.2431; T=303.15; case 2 m=0.6; V=0.2448; T=313.15; case 3 m=0.6; V=0.2491; T=323.15; end r2=-(k(1)*exp(-k(2)/8.314/T)*x(1)/V*x(2)/V-k(3)*exp(-k(4)/8.314/T)*x(3)/V); r3=-r2; r4=k(5)*exp(-k(6)/8.314/T)*x(1)/V*x(1)/V; r1=r2-2*r4; dxdt=[m*r1;m*r2;m*r3;m*r4]; %因?yàn)檫@里r的單位是mol/(g min L)所以需要乘以質(zhì)量 |
回復(fù)此樓» 猜你喜歡
西安電子科學(xué)大學(xué)杭州研究院劉麗香教授招收智能多模態(tài)傳感器和微型儲(chǔ)能器件方向博士
已經(jīng)有13人回復(fù)
-
重慶交通大學(xué)光子學(xué)微結(jié)構(gòu)與器件課題組招收2026年碩士研究生信息
已經(jīng)有1人回復(fù)
-
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)?
已經(jīng)有300人回復(fù)
-
一志愿鄭大材料學(xué)碩298分，求調(diào)劑
已經(jīng)有6人回復(fù)
-
尋合作：應(yīng)力腐蝕多尺度模擬
已經(jīng)有3人回復(fù)
-
考研交流
已經(jīng)有0人回復(fù)
-
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有0人回復(fù)
-
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有0人回復(fù)
dingd
鐵桿木蟲 (職業(yè)作家)
- 計(jì)算強(qiáng)帖: 4
- 應(yīng)助: 1641 (講師)
- 金幣: 15037.3
- 散金: 101
- 紅花: 234
- 帖子: 3410
- 在線: 1223.7小時(shí)
- 蟲號(hào): 291104
- 注冊(cè): 2006-10-28
修竹好衰2樓2022-03-05 19:23:00
修竹好衰
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 7.1
- 帖子: 28
- 在線: 22.6小時(shí)
- 蟲號(hào): 21829020
- 注冊(cè): 2020-04-14
- 專業(yè): 能源化工
3樓2022-03-06 13:56:17
4樓2022-03-06 20:23:03
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 修竹好衰 的主題更新
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 286 +3 | ksncj 2026-03-04 | 3/150 |
|
|---|---|---|---|---|
|
[考研]
 一志愿中南大學(xué)理學(xué)化學(xué)
+6
|
15779376950 2026-03-01 | 8/400 |
|
|
|
[考研] 347求調(diào)劑 +5 | 啊歐歐歐 2026-03-03 | 5/250 |
|
|
|
[考博] 26申博-目前4篇SCI一作 +3 | chen_2024 2026-03-02 | 3/150 |
|
|
|
[考研] 0703 總分319求調(diào)劑 +4 | Xinyuu 2026-03-02 | 4/200 |
|
|
|
[考研] 276求調(diào)劑 +8 | 路lyh123 2026-02-28 | 10/500 |
|
|
|
[考研] 085602化學(xué)工程350，調(diào)劑，有沒有211的 +5 | 利好利好. 2026-03-02 | 9/450 |
|
|
|
[考研] 0805總分292，求調(diào)劑 +12 | 幻想之殤 2026-03-01 | 12/600 |
|
|
|
[考研] 284求調(diào)劑 +6 | 天下熯 2026-03-02 | 6/300 |
|
|
|
[考研] 課題組接收材料類調(diào)劑研究生 +6 | gaoxiaoniuma 2026-02-28 | 9/450 |
|
|
|
[論文投稿]
通訊作者寫誰，問題是你意想不到的問題
 15+3
|
阿爾法啊 2026-03-01 | 3/150 |
|
|
|
[考研] 290分材料工程085601求調(diào)劑 數(shù)二英一 +8 | llx0610 2026-03-02 | 9/450 |
|
|
|
[考研] 求調(diào)劑 +7 | repeatt?t 2026-02-28 | 7/350 |
|
|
|
[考研] 285求調(diào)劑 +9 | 滿頭大汗的學(xué)生 2026-02-28 | 9/450 |
|
|
|
[考研] 0856求調(diào)劑285 +11 | 呂仔龍 2026-02-28 | 11/550 |
|
|
|
[考研] 281求調(diào)劑 +5 | 2026計(jì)算機(jī)_誠(chéng)心 2026-03-01 | 8/400 |
|
|
|
[碩博家園] 博士自薦 +7 | 科研狗111 2026-02-26 | 11/550 |
|
|
|
[考研] 272求調(diào)劑 +6 | 田智友 2026-02-28 | 6/300 |
|
|
|
[考研] 311求調(diào)劑 +6 | 亭亭亭01 2026-03-01 | 6/300 |
|
|
|
[考研] 295復(fù)試調(diào)劑 +3 | 簡(jiǎn)木ChuFront 2026-03-01 | 3/150 |
|