矩阵理论分析matlab实验..docxVIP

矩阵理论实验专业：控制科学与工程姓名：刘思学号：211508010020矩阵的LU分解原理定义1.1设，若A可以表示成一个下三角矩阵与一个上三角矩阵的乘积（1.1）则称其为矩阵的分解（三角分解）.定理1.1设，如果的顺序主子式，，…，（1.2）则存在唯一的主对角线上元素全为一的下三角矩阵与唯一的上三角矩阵，使得.算法此算法较为简单，直接叙述如下：由公式（1.3），可以按照从的第一行求起，再求的第一列，再求的下一行，再求的下一列，循环下去，直到求出.流程图 程序使用MATLAB的自带程序进行LU分解A=[2 1 1;4 1 0;-2 2 1];[L,U,P] = lu(A);运行结果A=[2 1 1;4 1 0;-2 2 1]; [L,U,P]=lu(A)L = 1.0000 0 0 -0.5000 1.0000 0 0.5000 0.2000 1.0000U = 4.0000 1.0000 0 0 2.5000 1.0000 0 0 0.8000P = 0 1 0 0 0 1 1 0 0自编程序进行LU分解function [ L,U ] = LSLU( A )%UNTITLED 刘思编写的矩阵的LU分解程序[n,n]=size(A); %获取矩阵A的阶L=zeros(n,n); %将矩阵LU置零U=zeros(n,n); for i=1:n %对L矩阵的对角线元素赋值1L(i,i)=1; end for k=1:n %根据算法求L,U矩阵for j=k:n U(k,j)=A(k,j)-sum(L(k,1:k-1).*U(1:k-1,j)); end for i=k+1:n L(i,k)=(A(i,k)-sum(L(i,1:k-1).*U(1:k-1,k)))/U(k,k); endendend运行结果： A=[2 1 1;4 1 0;-2 2 1]A = 2 1 1 4 1 0 -2 2 1 A=[2 1 1;4 1 0;-2 2 1]A = 2 1 1 4 1 0 -2 2 1 [L,U]=LSLU(A)L = 1 0 0 2 1 0 -1 -3 1U = 2 1 1 0 -1 -2 0 0 -4与MATLAB 自带分解分解程序一致。矩阵的QR分解原理定理2.1设，且，则比存在非奇异上三角矩阵及矩阵，，使得（2.1）则称其为矩阵的分解.算法1 利用 Schmidt正交化求矩阵的 分解，Schmidt正交化方法是矩阵的 QR 分解最常用的方法. 主要依据下面的两个结论 :结论 1 设 是 n阶实非奇异矩阵 ,则存在正交矩阵 和实非奇异上三角矩阵 使 有 分解 ;且除去相差一个对角元素的绝对值（模）全等于 1的对角矩阵因子外 ,分解是唯一的. 2设 是实矩阵 ,且其个列向量线性无关 ,则 有分解 ,其中 是 实矩阵 ,且满足 , 是阶实非奇异上三角矩阵该分解除去相差一个对角元素的绝对值 （模） 全等于1的对角矩阵因子外是唯一的。步骤 : 1、写出矩阵的列向量；2、把列向量按照 Schmidt正交化方法进行正交； 3、得出矩阵的 ；4、对的列向量单位化得到，并在的每一行乘以每一列模得到.流程图程序1.MATLAB自带函数进行QR分解B=[1 1;0 1;1 1];[Q,R]=qr(B)[Q,R,E]=qr(B)运行结果B=[1 1;0 1;1 1];[Q,R]=qr(B)[Q,R,E]=qr(B)Q = -0.7071 0 -0.7071 0 1.0000 0 -0.7071 0 0.7071R = -1.4142 -1.4142 0 1.0000 0 0Q = -0.5774 0.4082 -0.7071 -0.5774 -0.8165 0 -0.5774 0.4082 0.7071R = -1.7321 -1.1547 0 0.8165 0 0E =10自编算法进行QR分解function [ Q,R ] = LSQR( A )%LSQR 刘思编写的矩阵的QR分解程序A=[1,1;0,1;1,1];[m,n]=size(A); %得到A的维数if (m=n)(n==rank(A)) R