北京科技大学计算方法大作业.doc

计算方法大作业 机械电子工程系 老师：廖福成 注：本文本只有程序题，证明题全部在手写已交到理化楼204了。 2. 证明方程 在上有一实根，并用二分法求这个根。要求。请给出程序和运行结果。 证明： 设f(x)=x3-x-1 则f(1)= -1， f(2)= 5，f(1)*f(2)= -50 因此，方程在[1,2]上必有一实根。 二分法求解程序： %预先定义homework2.m文件如下： function lc=homework2(x) lc=x^3-x-1; 在MALAB窗口运行： clear a=1;b=2;tol=10^(-3);N=10000; k=0; fa=homework2(a); % f 需事先定义 for k=1:N p=(a+b)/2; fp=homework2(p); if( fp==0 || (b-a)/2tol) break end if fa*fp0 b=p; else a=p; end end k,p 程序运行结果： k = 10 p = 1.325195312500000 3. 用Newton迭代法求方程 的一个正根，计算结果精确到7位有效数字. 要求给出程序和运行结果. 解： 取迭代初值 ,并设,则． 牛顿迭代函数为 牛顿迭格式为： Matlab程序如下： %定义zuoye3.m文件 function x=zuoye3(fname,dfname,x0,e,N) if nargin5,N=500;end if nargin4,e=1e-7;end x=x0;x0=x+2*e;k=0; while abs(x0-x)ekN, k=k+1; x0=x;x=x0-feval(fname,x0)/feval(dfname,x0); disp(x) end if k==N,warning(已达上限次数);end 在Matlab窗口中执行： zuoye3(inline(x^3+2*x^2+10*x-20),inline(3*x^2+4*x+10),1,1e-7) 结果如下： 1.41176470588235 1.36933647058824 1.36880818861753 1.36880810782137 ans = 1.36880810782137 4. 用牛顿迭代法求方程在附近的根. 要求给出程序和运行结果. 解：令：,则． 牛顿迭代函数为 牛顿迭格式为： Matalb程序如下： %定义zuoye4.m文件 function x=zuoye4(fname,dfname,x0,e,N) if nargin5,N=500;end if nargin4,e=1e-7;end x=x0;x0=x+2*e;k=0; while abs(x0-x)ekN, k=k+1; x0=x;x=x0-feval(fname,x0)/feval(dfname,x0); disp(x) end if k==N,warning(已达上限次数);end 在Matlab窗口执行： zuoye4(inline(x^3-x-1),inline(3*x^2-1),1,1e-7) 结果如下： 1.50000000000000 1.34782608695652 1.32520039895091 1.32471817399905 1.32471795724479 1.32471795724475 ans = 1.32471795724475 6. 编写用全主元Gauss消去法解线性方程组的程序，并求解 解： Matlab程序如下： A=[2 -1 4 -3 1;-1 1 2 1 3;4 2 3 3 -1;-3 1 3 2 4;1 3 -1 4 4] b=[11 14 4 16 18] function x=zuoye6(A,b) [n,v]=size(b); D=[A,b;eye(n),zeros(n,v)],[s1,m]=size(D); for k=1:(n-1) s=abs(A(k,k));p=k;q=k; for i=k:n for j=k:n if abs(A(i,j))s s=abs(A(i,j));p=i;q=j; end end end if pk t=D(k,:); D(k,:)=D(p,:); D(p,:)=t; end if qk t1=D(:,k); D(:,k)=D(:,q); D(:,q)=t1; end h=D(k+1:n,k)/D(k,k); D(k+1:n,