(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