数值计算试卷.docVIP

一 编写 Gauss 列主元消元法求解线性方程组的程序，要求附有算例。 程序： % 输出的量：系数矩阵和增广矩阵的秩RA,RB, 方程组中未知量的个数n和有关方程组解及其解X的信息. function [RA,RB,n,X]=gaus(A,b) B=[A b]; n=length(b); RA=rank(A); RB=rank(B); zhica=RB-RA; %求矩阵A、B的秩，并作比较 if zhica0, disp(RA~=RB，此方程组无解.) return end if RA==RB if RA==n disp(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))); %返回行向量Y和j，y向量记录abs(B(p:n,p)) 的每列的最大值，j向量记录每列最大值的行号 C=B(p,:); %把矩阵B中第p行所有列的值全都赋给矩阵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); end end 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 %从n-1循环到1 X(q)=(b(q)-sum(A(q,q+1:n)*X(q+1:n)))/A(q,q); %回代计算 end else disp(RA=RBn，此方程组有无穷多解.) end End 算例：求解线性方程组 在命令窗口输入： A=[5 2 1;2 8 -3;1 -3 -6]; b=[8;21;1]; [RA,RB,n,X] =gaus (A,b) RA=RB=n，此方程组有唯一解. RA = 3 RB = 3 n = 3 X = 1 2 -1 二 编写多项式插值的程序，要求附有算例 Lagrange差值法 程序 function yy=lagrange(x,y,x0) //定义拉格朗日差值函数 xlen=length(x); //X向量的维数 ylen=length(y); //y向量的维数 if xlen~=ylen warning(The length of x and y is not equal); //提示x.y维数不等 end n=length(x0); for i=1:n //循环。从1到n z=x0(i); s=0.0; for k=1:xlen p=1.0; for j=1:ylen if j~=k p=p*(z-x(j))/(x(k)-x(j)); end end s=p*y(k)+s; end yy(i)=s; end 算例：已知数据如下表，试用lagrange差值多项式求x分别为0.5626;0.5635;0.5645时函数的近似值。 0.56160 0.56280 0.56401 0.56521 0.82741 0.82659 0.82577 0.81495 解：在命令窗口输入执行命令，可得结果： x=[0.56160;0.56280;0.56401;0.56521]; y=[0.82741;0.82659;0.82577;0.81495]; x0=[0.5626;0.5635;0.5645]; y0=lagrange(x,y,x0) y0 = 0.8265 0.8268 0.8231 plot(x,y,o,x0,y0,k^) 得到图形 三 编写 Gauss 积分方法的程