1、实验一、误差分析一、实验目的1通过上机编程，复习巩固以前所学程序设计语言及上机操作指令；2通过上机计算，了解误差、绝对误差、误差界、相对误差界的有关概念；3通过上机计算，了解舍入误差所引起的数值不稳定性。二实验原理误差问题是数值分析的基础，又是数值分析中一个困难的课题。在实际计算中，如果选用了不同的算法，由于舍入误差的影响，将会得到截然不同的结果。因此，选取算法时注重分析舍入误差的影响，在实际计算中是十分重要的。同时，由于在数值求解过程中用有限的过程代替无限的过程会产生截断误差，因此算法的好坏会影响到数值结果的精度。三实验内容对 ，计算定积分20,1n.105dxynn算法 1：利用递推公式,
2、 ,15nny2,1n取 .0 83.05l6dx算法 2：利用递推公式.nnyy5111,920注意到,101022102 556dxxd取 .:873.)16(2y四.实验程序及运行结果程序一： t=log(6)-log(5);n=1;y(1)=t;for k=2:1:20y(k)=1/k-5*y(k-1);n=n+1;endyy =0.0884y =0.0581y =0.0431y =0.0346y =0.0271y =0.0313y =-0.0134y =0.1920y =-0.8487y =4.3436y =-21.6268y =108.2176y =-541.0110y =2.70
3、51e+003y =-1.3526e+004y =6.7628e+004y =-3.3814e+005y =1.6907e+006y =-8.4535e+006y =4.2267e+007程序 2：y=zeros(20,1);n=1;y1=(1/105+1/126)/2;y(20)=y1;for k=20:-1:2y(k-1)=1/(5*k)-(1/5)*y(k);n=n+1;end运行结果：y =0.08840.05800.04310.03430.02850.02120.01880.01690.01540.01410.01300.01200.01120.01050.00990.00930.0
4、089实验二、插值法一、实验目的1、理解插值的基本概念，掌握各种插值方法，包括拉格朗日插值和牛顿插值等，注意其不同特点；2、通过实验进一步理解并掌握各种插值的基本算法。二实验原理插值法是函数逼近的一种重要方法，它是数值积分、微分方程数值解等数值计算的基础与工具，其中多项式插值是最常用和最基本的方法。拉格朗日插值多项式的优点是表达式简单明确，形式对称，便于记忆，它的缺点是如果想要增加插值节点，公式必须整个改变，这就增加了计算工作量。而牛顿插值多项式对此做了改进，当增加一个节点时只需在原牛顿插值多项式基础上增加一项，此时原有的项无需改变，从而达到节省计算次数、节约存储单元、应用较少节点达到应有精度
5、的目的。三实验内容1、 已知函数表0.56160 0.56280 0.56401 0.56521ix0.82741 0.82659 0.82577 0.82495iy用二次拉格朗日插值多项式求 时的函数近似值。563.0x2、 已知函数表0.4 0.55 0.65 0.8 0.9ix0.41075 0.57815 0.69675 0.88811 1.02652iy用牛顿插值多项式求 和 。)596.0(3N)895.0(4四实验程序及运行结果1.程序如下：function y,R=lagranzi(X,Y,x,M)n=length(X); m=length(x);for i=1:mz=x(i)
6、;s=0.0;for k=1:np=1.0; q1=1.0; c1=1.0;for j=1:nif j=kp=p*(z-X(j)/(X(k)-X(j);endq1=abs(q1*(z-X(j);c1=c1*j;ends=p*Y(k)+s;endy(i)=s;endR=M*q1/c1;在 MATLAB 工作窗口输入程序： x=0.5635; M=1; X=0.56280,0.56401; Y=0.82659,0.82577; y,R=lagranzi(X,Y,x,M)实验结果如下：y =0.8261R =1.7850e-0072.程序如下：function y,R= newcz(X,Y,x,M)
7、n=length(X); m=length(x);for t=1:mz=x(t); A=zeros(n,n);A(:,1)=Y; s=0.0; p=1.0; q1=1.0; c1=1.0;for j=2:nfor i=j:n A(i,j)=(A(i,j-1)-A(i-1,j-1)/(X(i)-X(i-j+1);endq1=abs(q1*(z-X(j-1);c1=c1*j;endC=A(n,n);q1=abs(q1*(z-X(n);for k=(n-1):-1:1C=conv(C,poly(X(k);d=length(C); C(d)=C(d)+A(k,k);endy(k)= polyval(C
8、, z);endR=M*q1/c1;在 MATLAB 工作窗口输入程序：三次牛顿差值程序如下： X=0.4,0.55,0.6; Y=0.41075,0.57815,0.69675; y,R= newcz(X,Y,x,M)实验结果如下：y =0.6861R =6.0107e-006四次牛顿差值程序如下： x=0.895; X=0.55,0.65,0.8,0.9; Y=0.57815,0.69675,0.88811,1.02652;实验结果如下：y =1.0194R =1.6729e-006实验三、解线性方程组的直接法一、实验目的1了解求线性方程组的直接法的有关理论和方法；2会编制列主元消去法、L
9、U 分解法的程序；二实验原理解线性方程组的直接法是指经过有限步运算后能求得方程组精确解的方法。但由于实际计算中舍入误差是客观存在的，因而使用这类方法也只能得到近似解。目前较实用的直接法是古老的高斯消去法的变形，即主元素消去法及矩阵的三角分解法。引进选主元的技巧是为了控制计算过程中舍入误差的增长，减少舍入误差的影响。一般说来，列主元消去法及列主元三角分解法是数值稳定的算法，它具有精确度较高、计算量不大和算法组织容易等优点，是目前计算机上解中、小型稠密矩阵方程组可靠而有效的常用方法。三实验内容1. 用列主元高斯消去法求解方程组.03564.5072.1835.217248.3 xx2用矩阵直接三角
10、分解法求解方程组 ，其中bA， .3812659734549127b四实验程序及运行结果1.主程序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=RBif RA=ndisp(请注意：因为 RA=RB=n，所以此方程组有唯一解.) X=zeros(n,1); C=zeros(1,n+1);for p= 1:n-1Y,j=max(abs(B(p:n,p); C=B(p,:)
11、;B(p,:)= B(j+p-1,:); B(j+p-1,:)=C;for k=p+1:nm= B(k,p)/ B(p,p);B(k,p:n+1)= B(k,p:n+1)-m* B(p,p:n+1);endendb=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);endelsedisp(请注意：因为 RA=RB A=0.101 2.304 3.555;-1.347 3.712 4.623;-2.835 1.072 5.643;b=1.183;2.
12、137;3.035;RA,RB,n,x=liezhu(A,b)运行结果 RA =3RB =3n =3x =-0.39820.01380.33512.主程序function X=LUjfcz(A,b)n,n =size(A); X=zeros(n,1); Y=zeros(n,1); C=zeros(1,n);r=1:n;for p=1:n-1max1,j=max(abs(A(p:n,p);C=A(p,:); A(p,:)= A(j+p-1,:); A(j+p-1,:)=C;g=r(p); r(p)= r(j+p-1);r(j+p-1)=g;if A(p,p)=0disp(A 是奇异阵,方程组无唯
13、一解);break;endfor k=p+1:nH= A(k,p)/A(p,p); A(k,p) = H; A(k,p+1:n)=A(k,p+1:n)- H* A(p,p+1:n);endendY(1)=b(r(1);for k=2:nY(k)= b(r(k)- A(k,1:k-1)* Y(1:k-1);endX(n)= Y(n)/ A(n,n);for i=n-1:-1:1X(i)= (Y(i)- A(i, i+1:n) * X (i+1:n)/ A(i,i);end在 MATLAB 工作窗口输入程序： A=1 2 -12 8;5 4 7 -2;-3 7 9 5;6 -12 -8 3;b=2
14、7;4;11;49;X=LUjfcz(A,b)运行结果X =3.0000-2.00001.00005.0000实验四、解线性方程组的迭代法一、实验目的 1、熟悉迭代法的有关理论和方法；2、会编制雅可比迭代法、高斯-塞德尔迭代法的程序；3、注意所用方法的收敛性及其收敛速度问题。二实验原理解线性方程组的迭代法是用某种极限过程去逐步逼近线性方程组精确解的方法，即是从一个初始向量 出发，按照一定的迭代格式产生一个向量序列 ，使其收敛到方程组 的解。迭代法的优点是所需计算机存储单元少，程序设计简单，原始系数矩阵在计算过程中始终不变等。但迭代法存在收敛性及收敛速度问题。迭代法是解大型稀疏矩阵方程组的重要方
15、法。三.实验内容1.用改进欧拉法解 .2)0(5yx2用龙格-库塔方法求解 )(0y四.实验程序及运行结果1.主程序function X=jacdd(A,b,X0,P,wucha,max1)n m=size(A);for j=1:ma(j)=sum(abs(A(:,j)-2*(abs(A(j,j);endfor k=1:max1kfor j=1:mX(j)=(b(j)-A(j,1:j-1,j+1:m)*X0(1:j-1,j+1:m)/A(j,j);endX,djwcX=norm(X-X0,P); xdwcX=djwcX/(norm(X,P)+eps); X0=X;X1=Ab;if (djwcX
16、wucha)jX=X1,在 MATLAB 工作窗口输入程序： A=1 2 -2;1 1 1;2 2 1;b=7;2;5;X=jacdd(A,b,X0,P,wucha,max1)运行结果X =1 2 -12.主程序function X=gsdddy(A,b,X0,P,wucha,max1)D=diag(diag(A);U=-triu(A,1);L=-tril(A,-1); dD=det(D);if dD=0disp(请注意：因为对角矩阵 D 奇异，所以此方程组无解.)elsedisp(请注意：因为对角矩阵 D 非奇异，所以此方程组有解.)iD=inv(D-L); B2=iD*U;f2=iD*b;
17、jX=Ab; X=X0; n m=size(A);for k=1:max1X1= B2*X+f2; djwcX=norm(X1-X,P);xdwcX=djwcX/(norm(X,P)+eps);if (djwcX A=1 0.9 0.9;0.9 1 0.9;0.9 0.9 1;b=1.9;2.0;1.7;X0=0 0 0;X=gsdddy(A,b,X0,inf, 0.001,100)运行结果k =1ans =1.9000 0.2900 -0.2710k =2ans =1.8829 0.5493 -0.4890k =3ans =1.8457 0.7789 -0.6622k =4ans =1.79
18、49 0.9805 -0.7979k =5ans =1.7356 1.1560 -0.9025k =6ans =1.6718 1.3076 -0.9815k =7ans =1.6065 1.4375 -1.0396k =8ans =1.5419 1.5479 -1.0808k =9ans =1.4796 1.6411 -1.1086k =10ans =1.4208 1.7191 -1.1259k =11ans =1.3661 1.7838 -1.1349k =12ans =1.3160 1.8370 -1.1377k =13ans =1.2706 1.8804 -1.1359k =14ans
19、=1.2300 1.9153 -1.1308k =15ans =1.1939 1.9432 -1.1234k =16ans =1.1622 1.9651 -1.1145k =17ans =1.1345 1.9820 -1.1049k =18ans =1.1106 1.9949 -1.0949k =19ans =1.0900 2.0044 -1.0850k =20ans =1.0725 2.0112 -1.0754k =21ans =1.0577 2.0159 -1.0662k =22ans =1.0453 2.0188 -1.0577k =23ans =1.0350 2.0204 -1.049
20、9k =24ans =1.0265 2.0210 -1.0428k =25ans =1.0196 2.0209 -1.0364k =26ans =1.0140 2.0202 -1.0308k =27ans =1.0095 2.0191 -1.0258k =28ans =1.0060 2.0178 -1.0214k =29ans =1.0032 2.0164 -1.0176k =30ans =1.0012 2.0148 -1.0144X =0.99852.0118-1.0093实验心得通过本次实验，我学会了将在课堂上学到的一些基本的数值分析方法运用到实践之中进行数值分析，感受到了利用计算机这个工具解决实际计算的优越性，由此也明白了为了更简单、更快速地得到一个更准确的结果，选择一个好的算法的重要性。因此，我们在今后的学习工作中应该更加重视算法的设计。另外，由于自己之前没有学习 Matlab 语言，所以在实验时遇到了很多问题，这次实验同时也明白了要学好 Matlab 语言的重要性。今后一定要利用自己的课余时间补上这一课，为将来工作中的运用打下坚实的基础。
