Kaczmarz - Ansorge implementation method

1 次查看(过去 30 天)
Hi there, I have to implement the Kaczmarz-Ansorge program, but they way I make it, I am sure that it is the same as original Kaczmarz.
Basics:
- we take a random value between 1 and M ( this is ik )
- then we make the program on only 1 projection, that random value.
Formula:
Here is what I made, am I on the right path ?
%%%%%Pas 2 Algoritm - K-A %%%%%
% for ik = 1:M
% nlin(ik) = A(ik,:) * A(ik,:)';
% end
%
% for cont = 1:NITER
%
% for jj = 1:M
% ik = M -jj+1;
% x = x - (A(ik,:) * x - b(ik))/nlin(ik) * A(ik,:)';
% end
%
% end
%%%%%End Pas 2 Algoritm - K-A %%%%%
Where NITER are the number of iterations that I have declared in the beginning of program.
Thank you!
  2 个评论
Cristian
Cristian 2012-11-23
My question is if I am on the right path on this. This is about reconstruction of images.

请先登录,再进行评论。

回答(0 个)

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by