matlab高斯法实验题答案

高斯-若当全主元消去法解线性方程组我们知道对于一个方程组而言，如果做相应的变换使系数矩阵变为单位阵的话那么就能得到它的解了，高斯-若当法就是这样的，但是为了减少误差，使用了全选主元法,下面是程序。保存为gauss.m在matlab中运行functionans=gauss(a,n,b)indexC=[];%用于纪录主元indexR=[];ipiv=[];m=1;forj=1:nipiv(j)=0;endfori=1:nbig=0;forj=1:nif(ipiv(j)~=1)fork=1:nifabs(a(j,k))>big%寻找最大的主元big=abs(a(j,k));irow=j;icol=k;endendendendipiv(icol)=ipiv(icol)+1;if(irow~=icol)%需要调换forl=1:ntemp=a(irow,l);a(irow,l)=a(icol,l);a(icol,l)=temp;endforl=1:mtemp=b(irow,l);b(irow,l)=b(icol,l);b(icol,l)=temp;endendindexR(i)=irow;indexC(i)=icol;if(a(icol,icol)==0)disp('error');endpivinv=1/a(icol,icol);forl=1:na(icol,l)=a(icol,l)*pivinv;endforl=1:mb(icol,l)=b(icol,l)*pivinv;endforll=1:nif(ll~=icol)dum=a(ll,icol);a(ll,icol)=0.0;forl=1:na(ll,l)=a(ll,l)-a(icol,l)*dum;endforl=1:mb(ll,l)=b(ll,l)-b(icol,l)*dum;endendendendforl=n:-1:1%已经进行过交换，现在交换回来if(indexR(l)~=indexC(l))fork=1:n104060104060temp=a(k,indexR(l));a(k,indexR(l))=a(k,indexC(l));a(k,indexC(l))=a(k,indexR(l));endendendans=b;面解这个简单的方程组>>a=[1000;0200;030;0004]面解这个简单的方程组>>a=[1000;0200;030;0004]>>>>b=[10;40;60;80]8080>>gauss(a,4,b)ans=10202020昨天才回答过这个问题..你可以再搜搜的Gauss消去法的分析。其包括两个过程：消去过程：把方程组系数矩阵A化为同解的上三角矩阵；回代过程：按相反的顺序，从xn至x1逐个求解上三角方程组。%高斯消去法的MATLAB程序functionx=gauss(a,b);%编写高斯消去法函数%a表示方程组的系数矩阵，b表示方程组的值%X表示最终的输出结果，即方程组的解n=length(b);%计算方程组的维数%下面的程序在不断的消去，直到变成a变成上三角矩阵未知fork=1:n-1fori=k+1:na(i,k)=a(i,k)/a(k,k);forj=k+1:na(i,j)=a(i,j)-a(i,k)*a(k,j);endb(i)=b(i)-a(i,k)*b(k);endend%表示高斯消去法的回带过程x=zeros(n,1);x(n)=b(n)/a(n,n);fork=n-1:-1:1s=b(k);forj=k+1:ns=s-a(k,j)*x(j);endx(k)=s/a(k,k);end实例验证：%调用编好的消去法函数A=[1,2,3;2,2,3;-1,-3,10];B=[0,3,2];gauss(A,B)ans=3.0000-1.55170.0345A=[1,2,3;2,2,3;-1,-3,10];B=[0,3,2];x=gauss(A,B)x=3.0000-1.55170.0345A*x%反代求解进行比较ans=0.00003.00002.0000matlab程序报错意思释疑（转自他处）1・Subscriptindicesmusteitherberealpositiveintegersorlogicals中文解释：下标索引必须是正整数类型或者逻辑类型出错原因：在访问矩阵（包括向量、二维矩阵、多维数组，下同）的过程中，下标索引要么从0开始，要么出现了负数。注：matlab的语法规定矩阵的索引从1开始，这与C等编程语言的习惯不一样。解决办法：自己调试一下程序，把下标为0或者负数的地方修正。2.Undefinedfunctionorvariable"a"中文解释：函数或变量a没有定义3・lnputargument"x"isundefined中文解释：输入变量x没有定义4・MatrixdimensionsmustagreeInnermatrixdimensionsmustagree中文解释：矩阵的维数必须一致出错原因：这是由于运算符(=+-/*等)两边的运算对象维数不匹配造成的，典型的出错原因是错用了矩阵运算符。matlab通过“来区分矩阵运算和元素运算Functiondefinitionsarenotpermittedatthepromptorinscripts中文解释：不能在命令窗口或者脚本文件中定义函数出错原因：一旦在命令窗口写functionc=myPlus(a,b),此错误就会出现，因为函数只能定义在m文件中1)Xmusthaveoneortwocolumns2・)Vectorsmustbethesamelengths中文解释：X必须是1或者2列向量长度必须一致7.Oneormoreoutputargumentsnotassignedduringcallto'・・・'中文解释：在调用…函数过程中，一个或多个输出变量没有被赋值8.???Errorusing==>mpowerMatrixmustbesquare中文解释：错误使用mpwoer函数，要求矩阵必须是方阵9・Explicitintegralcouldnotbefound・中文解释：显式解没有找到10」ndexexceedsmatrixdimensions.Attemptedtoaccessb(3,2);indexoutofboundsbecausesize(b)=[2,2]・中文解释：索引超出矩阵的范围附本人遇至附本人遇至U的问题：Attemptedtoaccessx(3);indexoutofboundsbecausenumel(x)=2.表示x矩阵元素不够，导致x(3)无法索引。11」nanassignment