偏微分的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 )