matlab求解pde方程的解析解,pde函数解一阶偏微分方程的问题

(1)形如 du/dt + a*du/dx = v*du2/dx2(二次偏导数) 的一阶偏微分方程可用pde函数求解

如下例：

0

编程如下：

function pdexvv

m = 0;

x = linspace(0,1,101);

t = linspace(0,0.00006,501);

sol = pdepe(m,@pdexvvpde,@pdexvvic,@pdexvvbc,x,t);

u = sol(:,:,1);

figure

plot(x,u(end,:),'-r')

title('Numerical solution and exact solution.')

xlabel('x')

ylabel('u(x,end)')

%----------------------------------------------

function [c,f,s] = pdexvvpde(x,t,u,DuDx)

c = 1;

f = DuDx;

s = -10000.*DuDx ;

%---------------------------------------------

function u0 = pdexvvic(x)

u0 = 0;

%----------------------------------------------

function [pl,ql,pr,qr] = pdexvvbc(xl,ul,xr,ur,t)

p