奇異值的計算可以通過以下步驟進行：
計算矩陣A的QR分解，表示爲A=QR，其中Q是m×m的正交矩陣，R是m×n的上三角矩陣。
對R進行LU分解，即R=LU，其中L是m×n的下三角矩陣，U是m×n的上三角矩陣。
定義對角矩陣S，其元素爲sqrt(diag(U'Q))的前n個元素。
計算左奇異矩陣U，公式爲U(1:n,:)=Q(1:n,:)/S(1:n)。
計算右奇異矩陣V，公式爲V(:,1:n)=(R(1:n,:)*S(1:n))/S(:,1:n)。
通過以上步驟，可以得到矩陣A的奇異值分解A=USV^T。需要注意的是，當矩陣A的秩爲r時，奇異值個數也爲r，其餘的奇異值爲零。