W11 常微分方程_2012.pdf

function Excer09_3 T=100:100:1000; Cp1=[29.38 29.88 30.42 30.98 31.57 32.15 32.73 33.29 33.82 34.31]; Cp2=[30.04 29.08 30.18 30.14 32.27 31.79 32.97 32.56 34.24 35.27]; Cp3=[28.52 29.79 31.41 31.18 31.16 32.81 32.38 34.26 34.72 33.69]; beta1=polyfit(T,Cp1,2); beta2=polyfit(T,Cp1,3); Tcal=100:10:1000; plot(T,Cp1, bo,Tcal,polyval(beta1,Tcal),r:,Tcal,polyval(beta2,Tcal),m- ) title(Polyfit) [beta3,r]=nlinfit(repmat(T,1,3),[Cp1,Cp2 ,Cp3],@fun1,beta1); [beta4,r2]=nlinfit(repmat(T,1,3),[Cp1,Cp2,Cp3],@fun2,beta2); [beta5,r3]=lsqcurvefit(@fun1,beta1,repmat(T,1,3),[Cp1,Cp2,Cp3]) norm1=norm(r) norm2=norm(r2) if norm1norm2 fprintf(Cp=%.4f+%.4fT+%.8fT^2\n,beta3(end:-1:1)) else fprintf(Cp=%.4f+%.4fT+%.4fT^2+%.4fT^3\n,beta4(end:-1:1)) end figure plot(repmat(T,1,3),[Cp1,Cp2,Cp3],bo,Tcal,fun1(beta3,Tcal),r:) hold on plot(Tcal,fun2(beta4,Tcal), m-) title(Nlinfit) function y=fun1(beta,T) y=beta(3)+beta(2)*T+beta(1)*T.^2; function y=fun2(beta,T) y=beta(4)+beta(3)*T+beta (2)*T.^2+beta(1)*T.^3; 第6章 第0章绪论 第6章微分方程数值解 • 常微分方程定义 • 初值问题数值解方法 • MATLAB求解初值问题函数：ode45 第6章 第0章绪论 微分方程在化工模型中的应用 • 间歇反应器的计算 • 活塞流反应器的计算 • 全混流反应器的动态模拟 • 定态一维热传导问题 • 逆流壁冷式固定床反应器一维模型 • 固定床反应器的分散模型 • 各种传递过程 第6章 第0章绪论 微分方程的定义 2 d s −mg 6 −1 2 ( ) dt 若设：t0 0 时 s 0, s 1 (6 −2) t 0时 s 0 ⎧0 或设： 6 −3 ⎨ ( ) t 1 s 0 ⎩1 时 •只有一个自变量的微分方程为常微分方程(ODE)； •方程中未知函数导数的最高阶数称为方程的阶； •方程中关于未知函数及其各阶导数均是一次的，则称为线性微 分方程。 •