matlab pdebc,偏微分的matlab数值解法

偏微分的matlab数值解法

偏微分的 MATLAB数值解法偏微分的 MATLAB数值解法•方法一： pdepe函数实现•方法二： pdetool实现•方法三：程序实现方法一： pdepe函数实现• MATLAB提供的 pdede函数，可以求解一般的偏微分方程组，其调用格式 ：• sol =pdepe(m,@pdefun,@pdeic,@pdebc,x,t)方法一： pdepe函数实现@pdefun(函数格式描述 )：[c,f,s]=pdefun(x,t,u,du)方法一： pdepe函数实现@pdebc( 边界条件描述)：[pa,qa,pb,qb]=pdebc(x,t,u,du)a:左边界， b:右边界方法一： pdepe函数实现• @pdeic( 初值条件描述 )：• 形式： u( x， t0) =u0• u0=pdeic( x)• m， x， t对应函数式中参数方法一： pdepe函数实现• 应用实例：方法一： pdepe函数实现• 应用实例：• @pdefun：• function [c,f,s]=pdefun(x,t,u,du)• c=1;• f=400*du;• s=0;• @pdebc：• function [pa,qa,pb,qb]=pdebc(xa,ua,xb,ub,t)• pa=ua;• qa=0;• pb=ub;• qb=0;方法一： pdepe函数实现• @pdeic：• function u0=pdeic(x)• if x<10 • u0=0;• elseif x<30• u0=1;• else• u0=0;• end方法一： pdepe函数实现• x=0:1:40;• t=0:0.01:0.2;• m=0;• sol=pdepe(m,@pdefun,@pdeic,@pdebc,x,t);• b=sol(20,:);• plot(x,b);• title( the solution of u )• xlabel( x )• ylabel( y )• zlabel( u )