2022年一维抛物线偏微分方程数值解法.docx

《2022年一维抛物线偏微分方程数值解法.docx》由会员分享，可在线阅读，更多相关《2022年一维抛物线偏微分方程数值解法.docx（18页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、名师归纳总结 精品学习资料 - - - - - - - - - - - - - - -一维抛物线偏微分方程数值解法（3）上一篇参看 一维抛物线偏微分方程数值解法（2）（附图及 matlab程序）解一维抛物线型方程 理论书籍可以参看孙志忠：偏微分方程数值解法）Ut-Uxx=0, 0x1,0t0 Ux,0=ex, 0=x=1, U0,t=et,U1,t=e1+t, 0t xlabelx;ylabelt;zlabele; title 误差曲面 细心整理归纳 精选学习资料 - - - - - - - - - - - - - - - 第 2 页,共 9 页 - - - - - - - - - 名师归纳总

2、结 精品学习资料 - - - - - - - - - - - - - - -plotx,e plott,e 细心整理归纳 精选学习资料 - - - - - - - - - - - - - - - 第 3 页,共 9 页 - - - - - - - - - 名师归纳总结 精品学习资料 - - - - - - - - - - - - - - -误差较向前欧拉法减小一半但是运行时间较长,约39 秒,而前两次运行只需l 秒左右；,运算量也提高五倍 第 4 页,共 9 页 u p e x t=CN0.01,0.01,100,100; 运行需三分钟左右,误差比前次提高五倍细心整理归纳 精选学习资料 - -

3、 - - - - - - - - - - - - - - - - - - - - - - 名师归纳总结 精品学习资料 - - - - - - - - - - - - - - -u p e x t=CN0.1,0.1,10,10;surfx,t,e u p e x t=CN0.1,0.2,10,5;surfx,t,e 运行需要 2秒；精度仍是挺高的；细心整理归纳 精选学习资料 - - - - - - - - - - - - - - - 第 5 页,共 9 页 - - - - - - - - - 名师归纳总结 精品学习资料 - - - - - - - - - - - - - - -误差仍可以接受

4、此种方法精度高,运算量较大二：用迭代法解线性方程组：Matlab 程序如下：function u e p x t k=CN1h1,h2,m,n,kmax,ep % 解抛物线型一维方程 C-N 格式 （ Ut-aUxx=fx,t,a0 %用g-s 高斯 - 赛德尔 迭代法解 %kmax为最大迭代次数%m,n 为x,t方向的网格数,例如（2-0 ）/0.01=200;%e为误差, p为精确解 syms temp ; u=zerosn+1,m+1; x=0+0:m*h1; t=0+0:n*h2; for i=1:n+1 ui,1=expti; ui,m+1=exp1+ti; end for i=1:

5、m+1 u1,i=expxi; endfori=1:n+1 第 6 页,共 9 页 细心整理归纳 精选学习资料 - - - - - - - - - - - - - - - - - - - - - - - - 名师归纳总结 精品学习资料 - - - - - - - - - - - - - - -for j=1:m+1 fi,j=0; end end a=zerosn,m-1;r=h2/h1*h1; %此处 r=a*h2/h1*h1）；a=1.fork=1:kmaxfori=1:nforj=2:m temp=r/2*ui,j-1+1-r*ui,j+r/2*ui,j+1+h2*fi,j+r/2*ui

6、+1,j-1+r/2*ui+1,j+1/1+r; ai+1,j=temp-ui+1,j*temp-ui+1,j; ui+1,j=temp;%此处留意是 ui+1,j,而不是 ui+1,j+1%end end ai+1,j=sqrtai+1,j;ifkkmax;breakend if maxmaxaep break ; end end for i=1:n+1 for j=1:m+1 pi,j=expxj+ti; ei,j=absui,j-pi,j; end end u e p x t k=CN10.1,0.005,10,200,10000,1e-10; 运行速度： 1秒 迭代次数 k = 81

7、surfx,t,e 细心整理归纳 精选学习资料 - - - - - - - - - - - - - - - 第 7 页,共 9 页 - - - - - - - - - 名师归纳总结 精品学习资料 - - - - - - - - - - - - - - -其次幅图为三角追逐法解方程作出的图,两者几乎一样；细心整理归纳 精选学习资料 - - - - - - - - - - - - - - - 第 8 页,共 9 页 - - - - - - - - - 名师归纳总结 精品学习资料 - - - - - - - - - - - - - - -由于迭代法速度很快,所以可以将区间分得更小u e p x t k=CN10.01,0.01,100,100,10000,1e-12;surfx,t,e;shading interp; k=6903细心整理归纳 精选学习资料 - - - - - - - - - - - - - - - 第 9 页,共 9 页 - - - - - - - - -