close all;
clear all;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%具体实例%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
a=-1;b=1;c=-1;d=1;
n=6;m=5;TOL=1e-10;
ITMAX=500;
%f=inline('x*exp(y)','x','y');
%ga=inline('0','x','y');gb=inline('2*exp(y)','x','y');
%gc=inline('x','x','y');gd=inline('exp(1)*x','x','y');
f=inline('2*pi^2*sin(pi*x)*sin(pi*y)','x','y');
ga=inline('0','x','y');gb=inline('0','x','y');
gc=inline('0','x','y');gd=inline('0','x','y');
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%程序部分%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
h=(b-a)/n;
k=(d-c)/m;
x=linspace(a,b,n+1);
x=x(2:n);
y=linspace(c,d,m+1);
y=y(2:m);
u=zeros(n-1,m-1);
lmd=h^2/k^2;
mu=2*(1+lmd);
xx=[];
yy=[];
uu=[];
ture=[];
l=1;
while (l<=ITMAX)
z=(-h*h*f(x(1),y(m-1))+ga(a,y(m-1))+lmd*gd(x(1),d)+...
lmd*u(1,m-2)+u(2,m-1))/mu;
NORM=abs(z-u(1,m-1));
u(1,m-1)=z;
for i=2:n-2
z=(-h*h*f(x(i),y(m-1))+lmd*gd(x(i),d)+u(i-1,m-1)+...
u(i+1,m-1)+lmd*u(i,m-2))/mu;
if(abs(u(i,m-1)-z)>NORM)
NORM=abs(u(i,m-1)-z);
end
u(i,m-1)=z;
end
z=(-h*h*f(x(n-1),y(m-1))+gb(b,y(m-1))+...
lmd*gd(x(n-1),d)+u(n-2,m-1)+lmd*u(n-1,m-2))/mu;
if(abs(u(n-1,m-1)-z)>NORM)
NORM=abs(u(n-1,m-1)-z);
end
u(n-1,m-1)=z;
for j=m-2:-1:2
z=(-h*h*f(x(1),y(i))+ga(a,y(j))+lmd*u(1,j+1)+...
lmd*u(1,j-1)+u(2,j))/mu;
if(abs(u(1,j)-z)>NORM)
NORM=abs(u(1,j)-z);
end
u(1,j)=z;
for i=2:n-2
z=(-h*h*f(x(i),y(j))+u(i-1,j)+lmd*u(i,j+1)+...
u(i+1,j)+lmd*u(i,j-1))/mu;
if( abs(u(i,j)-z)>NORM)
NORM=abs(u(i,j)-z);
end
u(i,j)=z;
end
z=(-h*h*f(x(n-1),y(j))+gb(b,y(j))+u(n-2,j)+...
lmd*u(n-1,j+1)+lmd*u(n-1,j-1))/mu;
if(abs(u(n-1,j)-z)>NORM)
NORM=abs(u(n-1,j)-z);
end
u(n-1,j)=z;
end
z=(-h*h*f(x(1),y(1))+ga(a,y(1))+lmd*gc(x(1),c)+...
lmd*u(1,2)+u(2,1))/mu;
if(abs(u(1,1)-z)>NORM)
NORM=abs(u(1,1)-z);
end
u(1,1)=z;
for i=2:n-2
z=(-h*h*f(x(i),y(1))+lmd*gc(x(i),c)+...
u(i-1,1)+lmd*u(i,2)+u(i+1,1))/mu;
if(abs(u(i,1)-z)>NORM)
NORM=abs(u(i,1)-z);
end
u(i,1)=z;
end
z=(-h*h*f(x(n-1),y(1))+gb(b,y(1))+lmd*gc(x(n-1),c)+...
u(n-2,1)+lmd*u(n-1,2))/mu;
if(abs(u(n-1,1)-z)>NORM)
NORM=abs(u(n-1,1)-z);
end
u(n-1,1)=z;
if(NORM<=TOL)
for i=1:n-1
for j=1:m-1
xx=[xx,x(i)];
yy=[yy,y(j)];
uu=[uu,u(i,j)];
%tur(i,j)=x(i)*exp(y(j));
tur(i,j)=-sin(pi*x(i))*sin(pi*y(j));
ture=[ture,tur(i,j)];
end
end
break;
end
l=l+1;
end
re=[xx' yy' uu' ture']
%%%%%%%%%%%%%%%%%%%%%%%问题%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
[Y,X]=meshgrid(y,x);
mesh(X,Y,u)
MATLAB源代码-自学Matlab必备的60个小程序代码.zip
版权申诉
104 浏览量
2023-08-27
23:34:02
上传
评论
收藏 32KB ZIP 举报
MATLAB源代码-自学Matlab必备的60个小程序代码.zip （60个子文件） 
Matlab 
matlab经典源码-stguass.m 885B matlab经典源码-改进平方根法.m 761B matlab经典源码-hermite.m 584B matlab经典源码-qgauss.m 920B matlab经典源码-cgls.m 666B matlab经典源码-sp.m 213B matlab经典源码-fun.m 54B matlab经典源码-fc.m 55B matlab经典源码-RK.m 435B matlab经典源码-wavequ.m 703B matlab经典源码-ft.m 33B matlab经典源码-erfen.m 444B matlab经典源码-gexian.m 204B matlab经典源码-fg.m 38B matlab经典源码-diedai.m 207B matlab经典源码-Euler1.m 656B matlab经典源码-Newton插值.m 341B matlab经典源码-giassl.m 787B matlab经典源码-Untitled9.m 324B matlab经典源码-dNeville.m 97B matlab经典源码-tranBB.m 358B matlab经典源码-dNewton.m 59B matlab经典源码-fitfun.m 2KB matlab经典源码-fh.m 30B matlab经典源码-LR.m 751B matlab经典源码-sanci.m 814B matlab经典源码-hothch.m 350B matlab经典源码-Taylor解常微分方程.m 416B matlab经典源码-追赶法求解三对角方程.m 352B matlab经典源码-sc.m 469B matlab经典源码-poisson.m 3KB matlab经典源码-Euler解常微分方程.m 655B matlab经典源码-newton.m 154B matlab经典源码-Untitled.m 2KB matlab经典源码-xycf.m 1KB matlab经典源码-Langrange.m 424B matlab经典源码-CN.m 684B matlab经典源码-hotqch.m 491B matlab经典源码-poiss.m 3KB matlab经典源码-gauss.m 1KB matlab经典源码-df.m 37B matlab经典源码-平方根法.m 729B matlab经典源码-seidel.m 188B matlab经典源码-nihe1.m 348B matlab经典源码-fun1.m 54B matlab经典源码-trisys.m 211B matlab经典源码-lbg.m 430B matlab经典源码-ytiao.m 399B matlab经典源码-fx.m 432B matlab经典源码-sanci2.m 866B matlab经典源码-jtf.m 366B matlab经典源码-fitpt.m 2KB matlab经典源码-zjf.m 1KB matlab经典源码-nihe.m 334B matlab经典源码-euler.m 364B matlab经典源码-Neville插值.m 242B matlab经典源码-trapr1.m 140B matlab经典源码-yangtiao.m 427B matlab经典源码-dEuler.m 81B matlab经典源码-jacobi.m 177B资源评论
