# ADD YOUR CODE HERE

# calculate P$\_$s$($k$,$l$)$ from P$\_$g$($k$,$l$)$ as described in the lecture slides

noise$\_$var $=$ noise$\_$Std$**2$

kernel $=$ np$.$ones$($size$\_$Blur$)/$ np$.$prod$($size$\_$Blur$)$

H $=$ ft$.$fft$2($kernel$)$

H$\_$squared $=$ np$.$abs$($H$)**2$

P$\_$s $=($P$\_$g $-$ noise$\_$var$)/$ H$\_$squared

# obtain the restoration filter frequency response Phi$($k$,$l$)$ as described in the lecture slides

# take the inverse DFT of Phi$($k$,$l$)$ to obtain the impulse response of the restoration filter

#Phi$\_$DFT$=$ ft$.$ifft$2($Phi$)$ Could you please do the parts in Step $4$ one by one and explain how H$($k$,$l$)$ is output with code?