Kaczmarz - Ansorge implementation method
2 views (last 30 days)
Show older comments
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:
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/149755/image.png)
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 Comments
Answers (0)
See Also
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!