数学建模社区-数学中国

标题: matlab使用pde解三维偏微分方程的程序 [打印本页]

作者: 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[/math]

写的程序如下，不知道怎么改：

%pdefun1

function [c,f,s]=pdefun1(x,y,z,dt)

c=[1];

f=[-dt,-dt];

s=0

%pdebc1

function [pa,qa,ka,pb,qb,kb]=pdebc1(xa,ya,ta,xb,yb,tb,z)

pa=[0];

qa=[0];

ka=[300];

pb=[100];

qb=[120];

kb=[300];

%pdeic1

function t0=pdeic1(x,y)

u0=[300 300];

主调动程序：

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')

请求高手解答.[math](∂^2 t)/(∂x^2 )+(∂^2 t)/(∂y^2 )+(∂^2 t)/(∂z^2 )=0

欢迎光临数学建模社区-数学中国 (http://www.madio.net/)