清华大学数值分析A第五次作业.docxVIP

清华大学数值分析A第五次作业 清华大学数值分析A第五次作业 PAGE 清华大学数值分析A第五次作业 定义f∶R2→R如下：fx1,x2 证明：? 同理，? 而沿着x1=x ∴f在点(0,0)处不可导. 2.给出一个线性映射B:R2→R2 解：令Bx=y:y=Ax，其中A 易知A1=，则存在m,n 则 B ∴B:R 而对任意u,v∈R ∴B:R2→ 证明由Gx=ln1+e 证明：①压缩性： 要证明“G:R→R，在任何闭区间a 而Gx-G ∴只需取L= 即可满足在任何闭区间a,b 压缩性得证 ②没有不动点： 假设Gx有不动点 则Gx*=x*，即ln1+e ∴原假设不成立，即Gx 9. 给定非线性方程组 x (1)应用压缩映射定理证明G=g1 (2)用不动点迭代法求方程组的解，使xk 解：（1）证明：依题意知需证明G=g 映内性： 对 D 显然，0＜0.7sinx1 0.7cosx1 又 0.7cosx ∴0＜ 即g 映内性得证. 压缩性： 假设m=m1 n=n1 其中0 则只需证明Gm- 即0.7sin =0.7sinm 这里取1范数进行验证， Gm-Gn1=0.7sinm 不妨假设m1＞ 则sin 假设f 则 sin =fξ cosξ+sinξ ∴sin 易知无论m1＞n1或m1 均有0.7sinm 同理易得0.2cosm ∴0.7sinm ≤0.72|m1-n 又m1- 取L=0.99，则有 则由压缩映射原理可得G=g1 （2）使用matlab进行计算， 程序： x=[;]; y=[*sin(x(1,1))+*cos(x(2,1));*cos(x(1,1))*sin(x(2,1))]; k=1; f=abs(x(1,1)-y(1,1))+abs(x(2,1)-y(2,1)); while f k=k+1; x=y; y=[*sin(x(1,1))+*cos(x(2,1));*cos(x(1,1))*sin(x(2,1))]; f=abs(x(1,1)-y(1,1))+abs(x(2,1)-y(2,1)); end disp(k); disp(x); 运行结果： k=9，即共迭代9次 x=0.5258 10(2). 用牛顿法求非线性方程组 x12+x2 xk 解：记f1x=x12 Fx= Fx 牛顿方程为： Fx 即2x1 则? 则x 用matlab求解，程序如下： x=[;]; Fd=[2*x(1,1),2*x(2,1);2*x(1,1),-2*x(2,1)]; F=-[x(1,1)*x(1,1)+x(2,1)*x(2,1)-4;x(1,1)*x(1,1)-x(2,1)*x(2,1)-1]; y=x+inv(Fd)*F; f=abs(x(1,1)-y(1,1))+abs(x(2,1)-y(2,1)); k=1; while f x=y; k=k+1; Fd=[2*x(1,1),2*x(2,1);2*x(1,1),-2*x(2,1)]; F=-[x(1,1)*x(1,1)+x(2,1)*x(2,1)-4;x(1,1)*x(1,1)-x(2,1)*x(2,1)-1]; y=x+inv(Fd)*F; f=abs(x(1,1)-y(1,1))+abs(x(2,1)-y(2,1)); end disp(k); disp(y); 迭代结果如下： 迭代次数：3 x=