• 设为首页
  • 点击收藏
  • 手机版
    手机扫一扫访问
    迪恩网络手机版
  • 关注官方公众号
    微信扫一扫关注
    迪恩网络公众号

利用Matlab求解Laplace方程

原作者: [db:作者] 来自: [db:来源] 收藏 邀请

相关链接:http://wenku.baidu.com/view/f17e57e8856a561252d36f5e.html

问题描述:

求解泊松方程 
                -Δu=1     
求解区域为单位圆盘,边界条件为在圆盘边界上u=0 

下面求它的数值解,编写程序如下:

%(1)问题定义 
g=\'circleg\';        %单位圆 
b=\'circleb1\';       %边界上为零条件 
c=1;a=0;f=1; 
%(2)产生初始的三角形网格 
[p,e,t]=initmesh(g);  
%(3)迭代直至得到误差允许范围内的合格解 
error=[]; err=1; 
while err > 0.01, 
[p,e,t]=refinemesh(g,p,e,t); 
u=assempde(b,p,e,t,c,a,f);  %求得数值解 
exact=(1-p(1,:).^2-p(2,:).^2)/4; 
err=norm(u-exact\',inf); 
error=[error err]; 
end 
%结果显示 
subplot(2,2,1),pdemesh(p,e,t); 
subplot(2,2,2),pdesurf(p,t,u) 
subplot(2,2,3),pdesurf(p,t,u-exact\')

显示结果如下:

 

 

 

 

 

 


鲜花

握手

雷人

路过

鸡蛋
该文章已有0人参与评论

请发表评论

全部评论

专题导读
上一篇:
delphi线程使用发布时间:2022-07-22
下一篇:
[delphi]SetWindowsHookExA函数入口处修改发布时间:2022-07-22
热门推荐
阅读排行榜

扫描微信二维码

查看手机版网站

随时了解更新最新资讯

139-2527-9053

在线客服(服务时间 9:00~18:00)

在线QQ客服
地址:深圳市南山区西丽大学城创智工业园
电邮:jeky_zhao#qq.com
移动电话:139-2527-9053

Powered by 互联科技 X3.4© 2001-2213 极客世界.|Sitemap