矩阵理论及应用课程设计报告
目录
1、第一题 QR分解········································3
2、第二题 奇异值分解·····································5
3、第三题 图像压缩处理···································9
4、附录代码··············································12
1、对以下矩阵A、B作QR分解。
1.1QR分解的原理
Schmidt正交化方法是矩阵的QR分解最常用的方法
用Schmidt方法对矩阵进行QR分解时,矩阵必须是列满秩矩阵。
1.2QR分解的算法
1、写出矩阵的列向量
2、列向量按照Schmidt正交化正交
3、得到Q’,单位化得到Q
4、根据R=QTA可计算出R
1.3QR分解流程图:
1.4MATLAB代码
见附件三
1.5解题运行与分析
(1)
运用附件三代码
A=[2 3 4;8 4 9;5 3 1];
[Q1,R1]=QRDecmsch1(A)
结果如下:
Q1 =
0.2074 0.9548 0.2129
0.8296 -0.2870 0.4790
0.5185 0.0773 -0.8516
R1 =
9.6437 5.4958 8.8141
0 1.9483 1.3136
0 0 4.3112
运用MATLAB自带的qr(A)进行验证
A=[2 3 4;8 4 9;5 3 1];
[Q1,R1]=qr(A)
Q1 =
-0.2074 0.9548 -0.2129
-0.8296 -0.2870 -0.4790
-0.5185 0.0773 0.8516
R1 =
-9.6437 -5.4958 -8.8141
0 1.9483 1.3136
0 0 -4.3112
(2)
运用附件三代码
A=[1 2 3 4;3 5 6 2;3 6 9 12];
[Q1,R1]=QRDecmsch1(A)
结果如下:
Q1 =
0.2294 0