五点差分法matlab解椭圆型偏微分方程图文.pdfVIP

五点差分法matlab解椭圆型偏微分方程 五点差分法matlab解椭圆型偏微分方程 用差分法解椭圆型偏微分方程 用差分法解椭圆型偏微分方程 -(U + y ) (p *pi-1) ^xsi ( *y) 0 ; 0y1 -(U + y ) (p *pi-1) ^xsi ( *y) 0 ; 0y1 ｘｘ Ｕ ｙ ＝ ｉ ｅ ｎ ｐｉ ｘ ２ ｘｘ Ｕ ｙ ＝ ｉ ｅ ｎ ｐｉ ｘ ２ U(0,y)=sin(pi*y),U(2,y) ^2 n(pi*y); 0=y=1 U(0,y)=sin(pi*y),U(2,y) ^2 n(pi*y); 0=y=1 ＝ｅ ｓｉ ＝ｅ ｓｉ U(x, )=0, U(x,1)=0; 0= x =2 U(x, )=0, U(x,1)=0; 0= x =2 ０ ＜ ＜ ０ ＜ ＜ 先自己去瞧一下关于五点差分法的理论书籍 先自己去瞧一下关于五点差分法的理论书籍 a b : a b : Ｍ ｔｌａ 程序 Ｍ ｔｌａ 程序 unct on p e x y k]=wud ancha enfa (h, ,n,k ax,e ) unct on p e x y k]=wud ancha enfa (h, ,n,k ax,e ) ｉ ［ ｕ ｉ ｆ ｐ ｉ ［ ｕ ｉ ｆ ｐ % g-s % g-s迭代法解五点差分法问题 迭代法解五点差分法问题 % a %ｋ aｘ为最大迭代次数 ｋ ｘ为最大迭代次数 % , x,y , (2-0)/0 1= 00; % , x,y , (2-0)/0 1= 00; ｎ为 方向的网格数 例如 、 ０ ２ ｎ为 方向的网格数 例如 、 ０ ２ % ,p % ｅ为误差,p为精确解 ｅ为误差 为精确解 sy s te p; sy s te p; ｍ ｍ =z ro (n+1, 1); =z ro (n+1, 1); ｕ ｅ ｓ ＋ ｕ ｅ ｓ ＋ =0 ( : )* ; =0 ( : )* ; ｘ ＋ ０ ｍ ｈ ｘ ＋ ０ ｍ ｈ y 0+ (0:n) h; y 0+ (0:n) h; ＝ ＊ ＝ ＊ for ( =1:n+1) for ( =1:n+1) ｉ ｉ u ( ,1)= i (pi*y (i)); u ( ,1)= i (pi*y (i)); ｉ ｓ ｎ ｉ ｓ ｎ u (i, + )=exp ( ) exp (1)*sin (pi y ( )); u (i, + )=exp ( ) exp (1)*sin (pi y ( )); １ １ ＊