pleas any one can help me .....noisy image for heart by CT cardiac gating

I have a noisy CT image for the heart using cardicgating by motion artifacts and i need to restore it using image restoration , and it have a model r(t)=(0.05*sqrt(t)/T)x+(0.04*sqrt(t)/T)y so i Calculate the degradation function H , i do it but is my code right ?