ImageVerifierCode 换一换
格式:DOCX , 页数:9 ,大小:97.20KB ,
资源ID:1391103      下载积分:5 文钱
快捷下载
登录下载
邮箱/手机:
温馨提示:
快捷下载时,用户名和密码都是您填写的邮箱或者手机号,方便查询和重复下载(系统自动生成)。 如填写123,账号就是123,密码也是123。
特别说明:
请自助下载,系统不会自动发送文件的哦; 如果您已付费,想二次下载,请登录后访问:我的下载记录
支付方式: 支付宝    微信支付   
验证码:   换一换

加入VIP,省得不是一点点
 

温馨提示:由于个人手机设置不同,如果发现不能下载,请复制以下地址【https://www.wenke99.com/d-1391103.html】到电脑端继续下载(重复下载不扣费)。

已注册用户请登录:
账号:
密码:
验证码:   换一换
  忘记密码?
三方登录: QQ登录   微博登录 

下载须知

1: 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。
2: 试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓。
3: 文件的所有权益归上传用户所有。
4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
5. 本站仅提供交流平台,并不能对任何下载内容负责。
6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

版权提示 | 免责声明

本文(数值上机作业.docx)为本站会员(h****)主动上传,文客久久仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对上载内容本身不做任何修改或编辑。 若此文所含内容侵犯了您的版权或隐私,请立即通知文客久久(发送邮件至hr@wenke99.com或直接QQ联系客服),我们立即给予删除!

数值上机作业.docx

1、 数值实验题一 班级:研开发 14-1 班 姓名:马高强 学号: 2014212090 4.用有限差分方法(五点差分格式)求解正方形域上的 Poisson 方程边值问题 2222 ,3 0 , 10 , 1 , 0 , x , 0 x , 1 1uu f x yxy xyu y u y u u x x 用 MATLAB 语言编写求解线性方程组 Au f 的算法程序,采用下列方法计算,比较计算结果和算法性能,对计算结果给出结论。 ( 1)用 SOR 迭代法求解线性方程组 Au f , 用试算法确定最佳松弛因子 。 ( 2)用块 SOR 迭代法求解线性方程组 Au f , 用试算法确定最佳松弛因子

2、 。 ( 3)用(预条件)共轭斜向量法求解线性方程组 Au f 。 解:首先,剖分求解区域。用平行于坐标轴的直线 1, , , 0 , 1 , , , 11ijx x i h y y j h h i j n nn 然后 用数值微分公式对微分算子进行离散 。 211 222211 2, 2 , , 2 , ,i j i j i jiji j i j i jiju x y u x y u x yux y hxhu x y u x y u x yux y hyh 既有 21 , , 1 ,222, 1 , , 1222,2,i j i j i jiji j i j i jiju u uu xyxhu

3、 u uu xyyh 其中 ,iju 表示 ,iju x y的近似值 。 得到在每个点 ,ijxy上的有限差分方程为 : 2, 1 , 1 , , 1 , 14 = 1 ,i j i j i j i j i j iju u u u u h f i j n 又称为五点差分格式 .其中 =,ij i jf f x y。 在边界上有 0 , + 1 , , 0 , 1= 0 , = 1 1 1 1 ,j n j i i nu u u u i h i h i j n 可得差分格式为 : 2, 1 , 1 , , 1 , 10 , + 1 , , 0 , 14= 1,= 0 , = 1 1 1i j i

4、 j i j i j i j ijj n j i i nu u u u u h f i j nu u u u i h i h 将差分方程组写成矩阵形式为 Au f , 其中 ,n n n n nnA R u f R 。 如果把每一条线上的网格点看作一个组 , 如 1 1 , 1 2 , 1 n , 1 2 1 , 2 2 , 2 n , 2 1 , n 2 , n n , n2 2 21 1 , 1 2 , 1 n , 1 2 1 , 2 2 , 2 n , 2 1 , n 2 , n n , n, , , , , , , , , , , , , , , , , , , , , , ,T T

5、TnT T Tnv u u u v u u u v u u ub h f f f b h f f f b h f f f 对矩阵也相应地分块 1 1 1 12 2 2 2,n n n nA I v bI A v bA v bII A v b 其中 114114 ,114nniiA A R则差分方程等价于 Av b 。 这样就可以用块迭代法进行求解 。 ( 1)用 SOR 迭代法求解线性方程组 Au f , 用试算法确定最佳松弛因子 。 SOR 迭代法: function u,k,t=SOR(n,w) %用 SOR迭代法求解 Poisson方程 %输出结果 : u, k, t 输入 : n, w

6、 h=1/(n+1); f(2:n+1,2:n+1)=3*h2;%定义右端项 u=zeros(n+2,n+2); %对边界进行差分 for i=1:n+2 u(i,1)=(i-1)*h*(1-(i-1)*h); u(i,n+2)=(i-1)*h*(1-(i-1)*h); end e=0.000000001; tic%程序开始的时间 %SOR迭代法过程 for k=1:10000 er=0; for i=2:n+1 for j=2:n+1 ub=u(i,j); u(i,j)=w*(u(i-1,j)+u(i+1,j)+u(i,j+1)+u(i,j-1)+f(i,j)/4)+(1-w)*u(i,j)

7、; er=er+abs(ub-u(i,j); end end if er/n2 goodw,goodk=Bestw(9) goodw = 1.5330 goodk = 33 u,k,t=SOR(9,goodw) k = 33 t = 2.3291e-004 u 的结果如下表 : 将 u 画成三维图 surf(u) ( 2)用块 SOR 迭代法求解线性方程组 Au f , 用试算法确定最佳松弛因子 。 块 SOR 迭代法: function u,k,t=BSOR(n,w) %用块 SOR求解 Poisson问题 ,其中用到追赶法 %输出结果 : u, k, t 输入 : n, w h=1/(n+

8、1); f(2:n+1,2:n+1)=3*h2; u=zeros(n+2,n+2); a=-1*ones(1,n);b=4*ones(1,n);c=-1*ones(1,n); for i=1:n+2%对边界条件进行差分 u(i,1)=(i-1)*h*(1-(i-1)*h); u(i,n+2)=(i-1)*h*(1-(i-1)*h); end tic;%块 SOR迭代法过程 e=0.000000001; for k=1:1000 er=0; for j=2:n+1 ub=u(:,j); d=f(2:n+1,j)+u(2:n+1,j-1)+u(2:n+1,j+1); x=zg(a,b,c,d);%

9、调用追赶法求解 for i=1:n u(i+1,j)=w*x(i)+(1-w)*u(i+1,j); end er=er+norm(ub-u(:,j),1); end if er/n2 goodw,goodk=BBestw(9) goodw = 1.4200 goodk = 25 u,k,t=BSOR(9,goodw) k = 25 t = 0.0156 u 的结果如下表 : 将 u 画成三维图 surf(u) ( 3)用(预条件)共轭斜向量法求解线性方程组 Au f 。 共轭斜量法: function u,k,t=CG(n) %共轭 斜量法求解 Poisson问题 A=zeros(n2,n2)

10、; for i=1:n2%系数 A(i,i)=4; if i+n u,k,t=CG(9) k = 13 t = 0.0019 u 的结果如下表 : 将 u 画成三维图 surf(u) 将三个程序结果进行比较: 对于同一问题 , SOR 迭代方法迭代次数为 33 次,迭代时间为 0.00219s;块 SOR 迭代方法迭代次数为 25 次,迭代时间为 0.0156s;共轭斜量法迭代次数为 13 次,迭代时间为: 0.0019s,由于在计算迭代时间时 SOR 方法和块 SOR 方法没有算上求解最佳松弛因子的时间,所以用迭代次数作为评价三种方法好坏的标准,显然 ,共轭斜量 法优于块 SOR 迭代法,块 SOR 迭代方法优于 SOR 迭代方法,即 CGBSORSOR。

Copyright © 2018-2021 Wenke99.com All rights reserved

工信部备案号浙ICP备20026746号-2  

公安局备案号:浙公网安备33038302330469号

本站为C2C交文档易平台,即用户上传的文档直接卖给下载用户,本站只是网络服务中间平台,所有原创文档下载所得归上传人所有,若您发现上传作品侵犯了您的权利,请立刻联系网站客服并提供证据,平台将在3个工作日内予以改正。