matlab之列主元素高斯消去法,高斯消元法,列主元素消元法及lu分解法的matlab程序.doc...

高斯消元法,列主元素消元法及lu分解法的matlab程序.doc

还剩

5页未读，

继续阅读

下载文档到电脑，马上远离加班熬夜！

亲，喜欢就下载吧，价低环保！

内容要点：

1. 在 MATLAB 工作窗口输入程序>> A=[0 -1 -1 1；1 -1 1 -3；2 -2 -4 6；1 -2 -4 1]； b=[0；1；-1；-1]； [RA,RB,n,X]=liezhu(A,b)运行后输出结果请注意：因为 RA=RB=n，所以此方程组有唯一解.RA = 4，RB = 4，n = 4，X =[0 -0.5 0.5 0]’例 2.3 将矩阵 直接进行 LU 分解，然后解矩阵方程 。AbAX?，???????30142???????512解 1. 在 MATLAB 工作窗口输入程序>> A=[1 0 2 0；0 1 0 1；1 2 4 3；0 1 0 3]； hl=zhjLU(A)运行后输出结果请注意：因为 A 的各阶主子式都不等于零，所以 A 能进行 LU 分解.A 的秩 RA和各阶顺序主子式值 hl 依次如下：RA = 4 U = 1 0 2 0 0 1 0 1 0 0 2 1 0 0 0 2 2. 再在工作窗口输入>> U=[1 0 2 0；0 1 0 1；0 0 2 1；0 0 0 2]； L=[1 0 0 0；0 1 0 0；1 2 1 0；0 1 0 1]；b=[1；2；-1；5]；U1=inv(U)； L1=inv(L)； X=U1*L1*b,x=A\b运行后输出方程组的解X = 8.500000000000000.50000000000000-3.750000000000001.50000000000000L = 1 0 0 0 0 1 0 0 1 2 1 0 0 1 0 1 hl = 1 1 2 4 0 0.9300 -3.97110 0 -0.0837]；>> L= [1.0000 0 00.5896 1.0000 01.0677 1.5696 1.0000]；>> b=[0.05；1.03；-0.53]；U1=inv(U)； L1=inv(L)； X=U1*L1*b,x=A\bX =-111.8440110.953125.7324x =1.4531-1.5892-0.2749例 2.1： 用高斯消元法求解下面的非齐次线性方程组。 ????????.142,60,3431421xx解 1. 在 MATLAB 工作窗口输入程序>> A=[1 -1 1 -3； 0 -1 -1 1；2 -2 -4 6；1 -2 -4 1]； b=[1；0； -1；-1]； [RA,RB,n,X] =gaus (A,b)运行后输出结果请注意：因为 RA=RB=n，所以此方程组有唯一解.RA = 4，RB = 4，n = 4，X =[0 -0.5 0.5 0]’例 2.2 用列主元素消元法解线性方程组。.????????142,63,04142xx解 L(1,1)=1；L(i,i)=1；if i>jL(1,1)=1；L(2,1)=A(2,1)/U(1,1)； L(i,1)=A(i,1)/U(1,1)；L(i,k)=(A(i,k)- L(i,1:k-1)*U(1:k-1,k))/U(k,k)；elseU(k,j)=A(k,j)-L(k,1:k-1)*U(1:k-1,j)；endendendendhl；RA,U,Lendend运行命令及结果a=[2.51 1.48 4.53；1.48 0.93 -1.30 ；2.68 3.04 -1.48]；hl=zhjLU(A)请注意：因为 A 的各阶主子式都不等于零，所以 A

发表评论

暂无评论，赶快抢占沙发吧。