您好,欢迎光临本网站![请登录][注册会员]  
文件名称: 用SOR迭代法解法方程①
  所属分类: 其它
  开发工具:
  文件大小: 32kb
  下载次数: 0
  上传时间: 2009-12-19
  提 供 者: gy111******
 详细说明: 用SOR迭代法解法方程① SOR.m文件为 function [x,n]=SOR(a,b,E,p) temp1=size(a);temp2=size(b); if temp1(1)~=temp1(2)|temp1(1)~=temp2(1)|temp2(2)~=1|det(a)==0 error('矩阵或向量的大小不对应或矩阵为奇异矩阵') return; end %%化为三角矩阵过程 leng=temp1(1);N=-tril(a,-1)-triu(a,1);M=a+N;N=-a; t=1; for i=1:leng t=t*M(i,i); end if t==0 error('对角元素都为零!') end x=zeros(leng,1); t=0;x1=zeros(leng,1); for i=1:leng for j=1:leng if j=i x(i)=N(i,j)*x1(j)+x(i); end end x(i)=x1(i)+(x(i)+b(i))/M(i,i)*p; end for i=1:leng t=t+(x1(i)-x(i))^2; end e=E^2; n=1; while t>e t=0; x1=x; x=zeros(leng,1); for i=1:leng temp=0; for j=1:leng if j=i x(i)=N(i,j)*x1(j)+x(i); end end x(i)=x1(i)+(x(i)+b(i))/M(i,i)*p; t=t+(x1(i)-x(i))^2; end n=n+1; if t>1000 break; error('不收敛') end end 控制命令窗口命令为: 当松弛因子p=1.3时 [x,n]=SOR(a,b,0.001,1.3) x = 3.0000 2.0000 1.0001 n = 11 当松弛因子p=1时也即是高斯-塞德尔方法 >> [x,n]=SOR(a,b,0.001,1) x = 2.9998 2.0001 1.0001 n = 5 当松弛因子p=0.9时 [x,n]=SOR(a,b,0.001,0.9) x = 3.0001 1.9999 0.9999 n = 6 所以当p=1时松弛因子比较好 通过比较发现 高斯-塞德尔迭代法收敛的速度明显比雅各比收敛速度快 ...展开收缩
(系统自动生成,下载前可以参看下载内容)

下载文件列表

相关说明

  • 本站资源为会员上传分享交流与学习,如有侵犯您的权益,请联系我们删除.
  • 本站是交换下载平台,提供交流渠道,下载内容来自于网络,除下载问题外,其它问题请自行百度
  • 本站已设置防盗链,请勿用迅雷、QQ旋风等多线程下载软件下载资源,下载后用WinRAR最新版进行解压.
  • 如果您发现内容无法下载,请稍后再次尝试;或者到消费记录里找到下载记录反馈给我们.
  • 下载后发现下载的内容跟说明不相乎,请到消费记录里找到下载记录反馈给我们,经确认后退回积分.
  • 如下载前有疑问,可以通过点击"提供者"的名字,查看对方的联系方式,联系对方咨询.
 相关搜索: sor 迭代
 输入关键字,在本站1000多万海量源码库中尽情搜索: