Abstract

In wireless communication, the travelling electromagnetic (EM) wave nature of the signal is not reported comprehensively in the literature. To date, the majority of reported work analyzed the nature of the signal propagating through the channel in time and frequency domains only, whereas the travelling EM wave which is actually propagating through the channel, has both time and spatial dependencies. Very little work has been reported in which this aspect of EM wave is being analyzed, but the time dependency of frequency and spatial dependency of propagation constant associated with the EM wave have not been included, which will affect the performance of the channel in practical conditions. In this paper, all these aspects of EM wave and channel are considered to define a new method (DDT-FRFT) to analyze the travelling EM wave for wireless communication.

1. Introduction

The fractional Fourier transform (FRFT) is an integral transform defined in time-frequency plane, that is, it can be considered as a generalization to the Fourier transform (FT) with an order (or power) parameter π‘Ž. The π‘Žth-order fractional Fourier transform operator is the π‘Žth power of the ordinary Fourier transform operator. The FRFT has already found many applications in the areas of signal processing and optics [1–9]. While originally formulated in an optical context, the same result equally applies to electromagnetic waves satisfying the linear wave equation. One of the central results of diffraction theory is that the far-field diffraction pattern is the Fourier transform of the diffracting object. It is possible to generalize this result by showing that the field patterns at closer distances are the fractional Fourier transforms of the diffracting object [10–12]. As the wave field propagates, its distribution evolves through fractional transforms of increasing orders. So the application of FRFT in wave propagation needs to be investigated more.

Later on, the travelling wave has been analyzed with the help of dual-domain transform (DDT-FT) [13], by transforming the travelling wave expression twice, first from time domain to frequency domain and second in spatial coordinate domain to propagation constant domain, by considering transform method as Fourier transform in both cases. The results of this transform are applied to estimate the transfer function of a wireless communication channel with multiple phase shifts, and, has been established that the channel acts as a band pass filter in both the frequency domain and the propagation constant domain.

But the time dependency of frequency components and spatial variable dependency of propagation constant have not been considered in [13], namely, the case when the electromagnetic wave propagating in a medium has different permittivity for different portions of the medium (e.g., in satellite communication, the wave has to travel from earth station to satellite through different layers of earth’s atmosphere and space) and also suffers Doppler shift in frequency due to motion of transmitter and/or receiver. This enforce a constraint on the analysis of electromagnetic wave that analysis should be carry out in time-frequency plane and spatial coordinate-propagation constant plane for time and spatial dependencies respectively, and the solution is FRFT. Hence, one another type of dual-domain transform is needed to be defined, which is based on the FRFT.

Some other works were reported in literature on synthetic aperture radar which has applied FRFT for estimating some of the parameters associated with the received signal. In Miaohong et al.’s work [14], the mechanism of radar target micro-Doppler effect has been applied to analyze the passive positioning system, and there it is demonstrated that the micro-Doppler frequency shifts on the received signal are relevant to both terminal motion and multipath propagation characteristics. And micro-Doppler parameters like terminal velocity and acceleration were estimated by FRFT. Yi et al. [15] gave a method for imaging and locating multiple ground moving targets which is based on the first-order Keystone transform and fractional Fourier transform. With this method, the focused image of moving targets can be obtained, and estimates of moving parameters such as radial velocity, azimuth speed, and radial acceleration can be determined. Similarly, Gao and Su [16] analyzed inverse synthetic aperture radar (ISAR) imaging by using FRFT. Here, the angular velocity can be calculated by knowing the chirp rate of chirp signal. In all these works, the effect of motion of transmitter/receiver on frequency is considered only while the EM wave nature of transmitted signal is not considered, and also the effect of the variation in the permittivity and permeability values for different portion of the overall media is neglected. However, in this paper, a new dual-domain transform based on the FRFT (DDT-FRFT) has been proposed, and with the help of some simulation results, the utility of this transform method in the above- mentioned case has been shown.

The paper is organized as follows. In Section 2, a brief introduction of FRFT is given along with the dual-domain transform defined for FT. The derivation of the proposed dual-domain transform based on FRFT along with its application in channel modeling is included in Section 3. It has been shown in Section 4 that the simulation results are in conformity with the proposed theory. The conclusive remarks are made in Section 5.

2. FRFT and DDT-FT

For an electromagnetic wave, the electric field intensity expression for transverse electromagnetic (TEM) mode, which is the function of both time and spatial domain variable, is given as𝑓(π‘₯,𝑑)=cos(π‘˜π‘₯βˆ’πœ”π‘‘).(2.1)

The dual-domain transform based on Fourier transform maps the signal defined in time and spatial β€œπ‘₯” domain to frequency and propagation constant β€œπ‘˜β€ domain, respectively. DDT-FT of the function 𝑓(π‘₯,𝑑) is defined as [13]𝐹(π‘˜,πœ”)=βˆžβˆ’βˆžπ‘“(π‘₯,𝑑)𝑒𝑗(π‘˜π‘₯βˆ’πœ”π‘‘)𝑑π‘₯𝑑𝑑.(2.2)

In DDT-FT, propagation constant should be independent of the spatial variable, and frequency should be independent of time variable. This constraint is violated when trans receiver is moving and channel changes its characteristics with distance. In this scenario, FRFT can play an important role. FRFT was introduced wayback in 1920s but remained largely unknown until the work of Namias in 1980 [17]. The FRFT can be viewed as the chirp-basis expansion directly from its definition, but essentially it can be interpreted as a rotation in the time-frequency plane, that is, the unified time-frequency transform, that is, FRFT is a generalization of Fourier transform (FT). The FRFT, represented by β€œβ„‘β€, of a signal π‘₯(𝑑) with angle parameter β€œπ›Όβ€, represented by 𝑋𝛼(𝑒), is defined for the entire time-frequency plane as [7]β„‘[]π‘₯(𝑑)=π‘‹π›Όξ€œ(𝑒)=βˆžβˆ’βˆžπ‘₯(𝑑)𝐾𝛼(𝑑,𝑒)𝑑𝑑,(2.3) where 𝐾𝛼(𝑑,𝑒) is the kernel of FRFT and is given asπΎπ›ΌβŽ§βŽͺβŽͺ⎨βŽͺβŽͺβŽ©ξ‚™(𝑑,𝑒)=1βˆ’π‘—cot𝛼𝑒2πœ‹π‘—/2{(𝑑2+𝑒2)cotπ›Όβˆ’2𝑑𝑒csc𝛼},ifπ›Όβ‰ π‘›πœ‹,𝛿(π‘‘βˆ’π‘’),if𝛼=2π‘›πœ‹,𝛿(𝑑+𝑒),if𝛼=(2𝑛+1)πœ‹.(2.4)

Many properties of this transform have been documented in the recent past including multiplication, mixed product, differentiation, integration, shift, similarity, and modulation in [7, 8], the shifting property will be utilized in formulating the channel response by DDT-FRFT. The shifting property of FRFT is given as [9]β„‘[]π‘₯(π‘‘βˆ’πœ)=𝑋𝛼(π‘’βˆ’πœcos𝛼)𝑒𝑗(𝜏2/2sin𝛼cosπ›Όβˆ’π‘’πœsin𝛼).(2.5)

3. DDT-FRFT and Its Application in Channel Modeling

The dual-domain transform, represented as β€œβ„‘β€, based on FRFT, is a technique which maps the signal defined in time and spatial β€œπ‘₯” domain to anywhere in time (spatial coordinate)-frequency (propagation constant) plane. The proposed DDT-FRFT of the signal 𝑓(π‘₯,𝑑) is defined as β„‘[]𝑓(π‘₯,𝑑)=𝐹𝛼(𝛽,𝑒)=βˆžβˆ’βˆžπ‘“(π‘₯,𝑑)𝐾𝛼(𝑑,𝑒)𝐾𝛼(π‘₯,𝛽)𝑑π‘₯𝑑𝑑.(3.1)

Based on the definition of DDT-FRFT, the shifting property of this transform is defined asβ„‘[]𝑓(π‘₯βˆ’π›Ώ,π‘‘βˆ’πœ)=𝑒𝑗sin𝛼((𝜏2/2)cosπ›Όβˆ’π‘’πœ)𝑒𝑗sin𝛼((𝛿2/2)cosπ›Όβˆ’π›½π›Ώ)𝐹𝛼(π›½βˆ’π›Ώcos𝛼,π‘’βˆ’πœcos𝛼).(3.2)

This is modified asβ„‘[]𝑓(π‘₯βˆ’π›Ώ,π‘‘βˆ’πœ)=π‘’π‘—πœ™πœπ‘’π‘—πœ™π›ΏπΉπ›Ό(π›½βˆ’π›Ώcos𝛼,π‘’βˆ’πœcos𝛼),(3.3) whereπœ™πœξ‚΅πœ=sin𝛼22ξ‚Ά,πœ™cosπ›Όβˆ’π‘’πœπ›Ώξ‚΅π›Ώ=sin𝛼22ξ‚Ά.cosπ›Όβˆ’π›½π›Ώ(3.4) Assume thatπœ™π‘›=πœ™πœ+πœ™π›Ώ.(3.5)

The output of channel having β€œπ‘β€ multipath, each has an associated delay β€œπœβ€ and associated path β€œπ›Ώβ€, will be given as𝑂𝑃=π‘βˆ’1𝑛=0π΄π‘›π‘’π‘—πœ™π‘›πΉπ›Όξ€·π›½βˆ’π›Ώπ‘›cos𝛼,π‘’βˆ’πœπ‘›ξ€Έ,cos𝛼(3.6) where β€œπ΄β€ is the amplitude of phase-shifted component. From this expression, it is viewable that the channel imposes the following parameters on the signal propagating through it:(i)introduces a gain β€œπ΄β€,(ii)a time delay 𝜏 and path delay 𝛿 when transformation is not carried,(iii)a phase shift β€œπœ™π‘›β€ in FRFT domain, which becomes equal to the phase shift in FT domain when angle parameter 𝛼=πœ‹/2,(iv)a delay (𝛿𝑛cos𝛼,πœπ‘›cos𝛼) in FRFT domain, which becomes equal to zero at 𝛼=πœ‹/2.

Since FRFT converts into FT at angle parameter 𝛼=πœ‹/2, so the output of channel is derived for this value of 𝛼 and is given as𝑂𝑃=π‘βˆ’1𝑛=0π΄π‘›π‘’π‘—πœ™π‘›,𝛼=πœ‹/2𝐹(π‘˜,πœ”).(3.7)

This expression resembles the output expression given for DDT-FT in [7]. The effects introduced by the channel to the propagating signal is actually:

(i) introduces a gain term β€œπ΄β€,

(ii) introduces a phase shift term in Fourier domain.

Here, the output signal which encounters multiple phase shifts in a wireless communication channel can be represented by𝐴0π‘’π‘—πœ™0𝐹(π‘˜,πœ”)+𝐴1π‘’π‘—πœ™1𝐹(π‘˜,πœ”)+𝐴2π‘’π‘—πœ™2𝐹(π‘˜,πœ”)+β‹―,(3.8) where πœ™π‘› is the various phase shifts in FT domain. So the transfer function (TF) of the channel will be calculated asTF=Channelπ‘œ/𝑝=𝐴Channel𝑖/𝑝0π‘’π‘—πœ™0𝐹(π‘˜,πœ”)+𝐴1π‘’π‘—πœ™1𝐹(π‘˜,πœ”)+β‹―.𝐹(π‘˜,πœ”)(3.9)

By expanding each exponential term as a power series in (3.9) and rearranging the expression𝐴TF=0+𝐴1𝐴+β‹―+𝑗0πœ™0+𝐴1πœ™1ξ€Έ+𝑗+β‹―22𝐴0πœ™02+𝐴1πœ™12ξ€Έ+β‹―+β‹―,(3.10)

and comparing it with the Fourier transform of a rectangular pulse in time domain [9], it can be easily predicted that the TF of channel is the transform of pulse-like characteristic in time domain and pulse-like characteristic in spatial domain.

As the FRFT of a rectangular function is calculated from 𝛼=0 to 𝛼=πœ‹/2, the rectangular function changes its shape from rectangular to sync function. Similarly, the channel response to a signal is analyzed in FRFT domain from 𝛼=0 to 𝛼=πœ‹/2, and it is clear from the discussion above that at 𝛼=0, channel imposes time and path delay to the signal, and at 𝛼=πœ‹/2, channel imposes phase shift term in frequency and propagation constant domain, that is, filtering action will be performed in these domains, and for 0<𝛼<πœ‹/2, channel will introduce not only phase term but also a delay term in that domain.

4. Utility of DDT-FRFT by Simulation

For establishing the utility of DDT-FRFT in the environment where frequency components of electromagnetic wave become time variant, and propagation constant associated with wave depends on spatial coordinate, a case has been studied with the help of simulation. Here, for simplicity, only the linear variation of propagation constant with spatial variable is assumed. In this example, the propagation constant β€œπ‘˜β€ and instantaneous frequency β€œπœ”β€, as shown in (2.1), having respective dependencies are given asξ€·π‘“πœ”=2πœ‹0ξ€Έ,ξ€·πœ‡+π‘Žπ‘‘π‘˜=2πœ‹0ξ€Έ.+𝑏π‘₯(4.1) And for these dependencies of β€œπœ”β€ and β€œπ‘˜β€ on time and spatial variable, the travelling electromagnetic wave has the expressionξ€·ξ€·πœ‡π‘“(π‘₯,𝑑)=cos2πœ‹0𝑓+𝑏π‘₯π‘₯βˆ’2πœ‹0𝑑.+π‘Žπ‘‘(4.2)

The values of𝑓0, π‘Ž, πœ‡0, and 𝑏 are chosen as 2, 10, 1.5, and 20, respectively. And the DDT-FRFT has been observed for a duration double of the time period of EM wave in time domain and double of the spatial wavelength in the spatial domain. First, the DDT-FRFT of travelling wave expression is evaluated for 𝛼=πœ‹/2, as it is the same as DDT-FT (because DDT-FRFT is converted to DDT-FT at 𝛼=πœ‹/2), as shown in Figure 1. Next, the DDT-FRFT is simulated for 𝛼=0.2, as shown in Figure 2. The results are shown in Figures 1 and 2

5. Conclusion

In this paper, a new dual-domain transform based on fractional Fourier transform (DDT-FRFT) has been proposed. The proposed technique has the application where EM wave encounters a Doppler shift due to the motion of transmitter and/or receiver along with having different phase velocity for the different portion of the medium in which the wave is propagating. Simulation results shown above are in the conformity of this statement, because at angle parameter β€œπ›Όβ€ associated with FRFT having value of β€œπœ‹/2”, that is, the case of Fourier transform, the transformed quantity does not have a clear peak (Figure 1), while for 𝛼=0.2, a clear peak is observable in the magnitude response of DDT-FRFT of travelling wave (Figure 2). It shows that in noisy conditions, due to a clear peak, the analysis of travelling wave with DDT-FRFT corresponding to 𝛼=0.2 gives better results than DDT-FT.

Also the application of DDT-FRFT is discussed for the channel modeling. And it has been established that a multipath wireless channel exhibits pulse-like characteristic not only in time domain but also in spatial domain, that is, for 𝛼=0. Similarly, the channel acts as a band pass filter in both the frequency domain and propagation constant domain, that is, for 𝛼=πœ‹/2. And for the values in between these two limits, the channel introduces a delay in that domain along with the filtering action.