MATLAB编程求解二维泊松方程
![MATLAB编程求解二维泊松方程_第1页](https://file7.zhuangpeitu.com/fileroot7/2021-7/25/04915471-da3c-4b36-9800-b5a4c091799c/04915471-da3c-4b36-9800-b5a4c091799c1.gif)
![MATLAB编程求解二维泊松方程_第2页](/images/s.gif)
![MATLAB编程求解二维泊松方程_第3页](/images/s.gif)
《MATLAB编程求解二维泊松方程》由会员分享,可在线阅读,更多相关《MATLAB编程求解二维泊松方程(3页珍藏版)》请在装配图网上搜索。
1、% 真解 u=sin(pi*x)*sin(pi*y) % 方程 -Laplace(u)=f % f=2*pi2*sin(pi*x)*sin(pi*y) %difference code for elliptic equations with constant coefficient %clear all%clcN=20;h=1/N;S=h2;x=0:h:1;y=0:h:1;% Stiff matrix A=zeros(N-1)2,(N-1)2);for i=1 A(i,i)=4/h2; A(i,i+1)=-1/h2; A(i,i+(N-1)=-1/h2;endfor i=N-1 A(i,i-1
2、)=-1/h2; A(i,i)=4/h2; A(i,2*i)=-1/h2; %A(i,i+(N-1)=-1/h2endfor i=(N-2)*(N-1)+1 A(i,i-(N-1)=-1/h2; A(i,i)=4/h2; A(i,i+1)=-1/h2;endfor i=(N-1)2 A(i,i-(N-1)=-1/h2; A(i,i)=4/h2; A(i,i-1)=-1/h2;endfor n=2:N-2 i=(N-2)*(N-1)+n; A(i,i-(N-1)=-1/h2; A(i,i-1)=-1/h2; A(i,i)=4/h2; A(i,i+1)=-1/h2;endfor i=2:N-2 A
3、(i,i-1)=-1/h2; A(i,i)=4/h2; A(i,i+1)=-1/h2; A(i,i+(N-1)=-1/h2;endfor m=1:N-3 i=m*(N-1)+1; A(i,i-(N-1)=-1/h2; A(i,i)=4/h2; A(i,i+1)=-1/h2; A(i,i+(N-1)=-1/h2;endfor m=2:N-2 i=m*(N-1); A(i,i-(N-1)=-1/h2; A(i,i-1)=-1/h2; A(i,i)=4/h2; A(i,i+(N-1)=-1/h2;end% for m=1:N-3% i=m*(N-1)+(N-1);% A(i,i-(N-1)=-1/h
4、2;% A(i,i-1)=-1/h2;% A(i,i)=4/h2;% A(i,i+(N-1)=-1/h2;% endfor m=1:N-3 for n=2:N-2 i=m*(N-1)+n; A(i,i-(N-1)=-1/h2; A(i,i-1)=-1/h2; A(i,i)=4/h2; A(i,i+1)=-1/h2; A(i,i+(N-1)=-1/h2; endend% Right termF=zeros(N-1)2,1);for m=0:N-2 for n=1:N-1 i=m*(N-1)+n; F(i)=2*pi2*sin(pi*n*h)*sin(pi*(m+1)*h); endend%U=z
5、eros(N-1)2,1);u1=AF; u=zeros(N+1)2,1);for m=1:N-1 u(m*(N+1)+2:m*(N+1)+N)=u1(m-1)*(N-1)+1:m*(N-1);endU=zeros(N+1,N+1);for m=1:N+1 U(m,:)=u(m-1)*(N+1)+1:m*(N+1);endsurf(x,y,U)u_exact=zeros(N+1)2,1);for m=0:N for n=1:N+1 i=m*(N+1)+n; u_exact(i)=sin(pi*n*h)*sin(pi*m*h); endendU_exact=reshape(u_exact,N+1,N+1);subplot(1,2,)err=max(abs(u-u_exact);l2_err=norm(u-u_exact)*h;errl2_err
- 温馨提示:
1: 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
2: 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
3.本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
5. 装配图网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 新人版英语八年级下册Unit5总复习ppt课件
- 新人教部编版一年级语文上第五单元ppt课件(全套)
- 高鸿业经济学基础第十五章-总需求-总供给模型-授课-河北工大宋建林课件
- 新人教版高中数学《等差数列前n项和》课件
- 新人教部编版五年级语文上册第六单元测试卷课件
- 高鸿业微观经济学课件第4章生产论
- 高鸿业--微观经济学-第一章课件
- 新人教版部编本五年级下册语文13 人物描写一组 ppt课件
- 新人教版高中化学必修第一册——电解质的电离ppt课件
- 新人教版部编教材二年级下册第一单元3《贝的故事》优质课教学ppt课件
- 高风险作业培训讲义_002
- 新人教版语文三年级下册第五单元全套ppt课件部编版
- 新人教版英语八年级上册第二单元全部ppt课件
- 《走一步再走一步》重点课件
- 新人教版语文一年级上册:识字1《天地人》课件