spaugment--生成最小二乘增广矩阵

【功能简介】生成最小二乘增广矩阵。

【语法格式】

S=spaugment(A,c) 

S返回一个稀疏、对称的非正定方阵，S=[c*I,A;A',0]，c为残留比例因子，I为单位矩阵。函数的功能是求解一个最小二乘问题：

min(norm(b-A*x))

norm函数用于求第二范式，如向量[1,1]的第二范式为sqrt(12+12)=1.4142。最小二乘问题可以表述为：

r=b-A*x

S*[r/c;x]=[b;0]

【实例3.20】生成一个最小二乘增广矩阵。

>> a=[1,2,3;4,5,6];  
>> b=spaugment(a,1);        %生成最小二乘增广矩阵  
>> full(b)              %将稀疏矩阵转化为普通矩阵  
ans =  
     1     0     1     2     3  
     0     1     4     5     6  
     1     4     0     0     0  
     2     5     0     0     0  
     3     6     0     0     0