matlab使用pde解三维偏微分方程的程序
t(x,y,0)=300t(x,y,100)=300
t(0,y,z)=t(120,y,z)
t(x,0,z)=t(x,100,z)
0≤ x≤120
0≤ y≤100
0≤ z≤10
写的程序如下,不知道怎么改:
%pdefun1function =pdefun1(x,y,z,dt) c=; f=[-dt,-dt]; s=0
%pdebc1function =pdebc1(xa,ya,ta,xb,yb,tb,z) pa=; qa=; ka=; pb=; qb=; kb=;
%pdeic1function t0=pdeic1(x,y) u0=;
主调动程序:clear all;clc;x=0:1:100;y=0:1:120;z=0:1:100;m=0;sol=pdepe(m,@pdefun1,@pdeic1,@pdebc1,x,y,z);figure('pderunwindows')suft(x,y,z,sol(:,:,:,1)title('the solution of t')xlabel('x')ylabel('y')zlabel('z') 请求高手解答.(∂^2 t)/(∂x^2 )+(∂^2 t)/(∂y^2 )+(∂^2 t)/(∂z^2 )=0
页:
[1]