01应用实例课程.pptVIP

理论分析 理论分析 理论分析 MATLAB求解方法 作图法（1） 作图法（2） 源程序（1）——作图法 %主程式 function cstr_nes clear clc % koCAO=exp(20); ER=10000; tau=0.25; H=250; To=450; % %check for solution % t=450:700; y1=(t-450)/250; a=.25*exp(20-10000./t); b=1./a+2; y2=.5*(b-sqrt(b.^2-4)); plot(t,y1,r,t,y2,c) xlabel(temperature) ylabel(conversion rate) title(mass balance vs. energy balance) % %find solutions:use fsolve(共有三组解) 作图法结果 源程序（2）——数值求解 %find solutions:use fsolve(共有三组解) % X1=fsolve(@fun_cstr, [0.9 650]); %0.9是指转化率x的初始猜测值; % 650是指出口温度T的初始猜测值 X2=fsolve(@fun_cstr, [0.3 500]); X3=fsolve(@fun_cstr,[0.1 450]); % %print out % disp(sprintf(\n the possible solutions are:\n)) disp( conv. temp.) disp([X1]) disp([X2]) disp([X3]) %函数式:fun_cstr.m function y=fun_cstr(X) %传入的X有两个:x,T % global koCAO ER tau H To % %参数表: % koCAO= e^20 % ER=E/R % tau=V/v % H=(CAO*(-dH)/(rho*Cp)) 单位[k] % x=X(1); T=X(2); y(1)=x-koCAO*(exp(-ER/T))*((1-x).^2)*tau; y(2)=To-T+koCAO*(exp(-ER/T))*((1-x).^2)*tau*H; 数值计算结果 the possible solutions are: conv. temp. 0.8363 659.0833 0.2804 520.0955 0.0408 460.1964 作业 非线性方程组的同伦算法——简介 一、Newton迭代法的局限 求解各种方程需要用不同的方法，因没有万能的解方程方法。一般非线性方程组的求解方法（如Newton法），在初值选取不好时，往往得不到方程的解，而且在方程组有多组根时，通常只能得到一组解。 % 画第一条曲线 ezplot(x^2-y-1,[0 2]) hold on % 画第二条曲线 ezplot((x-1)^2+(y-1)^2-0.25,[0 2]) grid hold off Matlab求解 f1=inline([P(1)^2-P(2)-1;(P(1)-1)^2+(P(2)-1)^2-0.25],P); [x1,y1]=fsolve(f1,[0,0]) [x2,y2]=fsolve(f1,[5,5]) Optimization terminated: first-order optimality is less than options.TolFun. x1 = 1.2523 0.5683 y1 = 1.0e-008 * 0.0396 0.3343 Optimization terminated: first-order optimality is less than options.TolFun. x2 = 1.4717 1.1659 y2 = 1.0e-008 * 0.0602 0.5060 同伦算法 同伦算法的基本思想：从容易求解的方程组开始，逐步过渡到原方程组的求解，从而得到问题的解。 同伦算法 有多种方法可用来求取同伦方程组所要求的解，弧长法和微分方程法是其中的两种典型方法。 弧长法的基本思路：如果同伦函数连续可微，导数矩阵可逆，据隐函数的定理，可确保有一路径联接初值和期望的解。选取合适的函数，在从0到1变化时，使构造的同伦计算保持在路径上。 弧长法求取同伦方程组的基本步骤(1) 弧长法求取同伦方程组的基本步骤(2) 弧长法求取同伦方程组