Please just convert into Matlab file. It isa function only. Code is def Function(self,dz,lam, Angle, alpha): dx = abs(self.X[0,1]-self.X[0,0]) dy = abs(self.Y[1,0]-self.Y[0,0]) k0 = 2*pi/lam dv_x = 1/(self.Nx*dx) kappa_x = (arange(self.Nx)-self.Nx/2) * dv_x * 2*pi dv_y = 1/(self.Ny*dy) kappa_y = (arange(self.Ny)-self.Ny/2) * dv_y * 2*pi Kappa_x,Kappa_y = meshgrid(kappa_x,kappa_y) kx = Kappa_x + k0 * np.sin(Angle) ky = Kappa_y.copy() kz = sqrt(k0**2 - kx**2 - ky**2) kz_tilde = kz * np.cos(Angle) + kx * np.sin(Angle) R = dz/np.cos(Angle) guard = ones((self.Ny,self.Nx)) R_k = sqrt(Kappa_x**2 + Kappa_y**2) guard[R_k >= 2*pi/lam] = 0 filter = ones((self.Ny,self.Nx)) if alpha: vx_max = 1/(sqrt(1+(2*dz*dv_x)**2) * lam) vy_max = 1/(sqrt(1+(2*dz*dv_y)**2) * lam) Order = 20 filter = np.exp(-(abs(Kappa_x)/(vx_max*2*pi))**Order - (abs(Kappa_y)/(vy_max*2*pi))**Order) H0 = exp(-1j*kz_tilde * R) * guard * filter return RR(0)
Please just convert into Matlab file. It isa function only.
Code is
def Function(self,dz,lam, Angle, alpha):
dx = abs(self.X[0,1]-self.X[0,0])
dy = abs(self.Y[1,0]-self.Y[0,0])
k0 = 2*pi/lam
dv_x = 1/(self.Nx*dx)
kappa_x = (arange(self.Nx)-self.Nx/2) * dv_x * 2*pi
dv_y = 1/(self.Ny*dy)
kappa_y = (arange(self.Ny)-self.Ny/2) * dv_y * 2*pi
Kappa_x,Kappa_y = meshgrid(kappa_x,kappa_y)
kx = Kappa_x + k0 * np.sin(Angle)
ky = Kappa_y.copy()
kz = sqrt(k0**2 - kx**2 - ky**2)
kz_tilde = kz * np.cos(Angle) + kx * np.sin(Angle)
R = dz/np.cos(Angle)
guard = ones((self.Ny,self.Nx))
R_k = sqrt(Kappa_x**2 + Kappa_y**2)
guard[R_k >= 2*pi/lam] = 0
filter = ones((self.Ny,self.Nx))
if alpha:
vx_max = 1/(sqrt(1+(2*dz*dv_x)**2) * lam)
vy_max = 1/(sqrt(1+(2*dz*dv_y)**2) * lam)
Order = 20
filter = np.exp(-(abs(Kappa_x)/(vx_max*2*pi))**Order - (abs(Kappa_y)/(vy_max*2*pi))**Order)
H0 = exp(-1j*kz_tilde * R) * guard * filter
return RR(0)
Step by step
Solved in 2 steps
