功能: 使用QR方法求解线性方程组
 
格式:
[x, er] = SparseSolveByQR(A, b, n)
[x, er] = SparseSolveByQR(A, b)

A : 一个列数为3的系数矩阵
b : 一个行数与A一致的常数项矩阵
n : 数值,尝试求解次数,默认为16

x : 求解成功返回的解
er: 求解成功返回的误差平方和

说明:
1. 本函数本质是求解如下方程组
$$
Ax=b
$$

2. 本函数主要求解稀疏矩阵,其中要求A的第1行为非0元素行索引,第2行为非0元素列索引,第3行为对应元素值。

3. A当中存在非数字,或者索引值异常,或者存在重复索引,会返回失败。

4. A当中的行列索引,均从0开始。

5. 如果b有多列，则x也有多列。每一列单独使用1里的方程进行求解。

6. 本函数采用QR方法求解，此方法适合：适合A不是方阵的情况。

例子:
A =
[ 40.0000000000000    5.00000000000000    0.00000000000000    0.00000000000000
  5.00000000000000    15.0000000000000    2.00000000000000    0.00000000000000
  0.00000000000000    2.00000000000000    10.0000000000000    0.00000000000000
  0.00000000000000    0.00000000000000    0.00000000000000    6.00000000000000 ]
b =
[ 4.19307872848263
  9.59621052704575
  2.47985608991229
  1.13943188038628 ]
//假如系统存在如上形式,求A*x=b的解x
A = [0,0,40
     0,1,5
     1,0,5
     1,1,15
     1,2,2
     2,1,2
     2,2,10
     3,3,6];
b = [ 4.19307872848263
      9.59621052704575
      2.47985608991229
      1.13943188038628 ];

[x, er] = SparseSolveByQR(A, b)//回车得到如下结果
x =
[ 0.02811762592223
  0.61367473831863
  0.12525066132750
  0.18990531339771 ]
er =
[ 4.3387349787E-30 ]