多种矩阵分解的matlab程序

中国科学院大学矩阵分析与应用的课程中所涉及的各种矩阵分解（如lu分解）的matlab源程序，也是这门课程的大作业。对学习这门课程有所帮助，希望采纳！！！

function [gP gT]=givens（A）
Q=A;
m=size（A1）;%A的行数
gP=eye（m）;%初始化P，为防止与其余分解混淆，用gP
P1=eye（m）;%初始化每步的旋转矩阵
for i=1:m
    for j=i+1:m   %当mn时，将第n行及以下全化简至0
        if （Q（ii）^2+Q（ji）^2）==0
            continue
        end                      %若第ij位置上的元素均为0，终止此次化简
        c=Q（ii）/sqrt（Q（ii）^2+Q（ji）^2）;
        s=Q（ji）/sqrt（Q（ii）^2+Q（ji）^2）;        
        P1（ii）=c;
        P1（jj）=c;
        P1（ij）=s;
        P1（ji）=-s;%求得旋转矩阵Pi
        gP=P1*gP;    %循环求得最终旋转的P
        Q=P1*Q;    %循环求得最终的T
        P1=eye（m）;%初始化P1
    end
end
gT=gP*A;%求得T，为防止与其余分解混淆，用gT

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2017-12-02 21:12  JZFXDZY\
     文件         721  2017-12-02 21:05  JZFXDZY\givens.m
     文件         728  2017-12-02 21:03  JZFXDZY\householder.m
     文件        1280  2017-12-02 21:00  JZFXDZY\lufenjie.m
     文件         435  2017-12-02 20:51  JZFXDZY\main.m
     文件         614  2017-12-02 21:01  JZFXDZY\qrfenjie.m
     文件       17358  2017-12-02 20:56  JZFXDZY\程序清单.docx
     文件       13607  2017-12-02 21:12  JZFXDZY\程序说明.docx