Matlab 求解化工常微分方程和偏微分方程.pptVIP

调用方法 sol=pdepe(m,@pdefun, ,@iCfun @BCfun, xspan,tspan,option,p1,p2..) @Bcfun:为定义边界条件方程 @iCfun:为定义初始条件方程。 其他符号意义同上，注意边界条件必须写成以上形式： 应用策略 Pdepe 内部函数具体应用时，需将实际的偏微分方程对照标准模型，确定c、f、s函数的具体形式及边界条件p、q的具体形式及m 值。 蒸汽入口 流体入口，u，t0 冷凝液出口 流体出口，u，tL 套管动态传热问题 原方程： 代入具体数据，并对长度归一化无量纲处理后，得到以下方程： 对照标准模型，T就是标准模型中的函数u，m=0, a=0,b=1,t0=0, t f = 1 标准模型： 初始条件：零时刻所有位置温度为30，即u0=30 边界条件：已知在零位置处任意时间温度为30，在1位置处，偏导为0 pa=u-30,qa=0;pb=0,qb=1 程序清单 function l51 clc clear all global ua ub alpha = 1.0; % cm/s ua = 30; ub = 30; m = 0; a = 0; b = 1; t0 = 0; tf =1 x = linspace(a,b,11); t = linspace(t0,tf,101); sol = pdepe(m,@PDEfun,@ICfun,@BCfun,x,t); u = sol(:,:,1) % surface plot of the solution figure; surf(x,t,u); title(Numerical solution computed with 11 mesh points.); xlabel(Disuance x); ylabel(Time t); ppp=u % solution profile at t=1 figure; subplot(2,2,1); plot(x,u(1,:),o-); title(Solutions at t = 0.); xlabel(Distance x); ylabel(u(t=0,)); grid on; subplot(2,2,2); 程序清单 plot(x,u(50,:),o-); title(Solutions at t = 0.5.); xlabel(Distance x); ylabel(u(t=0.5)); grid on; subplot(2,2,3); plot(x,u(end,:),o); title(Solutions at t = 1.); xlabel(Distance x); ylabel(u(t=1,)); grid on; % ------------------------------------------------------------------ function [c1,f,s] = PDEfun(x,t,u,Du) c1 = 1 f = 0.001*Du; s = 2*(150-t)-3*Du % ------------------------------------------------------------------ function u0 = ICfun(x) u0 = 30; % ------------------------------------------------------------------ function [pa,qa,pb,qb] = BCfun(xl,ul,xr,ur,t) global ua ub pa = ul-ua; qa = 0; pb = 0 qb =1 计算结果 分析讨论 Matlab 求解化工常微分方程和偏微分方程 方利国 Matlab 求解化工常微分方程和偏微分方程 1、常微分方程（组）求解 1.1 问题描述及Matlab调用命令 1.2 初值问题求解 1.3 边值问题求解 1.4 加权问题求解（自学内容） 2、偏微分方程（组）求解 2.1 问题描述及一维动态PDE方程求解 2 .2 二维求解 1、常微分方程（组）求解 1.1 问题描述及Matlab调用命令 常微分方程：（初值问题） 常微分方程：（两点边值问题） 1、常微分方程（组）求解 1.1 问题描述及Matlab调用命令 微分方程组： 1、常微分方程（组）求解 1.1 问题描述及Matlab调用命令 高级微分方程： 1、常微分方程（组）求解 1