老鸭窝在线观看免费观看,日日噜噜夜夜狠狠免费,91亚洲国产成人久久蜜臀

版塊導航: 正在加載中...

客戶端APP下載

論文輔導

調劑小程序

登錄注冊

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

24小時熱門版塊排行榜

返回列表

justin534

新蟲 (初入文壇)

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

???????:

10?: Originally posted by ????? at 2014-03-29 15:11:44
?????????????????????????????????
?????????????dingd ??1stopt???????????????????...

???????,???!????????}??!??????????????????W?,?F????????3??,??6??k1 k2 k3 k4 k5 k6!
R1=k(1)*DEC*PA-(k(1)/k(2))*EPC*EA
R2=k(3)*EPC*PA-(k(3)/k(4))*DPC*EA
R3=k(5)*EPC.^2-(k(5)/k(6))*DPC*DEC
???????????????a!???????RUN,??????????!

function ODE_parafit2
%-------data
yexp=[5.14  12.76  45.41  21.52  3.31;
5.36  13.20  46.94  20.90  3.17;
5.51  13.47  48.68  20.01  3.08;
5.78  14.03  48.16  19.34  3.04;
6.07  14.33  47.80  18.65  2.94]';
x0=[93.7 6.3 0 0 0];
k0 = [1 1 1 1 1 1];
lb = [0 0 0 0 0 0];
ub = [1 1 1 1 1 1]*1e5;
tspan = [0 38.46 50 62.5 100 166.666 250 500];
%-----------------
[k,fval,flag] =fmincon(@ObjFuncFmincon,k0,[],[],[],[],lb,ub,[],[],tspan,x0,yexp);
fprintf('\n\n??ú???fmincon()?????????????:\n')
fprintf('\tk1 = %.16f \n',k(1))
fprintf('\tk2 = %.16f \n',k(2))
fprintf('\tk3 = %.16f \n',k(3))
fprintf('\tk4 = %.16f \n',k(4))
fprintf('\tk5 = %.16f \n',k(5))
fprintf('\tk6 = %16f \n',k(6))
fprintf('The sum of squares is: %.16f \n\n',fval)

[k,resnorm,residual,exitflag,output,lambda,jacobian] =lsqnonlin(@ObjFunc,k0,lb,ub,[],tspan,x0,yexp);
ci = nlparci(k,residual,jacobian);
fprintf('\n\n??ú???lsqnonlin()?????????????:\n')
fprintf('\tk1 = %.16f \n',k(1))
fprintf('\tk2 = %.16f \n',k(2))
fprintf('\tk3 = %.16f \n',k(3))
fprintf('\tk4 = %.16f \n',k(4))
fprintf('\tk5 = %.16f \n',k(5))
fprintf('\tk6 = %16f \n',k(6))
fprintf('The sum of squares is: %.16f \n\n',resnorm)

%%  plot figure
tspan = [0 38.46 50 62.5 100 166.666 250 500];
for i=1:5
[t x] = ode45(@KineticsEqs,tspan,x0,[],k);
fy(i,

=x(end,

;
end
Ypre=fy(

;
Yexp=yexp(

;
plot(Ypre,'or',Yexp,'b.')
legend('prediction','experiment')
%----------fmincon fun-------------------
function f = ObjFuncFmincon(k,x0,yexp)
tspan = [0 38.46 50 62.5 100 166.666 250 500];
for i=1:5
[t x] = ode45(@KineticsEqs,tspan,x0,[],k);
fy(i)=sum((yexp(i,

-x(end,

).^2);
end
f = sum(fy) ;
% -------------------lsqnonlin fun-------------------
function f = ObjFunc(k,x0,yexp)
tspan = [0 38.46 50 62.5 100 166.666];
for i=1:5
[t x] = ode45(@KineticsEqs,tspan,x0,[],k);
fy(i,

=yexp(i,

-x(end,

;
end
f =fy(

;
%-----------odefun--------------------------------
function dYdt = KineticsEqs(t,y,k)
k1=k(1);k2=k(2);k3=k(3);k4=k(4);k5=k(5);k6=k(6);
DEC=y(1);PA=y(2);DPC=y(3);EA=y(4);EPC=y(5);
R1=k(1)*DEC*PA-(k(1)/k(2))*EPC*EA;
R2=k(3)*EPC*PA-(k(3)/k(4))*DPC*EA;
R3=k(5)*EPC.^2-(k(5)/k(6))*DPC*DEC;
dy1dt=-R1+R3;
dy2dt=-R1-R2;
dy3dt=R2+R3;
dy4dt=R1+R2;
dy5dt=R1-R2-R3;
dYdt=[dy1dt;dy2dt;dy3dt;dy4dt;dy5dt];

贊一下

回復此樓

11樓2014-04-24 14:50:42

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

justin534

新蟲 (初入文壇)

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

引用回帖:

10樓: Originally posted by 月只藍 at 2014-03-29 15:11:44
可能是方程的問題，也可能是實驗數(shù)據(jù)的問題。
你可以請實習版主dingd 用1stopt軟件做一下，來驗證一下。...

sorry,剛剛的編碼有誤,這個是OK的編碼!還有data是隨意打的,所以只是請你幫我看可不可以RUN!

function ODE_parafit2
%-------data
yexp=[5.14  12.76  45.41  21.52  3.31 1 1;
5.36  13.20  46.94  20.90  3.17 1 1;
5.51  13.47  48.68  20.01  3.08 1 1;
5.78  14.03  48.16  19.34  3.04 1 1;
6.07  14.33  47.80  18.65  2.94 1 1]';
x0=[93.7 6.3 0 0 0];
k0 = [1 1 1 1 1 1];
lb = [0 0 0 0 0 0];
ub = [1 1 1 1 1 1]*1e5;
tspan = [0 38.46 50 62.5 100 166.666 250 500];
%-----------------
[k,fval,flag] =fmincon(@ObjFuncFmincon,k0,[],[],[],[],lb,ub,[],[],tspan,x0,yexp);
fprintf('\n\n使用函數(shù)fmincon()估計得到的參數(shù)值為:\n')
fprintf('\tk1 = %.16f \n',k(1))
fprintf('\tk2 = %.16f \n',k(2))
fprintf('\tk3 = %.16f \n',k(3))
fprintf('\tk4 = %.16f \n',k(4))
fprintf('\tk5 = %.16f \n',k(5))
fprintf('\tk6 = %16f \n',k(6))
fprintf('The sum of squares is: %.16f \n\n',fval)

[k,resnorm,residual,exitflag,output,lambda,jacobian] =lsqnonlin(@ObjFunc,k0,lb,ub,[],tspan,x0,yexp);
ci = nlparci(k,residual,jacobian);
fprintf('\n\n使用函數(shù)lsqnonlin()估計得到的參數(shù)值為:\n')
fprintf('\tk1 = %.16f \n',k(1))
fprintf('\tk2 = %.16f \n',k(2))
fprintf('\tk3 = %.16f \n',k(3))
fprintf('\tk4 = %.16f \n',k(4))
fprintf('\tk5 = %.16f \n',k(5))
fprintf('\tk6 = %16f \n',k(6))
fprintf('The sum of squares is: %.16f \n\n',resnorm)

%%  plot figure
tspan = [0 38.46 50 62.5 100 166.666 250 500];
for i=1:8
[t x] = ode45(@KineticsEqs,tspan,x0,[],k);
fy(i,

=x(end,

;
end
Ypre=fy(

;
Yexp=yexp(

;
plot(Ypre,'or',Yexp,'b.')
legend('prediction','experiment')
%----------fmincon fun-------------------
function f = ObjFuncFmincon(k,x0,yexp)
tspan = [0 38.46 50 62.5 100 166.666 250 500];
for i=1:8
[t x] = ode45(@KineticsEqs,tspan,x0,[],k);
fy(i)=sum((yexp(i,

-x(end,

).^2);
end
f = sum(fy) ;
% -------------------lsqnonlin fun-------------------
function f = ObjFunc(k,x0,yexp)
tspan = [0 38.46 50 62.5 100 166.666];
for i=1:8
[t x] = ode45(@KineticsEqs,tspan,x0,[],k);
fy(i,

=yexp(i,

-x(end,

;
end
f =fy(

;
%-----------odefun--------------------------------

function dydt = KineticsEqs(t,y,k)
k1=k(1);k2=k(2);k3=k(3);k4=k(4);k5=k(5);k6=k(6);
DEC=y(1);PA=y(2);DPC=y(3);EA=y(4);EPC=y(5);
R1=k(1)*DEC*PA-(k(1)/k(2))*EPC*EA;
R2=k(3)*EPC*PA-(k(3)/k(4))*DPC*EA;
R3=k(5)*EPC.^2-(k(5)/k(6))*DPC*DEC;
dy1dt=-R1+R3;
dy2dt=-R1-R2;
dy3dt=R2+R3;
dy4dt=R1+R2;
dy5dt=R1-R2-R3;
dydt=[dy1dt;dy2dt;dy3dt;dy4dt;dy5dt];

贊一下

回復此樓

12樓2014-04-24 15:11:47

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

WULONGLOVE

鐵蟲 (初入文壇)

應助: 0 (幼兒園)
金幣: 111
帖子: 3
在線: 10.9小時
蟲號: 2352853
注冊: 2013-03-16
專業(yè): 建筑學

引用回帖:

2樓: Originally posted by 月只藍 at 2014-03-26 20:00:43
比較粗糙的辦法是求數(shù)值導數(shù)，根據(jù)t CA的數(shù)據(jù)即可求出，如果需要MATLAB代碼參見：http://www.gaoyang168.com/bbs/viewthread.php?tid=7020495；
精確點的方法是先擬合出動力學方程中的各參數(shù)，以第一組數(shù)據(jù)（即T=325.15K這組 ...

大神， lsqnonlin(@ObjFunc,k0,lb,ub,[],tspan,y0,yexp); 中的y0代表啥意思啊

贊一下

回復此樓

13樓2014-12-03 17:15:12

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

月只藍

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

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

引用回帖:

13樓: Originally posted by WULONGLOVE at 2014-12-03 17:15:12
大神， lsqnonlin(@ObjFunc,k0,lb,ub,[],tspan,y0,yexp); 中的y0代表啥意思啊...

預估的初值

回復此樓

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

14樓2014-12-03 18:09:47

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

相關版塊跳轉我要訂閱樓主 liqianmelody 的主題更新

返回列表

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

[考研] 276求調劑。有半年電池和半年高分子實習經(jīng)歷 +3	材料學257求調劑 2026-03-23	4/200	2026-03-23 09:25 by 材料學257求調劑
[考研] 317求調劑 +12	申子申申 2026-03-19	18/900	2026-03-22 22:23 by luoyongfeng
[考研] 293求調劑 +3	濤濤Wjt 2026-03-22	5/250	2026-03-22 22:21 by jiangpengfei
[考研] 328求調劑，英語六級551，有科研經(jīng)歷 +6	生物工程調劑 2026-03-17	10/500	2026-03-22 20:22 by edmund7
[考研] 311求調劑 +6	冬十三 2026-03-18	6/300	2026-03-22 20:18 by edmund7
[考研] 一志愿華中農業(yè)071010，總分320求調劑 +5	困困困困坤坤 2026-03-20	6/300	2026-03-22 17:41 by hxsm
[考研] 求調劑 +6	十三加油 2026-03-21	6/300	2026-03-22 17:00 by i_cooler
[考研] 384求調劑 +3	子系博 2026-03-22	4/200	2026-03-22 11:04 by 搏擊518
[考研] 085600材料與化工306 +4	z1z2z3879 2026-03-21	4/200	2026-03-21 23:44 by ms629
[考研] 306求0703調劑一志愿華中師范 +5	紙魚ly 2026-03-21	5/250	2026-03-21 17:11 by 學員8dgXkO
[考研] 316求調劑 +6	梁茜雯 2026-03-19	6/300	2026-03-21 06:32 by Ecowxq666！
[考研] 華東師范大學-071000生物學-293分-求調劑 +3	研究生何瑤明 2026-03-18	3/150	2026-03-21 01:30 by JourneyLucky
[考研] 288求調劑 +16	于海海海海 2026-03-19	16/800	2026-03-20 22:28 by JourneyLucky
[考研] 329求調劑 +9	想上學吖吖 2026-03-19	9/450	2026-03-20 22:01 by luoyongfeng
[考研] 中南大學化學學碩337求調劑 +3	niko- 2026-03-19	6/300	2026-03-20 21:58 by luoyongfeng
[考研] 261求B區(qū)調劑，科研經(jīng)歷豐富 +3	牛奶很忙 2026-03-20	4/200	2026-03-20 19:34 by JourneyLucky
[考研] 材料學碩318求調劑 +5	February_Feb 2026-03-19	5/250	2026-03-19 23:51 by 23Postgrad
[考研] 0703化學調劑 +5	pupcoco 2026-03-17	8/400	2026-03-19 13:58 by houyaoxu
[考研] 0703化學調劑 +3	妮妮ninicgb 2026-03-17	3/150	2026-03-18 10:29 by macy2011
[考研] 有沒有道鐵/土木的想調劑南林，給自己招師弟中～ +3	TqlXswl 2026-03-16	7/350	2026-03-17 15:23 by TqlXswl

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

24小時熱門版塊排行榜

justin534

justin534

WULONGLOVE

月只藍

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