matlab 偏微分,matlab偏微分方程求解

该楼层疑似违规已被系统折叠 隐藏此楼查看此楼

用matlab求解偏微分方程不会啊。小白一枚，求大神指导！十分感谢！！

完全不会，参考别人的代码改的:

function ABC

x=0:1:90;

t=0:1:5;

m=0;

sol=pdepe(m,@pdefun,@pdeic,@pdebc,x,t);

figure(*numbertitle*,*off*,*name*,*PDE Demo*)

u=sol(:,:,1);

surf(x,t,u)

title(*Solution*)

xlabel(*X*)

ylabel(*T*)

zlabel(*U*)

figure(2);

M=length(t);

xout=linspace(0,90,91);

[uout,du]=pdeval(m,x,u(M,:),xout);

plot(xout,uout);

title(*ss*)

xlabel(*x*)

ylabel(*y*)

%问题描述函数

function [c,f,s]=pdefun (x,t,u,du)

c=1;

f=-u;

d=(3.788*exp(-((x-88.45)/0.8369)^2)+7.415*exp(-((x-82.75)/4.294)^2)+8.875*exp(-((x-80.25)/0.2927)^2)+617700000*exp(-((x-412.1)/82.86)^2));

s=-u*d;

%边界条件

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

kr=2.6026*sin(0.04757*x)-1.5332*cos(0.04757*x)+1.8762*sin(0.09514*x)-0.2125*cos(0.09514*x)+1.2704*sin(0.14271*x)+1.015*cos(0.14271*x)+0.2373*sin...

(0.19028*x)+1.1607*cos(0.19028*x)-0.4995*sin(0.23785*x)+0.7164*cos(0.23785*x)-0.526*sin(0.28542*x)+0.2621*cos(0.28542*x)-0.105*sin(0.33299*x)-0.0035*cos...

(0.33299*x)+0.0062*sin(0.38056*x)-0.1477*cos(0.38056*x);

jifen=int(kr*u*exp(-((x-15)/2))*(x-15)^5,15,49); %积分这里有问题。

p=1.46/3840*jifen;

pl=ul-p;

ql=0;

pr=ur;

qr=0;

%初值条件

function u0=pdeic(x)

u0=2.6026*sin(0.04757*x)-1.5332*cos(0.04757*x)+1.8762*sin(0.09514*x)-0.2125*cos(0.09514*x)+1.2704*sin(0.14271*x)+1.015*cos(0.14271*x)+0.2373*sin...

(0.19028*x)+1.1607*cos(0.19028*x)-0.4995*sin(0.23785*x)+0.7164*cos(0.23785*x)-0.526*sin(0.28542*x)+0.2621*cos(0.28542*x)-0.105*sin(0.33299*x)-0.0035*cos...

(0.33299*x)+0.0062*sin(0.38056*x)-0.1477*cos(0.38056*x);

求大神啊！！！谢谢！！