MATLAB使用非负最小二乘法NNLS求解线性方程

function [o_x] = NNLS(i_A,i_b,i_iter) %参考文献《基于反卷积的近场声源稳健高分辨被动定位方法》 %详细介绍请阅读文献p29-30 A=i_A; b=i_b; iter=i_iter; x=zeros(length(b),iter); r=zeros(length(b),iter); w=zeros(length(b),iter); w_1=zeros(length(b),iter); lamda=zeros(iter); for i=1:iter r(:,i)=A*x(:,i)-b; w(:,i)=-A.'* r(:,i); for l=1:length(b) if(w(l,i)<0 && x(l,i)==0) w_1(l,i)=0; else w_1(l,i)=w(l,i); end end g(:,i)=A*w_1(:,i); lamda(i)=-(g(:,i).'*r(:,i))/(g(:,i).'*g(:,i)); for l=1:length(b) x(l,i+1)=max(x(l,i)+lamda(i)* w_1(l,i),0); end end o_x=x(:,iter+1); end