天津大学化工数学课件四常微分方程数值解.PPTVIP

4.1 引言 ;问题1 在间歇蒸馏釜中，最初含有25mol的正辛烷和75mol的正庚烷，如果在1atm下进行操作，当最后釜中剩有溶液的总摩尔数为10mol时，试求正庚烷的摩尔分率。;总的物料平衡 ;问题2 在一个与外界没有质量交换的封闭系统中，三个组分A,B,C浓度分别为C1、C2、C3。当系统在特定频率光照下发生反应，求每种物质浓度随时间变化规律。;;问题3 一均匀细杆稳态导热问题。设其杆的一端恒温为T1，另一端绝热，试确定杆温度沿杆的分布情况。 ;对于一阶初值问题：;数值解就是求取在节点处的函数近似值，;1、显式尤拉法;分析解为 ;（2）误差分析;（3）收敛性和稳定性;一般以下列初值问题讨论方法的稳定性;2、隐式尤拉法;4、改进尤拉法;4.2.2 龙格库塔法;例2：用四阶龙格库塔法解例4-1; 例4-2 在一个与外界没有质量交换的封闭系统中，三个组分A,B,C浓度分别为C1、C2、C3。当系统在特定???率光照下发生反应，求每种物质浓度随时间变化规律。;;…………………;回顾用四阶龙格－库塔法解常微分方程初值问题：;常微分方程组初值问题求解公式：;取h＝0.05，求y1(0.1)，由y2(0.1)的值，计算过程保留4位小数。;设 Cji≈uji;例4－3 用四阶龙格—库塔法解 常微分方程组初值问题 ;4.2.4 高阶常微分方程求解;将