Olivier1 发表于 2015-1-27 22:58

matlab使用pde解三维偏微分方程的程序

t(x,y,0)=300
t(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]
查看完整版本: matlab使用pde解三维偏微分方程的程序