《作业1-线性方程组求解》由会员分享,可在线阅读,更多相关《作业1-线性方程组求解(5页珍藏版)》请在金锄头文库上搜索。
1、1 分别用矩阵求逆、矩阵除法以及矩阵分解求线性方程的解。2 下面是一个线性病态方程组:(1) 求方程的解。(2) 将方程右边向量元素b3改为0.5:0.01:0.6,再求解,并比较b3的变化和解的相对变化。(3) 计算系数矩阵A和条件数并分析结论。解: 1-1 A=2,3,5;3,7,4;1,-7,1;B=10,3,5X=AB.B = 10 3 5X = -1.8060 -0.5373 3.0448 A=2,3,5;3,7,4;1,-7,1;B=10,3,5B = 10 3 5 L,U=lu(A);X=U(LB)X = -1.8060 -0.5373 3.04481-2 clear all A
2、=6,5,-2,5;9,-1,4,-1;3,4,2,-2;3,-9,0,2A = 6 5 -2 5 9 -1 4 -1 3 4 2 -2 3 -9 0 2 b=-4,13,1,-1b = -4 13 1 -1 x=Abx = -5.3333 1.0000 18.5000 12.0000 L,U=lu(A);x=U(Lb)x = -5.3333 1.0000 18.5000 12.00002-1 clear all A=1/2,1/3,1/4;1/3,1/4,1/5;1/4,1/5,1/6; b=0.95;0.67;0.52b= 0.9500 0.6700 0.5200 B=sym(b) B =
3、 19/20 67/100 13/25 x=AB x = 6/5 3/5 3/52-2 A=1/2,1/3,1/4;1/3,1/4,1/5;1/4,1/5,1/6; b=0.95;0.67;0.52B=sym(b)b = 0.9500 0.6700 0.5200 B = 19/20 67/100 13/25 x=AB x = 6/5 3/5 3/5 b3=0.5:0.01:0.6b1=0.95b2=0.67for i=1:length(b3)b(1)=b1;b(2)=b2;b(3)=b3(i);x0=Ab;x1(i)=x0(1);x2(i)=x0(2);x3(i)=x0(3);endb3 = 0.5000 0.5100 0.5200 0.5300 0.5400 0.5500 0.5600 0.5700 0.5800 0.5900 0.6000b1 = 0.9500b2 = 0.6700 plot(b3,x3)在高温或低温情况下进行的高处作业。高温是指作业地点具有生产性热源,其气温高于本地区夏季室外通风设计计算温度的气温2及以上时的温度。低温是指作业地点的气温低于5。