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

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

24小時(shí)熱門版塊排行榜

返回列表

下一頁(yè)

liqianmelody

至尊木蟲 (著名寫手)

應(yīng)助: 0 (幼兒園)
金幣: 11387.8
散金: 2
紅花: 4
帖子: 1046
在線: 73.9小時(shí)
蟲號(hào): 1334159
注冊(cè): 2011-06-29
性別: MM
專業(yè): 分離過(guò)程

[求助] 求助用matlab三次樣條法算一個(gè)反應(yīng)的反應(yīng)速率已有1人參與

本人不會(huì)用matlab這個(gè)軟件，我有實(shí)驗(yàn)數(shù)據(jù)，請(qǐng)各位精通的人幫我算一下表格標(biāo)黃的那列反應(yīng)速率，謝謝了！數(shù)據(jù)在附件里了，PS：附件中上傳時(shí)間較晚的那個(gè)是正確文件

回復(fù)此樓

» 本帖附件資源列表

歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺(tái)，不對(duì)該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識(shí)產(chǎn)權(quán)問(wèn)題，其責(zé)任在于用戶本人，如對(duì)版權(quán)有異議，請(qǐng)聯(lián)系郵箱：xiaomuchong@tal.com
附件 1 : 動(dòng)力學(xué)實(shí)驗(yàn)數(shù)據(jù).doc

2014-03-26 14:46:39, 79.5 K

附件 2 : 動(dòng)力學(xué)實(shí)驗(yàn)數(shù)據(jù).doc

2014-03-26 14:48:30, 79.5 K

» 猜你喜歡

求調(diào)劑已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有0人回復(fù)
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有100人回復(fù)
請(qǐng)問(wèn)還有沒有用Latex寫文章的小伙伴們？已經(jīng)有0人回復(fù)
光學(xué)工程學(xué)碩調(diào)劑信息已經(jīng)有32人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
散金幣，求好運(yùn)，祝面上順利！已經(jīng)有76人回復(fù)
華東師范大學(xué)芯片設(shè)計(jì)徐瓏真實(shí)水平咋樣已經(jīng)有1人回復(fù)
調(diào)劑已經(jīng)有0人回復(fù)
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有0人回復(fù)
中德博士后2026項(xiàng)目5月30日截止- 國(guó)內(nèi)博士后入站、德國(guó)兩年全職已經(jīng)有0人回復(fù)

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

反應(yīng)速率常數(shù)k的置信區(qū)間如何算？已經(jīng)有5人回復(fù)
求gaussian計(jì)算反應(yīng)速率常數(shù)，指前因子方法。已經(jīng)有25人回復(fù)
幫我看看聚合反應(yīng)三個(gè)速率常數(shù)怎么算？已經(jīng)有9人回復(fù)
【求助】催化反應(yīng)動(dòng)力學(xué)matlab計(jì)算各基元反應(yīng)的速率常數(shù)時(shí)，該如何避免較小量被忽略？已經(jīng)有3人回復(fù)
【求助】反應(yīng)速率常數(shù)怎么算已經(jīng)有5人回復(fù)

1樓 2014-03-26 14:46:44

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

月只藍(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
fegg7502: 金幣+10, 3ks 2014-03-29 08:41:56

比較粗糙的辦法是求數(shù)值導(dǎo)數(shù)，根據(jù)t CA的數(shù)據(jù)即可求出，如果需要MATLAB代碼參見：http://www.gaoyang168.com/bbs/viewthread.php?tid=7020495；
精確點(diǎn)的方法是先擬合出動(dòng)力學(xué)方程中的各參數(shù)，以第一組數(shù)據(jù)（即T=325.15K這組）為例，樓主沒給出反應(yīng)方程式，根據(jù)三個(gè)濃度隨時(shí)間的增幅，推測(cè)反應(yīng)式為：
                              A +    B    =    C
C為生成物，初始濃度為0。忽略反應(yīng)前后體系體積變化，則CA0=0.1530+0.0206=0.1736；CB0=3.4505+0.0206=3.4711；設(shè)物質(zhì)A的瞬時(shí)濃度為CA，有：
                              A       +       B          =          C
初始濃度                   CA0             CB0                      0
瞬時(shí)濃度                   CA             [CB0-(CA0-CA)]       CA0-CA

-(dCA/dt)=kCAmCBn=k0*exp(-Ea/RT)* CA^m*CB^n
得dCA/dt=-k0*exp(-Ea/RT)* CA^m* [CB0-(CA0-CA)]^n    （1）
該微分方程無(wú)法積分求得關(guān)于CA的代數(shù)式解析解，故直接對(duì)式（1）作常微分方程的擬合，求出指前因子k0、活化能Ea、兩個(gè)反應(yīng)級(jí)數(shù)m和n，這樣把不同時(shí)間t對(duì)應(yīng)的CA代入式（1）即得dCA/dt，取負(fù)得-dCA/dt。
MATLAB代碼如下：

CODE:

function ode3333

clear all;clc

format long

global T

T=323.15;

tspan=[0 

    30

60

90

120

180

240

300

360

420

480

540

600

];

t1=length(tspan);

yexp=[0.1530

0.1222

0.1036

0.0874

0.0619

0.0331

0.0305

0.0241

0.0185

0.0120

0.0080

0.0050

]';

 t2=length(yexp);

y0=0.1736;

k0=[0.1 10000 2 2];

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

ub=[1 1 1 1]*1e9;        

yy=[y0 yexp];

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

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

ci = nlparci(k,residual,jacobian);

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

fprintf('\t k0 = %.4f\n',k(1))

fprintf('\t Ea = %.4f\n',k(2))

fprintf('\t m = %.4f\n',k(3))

fprintf('\t n = %.4f\n',k(4))

%fprintf('\t反應(yīng)級(jí)數(shù) n = %.4f\n',k(2))

fprintf('  The sum of the squares is: %.1e\n\n',resnorm)

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);

k0=k(1);

Ea=k(2);

m=k(3);

n=k(4);

y=yexp';

dydt_a=-k0*exp(-Ea/8.314/T)*y.^m.*(3.2975+y).^n;

[-dydt_a]

figure(1)

plot(ts,ys,'b',tspan,yy,'or'),legend('計(jì)算值','實(shí)驗(yàn)值','Location','best'),

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

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

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

Xsim;

ysim = Xsim(2:end);

f=ysim-yexp';

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

function dydt = KineticsEqs(t,y,k)

global T

k0=k(1);

Ea=k(2);

m=k(3);

n=k(4);

dydt=-k0*exp(-Ea/8.314/T)*y.^m.*(3.2975+y).^n;

計(jì)算結(jié)果顯示：
k0 = 923.7021
   Ea = 9079.9827
   m = 1.1387
   n = -6.7529
相關(guān)系數(shù)之平方R^2 = 0.994517

-dCA/dt：

1.0e-003 *

0.865245002514304
0.711666673454071
0.611807848360504
0.520628349512815
0.369919263546430
0.192217346733282
0.176046594470231
0.136397181433893
0.102089651326710
0.063193844805022
0.040152229316972
0.023656194149690
擬合結(jié)果見附圖1。
其他T下面的數(shù)據(jù)，修改代碼中T、tspan、yexp對(duì)應(yīng)的數(shù)據(jù)即可。

如果樓主嫌麻煩，直接用前面說(shuō)的求數(shù)值導(dǎo)數(shù)的方法也可以。

求助用matlab三次樣條法算一個(gè)反應(yīng)的反應(yīng)速率

附圖1.jpg

贊一下

回復(fù)此樓

MATLAB、MS小問(wèn)題、普通問(wèn)題請(qǐng)發(fā)帖求助！時(shí)間精力有限，恕不接受無(wú)償私信求助。

2樓2014-03-26 20:00:43

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

月只藍(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)助回帖

★ ★
fegg7502: 金幣+2, 應(yīng)助指數(shù)+1, 3ks 2014-03-29 08:42:05

注意到你要求用三次樣條法，三次B樣條法求數(shù)值導(dǎo)數(shù)的例子參見下列鏈接的第12樓
http://www.gaoyang168.com/bbs/viewthread.php?tid=6744599&fpage=1&target=self&page=2
涉及代碼為”

CODE:

% 動(dòng)力學(xué)數(shù)據(jù)

t = [0  20  40  60  120  180  300];

CAm = [10  8  6  5  3  2  1];

% 用最小二乘樣條擬合法計(jì)算微分dCA/dt--使用不經(jīng)過(guò)實(shí)驗(yàn)點(diǎn)的B樣條插值函數(shù)

knots = 3;

K = 3;                  % 三次B樣條

sp = spap2(knots,K,t,CAm);

pp = fnder(sp);         % 計(jì)算B樣條函數(shù)的導(dǎo)函數(shù)

dCAdt = fnval(pp,t)    % 計(jì)算t處的導(dǎo)函數(shù)值

rAm = dCAdt;

“
該代碼轉(zhuǎn)自原鏈接第12樓。

贊一下

回復(fù)此樓

MATLAB、MS小問(wèn)題、普通問(wèn)題請(qǐng)發(fā)帖求助！時(shí)間精力有限，恕不接受無(wú)償私信求助。

3樓2014-03-26 20:19:21

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

liqianmelody

至尊木蟲 (著名寫手)

應(yīng)助: 0 (幼兒園)
金幣: 11387.8
散金: 2
紅花: 4
帖子: 1046
在線: 73.9小時(shí)
蟲號(hào): 1334159
注冊(cè): 2011-06-29
性別: MM
專業(yè): 分離過(guò)程

★
fegg7502: 金幣+1, 鼓勵(lì)交流 2014-03-29 08:42:18

引用回帖:

3樓: Originally posted by 月只藍(lán) at 2014-03-26 20:19:21
注意到你要求用三次樣條法，三次B樣條法求數(shù)值導(dǎo)數(shù)的例子參見下列鏈接的第12樓
http://www.gaoyang168.com/bbs/viewthread.php?tid=6744599&fpage=1&target=self&page=2
涉及代碼為”
% 動(dòng)力學(xué)數(shù)據(jù)
t = ;
C ...

謝謝你幫我解釋那么多，您樓上做出了動(dòng)力學(xué)方程的結(jié)果和我之前找人做的不太一樣，您做出的結(jié)果不太跟我的實(shí)際一樣，之前是通過(guò)1stopt做出的，沒有做出反應(yīng)速率，所以這次就想做出這個(gè)，我是做分離的，對(duì)matlab一點(diǎn)也不懂，所以插入什么代碼的也不會(huì)，就想請(qǐng)人幫我做出結(jié)果，我看您對(duì)此比較精通，如果您可以幫我做出，感激涕零，順便可以追加金幣，謝謝了�。�！

贊一下

回復(fù)此樓

4樓2014-03-26 21:21:56

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

月只藍(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)助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
liqianmelody: 金幣+50, ★★★★★最佳答案 2014-03-27 11:15:20
fegg7502: 金幣+2, 應(yīng)助指數(shù)+1, 3ks 2014-03-29 08:42:30

引用回帖:

4樓: Originally posted by liqianmelody at 2014-03-26 21:21:56
謝謝你幫我解釋那么多，您樓上做出了動(dòng)力學(xué)方程的結(jié)果和我之前找人做的不太一樣，您做出的結(jié)果不太跟我的實(shí)際一樣，之前是通過(guò)1stopt做出的，沒有做出反應(yīng)速率，所以這次就想做出這個(gè)，我是做分離的，對(duì)matlab一點(diǎn) ...

你說(shuō)的這個(gè)帖子么：
http://www.gaoyang168.com/bbs/viewthread.php?tid=6727491
實(shí)習(xí)版主dingd已經(jīng)把各個(gè)參數(shù)求出來(lái)了：
k0       870504.566962617
ea       53899.8216594799
m       1.02066293612806
n       1.08210287262186
再把CA代回-k0*exp(-Ea/(R*T1))*CA^m*CB^n不就得到想要的-dCA/dt了么？

贊一下

回復(fù)此樓

MATLAB、MS小問(wèn)題、普通問(wèn)題請(qǐng)發(fā)帖求助！時(shí)間精力有限，恕不接受無(wú)償私信求助。

5樓2014-03-26 21:29:35

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

justin534

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 175.5
帖子: 22
在線: 24小時(shí)
蟲號(hào): 3005267
注冊(cè): 2014-02-28
性別: GG
專業(yè): 化學(xué)反應(yīng)工程

★
fegg7502: 金幣+1, 鼓勵(lì)交流 2014-03-29 08:42:40

引用回帖:

5樓: Originally posted by 月只藍(lán) at 2014-03-26 21:29:35
你說(shuō)的這個(gè)帖子么：
http://www.gaoyang168.com/bbs/viewthread.php?tid=6727491
實(shí)習(xí)版主dingd已經(jīng)把各個(gè)參數(shù)求出來(lái)了：
k0       870504.566962617
ea       53899.8216594799
m       1.02066293612806
n    ...

不好意思版大,我修改了一下你的編碼,但是差不多只改數(shù)據(jù)而已,但不知會(huì)和不能run,可否幫我看一下,謝謝!

function ode3333
clear all;clc

format long

tspan=[50
62.5
80
100
125
166.666
250
333.3333
500
];

t1=length(tspan);
yexp=[1.414
1.1369
0.91816
0.81575
0.7424811
0.6323514
0.6139
0.6109
0.60789
]';

t2=length(yexp);

y0=1.9880715;

k0=[0.1 1];
lb=[0 0 ];
ub=[1 1 ]*1e9;

yy=[y0 yexp];

[k,resnorm,residual,exitflag,output,lambda,jacobian] = ...
lsqnonlin(@ObjFunc,k0,lb,ub,[],tspan,y0,yexp);
ci = nlparci(k,residual,jacobian);
fprintf('\n\n使用函?lsqnonlin()估?得到的??值?:\n')
fprintf('\t kf1 = %.4f\n',k(1))
fprintf('\t kf2 = %.4f\n',k(2))
fprintf('  The sum of the squares is: %.1e\n\n',resnorm)
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相?系?之平方R^2 = %.6f',R2);
kf1=k(1);
kf2=k(2);
y=yexp';
dydt_a=-(k(1)*y*((y+3.97614).^2)-(k(1)/k(2)).*(y-1.9880715).*((y-1.9880715).^2));
[-dydt_a]

figure(1)
plot(ts,ys,'b',tspan,yy,'or'),legend('?算值','??值','Location','best'),

%---------------------------------------------------------
function f = ObjFunc(k,tspan,y0,yexp)          % 目?函?
[t Xsim] = ode45(@KineticsEqs,tspan,y0,[],k) ;
Xsim;
ysim = Xsim(2:end);
f=ysim-yexp';
%----------------------------------------------------------

function dydt = KineticsEqs(t,y,k)
kf1=k(1);
kf2=k(2);
dydt=-(kf1*y*((y+3.97614).^2)-(kf1/kf2)*(y-1.9880715)*((y-1.9880715).^2));

贊一下

回復(fù)此樓

6樓2014-03-28 12:34:46

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

月只藍(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ū)

★
fegg7502: 金幣+1, 3ks 2014-03-29 08:42:52

引用回帖:

6樓: Originally posted by justin534 at 2014-03-28 12:34:46
不好意思版大,我修改了一下你的編碼,但是差不多只改數(shù)據(jù)而已,但不知會(huì)和不能run,可否幫我看一下,謝謝!

function ode3333
clear all;clc

format long

tspan=;

t1=length(tspan);
yexp=';

t2=len ...

t=0時(shí)，y=1.9880715么？
以下代碼按上述設(shè)定編寫：

CODE:

function ode3333

clear all;clc

format long

tspan=[0

    50

62.5

80

100

125

166.666

250

333.3333

500

];

t1=length(tspan);

yexp=[1.414

1.1369

0.91816

0.81575

0.7424811

0.6323514

0.6139

0.6109

0.60789

]';

t2=length(yexp);

y0=1.9880715;

k0=[1e-6 -0.1];

lb=[1 1]*-1e9;

ub=[1 0];        

yy=[y0 yexp];

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

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

ci = nlparci(k,residual,jacobian);

fprintf('\n\n使用函數(shù)lsqnonlin()估計(jì)得到:\n')

fprintf('\t kf1 = %.9f\n',k(1))

fprintf('\t kf2 = %.9f\n',k(2))

fprintf('  The sum of the squares is: %.1e\n\n',resnorm)

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);

kf1=k(1);

kf2=k(2);

y=yexp';

%dydt_a=-(k(1)*y*((y+3.97614).^2)-(k(1)/k(2)).*(y-1.9880715).*((y-1.9880715).^2));

%[-dydt_a]

figure(1)

plot(ts,ys,'b',tspan,yy,'or'),legend('計(jì)算值','實(shí)驗(yàn)值','Location','best'),

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

function f = ObjFunc(k,tspan,y0,yexp)           % 目?函?

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

Xsim;

ysim = Xsim(2:end);

f=ysim-yexp';

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

function dydt = KineticsEqs(t,y,k)

kf1=k(1);

kf2=k(2);

dydt=-(kf1*y*((y+3.97614).^2)-(kf1/kf2)*(y-1.9880715)*((y-1.9880715).^2));

計(jì)算結(jié)果：
使用函數(shù)lsqnonlin()估計(jì)得到:
   kf1 = 0.000315091
   kf2 = -0.232971187
  The sum of the squares is: 4.9e-002

相關(guān)系數(shù)之平方R^2 = 0.923278>>

求助用matlab三次樣條法算一個(gè)反應(yīng)的反應(yīng)速率-1

附圖1.png

贊一下

回復(fù)此樓

MATLAB、MS小問(wèn)題、普通問(wèn)題請(qǐng)發(fā)帖求助！時(shí)間精力有限，恕不接受無(wú)償私信求助。

7樓2014-03-28 16:31:03

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

pdl9527

專家顧問(wèn) (小有名氣)

專家經(jīng)驗(yàn): +8
應(yīng)助: 100 (初中生)
金幣: 2109.3
紅花: 23
帖子: 282
在線: 112小時(shí)
蟲號(hào): 1227333
注冊(cè): 2011-03-09
性別: GG
專業(yè): 分離過(guò)程
管轄: 計(jì)算模擬

★
fegg7502: 金幣+1, 鼓勵(lì)交流 2014-03-29 08:42:59

如果其他溫度求出來(lái)的參數(shù)差別很大怎么辦，能不能求得不同溫度的參數(shù)然后求平均值做呢？

贊一下

回復(fù)此樓

8樓2014-03-28 21:42:17

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

justin534

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 175.5
帖子: 22
在線: 24小時(shí)
蟲號(hào): 3005267
注冊(cè): 2014-02-28
性別: GG
專業(yè): 化學(xué)反應(yīng)工程

???????:

7?: Originally posted by ????? at 2014-03-28 16:31:03
t=0???y=1.9880715???
????????????趨??д??

function ode3333
clear all;clc

format long

tspan=;

t1=length(tspan);
yexp=';

t2=length(yexp);

y0=1.9880715;

k0=;
lb=*-1 ...

?x?x???,????kf2????????x???????,?????????N??????kf2?????,?M??Y??????e??,???????????????????},??????lb ub??^?g?O???????r??,?M??Y??????!

贊一下

回復(fù)此樓

9樓2014-03-29 12:28:35

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

月只藍(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ū)

???????:

9?: Originally posted by justin534 at 2014-03-29 12:28:35
?x?x???,????kf2????????x???????,?????????N??????kf2?????,?M??Y??????e??,???????????????????},??????lb ub??^?g?O???????r??,?M??Y??????!...

?????????????????????????????????
?????????????dingd ??1stopt???????????????????

贊一下

回復(fù)此樓

MATLAB、MS小問(wèn)題、普通問(wèn)題請(qǐng)發(fā)帖求助！時(shí)間精力有限，恕不接受無(wú)償私信求助。

10樓2014-03-29 15:11:44

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

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

返回列表

下一頁(yè)

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

[考研] 材料與化工考研調(diào)劑 +3	孅華 2026-03-22	3/150	2026-03-23 11:20 by barlinike
[考研] 一志愿中南大學(xué)化學(xué)學(xué)碩0703總分337求調(diào)劑 +4	niko- 2026-03-22	4/200	2026-03-23 07:56 by Iveryant
[考研] 280分求調(diào)劑一志愿085802 +4	PUMPT 2026-03-22	7/350	2026-03-22 22:13 by 星空星月
[考研] 08工科 320總分求調(diào)劑 +11	梨花珞晚風(fēng) 2026-03-17	11/550	2026-03-22 17:42 by luoyongfeng
[考研] 0805材料320求調(diào)劑 +3	深海物語(yǔ) 2026-03-20	3/150	2026-03-21 15:46 by 無(wú)際的草原
[考研] 求調(diào)劑 +6	Mqqqqqq 2026-03-19	6/300	2026-03-21 08:04 by JourneyLucky
[考研] 265求調(diào)劑 +3	Jack?k?y 2026-03-17	3/150	2026-03-21 03:17 by JourneyLucky
[考研] 二本跨考鄭大材料306英一數(shù)二 +3	z1z2z3879 2026-03-17	3/150	2026-03-21 02:29 by JourneyLucky
[考研] 材料 336 求調(diào)劑 +3	An@. 2026-03-18	4/200	2026-03-21 01:39 by JourneyLucky
[考研] 一志愿華中科技大學(xué)，080502，354分求調(diào)劑 +5	守候夕陽(yáng)CF 2026-03-18	5/250	2026-03-21 01:06 by JourneyLucky
[考研] 308求調(diào)劑 +3	阿姐阿姐家啊 2026-03-18	3/150	2026-03-20 23:24 by JourneyLucky
[考研] 288求調(diào)劑 +16	于海海海海 2026-03-19	16/800	2026-03-20 22:28 by JourneyLucky
[考研] 求調(diào)劑，一志愿:南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分 +4	@taotao 2026-03-19	4/200	2026-03-20 22:14 by JourneyLucky
[考研] 材料與化工 322求調(diào)劑 +4	然11 2026-03-19	4/200	2026-03-20 22:12 by luoyongfeng
[考研] 329求調(diào)劑 +9	想上學(xué)吖吖 2026-03-19	9/450	2026-03-20 22:01 by luoyongfeng
[考研] 290求調(diào)劑 +7	^O^乜 2026-03-19	7/350	2026-03-20 21:43 by JourneyLucky
[考研] 一志愿西南交通專碩材料355 本科雙非求調(diào)劑 +5	西南交通專材355 2026-03-19	5/250	2026-03-20 21:10 by JourneyLucky
[考研] 295復(fù)試調(diào)劑 +8	簡(jiǎn)木ChuFront 2026-03-19	8/400	2026-03-20 20:44 by zhukairuo
[考研] 一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑 +11	Ymlll 2026-03-18	15/750	2026-03-20 19:40 by 丁丁*
[考研] 一志愿福大288有機(jī)化學(xué)，求調(diào)劑 +3	小木蟲200408204 2026-03-18	3/150	2026-03-19 13:31 by houyaoxu

24小時(shí)熱門版塊排行榜

[求助] 求助用matlab三次樣條法算一個(gè)反應(yīng)的反應(yīng)速率 已有1人參與

» 本帖附件資源列表

» 猜你喜歡

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

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

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

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

[求助] 求助用matlab三次樣條法算一個(gè)反應(yīng)的反應(yīng)速率已有1人參與

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