﻿功能:对矩阵进行右极分解(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 =
[ 0.92015234614013    0.20444337865811    0.16283980737484   -0.29155747623717
 -0.18243018489557    0.97324931054763    0.04541428309339    0.13207024664420
  0.00764400550222   -0.10093879907229    0.88770633619969    0.44914406229253
  0.34638195205822   -0.02828503638451   -0.42825016134982    0.83415903717990 ]
u =
[ 90.1600096663750    59.6469572679696    32.5383853704088    22.3791819318738
  59.6469572679696    79.7602499219720    52.5635199143592    40.8190412412902
  32.5383853704088    52.5635199143592    70.3253422679695    45.1004130424659
  22.3791819318738    40.8190412412902    45.1004130424659    86.1282027383202 ]