﻿功能:对矩阵进行右极分解(Right Polar Decomposition)

格式:
[R,U]=RPD(A)
A:为方阵变量

说明:任何一个可逆方阵A都可以唯一的表示为A=R*U,其中R为一个正交矩阵,U为正定对称矩阵.

例子:

a =
[ 93.9290964947683    67.8490794579727    38.9889378282190    11.1702606133978
  46.0367843257435    74.5233182211049    54.3716235339509    49.0676413518692
  33.6044989217094    57.3996375116518    77.6278596732895    74.7707348665086
  34.2759118575025    29.9438441777340    17.2880263148286    59.1275361641904 ]

[r,u]=rpd(a)//回车得到如下结果
r =
[ 6.24113649081953    0.32410435452631    1.23994655255781    0.15836154717303
  3.05892279706310    0.35598614074302    1.72915475294376    0.69563530069241
  2.23285725406363    0.27418901795538    2.46876171411551    1.06002981188487
  2.27746941291922    0.14303702226752    0.54980283700074    0.83825511611250 ]
u =
[ 15.0499987675217   -2.2333056572E-14    1.2945025317E-22   -3.0653462829E-41
 -2.2333056572E-14    209.343313381690   -3.1288116663E-28   -6.2341109881E-27
  1.2945025317E-22   -3.1288116663E-28    31.4440471226691    1.9191110626E-17
 -3.0653462829E-41   -6.2341109881E-27    1.9191110626E-17    70.5364453227561 ]