版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進行實名認證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

返回列表

805592836_01

鐵蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 823
帖子: 10
在線: 23.4小時
蟲號: 2393210
注冊: 2013-04-01
性別: MM
專業(yè): 能源化工

[求助] 跪求matlab大神幫忙已有1人參與

matlab擬合動力學(xué)模型：
r0=dCA/dt=(k1*KA*KB*CA0*CB0^3)/(1+KA*CA0+KB^(1/3)*CB)^4
動力學(xué)數(shù)據(jù)(其中的三組)：
CB0為定值0.4
一、
t=0,2,5,7,10,15,20,30,40
CA=1.686,1.530,1.428,1.424,1.326,1.203,1.017,0.910,0.775
二、
t=0,2,3,5,10,15,20,30,40
CA=,1.742,1.444,1.371,1.276,1.123,0.991,0.888,0.665,0.521
三、
t=0,5,10,15,20,30,40，50
CA=,1.745,1.563,1.298,1.219,1.126,0.885,0.736,0.634

回復(fù)此樓

» 猜你喜歡

光學(xué)工程學(xué)碩調(diào)劑信息已經(jīng)有32人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費? 已經(jīng)有203人回復(fù)
散金幣，求好運，祝面上順利！已經(jīng)有79人回復(fù)
華東師范大學(xué)芯片設(shè)計徐瓏真實水平咋樣已經(jīng)有1人回復(fù)
調(diào)劑已經(jīng)有0人回復(fù)
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有0人回復(fù)
中德博士后2026項目5月30日截止- 國內(nèi)博士后入站、德國兩年全職已經(jīng)有0人回復(fù)
面向全球申請人的中德博士后2026項目5月30日截止- 國內(nèi)博后入站德國兩年全職已經(jīng)有0人回復(fù)
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有2人回復(fù)
重慶交大材料學(xué)院26年碩士招生復(fù)試實施細則已出。歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。已經(jīng)有0人回復(fù)

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

求在matlab中進行分類的程序已經(jīng)有5人回復(fù)
兩道材料力學(xué)題，還求大神幫助啊~~~~~ 已經(jīng)有5人回復(fù)
跪求大神幫看matlab小程序已經(jīng)有3人回復(fù)
CV作圖求助，大神們，急急急求已經(jīng)有14人回復(fù)
求大神用matlab給我處理一組數(shù)據(jù)，有人沒有��？已經(jīng)有15人回復(fù)
修改圖標問題，請各位大神幫忙！已經(jīng)有5人回復(fù)
求大神幫忙用1stOpt 擬合復(fù)數(shù)數(shù)據(jù) 已經(jīng)有13人回復(fù)
關(guān)于matlab進行傅里葉變換和逆變換的程序，求大神指導(dǎo)！已經(jīng)有9人回復(fù)
大神來幫忙解釋一下，上下兩個圖是什么意思已經(jīng)有10人回復(fù)
matlab解決重心法選址問題，運行不了啊，求助各位大神已經(jīng)有10人回復(fù)
求助大神解決問題呀~~ 已經(jīng)有4人回復(fù)
求大神告訴下模擬軟件已經(jīng)有4人回復(fù)
【求助MATLAB大神】MATLAB畫多行圖時，每行的高度不太一樣已經(jīng)有14人回復(fù)
MATLAB微分方程參數(shù)擬合問題，求大神已經(jīng)有7人回復(fù)
求大神幫我改一下下面的程序！matlab迭代法的不知道錯在哪里！已經(jīng)有7人回復(fù)
求大神幫忙解病態(tài)線性方程組已經(jīng)有5人回復(fù)
多元非線性回歸，求大神幫忙，，，，，，，，，，，，，，，，，，，，，，已經(jīng)有6人回復(fù)
求大神幫忙研究一下關(guān)于鋁合金的抗彎強度的問題！�。。。。�！已經(jīng)有13人回復(fù)

加油加油加油……

1樓 2014-02-23 10:35:36

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

月只藍

主管區(qū)長 (職業(yè)作家)

專家經(jīng)驗: +1059
計算強帖: 8
應(yīng)助: 1712 (講師)
貴賓: 8.888
金幣: 68127.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線: 3291.4小時
蟲號: 1122189
注冊: 2010-10-14
專業(yè): 宇宙學(xué)
管轄: 計算模擬區(qū)

【答案】應(yīng)助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
805592836_01: 金幣+10, ★★★很有幫助, 救人一命 2014-02-24 09:35:43

原方程式：dCA/dt=(k1*KA*KB*CA0*CB0^3)/(1+KA*CA0+KB^(1/3)*CB)^4中，沒有給出CB和CA的關(guān)系，在定義方程式時默認CB=CA。對于第一組數(shù)據(jù)，代碼如下，結(jié)果見附圖。
推薦你找版主用1stopt軟件幫你跑一下，MATLAB求解初值不好給出，結(jié)果一定是最優(yōu)。

CODE:

function ode_fit33333333333

clear all;clc

format long

tspan=[0,2,5,7,10,15,20,30,40];

t1=length(tspan)

yexp=[1.530,1.428,1.424,1.326,1.203,1.017,0.910,0.775]';

 t2=length(yexp)

y0=1.686;

k0=[5 -10 100];

lb=[-1 -1 -1]*1e9;

ub=[1 1 1]*1e9;        

yy=[y0 yexp'];

options = optimset('MaxFunEvals',1000);

[k,resnorm,residual,exitflag,output,lambda,jacobian] = ...

    lsqnonlin(@ObjFunc,k0,lb,ub,options,tspan,y0,yexp);

ci = nlparci(k,residual,jacobian);

fprintf('\n\n使用函數(shù)lsqnonlin()估計得到的參數(shù)值為:\n')

fprintf('\t反應(yīng)速率常數(shù) k1 = %.4f\n',k(1))

fprintf('\t反應(yīng)速率常數(shù) KA = %.4f\n',k(2))

fprintf('\t反應(yīng)速率常數(shù) KB = %.4f\n',k(3))

ts=0:1:max(tspan);

[ts ys]=ode45(@KineticsEqs,ts,y0,[],k);

[ttt XXsim] = ode45(@KineticsEqs,tspan,y0,[],k); 

y=XXsim(2:end);

xexp=yexp;

R2=1-sum((xexp-y).^2)./sum((xexp-mean(y)).^2);

fprintf('\n\t相關(guān)系數(shù)之平方R^2 = %.6f',R2);

plot(ts,ys,'b',tspan,yy,'or'),legend('計算值','實驗值','Location','best'),

%---------------------------------------------------------

function f = ObjFunc(k,tspan,y0,yexp)           % 目標函數(shù)

[t Xsim] = ode45(@KineticsEqs,tspan,y0,[],k) ;

Xsim;

ysim = Xsim(2:end);

f=ysim-yexp;

 %----------------------------------------------------------

function dydt = KineticsEqs(t,y,k)

y0=1.686;

CB=y;

CB0=0.4;

CA0=y0;

k1=k(1);

KA=k(2);

KB=k(3);

dydt=(k1*KA*KB*CA0*CB0^3)./(1+KA*CA0+KB.^(1/3)*CB).^4;

附圖.jpg

贊一下(1人)

回復(fù)此樓

MATLAB、MS小問題、普通問題請發(fā)帖求助！時間精力有限，恕不接受無償私信求助。

2樓2014-02-24 09:16:40

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 805592836_01 的主題更新

返回列表

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[材料工程] 一志愿C9材料與化工專業(yè)總分300求調(diào)劑 +3	曼111 2026-03-24	4/200	2026-03-24 14:45 by JourneyLucky
[考研] 08工學(xué)調(diào)劑 +7	用戶573181 2026-03-20	12/600	2026-03-24 11:17 by 用戶573181
[考研] 一志愿北京化工大學(xué) 070300 學(xué)碩 336分求調(diào)劑 +7	vv迷 2026-03-22	7/350	2026-03-23 23:44 by Txy@872106
[考研] 070300化學(xué)求調(diào)劑 +8	苑豆豆 2026-03-20	8/400	2026-03-23 20:57 by baobaoye
[考研] 化學(xué)308分求調(diào)劑 +3	你好明天你好 2026-03-23	3/150	2026-03-23 20:11 by macy2011
[考研] 336化工調(diào)劑 +4	王大坦1 2026-03-23	5/250	2026-03-23 18:32 by allen-yin
[考研] 291 求調(diào)劑 +4	化工2026屆畢業(yè)?/a> 2026-03-21	5/250	2026-03-23 16:46 by 化工2026屆畢業(yè)?/a>
[考研] 一志愿東華大學(xué)化學(xué)070300，求調(diào)劑 +7	2117205181 2026-03-21	8/400	2026-03-22 22:55 by chixmc
[考研] 280分求調(diào)劑一志愿085802 +4	PUMPT 2026-03-22	7/350	2026-03-22 22:13 by 星空星月
[考研] 315分，誠求調(diào)劑，材料與化工085600 +3	13756423260 2026-03-22	3/150	2026-03-22 20:11 by edmund7
[考研] 289材料與化工（085600）B區(qū)求調(diào)劑 +3	這么名字咋樣 2026-03-22	4/200	2026-03-22 17:56 by 云民大李老師
[考研] 289求調(diào)劑 +7	懷瑾握瑜l 2026-03-20	7/350	2026-03-22 15:57 by ColorlessPI
[考研] 085600材料與化工306 +4	z1z2z3879 2026-03-21	4/200	2026-03-21 23:44 by ms629
[考研] 求調(diào)劑 +4	要好好無聊 2026-03-21	4/200	2026-03-21 18:57 by 學(xué)員8dgXkO
[考研] 299求調(diào)劑 +5	shxchem 2026-03-20	7/350	2026-03-21 17:09 by ColorlessPI
[考研] 299求調(diào)劑 +6	△小透明* 2026-03-17	6/300	2026-03-21 02:42 by JourneyLucky
[考研] 一志愿武理材料305分求調(diào)劑 +6	想上岸的鯉魚 2026-03-18	7/350	2026-03-21 01:03 by JourneyLucky
[考研] 290求調(diào)劑 +7	^O^乜 2026-03-19	7/350	2026-03-20 21:43 by JourneyLucky
[考研] 求調(diào)劑 +3	@taotao 2026-03-20	3/150	2026-03-20 19:35 by JourneyLucky
[考研] 材料，紡織，生物（0856、0710），化學(xué)招生啦 +3	Eember. 2026-03-17	9/450	2026-03-18 10:28 by Eember.

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产 高清 中文字幕,99re热久久亚洲综合精品成人,熟妇 一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

24小時熱門版塊排行榜

805592836_01

[求助] 跪求matlab大神幫忙 已有1人參與

» 猜你喜歡

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

月只藍

【答案】應(yīng)助回帖

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产高清中文字幕,99re热久久亚洲综合精品成人,熟妇一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

[求助] 跪求matlab大神幫忙已有1人參與

» 本主題相關(guān)價值貼推薦，對您同樣有幫助: