顺序Gauss消去法解线性方程组.docx

上传人:p** 文档编号:810686 上传时间:2024-03-08 格式:DOCX 页数:6 大小:219.12KB
下载 相关 举报
顺序Gauss消去法解线性方程组.docx_第1页
第1页 / 共6页
顺序Gauss消去法解线性方程组.docx_第2页
第2页 / 共6页
顺序Gauss消去法解线性方程组.docx_第3页
第3页 / 共6页
顺序Gauss消去法解线性方程组.docx_第4页
第4页 / 共6页
顺序Gauss消去法解线性方程组.docx_第5页
第5页 / 共6页
顺序Gauss消去法解线性方程组.docx_第6页
第6页 / 共6页
亲,该文档总共6页,全部预览完了,如果喜欢就下载吧!
资源描述

《顺序Gauss消去法解线性方程组.docx》由会员分享,可在线阅读,更多相关《顺序Gauss消去法解线性方程组.docx(6页珍藏版)》请在第壹文秘上搜索。

1、第二章五、计算题xl+Ix2-x3=11.顺序Gauss消去法解线性方程组Xl-3jv2-3x3=-14.v1+2x2+2x3=3解:由线性方程组消去后两个方程的乐得,xl + 2x2 X3 = 1-5x2 - 2均=2-6x2 +6x3 = -1再消去最后一个方程的X2得回代得解:计并过程保留三位小数。r 1116、-18解12-331512-183-12 1r 1 12.用列卜:元索法求解线性方程纽 12-3-18 33-1T5)r-183-1-3315O-12.333116 ; 01.1670.944-15、55.167,-18O31.167-10.944-155.167r-i8O31.

2、167-10.944-155.167,回代得X3=3.000X2=2.000、O-12.3333.用直接三角分解法解方程组5,2-42k6-17O3丫须、8r20x3O=32122J3.1429.428Zx1=1.000IIXL=maxxi=ma(I2M-3,4=4:O0、2-13、1O042一1b、0O3)1.UM=N=l2l+-3l+l4l=9三f=ll=f1* Ii=I )0,5=22+(-3)2+42=29五、计算题IOx1-X2-2x3=72I.用JaCObi迭代法求解线性方程组-X1+IOx2-2x3=83,迭代两次。-X1-X2+5*3=42Xl=0.Ix22x3+7.2解8将方

3、程组改写成等价形式x2=0lxl02x3+8.3.x3=0.2x1+0.2x2+8.40.1+0.2+7.2Jacobi迭代法计算公式为.以+D =0.1 X产 +心 x,?=0.2x*+0.2设初始向量取x(0)= (0,l0, 0) 0.2峭+8.3 .产 +8.4第一次迭代X=(7.2,8.3,8.4)T:第二次迭代x;2)=0.18.3+0.28.4+7.2=9.71.=0.17.2+0.28.4+8.3=10.7.靖=0.27.2+0.2X8.3+8.4=11.50.IOx1-X2-Ix3=722.用Gauss-Seidel迭代法求解线性方程组一盯+1-2x3=83,迭代两次.-x1

4、-x2+5=42x1=0.1x2+0.zx3+7.2解S将方程组改写成等价形式,x2=OJx1+0.2x3+8.3,x3=0.2x1+0.2x2+8.4x1+=0.1蛾+0.2好+7.2得GaUSS-SeidCI迭代法计算公式xf+D=0.1x,+“+0.2Xf)+8.3铲)=0.2x;ft+n+0.2+D+8.4设初始向址取x(0)=(0,0,0)Yx;=7.2第一次迭代IXT=O.1X7.2+8.3=9.02/2)=O27.2+0.29.02+8.4=11,644x;2)=0.19.02+0.211.644+7.2=l0.4308第二次迭代J=0.1l0.4308+0.211.644+8.

5、3=11.6719Af)=O.2X10.4308+0.211.6719+8.4=l2.82053.设有方程组4x = b,其中A= 1/2 J/2解 由JaCobi迭代矩阵M =/ 。TA =1/21/2、11/2。求Jacobi迭代法解该/J程蛆时的收敛性。1/2L0-1/2-1/2、-1/20-1/2lk-l2-1/201/21/2 =( + l)(-l2)2=0, 1/2得到对应的特征方程I力r-MI=1/21/21/2所以特征值为4=N2=1/2,i=-.由谱半布的定义得迭代矩阵的p(M)=1.由迭代法收敛的充分必要条件知JaCobi迭代法发散。1 .已知/(一1)=2,/(1)=3.

6、/(2)=-4,求拉格朗日插值多项式L2(x)及/(1.5)的近似值,取5位小数。解:必)=2、昂第+3、需符-4监H(IXx-2)争+吐2)争+2)/(L5)W=0.046172 .-HI数据如卜代所小,分别用Lagrange插值法和NeWton插值法求/的:湫插值多项土并求了(2)的近似值。Xi1345解:(1) Lagrange :次插值多项式为f(xi)2654/2(X-3Xx-4)(X-5)(X_IXX_4)(x_5)(X-I)(X-3)(x-5)(X-I)(X-3)(x-4)7;-3X1-4X1-5)+6(3-iX3-4X3-5)5(4-iX4-34-5)4(5-iX5-35-4)

7、具体差商表如下:Xi%一阶差商二阶差商三阶差商1236245-1154-100.25因此M(x)=2+2(x-1)-(X-I)(X-3)+;(X-I)(X-3)(X-4)412M=5.53 .取打点Xo=O,X1=0.5,=L求函数/=”在区间0,1上的次Lagrange插值多项式L2(x),并估i*B储L(X)=JX(x-05rT)-0.5X(X-O)(XT)-1(AoXX-0.5)(-O.5X-l)(0.5-X.5-1)(i-Xl-0.5)=2(x-0.5)(X-1)+4e0Rx-1)+2e-,x(x-0.5),般!尸(步因为f(x)=ex,(x)=-e-x故截断误差R (小(M-G(XI

8、 =3!g(4M-5 心 Tl4.以100,,144为插值节点,用NeWlOn插值法计算Jnm的近似值,并利用余项估计误差.解:1)构造如下差分表:Xi一阶差商二阶差商I(K)121144IOH120.04761900.0434783-0.000094U36所以M(x)=/(XO)+小0,占卜-Xo)+/Lo,w,与卜-/)(XTl)=10+0.0476MK)O)-0.00009411362(X-1OOYx-121)令.E15,代入上式,11510.72276-_52)因y=/(%)=,所以/TH=。4/8结合W2W=(x-10)(a-121)(x-I44)得艮(叫=3!(115-10)(11

9、5-)(115-144)111OO156290.00163685.已知/(-1)=2/()=4./(1)=10./(2)=26,则2次NeWion插值多项式中/的系数为?3次NeWtOn插值多项式为?答:差商计第表如下X心)一阶差商二阶差商三阶差商-I2042I1062226165I因为22。)=()+Xo,(X-Xo)+/与,七,2(x-)(X-Xl)=2+2(x+D+2(x+D(X-O)=2x2+4x+42所以二次Ncwton插值多项式中X”的系数为2.3次Ncwton插值多项式即是在2次NCWlOn插值多项式基础上增加一项:V3(x)=N2(x)+fx0,xpx2,x3(x-x0)(x-

10、X1)(x-X2)=2x2+4x+4+1(x+1)(-0)(:-1)=2x2+4x+4+x3-X第6章函数翼近o,=x+2x+3x+4二、计算题1.已知数据表如卜所示,求其一次最小二乘拟合多项式。xi 12驰歹252.04.0 3.5 6.05.5Sx = 28SY = 24.00.52.52.04.03.56.05.5所以 6 (X) = 0.0714 + 0.8393-所以p2 (x)=14PiM7 1011=FX 10 73 fX0)p,2) = (x-4(-)+0)+ f第7章四、计算题1.求积公式为(1)求4、B使求积公式的代数精度尽量高,并求其代数褶度:(2)利用此公式求Z=J-d

11、x(保留四位小数)。解:(l)(x)=l,x,%2时,求积公式精确成立,即2A+2B22A+-B=-23IQ+ 4)得A=-,8=-,99求枳公式为f/(xZr(-l)+(l)+f2l/(X)=X3时.公式显然精确成立;考/(x)=X,时,左=一,右=一所以代数精度为3。532 1t=2x-3(2) / = dx =f.79 -1+397- 0.69291402,已知数值积分公式为,7(x)伙H勺/()+/()+2/,(0)一/,(创.试确定积分公式中的参数使其代数精确度尽量高,并指出其代数精确度。2解:F(X)=I显然精确成立:/(x)=X 时,Ja(O + )+2(1-1):fx)=x2时

12、,,2d=y=-(+A2)+2(-2)=y-2z=2=;f(x)=3时,J3dr=-=(+3)+-j2(-32):/(x)=x4时,j4d=-(+4)+-z2(-43)=-,所以,其代数精确尚为3。52126f21,3,用梯形公式和Simpson公式计算积分J()了心.解:/(O)=Ij=05J=0r212-018梯形公式:=(/(0)+7)=仃T0588(212-01152Simpson积分公式:|-tiv-(f(0)+4/(1)+f(2)=(1+2H)=1.0196jo+X6317514.=3,用复化梯形公式求也的近似值(取四位小数),并求有效数字的位数,/一1解:Tn=/(a)+)+2),LA=I_因此Iez4=ye0+e,+2(e,/3+e2z3)1.7342已知=-(b)/),由f。)=e,(x)=ex.广(X)=,TeUo.025 0.05 1232得1号。)1=吟则/倒铮=所以,至少有两:有效数字。五、综合题1 .数值求积公式-/(1)+/(2)是否为插值型求积公式?为什么?其代数精度是多少?解:1)是.因为/(X)在节点1,2处的插值多项式为P(X)=/(0f/(2)=I/(D+/(2)已知/(外公CP(X)心,所以(1)(2)2 )(x)=l,x时,求积公式精确成立:当/(x)=2时,左=9,右7.5,求积公式不成立。所以其代数精度为Lf21

展开阅读全文
相关资源
猜你喜欢
相关搜索

当前位置:首页 > 研究生考试 > 考研数学

copyright@ 2008-2023 1wenmi网站版权所有

经营许可证编号:宁ICP备2022001189号-1

本站为文档C2C交易模式,即用户上传的文档直接被用户下载,本站只是中间服务平台,本站所有文档下载所得的收益归上传人(含作者)所有。第壹文秘仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对上载内容本身不做任何修改或编辑。若文档所含内容侵犯了您的版权或隐私,请立即通知第壹文秘网,我们立即给予删除!