第二次电子作业 信息计算.docVIP

解：将点（1,6），（3,5），（7,2）分别代入抛物线方程y=a+bx+cx2中， 得， a + b + c =6 a + 3b + 9c =5 a + 7b + 49c =2 M文件： function [RA,RB,n,X]=liezhu(A,b)B=[A b]; n=length(b); RA=rank(A); RB=rank(B);zhica=RB-RA;if zhica0,disp(请注意：因为RA~=RB，所以此方程组无解.)returnendif RA==RB if RA==ndisp(请注意：因为RA=RB=n，所以此方程组有唯一解.)  X=zeros(n,1); C=zeros(1,n+1); for p= 1:n-1[Y,j]=max(abs(B(p:n,p))); C=B(p,:);B(p,:)= B(j+p-1,:); B(j+p-1,:)=C;for k=p+1:n m= B(k,p)/ B(p,p); B(k,p:n+1)= B(k,p:n+1)-m* B(p,p:n+1);endend b=B(1:n,n+1);A=B(1:n,1:n); X(n)=b(n)/A(n,n); for q=n-1:-1:1 X(q)=(b(q)-sum(A(q,q+1:n)*X(q+1:n)))/A(q,q); endelse disp(请注意：因为RA=RBn，所以此方程组有无穷多解.)endend MATLAB工作窗口输入程序： A=[1 1 1;1 3 9;1 7 49]; b=[6;5;2]; [RA,RB,n,X]=liezhu(A,b) 程序运行截图： 所以:抛物线方程为y=6.3750+（-0.3333）x+（-0.0417）x 2。 (模仿任玉杰书162页 例3.3.4) 解： M文件： function hl=zhjLU(A) [n n] =size(A); RA=rank(A); if RA~=n disp(请注意：因为A的n阶行列式hl等于零，所以A不能进行LU分解.A的秩RA如下:), RA,hl=det(A); return end if RA==n for p=1:n h(p)=det(A(1:p, 1:p)); end hl=h(1:n); for i=1:n if h(1,i)==0 disp(请注意：因为A的r阶主子式等于零，所以A不能进行LU分解.A的秩RA和各阶顺序主子式值hl依次如下：), hl;RA return end end if h(1,i)~=0 disp(请注意：因为A的各阶主子式都不等于零，所以A能进行LU分解.A的秩RA和各阶顺序主子式值hl依次如下：) for j=1:n U(1,j)=A(1,j); end for k=2:n for i=2:n for j=2:n L(1,1)=1;L(i,i)=1; if ij L(1,1)=1;L(2,1)=A(2,1)/U(1,1); L(i,1)=A(i,1)/U(1,1); L(i,k)=(A(i,k)- L(i,1:k-1)*U(1:k-1,k))/U(k,k); else U(k,j)=A(k,j)-L(k,1:k-1)*U(1:k-1,j); end end end end hl;RA,U,L end end MATLAB工作窗口输入程序： (1) A =[2 4 -6;1 5 3;1 3 2]; hl = zhjLU(A) (2) A =[1 1 6;-1 2 9;1 -2 3]; hl = zhjLU(A) 程序运行截图： (1) 请注意：因为A的各阶主子式都不等于零，所以A能进行LU分解，A的秩RA和各阶顺序主子式值h1依次如下： (2) 请注意：因为A的各阶主子式都不等于零，所以A能进行LU分解，A的秩RA和各阶顺序主子式值h1依次如下： (模仿任玉杰书168页 例3.4.3) MATLAB工作窗口输入程序： A=[10 7 8 7;7 5 6 5;8 6 10 9;7 5 9 10];b=[32;23;33;31];X=A\b A=[10 7 8 7;7 5 6 5;8 6 10 9;7 5 9 10];b=[31.9;23.1;32.9;31.1];X=A\b 程序运行截图： b和解X的相对误差限分别为 ≈0.0034 . ≤Cond(