برنامه نویسی با متلب

حل پروژه های ریاضی پیشرفته مهندسی شیمی - مهندسی فرایند    )با تشکر از مهندس مسافر(

استاد دکتر یعقوبی

^2T/dx^2+d^T2/dy^2=0

T(0,y)=0
T(1,y)=0
T(x,0)=100
T(x,1)=100

كد:

  % jzheats
  clear,clc
  nx=10;ny=10;hx=1/nx;hy=1/ny;
  x=0:hx:nx*hx;
  by0=100;
  byn=100;
  y=0:hy:ny*hy;
  bx0=0;
  bxn=0;
  [T,iter]=HeatS(nx,hx,ny,hy,bx0,bxn,by0,byn);
  T=rot90(T,1);
  iter
  % Graph
  T=flipud(T);
  [c,h]=contour(T);clabel(c,h),colorbar
  title('Temperature Distribution')
  xlabel('Node Number in x Direction')
  ylabel('Node Number in y Direction')
  [m,n]=size(T);
  set(gca,'xtick',1:n,'ytick',1:m)
  axis('equal')
  axis([1 n 1 m])

اين هم خروجيش:

 

ادامه نوشته