How to Use Harmonic Analysis to Determine g(x,y) in Different Cases?

Click For Summary

Homework Help Overview

The discussion revolves around using harmonic analysis to determine the function g(x,y) based on the complex amplitudes of a monochromatic wave at different planes. The original poster presents three specific cases of the function f(x,y) and seeks assistance in solving the problem.

Discussion Character

  • Exploratory, Conceptual clarification, Mathematical reasoning, Assumption checking

Approaches and Questions Raised

  • Participants discuss the application of Fourier transforms and the Fraunhofer approximation in the context of the problem. There are attempts to clarify the conditions under which the approximation is valid and the implications of the input functions extending over all space.

Discussion Status

The conversation includes various suggestions for approaches, such as starting with the Helmholtz equation and using the angular spectrum method. Some participants express a desire for specific solutions to the cases presented, while others emphasize the limitations of the Fraunhofer approximation in this context.

Contextual Notes

Participants note that the functions f(x,y) are defined for all values of x and y, which raises questions about the applicability of certain mathematical methods. There is also mention of previous related discussions that may provide additional context.

Amirjahi
Messages
7
Reaction score
0

Homework Statement

The complex amplitudes of a monochromatic wave of wavelength λ in the z=0 and z=d planes are f(x,y) and g(x,y), resprctively. Assume d=104λ, use harmonic analysis to determine g(x,y) in the following cases:

(a) f(x,y) = exp{(-jπ/λ)(x+y)}
(b) f(x,y) = cos(πx/2λ)
(c) f(x,y)=cos2(πy/2λ)

can anyone help me to solve this problem?
thank you
 
Last edited:
Physics news on Phys.org
The other two sections of the template you have deleted are compulsory before you can hope to receive help. What equations in harmonic analysis, and additionally about the propagation of light, do you already know?
 
Fourier Transform Equations:
F(νxy)=∫−∞f(x,y)exp(−i2π(νxx+νyy))dxdy

Transfer Function of Free Space (Fraunhofer Approximation):
g(x,y)=h0exp(iπ(x2+y2)/λd)F(x/λd,y/λd)
3. The Attempt at a Solution
f(x,y)=1
F(νxy)=∫−∞exp(−i2π(νxx+νyy)dxdy=δ(νx−0)δ(νy−0)
g(x,y)=∫−∞F(νxy)H(νxy)exp(+i2π(νxx+νyy))dνxdνy
g(x,y)=h0exp(iπ(x2+y2)/λd)δ(x/λd)δ(y/λd)
 
Last edited:
Good start. Unfortunately, the objects (the fields in the input plane) are defined for all spaces, i.e. both ##x## and ##y## run from ##-\infty## to ##\infty##. In this case, the condition necessary for Fraunhofer approximation cannot be satisfied. Instead start with the Helmholtz equation in free space
$$
\nabla^2 E(x,y,z) + k^2 E(x,y,z) = 0
$$
with ##k = 2\pi/\lambda##. Now apply Fourier transform on ##E(x,y,z)## to ##x## and ##y## only, so that in the end you will obtain a function ##\tilde{E}(k_x,k_y,z)##. Hint: use the rule ##\textrm{FT}[df(x)/dx] = ik_x\textrm{FT}[f(x)]##.
 
blue_leaf77 said:
Good start. Unfortunately, the objects (the fields in the input plane) are defined for all spaces, i.e. both ##x## and ##y## run from ##-\infty## to ##\infty##. In this case, the condition necessary for Fraunhofer approximation cannot be satisfied. Start with the Helmholtz equation in free space
$$
\nabla^2 E(x,y,z) + k^2 E(x,y,z) = 0
$$
with ##k = 2\pi/\lambda##. Now apply Fourier transform on ##E(x,y,z)## to ##x## and ##y## only, so that in the end you will obtain a function ##\tilde{E}(k_x,k_y,z)##. Hint: use the rule ##\textrm{FT}[df(x)/dx] = ik_x\textrm{FT}[f(x)]##.
yes i know that, but i want a solution of one of the cases above for better understand
 
Do you mean one of the (a), (b), and (c)? All of ##f(x,y)##'s given in each case extends from minus infinity to plus infinity in both ##x## and ##y## directions, Fraunhofer approximation is not applicable here as it requires that ##a^2/(z\lambda) << 0.1## with ##a## the size of the input field. For instance, look at your previous work, in the input plane you have plane wave ##f(x,y) = 1##, but at ##z=d## you got a point as your calculation resulted in an expression proportional a delta function. Does it look make sense to you?
 
Last edited:
The right way to solve it is by the use of angular spectrum method. In case you are not familiar with it, start with the suggestion I pointed out in post #4. I just want to stress again, that Fraunhofer approximation will not work here.
 
ok thank you for your help, I will try to solve it if possible
 
  • #10
blue_leaf77 said:
The right way to solve it is by the use of angular spectrum method. In case you are not familiar with it, start with the suggestion I pointed out in post #4. I just want to stress again, that Fraunhofer approximation will not work here.
ok thank you for your help, I will try to solve it if possible
 
  • #11
Amirjahi said:
ok thank you for your help, I will try to solve it if possible
hi

can i ask you the solution of FT[ f(x,y) = exp{(-jπ/λ)(x+y)}] ??
 

Similar threads

  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 2 ·
Replies
2
Views
2K
Replies
2
Views
4K
Replies
1
Views
2K
  • · Replies 61 ·
3
Replies
61
Views
12K
  • · Replies 23 ·
Replies
23
Views
3K
  • · Replies 102 ·
4
Replies
102
Views
11K
  • · Replies 33 ·
2
Replies
33
Views
9K
Replies
6
Views
2K
  • · Replies 42 ·
2
Replies
42
Views
11K