ex4.m

来自「包括使用修正Gram-Schmit算法实现QR分解」· M 代码 · 共 22 行

%     练习4：矩阵A的奇异值分解问题
% ——gogomx—— 2005.5

clc;
clear;

a=ones(50);
A=a+eye(50)*99;              % 产生矩阵A

%   a=[100,ones(1,49)];
%   A=toeplitz(aa);          % 另一种方法产生矩阵AA

%%%%%%%%%%%%%%  计算A的奇异值分解  %%%%%%%%%%%%%%
[U,S,V] = svd(A);            % 利用Matlab自带函数svd计算奇异值分解 
                             % A = U*S*V',其中S为奇异值组成的对角矩阵
                             % U,V为对应的正交矩阵     


disp('练习4运行完毕^_^');
disp('计算结果已存于Wokspace中，');
disp('请参阅本程序中的注释查看相关结果。');