当前位置:   article > 正文

一维线性元有限体积法MATLAB程序 (1D linear finite volume element method)

一维线性元有限体积法MATLAB程序 (1D linear finite volume element method)
  • We seek the solution between Ω = ( 0 , π ) \Omega=(0,\pi) Ω=(0,π);
  • Model Problem:
    − d d x ( p d u d x ) + r d u d x + q u = f -\frac{d}{dx}(p\frac{du}{dx})+r\frac{du}{dx}+qu=f dxd(pdxdu)+rdxdu+qu=f
    MATLAB code:
%一维一次有限体积法解两点边值精度阶测试
%修改版

iter=3;
pde.left=0;%左边界
pde.right=pi;%右边界
pde.Alpha=0;%左边值条件Dirichlet
pde.Beta=-1;%右边值条件Neumann
pde.p=@(x) x;
pde.r=@(x)  x;
pde.q=@(x) x;
pde.f=@(x) 2*x.*sin(x)-cos(x)+x.*cos(x);
pde.realu=@(x) sin(x);
pde.Drealu
  • 1
  • 2
  • 3
  • 4
  • 5
  • 6
  • 7
  • 8
  • 9
  • 10
  • 11
  • 12
  • 13
声明:本文内容由网友自发贡献,不代表【wpsshop博客】立场,版权归原作者所有,本站不承担相应法律责任。如您发现有侵权的内容,请联系我们。转载请注明出处:https://www.wpsshop.cn/w/繁依Fanyi0/article/detail/902516
推荐阅读
相关标签
  

闽ICP备14008679号