矩阵与数值分析_大连理工大学2011试卷

第三章 逐次逼近法1.1内容提要1、一元迭代法x n+1=φ(x n )收敛条件为：1）映内性x ∈[a,b]，φ(x) ∈[a,b] 2）压缩性∣φ(x) -φ(y)∣≤L ∣x-y ∣其中L ＜1，此时φ为压缩算子，在不断的迭代中，就可以得到最终的不动点集。

由微分中值定理，如果∣φ’∣≤L ＜1，显然它一定满足压缩性条件。

2、多元迭代法x n+1=φ(x n )收敛条件为：1）映内性x n ∈Ω，φ(x n ) ∈Ω 2）压缩性ρ（▽φ）＜1，其中▽φ为x n 处的梯度矩阵，此时φ为压缩算子，在不断的迭代中，就可以得到最终的不动点集。

3、当φ(x )= Bx+f 时，收敛条件为，ρ（B ）＜1，此时x n+1= Bx n +f ，在不断的迭代中，就可以得到线性方程组的解。

4、线性方程组的迭代解法，先作矩阵变换 U L D A --= Jacobi 迭代公式的矩阵形式 f Bx b D x U L D x n n n +=++=--+111)(Gauss-Seidel 迭代公式的矩阵形式 f Bx b L D Ux L D x n n n +=-+-=--+111)()( 超松弛迭代法公式的矩阵形式f Bx b L D x U D L D x k k k +=-++--=--+ωωωωω111)(])1[()(三种迭代方法当1)(<B ρ时都收敛。

5、线性方程组的迭代解法，如果A 严格对角占优，则Jacob 法和Gauss-Seidel 法都收敛。

6、线性方程组的迭代解法，如果A 不可约对角占优，则Gauss-Seidel 法收敛。

7、Newton 迭代法，单根为二阶收敛 2211'''21lim)(2)(lim---∞→+∞→--=-==--k k k k k k k k x x x x f f c x x ξξαα8、Newton 法迭代时，遇到重根，迭代变成线性收敛，如果知道重数m ， )()('1k k k k x f x f m x x -=+仍为二阶收敛 9、弦割法)()())((111--+---=k k k k k k k x f x f x x x f x x 的收敛阶为1.618，分半法的收敛速度为（b-a ）/2n-110、Aitken 加速公式11211112)(),(),(+----+-+--+---+---===k k k k k k k k k k k x x x x x x x x x x x ϕϕ1.2 典型例题分析1、证明如果A 严格对角占优，则Jacob 法和Gauss-Seidel 法都收敛。

-03cos 2lnlim 0=+=®xx （10分）四、解：(1)0)cos )((lim 00sin )(lim 00=-¢=÷øöçèæ-=®®x x g x x x g a x x （4分）(2)200sin )(lim )0()(lim )0(x xx g x f x f f x x-=-=¢®® =12)0(2sin )(lim 2cos )(lim 00=¢¢=+¢¢=-¢®®g x x g x x x g x x∴ ïîïíì=¹---¢=¢时时010,)sin )(()cos )(()(2x x x x x g x x g x x f （8分） （3）200)sin )(()cos )((lim )(lim x x x g x x g x x f x x ---¢=¢®® =xx x g x x g x x x g x 2)cos )(()sin )((cos )(lim 0-¢-+¢¢+-¢® =)0(12)0(f g ¢==¢¢，因此)(x f ¢在(-∞，+∞+∞))连续。

连续。

（10分）五、解五、解:: 设x x x f ln)(=，由2ln 1)('xxx f -=，可知,当e x >时)(x f 单调减少单调减少 （5分）若e a b >>，则有b b a a ln ln >，推出a b b a ln ln >，即有a b b a > 2011201220122011> （10分）分）所以六、解：2)()()(x x f x f x x x f -¢=¢÷øöçèæ（4分）分） 令)()()(x f x f x x g -¢=，)()(x f x x g ¢¢=¢，令0)(=¢x g ，得0=x （唯一驻点），当0<x 时，0)(<¢x g ，当0>x 时，0)(>¢x g ，故)0(g 为最小值，故0)0()0()(>-=³f g x g ，∴0)(>¢÷øöçèæx x f ，即x x f )(单调增加。

矩阵与数值分析学生：学号：任课老师：金光日教学班号：（2）班院系：电子信息与电气工程学部《矩阵与数值分析》课程数值实验题目1.给定n 阶方程组A x b =，其中6186186186A ⎛⎫ ⎪ ⎪⎪= ⎪ ⎪ ⎪⎝⎭，7151514b ⎛⎫ ⎪⎪ ⎪= ⎪ ⎪⎪⎝⎭则方程组有解(1,1,,1)T x = 。

对10n =和84n =，分别用Gauss 消去法和列主元消去法解方程组，并比较计算结果。

1答： 程序1. Gauss 消元法function x=DelGauss(A,b) % Gauss 消去法 [n,m]=size(A); det=1; %存储行列式值 x=zeros(n,1); for k=1:n-1 for i=k+1:n if A(k,k)==0 return endm=A(i,k)/A(k,k); for j=k+1:nA(i,j)=A(i,j)-m*A(k,j); endb(i)=b(i)-m*b(k); enddet=det*A(k,k); %计算行列式enddet=det*A(n,n);for k=n:-1:1 %回代求解for j=k+1:nb(k)=b(k)-A(k,j)*x(j);endx(k)=b(k)/A(k,k);end2. 列主元Gauss消去法：function x=detGauss(A,b)% Gauss列主元消去法[n,m]=size(A);nb=length(b);det=1; %存储行列式值x=zeros(n,1);for k=1:n-1amax=0; %选主元for i=k:nif abs(A(i,k))>amaxamax=abs(A(i,k));r=i;endendif amax<1e-10return;endif r>k %交换两行for j=k:nz=A(k,j);A(k,j)=A(r,j);A(r,j)=z;endz=b(k);b(k)=b(r);b(r)=z;det=-det;endfor i=k+1:n %进行消元m=A(i,k)/A(k,k);for j=k+1:nA(i,j)=A(i,j)-m*A(k,j);endb(i)=b(i)-m*b(k);enddet=det*A(k,k);enddet=det*A(n,n);for k=n:-1:1 %回代求解for j=k+1:nb(k)=b(k)-A(k,j)*x(j);endx(k)=b(k)/A(k,k);end矩阵A和b的构造clc;clear;n=10;%n=84;A=eye(n)*6+diag(ones(1,n-1)*8,-1)+diag(ones(1,n-1),1); b=[7,15*ones(1,n-2),14]';计算结果：(1)n=10时Gauss消元法>>x=DelGauss(A,b)x =1.00001.00001.00001.00001.00001.00001.00001.00001.00001.0000列主元Gauss消去法>>x=detGauss(A,b)x =1111111111(2) n=84时Gauss消元法>>x=DelGauss(A,b) x =1.0e+008 *0.00000.00000.00000.00000.00000.00000.00000.00000.00000.00000.00000.00000.00000.00000.00000.00000.00000.00000.00000.00000.00000.00000.00000.00000.00000.00000.00000.00000.00000.00000.00000.00000.00000.00000.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 -0.0000 0.0000 -0.0000 0.0000 -0.0000 0.0000 -0.0000 0.0000 -0.0000 0.0000 -0.0000 0.0000 -0.0001 0.0002 -0.0003 0.0007 -0.0013 0.0026 -0.0052 0.0105 -0.0209 0.0419 -0.0836 0.16650.6501-1.25822.3487-4.02635.3684列主元Gauss消去法>>x=detGauss(A,b) x =1.00001.00001.00001.00001.00001.00001.00001.00001.00001.00001.00001.00001.00001.00001.00001.00001.00001.00001.00001.00001.00001.00001.00001.00001.00001.00001.00001.00001.00001.00001.00001.00001.00001.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.00001.0000 1.0000 1.00001.00001.0000 1.0000结果分析由上述实验结果可知，对于n=10采用Gauss 消去法和Gauss 列主元消去法得到的实验结果是相同的，而对于n=84，Gauss 消去法所得到的结果是错误的，Gauss 列主元消去法得到的结果是正确的。

大连理工大学矩阵与数值分析上机作业课程名称：矩阵与数值分析研究生姓名：交作业日时间：2016 年12 月20日1.1程序：Clear all;n=input('请输入向量的长度n：')for i=1:n;v(i)=1/i;endY1=norm(v,1)Y2=norm(v,2)Y3=norm(v,inf)1.2结果n=10 Y1 =2.9290Y2 =1.2449Y3 =1n=100 Y1 =5.1874Y2 =1.2787Y3 =1n=1000 Y1 =7.4855Y2 =1.2822Y3 =1N=10000 Y1 =9.7876Y2 =1.2825Y3 =11.3 分析一范数逐渐递增，随着n的增加，范数的增加速度减小；二范数随着n的增加，逐渐趋于定值，无群范数都是1.2.1程序clear all;x(1)=-10^-15;dx=10^-18;L=2*10^3;for i=1:Ly1(i)=log(1+x(i))/x(i); d=1+x(i);if d == 1y2(i)=1;elsey2(i)=log(d)/(d-1);endx(i+1)=x(i)+dx;endx=x(1:length(x)-1);plot(x,y1,'r');hold onplot(x,y2);2.2 结果2.3 分析红色的曲线代表未考虑题中算法时的情况，如果考虑题中的算法则数值大小始终为1，这主要是由于大数加小数的原因。

第3题3.1 程序clear all;A=[1 -18 144 -672 2016 -4032 5376 -4608 2304 -512];x=1.95:0.005:2.05;for i=1:length(x);y1(i)=f(A,x(i));y2(i)=(x(i)-2)^9;endfigure(3);plot(x,y1);hold on;plot(x,y2,'r');F.m文件function y=f(A,x) y=A(1);for i=2:length(A); y=x*y+A(i); end;3.2 结果第4题4.1 程序clear all;n=input('请输入向量的长度n：')A=2*eye(n)-tril(ones(n,n),0);for i=1:nA(i,n)=1;endn=length(A);U=A;e=eye(n);for i=1:n-1[max_data,max_index]=max(abs(U(i:n,i))); e0=eye(n);max_index=max_index+i-1;U=e0*U;e1=eye(n);for j=i+1:ne1(j,i)=-U(j,i)/U(i,i);endU=e1*U;P{i}=e0;%把变换矩阵存到P中L{i}=e1;e=e1*e0*e;endfor k=1:n-2Ldot{k}=L{k};for i=k+1:n-1Ldot{k}=P{i}*Ldot{k}*P{i};endendLdot{n-1}=L{n-1};LL=eye(n);PP=eye(n);for i=1:n-1PP=P{i}*PP;LL=Ldot{i}*LL;endb=ones(n,2);b=e*b; %解方程x=zeros(n,1);x(n)=b(n)/U(n,n);for i=n-1:-1:1x(i)=(b(i)-U(i,:)*x)/U(i,i);endX=U^-1*e^-1*eye(n);%计算逆矩阵AN=X';result2{n-4,1}=AN;result1{n-4,1}=x;fprintf('%d:\n',n)fprintf('%d ',AN);4.2 结果n=51.0625 -0.875 -0.75 -0.5 -0.06250.0625 1.125 -0.75 -0.5 -0.06250.0625 0.125 1.25 -0.5 -0.06250.0625 0.125 0.25 1.5 -0.0625-0.0625 -0.125 -0.25 -0.5 0.0625n=101.0625 -0.875 -0.75 -0.5 -0.0625 1.0625 -0.875 -0.75 -0.5 -0.0625 0.0625 1.125 -0.75 -0.5 -0.0625 0.0625 1.125 -0.75 -0.5 -0.0625 0.0625 0.125 1.25 -0.5 -0.0625 0.0625 0.125 1.25 -0.5 -0.0625 0.0625 0.125 0.25 1.5 -0.0625 0.0625 0.125 0.25 1.5 -0.0625 -0.0625 -0.125 -0.25 -0.5 0.0625 -0.0625 -0.125 -0.25 -0.5 0.0625 1.0625 -0.875 -0.75 -0.5 -0.0625 1.0625 -0.875 -0.75 -0.5 -0.0625 0.0625 1.125 -0.75 -0.5 -0.0625 0.0625 1.125 -0.75 -0.5 -0.0625 0.0625 0.125 1.25 -0.5 -0.0625 0.0625 0.125 1.25 -0.5 -0.0625 0.0625 0.125 0.25 1.5 -0.0625 0.0625 0.125 0.25 1.5 -0.0625 -0.0625 -0.125 -0.25 -0.5 0.0625 -0.0625 -0.125 -0.25 -0.5 0.0625同样的方法可以算出n=20,n=30时的结果，这里就不罗列了。

大连理工数学分析试题及解答大连理工大学硕士生入学考试数学分析试题一. 从以下的1到8题中选答6题1. 证明:2()f x x =在区间[0,]M 内一致连续(M 为任意正数),但是在[0,)+∞不一致连续2. 证明:若()f x 在[,]a b 内连续,那么()f x 在[,]a b 内Riemann 可积.3. 证明:若1α>,那么广义积分1sin x dx α+∞收敛4. 证明:若()f x ,()g x 为区间(,)a b 上的连续函数,对任意的(,)(,)a b αβ?有:()()f x dx g x dx ββαα=??,那么, ()()f x g x ≡于(,)a b5. 证明:若1nn a∞=∑收敛,那么1nxn n a e∞-=∑在[0,)+∞一致收敛6. 已知：2,0()0,0x e x f x x -?≠?=?=??，求"(0)f7. 已知：()()1(,)()22x atx at x at x at u x t d aφφψαα+-++-=+. 其中, ψ和φ分别是可以求导一次和求导两次的已知函数,计算22222(,)(,)u x t u x t a t x ??-??8. 计算,半径为R 的球的表面积二. 从9到14题中选取6题9.已知: lim '()0x f x →∞=,求证: ()lim0x f x x→∞=10.证明: ()af x dx +∞收敛,且lim ()x f x λ→+∞=,那么0λ=11.计算曲面积分: 333SI x dydz y dzdx z dxdy =++??, 其中S 为旋转椭球面2222221x y z a b c++=的外侧12.设()[0,1]f x C ∈,(0)0f =,(1)1f =,0()1f x ≤<. 求证: ()()n n S x f x =对于任意小于1的正数δ,在区间(0,1]δ-一致收敛,但是不在(0,1)一致收敛13.设()[0,1]f x C ∈,(0)0f =,(1)1f =,0()1f x ≤<. 求证: 1lim ()0n n f x dx →∞=?14.证明:若()[,]n u x C a b ∈，1,2,...,...n =且1()n n u b ∞=∑发散,那么1()n n u x ∞=∑不在[,)a b 一致收敛大连理工大学2001年硕士生入学考试数学分析试题解答一.1. 证利用定义证明(1) 对于0ε?>,21M εδ?=+,12||x x δ?-<,那么12121212|()()||()()|2||2f x f x x x x x M x x M δε-=-+<-<<(2) 任取1ε=,0δ?>,1211,22x x δδδ==+, 1212121|()()||()()|1f x f x x x x x δδ-=-+>=,推出矛盾,从而命题得证■2. 证利用一致连续的定义和Riemann 可积的定义来做因为函数在闭区间内连续,所以一致连续. 根据一致连续的定义对0ε?>,δ?,12||x x δ?-<,12|()()|f x f x ε-<考虑可积的定义,对于一个[,]a b 分割112:...n a a a a b ?=<<<=,11max ||i i i na a λ+≤<=-下面证明:振幅函数121110,[,]1()limmax {()}()i i n i i x x a a i w x f x a a λ+-+→∈==-∑=0当λδ<时,12111110,[,]110()limmax {()}()()i i n n i i i i x x a a i i w x f x a a a a b λεε+--++→∈==≤=-≤-=∑∑.根据夹逼定理,不难得到()0w x =. 从而,命题得证■3. 证利用莱布尼兹交错级数:假设;n a n π=,1sin nn a n a s x dx α-=?考虑:111|||||sin ||sin |n nnn a a n n a a s s x dx x dx αα+-+-=-?1111[|sin ||sin |]n n n n x x dx xx dx ππππααππα--+-=+??1111[|sin |(2)|sin |]n n n n xx dx n x x dx ππππααππαπ--++=--??1111[(2)]|sin |0n n x n x x dx ππααπαπ--+=--<?11lim |||sin |||lim ||0nnn n a a n n a a n n s x dx dx n n s αααπππ--→∞→∞=≤=--?=??如此,不难看出1sin x dx α+∞是一个莱布尼兹交错级数,从而命题得证■4. 证不妨设:2a bc +=()()x c F x f t dt =?,那么()()F x G x =于(,)x a b ∈因为()f x ()g x 都是(,)x a b ∈上的连续函数,所以()'()()()f x F x G xg x ===■5. 证利用A-D 判别法做，也可以通过Abel 求和公式出发推导1nxn n a e∞-=∑中nxn b e-=,现在,根据原题:1n n a∞=∑收敛,1nxnb e -=≤一致有界所以,根据Abel 判别法,知该函数项级数在定义域一致收敛. ■6. 解题目有问题,在零点不连续■7. 解不断利用链式求导法则()()1(,)()22x atx at x at x at u x t d aφφψαα+-++-=+(,)()()()()()()()()()()22'()'()'()'()22u x t xx at x at x at x at x at x at x at x at x at x x at x x x ax at x at x at x at aφφψψφφψψ+?+?-?-?+?-++--?+??-=+++-+--=+22'()'()()()(,)()()()()22"()"()'()'()22x at x at x at x at u x t x at x at x at x at x ax at x at x at x at aφφψψφφψψ?+?-?+?-+-+?-?+?-=+++-+--=+同理:(,)()()()()()()()()()()22'()'()'()'()22u x t tx at x at x at x at x at x at x at x at x at t x at t t t aa x at a x at x at x at φφψψφφψψ+?+?-?-?+?-++--?+??-=++--++-=+222'()'()()()(,)()()()()"()"()'()'()22x at x at x at x at aa u x t x at x at x at x at x x at x at x at x at a aφφψψφφψψ?+?-?+?--++?-?+?-=+++-+--=+22222(,)(,)0u x t u x t a t x ??-=??■8. 解方法很多,此处介绍一种比较简单的假设:()V R 为半径R 为的球的体积2234()()3R R V R R x dx R ππ-=-=?假设: ()S R 为半径R 为的球的表面积20()()()'()4RV R S x dx S R V R R π=?==?■二9. 证L ’Hosptial 法则因为x →+∞，()'()lim lim lim '()0'x x x f x f x f x x x →∞→∞→∞===■10. 证反证法如果命题不成立,即0λ≠,那么,根据极限的定义,G ?,当x G >的时候, |()|||2f x λ>()Gf x dx +∞→∞?和收敛矛盾,从而命题得证■11. 解利用Gauss 定理加换元3332223()VSI x dydz y dzdx z dxdy x y z dxdydz =++=++换元sin cos sin sin ,[0,1],[0,2),[0,]cos x ar y br r z cr ?θ?θθπ?π?=??=∈∈∈??=?4222222223sin (sin sin sin cos cos )VI abc r a b c drd d ??θ?θ?θ?=++22322322200033(sin sin sin cos )cos sin 55abc abc a b d d c d πππ?θ?θ?θ=++332232203646()sin ()5555abc abc abc abc a b d a b πππ??=++=++?■12. 证首先由于在闭区间内连续,所以函数在闭区间内一致连续(1)(0,1]x δ?∈-,根据确界存在定理,存在上确界,且上确界不等于1,否则和题意矛盾不妨设:(0,1]sup ()1x f x m δ∈-=<根据定义,对于0ε?>,ln ln N mε=,当n N >,|()||()|n n n S x f x m ε=≤< 从而知一致收敛于0(2)首先,根据前半题,显然()n S x 于(0,1)x ?∈收敛于0由于(1)1f =,且函数一致收敛，存在一组数列：12...a a <<,1()1i fa n=- 如此,考虑11lim ()lim ()lim(1)0nnn n n n n n S a f a ne→∞→∞→∞==-=≠,从而不是一致收敛的. ■13. 证利用前一小题的结论因为()nf x 内闭一致收敛,对于0ε?>,2εδ?=,当n 足够大的时候:10()2n f x dx δε-<又1111|()|||2n f x dx dx δδε--<=所以,1111()()()n n n f x dx f x dx f x dx δδε--=+<?从而命题得证. ■14. 证反证法:假设命题不成立,那么1()n n u x ∞=∑在[,)a b 一致收敛.即0ε?>,N ?,,m n N ?>,(,)x a b ?∈,|()|m n nu x ε<∑因为|()|lim |()|m mn n x bnnu b u b ε→=≤∑∑,否则与()[,]n u x C a b ∈矛盾而1|()|n n u b ∞=∑发散,所以|()|n n Nu b ∞=∑发散,与|()|lim |()|m m n n x bnnu b u b ε→=≤∑∑矛盾从而命题得证. ■。

共享知识分享快乐大连理工大学矩阵与数值分析上机作业课程名称：矩阵与数值分析研究生姓名：12 交作业日时间：日20 月年2016卑微如蝼蚁、坚强似大象．共享知识分享快乐第1题1.1程序：Clear ;all n=input('请输入向量的长度n：') for i=1:n;v(i)=1/i;endY1=norm(v,1)Y2=norm(v,2)Y3=norm(v,inf)1.2结果n=10 Y1 =2.9290Y2 =1.2449Y3 =1n=100 Y1 =5.1874Y2 =1.2787Y3 =1n=1000 Y1 =7.4855Y2 =1.2822Y3 =1N=10000 Y1 =9.7876Y2 =1.2825Y3 =11.3 分析一范数逐渐递增，随着n的增加，范数的增加速度减小；二范数随着n的增加，逐渐趋于定值，无群范数都是1.卑微如蝼蚁、坚强似大象．共享知识分享快乐第2题2.1程序;clear all x(1)=-10^-15;dx=10^-18;L=2*10^3; i=1:L fory1(i)=log(1+x(i))/x(i); d=1+x(i); d == 1ify2(i)=1;elsey2(i)=log(d)/(d-1);endx(i+1)=x(i)+dx;end x=x(1:length(x)-1););'r'plot(x,y1,on holdplot(x,y2);卑微如蝼蚁、坚强似大象．共享知识分享快乐2.2 结果2.3 分析红色的曲线代表未考虑题中算法时的情况，如果考虑题中的算法则数值大小始终为1，这主要是由于大数加小数的原因。

第3题3.1 程序;clear all A=[1 -18 144 -672 2016 -4032 5376 -4608 2304 -512];x=1.95:0.005:2.05; i=1:length(x);for y1(i)=f(A,x(i)); y2(i)=(x(i)-2)^9;end figure(3);plot(x,y1);;on hold卑微如蝼蚁、坚强似大象．共享知识分享快乐);'r'plot(x,y2,F.m文件y=f(A,x)function y=A(1); i=2:length(A);for y=x*y+A(i);;end3.2 结果第4题卑微如蝼蚁、坚强似大象．共享知识分享快乐4.1 程序;clear all n=input('请输入向量的长度n：')A=2*eye(n)-tril(ones(n,n),0); i=1:n for A(i,n)=1;end n=length(A);U=A; e=eye(n);for i=1:n-1[max_data,max_index]=max(abs(U(i:n,i))); e0=eye(n);max_index=max_index+i-1; U=e0*U; e1=eye(n); j=i+1:n fore1(j,i)=-U(j,i)/U(i,i);endU=e1*U;中把变换矩阵存到P P{i}=e0;% L{i}=e1; e=e1*e0*e;endk=1:n-2for Ldot{k}=L{k}; i=k+1:n-1forLdot{k}=P{i}*Ldot{k}*P{i};endend Ldot{n-1}=L{n-1};LL=eye(n);PP=eye(n); i=1:n-1for PP=P{i}*PP;LL=Ldot{i}*LL;endb=ones(n,2);解方程 %b=e*b;x=zeros(n,1);x(n)=b(n)/U(n,n); i=n-1:-1:1for卑微如蝼蚁、坚强似大象．共享知识分享快乐x(i)=(b(i)-U(i,:)*x)/U(i,i);end计算逆矩阵%X=U^-1*e^-1*eye(n);AN=X'; result2{n-4,1}=AN;result1{n-4,1}=x;,n)'%d:\n'fprintf(fprintf('%d ',AN);4.2 结果n=51.0625 -0.875 -0.75 -0.5 -0.0625-0.0625 0.0625 -0.75 1.125 -0.5-0.0625 0.125 0.0625 1.25 -0.5-0.0625 0.1250.25 0.06251.50.0625-0.5-0.25-0.0625 -0.125n=101.0625 -0.875 -0.75 -0.5 -0.0625 1.0625 -0.875 -0.75 -0.5 -0.0625 -0.0625 1.125 0.0625 -0.75 -0.5 -0.5 0.0625 1.125 -0.75 -0.0625 -0.0625 0.0625 0.125 1.25 1.25 -0.0625 -0.5 0.0625 0.125 -0.5-0.0625 0.250.250.0625 0.1251.5 1.5 -0.0625 0.1250.06250.0625 -0.0625 -0.125 -0.25 0.0625 -0.5 -0.0625 -0.125 -0.25 -0.5 -0.0625 -0.75 1.0625 -0.5 -0.0625 -0.875 -0.5 -0.75 1.0625 -0.875 -0.0625 -0.5 0.0625 1.125 -0.5 0.0625 1.125 -0.75 -0.0625 -0.75 1.25 0.125 0.0625 -0.0625 -0.0625 -0.5 -0.5 0.0625 0.125 1.250.25-0.0625 -0.0625 1.50.1250.0625 0.0625 0.250.1251.5-0.0625 -0.125 -0.25 0.0625-0.5 0.0625 -0.0625 -0.125 -0.25-0.5同样的方法可以算出n=20,n=30时的结果，这里就不罗列了。