二维差分&紧致差分(matlab求解泊松方程)

二维差分&紧致差分(matlab求解泊松方程)

ID:36447702

大小:50.50 KB

页数:8页

时间:2019-05-10

二维差分&紧致差分(matlab求解泊松方程)_第1页
二维差分&紧致差分(matlab求解泊松方程)_第2页
二维差分&紧致差分(matlab求解泊松方程)_第3页
二维差分&紧致差分(matlab求解泊松方程)_第4页
二维差分&紧致差分(matlab求解泊松方程)_第5页
资源描述:

《二维差分&紧致差分(matlab求解泊松方程)》由会员上传分享,免费在线阅读,更多相关内容在应用文档-天天文库

1、实用文案%%%%真解u=sin(pi*x)*sin(pi*y)%%%%%%%方程-Laplace(u)=f%%%%%%%%%%f=2*pi^2*sin(pi*x)*sin(pi*y)%%%%%%%%%%differencecodeforellipticequationswithconstantcoefficient%%%%%clear%clcN=20;h=1/N;S=h^2;x=0:h:1;y=0:h:1;%%%StiffmatrixA=zeros((N-1)^2,(N-1)^2);fori=1A(i,i)=4/h^2;A(i,i+1)=-1/h^2;A(i,i

2、+(N-1))=-1/h^2;endfori=N-1A(i,i-1)=-1/h^2;A(i,i)=4/h^2;A(i,2*i)=-1/h^2;%A(i,i+(N-1))=-1/h^2endfori=(N-2)*(N-1)+1A(i,i-(N-1))=-1/h^2;A(i,i)=4/h^2;A(i,i+1)=-1/h^2;endfori=(N-1)^2A(i,i-(N-1))=-1/h^2;A(i,i)=4/h^2;A(i,i-1)=-1/h^2;endforn=2:N-2i=(N-2)*(N-1)+n;A(i,i-(N-1))=-1/h^2;A(i,i-1)=-

3、1/h^2;A(i,i)=4/h^2;A(i,i+1)=-1/h^2;endfori=2:N-2A(i,i-1)=-1/h^2;A(i,i)=4/h^2;A(i,i+1)=-1/h^2;标准文档实用文案A(i,i+(N-1))=-1/h^2;endform=1:N-3i=m*(N-1)+1;A(i,i-(N-1))=-1/h^2;A(i,i)=4/h^2;A(i,i+1)=-1/h^2;A(i,i+(N-1))=-1/h^2;endform=2:N-2i=m*(N-1);A(i,i-(N-1))=-1/h^2;A(i,i-1)=-1/h^2;A(i,i)=4/h

4、^2;A(i,i+(N-1))=-1/h^2;end%form=1:N-3%i=m*(N-1)+(N-1);%A(i,i-(N-1))=-1/h^2;%A(i,i-1)=-1/h^2;%A(i,i)=4/h^2;%A(i,i+(N-1))=-1/h^2;%endform=1:N-3forn=2:N-2i=m*(N-1)+n;A(i,i-(N-1))=-1/h^2;A(i,i-1)=-1/h^2;A(i,i)=4/h^2;A(i,i+1)=-1/h^2;A(i,i+(N-1))=-1/h^2;endend%%%RighttermF=zeros((N-1)^2,1)

5、;form=0:N-2forn=1:N-1i=m*(N-1)+n;F(i)=2*pi^2*sin(pi*n*h)*sin(pi*(m+1)*h);endEndU=zeros((N-1)^2,1);U=AF;U1=zeros((N+1)^2,1);标准文档实用文案form=1:N-1U1(m*(N+1)+2:m*(N+1)+N)=U((m-1)*(N-1)+1:m*(N-1));endplot(x,y,U1,’-k’)U-exact=inline(‘sin(xx(:,1).*sin(xx(:,2))’,’xx’);%真解plot(x,y,U-exact,’-r’

6、)Err=max(abs(U1-U-exact));%最大模L2-err=norm(U1-U-exact)*h;%L2errorErrL2-err标准文档实用文案%%%%真解u=sin(pi*x)*sin(pi*y)%%%%%%%方程-Laplace(u)=f%%%%%%%%%%f=2*pi^2*sin(pi*x)*sin(pi*y)%%%%%%%%%%compactdifferencecodeforellipticequationswithconstantcoefficient%%%%%clear%clcN=20;h=1/N;S=h^2;x=0:h:1;y=0

7、:h:1;%%%StiffmatrixA=zeros((N-1)^2,(N-1)^2);B=zeros((N-1)^2,(N-1)^2);F=zeros((N-1)^2,1);fori=1A(i,i)=10/(3*h^2);A(i,i+1)=-2/(3*h^2);A(i,i+(N-1))=-2/(3*h^2);A(i,i+(N-1)+1)=-1/(6*h^2);B(i,i)=25/36;B(i,i+1)=5/72;B(i,i+(N-1))=5/72;B(i,i+(N-1)+1)=1/144;endfori=N-1A(i,i-1)=-2/(3*h^2);A(i,i

8、)=10/(3*h^2)

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。