連串—平行及連串反應(yīng)的等溫優(yōu)化畢業(yè)論文-資料下載頁(yè)
2025-06-26 15:10本頁(yè)面
【正文】 yd=yc。 h=ba。 c=a+r2*h。 [z]=interp1(XP,FXP,c,39。spline39。)。%f(c)。 yc=z。 else a=c。 ya=yc。 c=d。 yc=yd。 h=ba。 d=a+r1*h。 [z]=interp1(XP,FXP,d,39。spline39。)。%f(d)。 yd=z。 end A(k)=a。B(k)=b。C(k)=c。D(k)=d。 end dp=abs(ba)。 dy=abs(ybya)。 p=a。 yp=ya。 if(ybya) p=b。 yp=yb。 end G=[A39。 C39。 D39。 B39。]。 S=[p,yp]。 E=[dp dy]。附計(jì)算程序2function CSTRTempOpt_lun%計(jì)算連串反應(yīng)166。ABC使得B的濃度最大的溫度分布global T dt Ci0%溫度，時(shí)間間隔，系統(tǒng)物質(zhì)初始濃度n=20%計(jì)算溫度，時(shí)間間隔的數(shù)目tspan =[:1/n:]。dt=1/nTemp0 =398:(398298)/n:298lb = ones(size(T))*298。ub = ones(size(T))*398。A = zeros(length(tspan),length(tspan))。%使得T1小于T2，T2小于T3…..Tn1小于Tn，即使得溫度由高%到低的%約束b = zeros(length(tspan),1)。 %使得T1小于T2，T2小于T3…..Tn1小于Tn，即使得溫度由高%到低的約束for i = 1:length(tspan)1 A(i,i) = 1。 A(i,i+1) = 1。endA=A。b=b。Ci0=[1 0 0][T,fval,exitflag] = fmincon(@ObjFunc,Temp0,A,b,[],[],lb,ub)plot(tspan,T,39。o39。)xlabel(39。time/s39。)ylabel(39。Temperature/161。230。39。)stairs(tspan,T,39。39。)TempCal=T%x_dist = x_distCal(TempCal)%figureplot(tspan,x_dist(:,2),39。o39。)xlabel(39。Time/s39。)ylabel(39。c_239。)T=TempCal(1)%[t,x] = ode45(@C1Equation,[0:dt/10:dt],Ci0)%% delta=。epsilon=。for i=1:length(Temp0) T=Temp0(i)。 [t,x] = ode45(@C1Equation,[0:1/10:1],Ci0)。 [S,E,G]=golden_point_opt_max(t(1),t(end),delta,epsilon,x(:,2),t)。 time_opt(i)=S(1)。 X2_max(i)=S(2)。endX2_max=X2_max[S,E,G]=golden_point_opt_max(Temp0(1),Temp0(end),delta,epsilon,X2_max,Temp0)。Temp_opt=S(1)x2_max=S(2)figureplot(Temp0,X2_max,39。o39。)hold onplot(Temp_opt,x2_max,39。rs39。)%c2最大時(shí)的溫度xlabel(39。Temperature/K39。)ylabel(39。c_2_m_a_x39。)grid on%function f = ObjFunc(Temp)% 目標(biāo)函數(shù)global T dt Ci0x0=Ci0。for i=1:length(Temp) T=Temp(i)。 [t,x] = ode45(@C1Equation,[0:dt/10:dt],x0)。 x0=x(end,:)。endf=x(end,2)。 % 目標(biāo)函數(shù)使得B的濃度最大%function x_dist = x_distCal(TempCal)% global T dt Ci0x0=Ci0。for i=1:length(TempCal) T=TempCal(i)。 [t,x] = ode45(@C1Equation,[0:dt/10:dt],x0)。 x0=x(end,:)。 x_dist(i,:)=x(end,:)。endfunction dxdt = C1Equation(t,x) %微分方程global T Ea1=2500。Ea2=5000。A1=4*10^3。A2=620*10^3。k1=A1*exp(Ea1./T)。k2=A2*exp(Ea2./T)。dxdt(1) =k1*x(1)*x(1)。dxdt(2) =k1*x(1)*x(1)k2*x(2)。dxdt(3) =k2*x(2)。dxdt = dxdt39。26
高考資料相關(guān)推薦
鄂ICP備17016276號(hào)-1