function?[xk]=Landweber（Abx0）
alfa=1/sum（diag（A*A‘））;
k=1;
L=200;
x=x0;
while?k????x1=x;
????x=x+alfa*A‘*（b-A*x）;