﻿功能:非奇异方阵的QR分解.

格式:
[Q,R]=QR(A)
A:为非奇异矩阵变量
Q:返回的正交矩阵
R:返回的上三角矩阵

说明:A=Q*R

例子:

a =
 [ 59.5926120223443     34.2398770312964     40.3722222616767    28.2503625975225
   11.736761271831      69.6454248715404     37.6082609582731    96.1153537948222
   80.2028119471869     6.89348467015358     12.8627523374105    32.1448729988862
   87.3011022747034     71.067518494589      60.0597060099522    9.06287506644748 ]

[q,r]=qr(a)//回车后得到如下解
q =
 [ 0.447382343876396    0.0253173625801201    0.694679063826229   0.562698025428526
   0.088111925104754    0.823167345938837    -0.40287872761467    0.390282641261342
   0.602110241131175   -0.47565906907231     -0.584350896084464   0.26409418353971
   0.655399561006074    0.309035599472391     0.116806503938702  -0.679194121241786]
r =
 [ 133.202869621535     72.1831670405959      68.4834555603685     46.4021838617169
   0                    76.8801476099097      44.4223129959617     67.3449960231626
   0                    0                     12.3931726747868    -36.8231785883689
   0                    0                     0                    55.7423999190258 ]