开发工具:
文件大小: 1kb
下载次数: 0
上传时间: 2010-11-15
详细说明: 高斯迭代法解矩阵方程 function [a,b] = gaussj(a,n,np,b,m,mp) for j=1:n ipiv(j)=0; end for i=1:n big=0; for j=1:n if(ipiv(j)~=1) for k=1:n if (ipiv(k)==0) if (abs(a(j,k))>=big) big=abs(a(j,k)); irow=j; icol=k; end elseif (ipiv(k)>1) pause, dips('gaussj') end end end end ipiv(icol)=ipiv(icol)+1; if (irow~=icol) for l=1:n dum=a(irow,l); a(irow,l)=a(icol,l); a(icol,l)=dum; end for l=1:m dum=b(irow,l); b(irow,l)=b(icol,l); b(icol,l)=dum; end end indxr(i)=irow; indxc(i)=icol; if (a(icol,icol)==0), pause, disp('in gaussj'), end pivinv=1/a(icol,icol); a(icol,icol)=1; for l=1:n a(icol,l)=a(icol,l)*pivinv; end for l=1:m b(icol,l)=b(icol,l)*pivinv; end for ll=1:n if(ll~=icol) dum=a(ll,icol); a(ll,icol)=0; for l=1:n a(ll,l)=a(ll,l)-a(icol,l)*dum; end for l=1:m b(ll,l)=b(ll,l)-b(icol,l)*dum; end end end end for l=n:-1:1 if(indxr(l)~=indxc(l)) for k=1:n dum=a(k,indxr(l)); a(k,indxr(l))=a(k,indxc(l)); a(k,indxc(l))=dum; end end end return ...展开收缩
(系统自动生成,下载前可以参看下载内容)
下载文件列表
相关说明
- 本站资源为会员上传分享交流与学习,如有侵犯您的权益,请联系我们删除.
- 本站是交换下载平台,提供交流渠道,下载内容来自于网络,除下载问题外,其它问题请自行百度。
- 本站已设置防盗链,请勿用迅雷、QQ旋风等多线程下载软件下载资源,下载后用WinRAR最新版进行解压.
- 如果您发现内容无法下载,请稍后再次尝试;或者到消费记录里找到下载记录反馈给我们.
- 下载后发现下载的内容跟说明不相乎,请到消费记录里找到下载记录反馈给我们,经确认后退回积分.
- 如下载前有疑问,可以通过点击"提供者"的名字,查看对方的联系方式,联系对方咨询.