数值分析作业.docVIP

数值分析作业及答案 Chap1 1、写出下列语句的运行结果。在MATLAB上执行它们以验证所得解答。 a=[1 2 3 ;4 5 6 ]’ b=[9;7;5;3;1] c=b(2:4) d=b(4:-1:1) e=sort(b) f=[3,b’] 解：a= b= c= d= e= F=[3 9 7 5 3 1] 3、给定一向量： a=[4 -1 2 -8 4 5 -3 -1 6 -7] 写一段程序计算a中正数的和。运行程序并显示结果。 解：a=[4 -1 2 -8 4 5 -3 -1 6 -7]; s=0; for i=1:length(a) if a(i)0 s=s+a(i); end end s 6、编写一个函数M文件fun_es(x)，计算如下函数： 其中参数可以为标量，也可以为向量。在MATLAB里键入如下命令检验此函数: fun_es(3) fun_es([1 2 3]) 解：function y=fun_es(x) y=0.5*exp(x/3)-x.^2.*sin(x); chap2 1、设，的相对误差为，求Lnx的误差。 解：Lnx-Lnx*=dLnx=dx/x= 2、设的相对误差为2%，求的相对误差。 解：dLnf(x)=xf’(x)/f(x)dLnx=4% 5、计算球体积要使相对误差限为1%，问度量半径R时允许的相对误差限是多少？ 解：dLnf(x)=xf’(x)/f(x)dLnx=3dLnx=1% dLnx=0.33% 9、正方形的边长大约为100cm，应怎样测量才能使其面积误差不超过1cm2? 解：s=x2 s-s*=2x(x-x*)=1 x-x*=1/(2x)=1/200=0.5*10-2 即测量边的误差不超过0.005cm 10、设，假定g是准确的，而对t的测量有±0.1秒的误差，证明当t增加时S的绝对误差增加，而相对误差却减少。 解：s的绝对误差为： gtdt 当t增加时，gtdt增大。 S的相对误差为：tgt/(gt2/2)dLnt=2/tdt 当t增加时，2/tdt变小。 Chap3 1.当 解：Φ(x)=l0(x)y0+l1(x)y1+l2(x)y2 =0+(x-1)(x-2)/((-1-1)(-1-2))*(-3)+(x+1)(x-1)/((2-1)(2+1)*4 =(5x2+9x-14)/6 2.求不大于3次的插值多项式，此多项式在处与的值相同 解：y0=1 y1=9 y2=36 y3=49 Φ(x)=l0(x)y0+l1(x)y1+l2(x)y2+l3(x)y3 =(x-3)(x-6)(x-7)/(-60)*1+(x-1)(x-6)(x-7)/24*9+(x-1)(x-3)(x-7)/(-15)*36+ (x-1)(x-3)(x-6)/24*49 15.给定数据 1 2 2 3 0 -1 试构造Hermite插值多项式并计算的值 解：a0(x)=(1+2*(x-1)/(2-1))((x-2)/(1-2))2=(2x-1)(x-2)2 a1(x)=(1+2*(x-2)/(1-2))((x-1)/(2-1))2=(2x-1)(x-2)2=(5-2x)(x-1)2 b0(x)=(x-1)((x-2)/(1-2))2=(x-1)(x-2)2 b1(x)=(x-2)((x-1)/(2-1))2=(x-2)(x-1)2 H3(x)=2a0(x)+3a1(x)-b1(x) Chap4 3．用复合梯形求积法计算下列积分，取n=2,4,8,16: 解：（a） I2=h/2(f1+2f2+f3)=pi/16(tan(0)+2tan(pi/8)+tan(pi/4))=0.35901 I4=h/2(f1+2f2+2f3+2f4+f5)=0.34976 I8=h/2(f1+2(f2+f3+f4+f5+f6+f7+f8)+f9)=0.34737 I16=h/2(f1+2(f2+f3+f4+f5+f6+f7+f8+f9+f10+f11+f12+f13+f14+f15+f16)+f17)=0.34677 (b) I2=h/2(f1+2f2+f3)=1/4(exp(0)+2*exp(0.5)+exp(1))=1.75393 I4=h/2(f1+2f2+2f3+2f4+f5)=1.72722 I8=h/2(f1+2(f2+f3+f4+f5+f6+f7+f8)+f9)=1.7205