偏微分方程在图一中,下面是编程的代码
function pde
m=0;n=101;x=linspace(0,0.34,n);t=linspace(0,200,n);
sol=pdepe(m,@klfunPDE,@klfunIC,@klfunBC,x,t);
u=sol(:,:,1);
surf(x,t,u)
xlabel('z'),ylabel('t'),zlabel('X')
%------------------------------------------
function [c,f,s]=klfunPDE(x,t,u,du)
c=457.5/0.0296;f=du;s=-1.6/0.0296*du+2.1/0.0296;
%------------------------------------------
function u0=klfunIC(x)
u0=0;
%------------------------------------------
function [pa,qa,pb,qb]=klfunBC(xa,ua,xb,ub,t)
pa=ua;qa=0;pb=0;qb=1;
结果如下在下面的图片
接下来怎么算Xe= 及过程的编程是什么?Qm=15L/h,m0=200g
![matlab的求解怎么求解Xe、Xe-t-z的图像及他们函数关系式]()
偏微分方程
![matlab的求解怎么求解Xe、Xe-t-z的图像及他们函数关系式-1]()
x-t-z 的图
![matlab的求解怎么求解Xe、Xe-t-z的图像及他们函数关系式-2]()
Xe 的表达式 |