Next Article in Journal
Hotel Comment Emotion Classification Based on the MF-DFA and Partial Differential Equation Classifier
Next Article in Special Issue
Parameter Estimation of LFM Signals Based on FOTD-CFRFT under Impulsive Noise
Previous Article in Journal
Continuous Dependence on the Initial Functions and Stability Properties in Hyers–Ulam–Rassias Sense for Neutral Fractional Systems with Distributed Delays
Previous Article in Special Issue
Sliding-Window TD-FrFT Algorithm for High-Precision Ranging of LFM Signals in the Presence of Impulse Noise
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Circuit of Quantum Fractional Fourier Transform

1
Information Science Teaching and Research Section, Northeastern University at Qinhuangdao, Qinhuangdao 066004, China
2
College of Marxism, Northeastern University at Qinhuangdao, Qinhuangdao 066004, China
*
Author to whom correspondence should be addressed.
Fractal Fract. 2023, 7(10), 743; https://doi.org/10.3390/fractalfract7100743
Submission received: 24 August 2023 / Revised: 28 September 2023 / Accepted: 7 October 2023 / Published: 8 October 2023
(This article belongs to the Special Issue Recent Advances in Fractional Fourier Transforms and Applications)

Abstract

:
In this paper, we first use the quantum Fourier transform (QFT) and quantum phase estimation (QPE) to realize the quantum fractional Fourier transform (QFrFT). As diverse definitions of the discrete fractional Fourier transform (DFrFT) exist, the relationship between the QFrFT and a classical algorithm is then established; that is, we determine the classical algorithm corresponding to the QFrFT. Second, we observe that many definitions of the multi-fractional Fourier transform (mFrFT) are flawed: when we attempt to propose a design scheme for the quantum mFrFT, we find that there are many invalid weighting terms in the definition of the mFrFT. This flaw may have very significant impacts on relevant algorithms for signal processing and image encryption. Finally, we analyze the circuit of the QFrFT and the reasons for the observed defects.

1. Introduction

In recent years, research into and the application of quantum information technology—represented by quantum computing, quantum communication, and quantum measurement—have accelerated globally, with many countries increasing investments and broadening project layouts. Quantum computing is expected to fundamentally improve the speed of information processing, quantum communication will greatly improve communication security, and quantum precision measurement and sensing technology will have extensive applications in the future digital age and the era of the Internet of Things. Quantum computing is a novel computation mode that follows the laws of quantum mechanics to regulate the calculation of quantum information units, utilizing the characteristics of quantum mechanics [1] (e.g., superposition and entangled states) for the storage, processing, and transmission of information. Since Feynman proposed the concept of quantum computing [2], researchers have shown great interest in the high parallelism and speed provided by such an approach. The discovery of Shor’s efficient decomposition algorithm [3] aroused interest in the quantum implementation of modular arithmetic operations, which are the basis of quantum decomposition circuits. The core of Shor’s algorithm is the quantum Fourier transform (QFT) [4]. In the field of quantum computing, involving topics such as quantum information processing [5], quantum machine learning [6], quantum neural networks [7], and so on, a series of significant theoretical improvements and encouraging experimental results have recently been achieved, in which the QFT plays an important role. Quantum image processing plays an important role in satisfying the demands of machine vision, and quantum image representation is the basis of quantum image processing. The representation of quantum images has been developed from the initial Qubit Lattice [8], Real Ket [9], and Entangled image [10] to the common representation forms FRQI [11] and NEQR [12]. Based on these characterizations, various methods for quantum image processing have also been proposed [13], such as local feature point extraction of quantum images [14], weighted filtering of quantum images in the spatial domain [15], and restoration of quantum image noise removal [16]. Among them, the QFT plays an important role in frequency domain filtering [17,18], and the quantum arithmetic operations associated with the QFT are equally important to consider [19]. It is clear that the QFT plays an important role in quantum algorithm design and information processing.
As an extended definition of the Fourier transform (FT), the fractional Fourier transform (FrFT) has become one of the most actively researched subjects in signal processing and has been widely used in optics, image and signal processing, and communication [20,21,22,23,24]. However, at present, there are few studies on the quantum fractional Fourier transform (QFrFT) [25,26,27,28]. The only existing studies have also been conducted from the perspective of quantum mechanics, in order to study the theory underlying the QFrFT, and have not presented circuits for realization of the QFrFT. In this paper, the QFrFT is designed by means of the QFT and quantum phase estimation (QPE) [1]. The idea is to retain the eigenvector of the QFT, realize the fractional power of the eigenvalue by introducing two phase gates, and finally obtain the QFrFT by combining the eigenvalue of the fractional power with the eigenvector.
The remainder of this paper is organized as follows. The QFrFT is proposed in Section 2. Section 3 explains the correspondence between the QFrFT and a classical algorithm. The flaws of the mFrFT are detailed in Section 4. In Section 5, circuits for the QFrFT and mFrFT are discussed. Finally, the conclusions are presented in Section 6.

2. Quantum Fractional Fourier Transform

At present, some related studies on the QFrFT have been proposed [25,26,27,28], but these studies have not presented specific quantum circuits. Parasa et al. have even stated that the QFrFT is impossible to achieve [25]. In this section, we use the QFT and QPE to design the QFrFT.
For a unitary matrix U , assuming that it has an eigenvector | u and the corresponding eigenvalue e 2 π i φ , then U | u = e 2 π i φ | u is satisfied. Therefore, we can calculate φ through the QPE. The circuit is shown in Figure 1. It is not difficult to find that the QFT is the key of phase estimation, and phase estimation is the key of many quantum algorithms.
We know the relationship between the eigenvalues and eigenvectors of the discrete Fourier transform (DFT), as shown in Equation (1):
F | u = D | u ,
where F denotes the Fourier transform, D is the eigenvalue, and | u is the eigenvector. The discrete fractional Fourier transform (DFrFT) has the same eigenvector | u as the DFT, and the correspondence between the eigenvalue and eigenvector is shown in Equation (2):
F α | u = D α | u ,
where F α denotes the FrFT, and D α is the eigenvalue. Thus, we can determine the eigenvalue of the QFrFT by means of the QPE, as depicted in Figure 2.
The process is as follows:
| ψ 1 = 1 2 ( | 0 + | 1 ) 1 2 ( | 0 + | 1 ) | u = 1 2 ( | 00 + | 01 + | 10 + | 11 ) | u .
Thus, we can obtain
| ψ 2 = 1 2 ( | 00 I | u + | 01 F | u + | 10 F 2 | u + | 11 F 3 | u ) .
Here, | u is the eigenvector of the FT. From Equation (1), we can obtain
| ψ 2 = 1 2 ( | 00 D 0 | u + | 01 D 1 | u + | 10 D 2 | u + | 11 D 3 | u ) .
The eigenvalue D can be expressed as:
D = ( λ 0 0 λ 1 0 λ n 1 ) ,
where λ j { 1 , i , 1 , i } ; j = 0 , 1 , , n 1 . The 2-qubit inverse quantum Fourier transform (IQFT) is used in the control register, and its matrix can be expressed as:
I Q F T = 1 4 ( 1 1 1 1 1 i 1 i 1 1 1 1 1 i 1 i ) .
After 2-qubit IQFT, we obtain
| ψ 3 = 1 4 ( Q | 00 + W | 01 + E | 10 + R | 11 ) | u
where
{ Q = D 0 + D 1 + D 2 + D 3 W = D 0 + i D 1 D 2 i D 3 E = D 0 D 1 + D 2 D 3 R = D 0 i D 1 D 2 + i D 3 .
Moreover, from Equations (6) and (9), we have
Q = D 0 + D 1 + D 2 + D 3 = ( λ 0 0 + λ 0 1 + λ 0 2 + λ 0 3 λ 1 0 + λ 1 1 + λ 1 2 + λ 1 3 λ n 1 0 + λ n 1 1 + λ n 1 2 + λ n 1 3 ) ,
W = D 0 + i D 1 D 2 i D 3 = ( λ 0 0 + i λ 0 1 λ 0 2 i λ 0 3 λ 1 0 + i λ 1 1 λ 1 2 i λ 1 3 λ n 1 0 + i λ n 1 1 λ n 1 2 i λ n 1 3 ) ,
E = D 0 D 1 + D 2 D 3 = ( λ 0 0 λ 0 1 + λ 0 2 λ 0 3 λ 1 0 λ 1 1 + λ 1 2 λ 1 3 λ n 1 0 λ n 1 1 + λ n 1 2 λ n 1 3 ) ,
and
R = D 0 i D 1 D 2 + i D 3 = ( λ 0 0 i λ 0 1 λ 0 2 + i λ 0 3 λ 1 0 i λ 1 1 λ 1 2 + i λ 1 3 λ n 1 0 i λ n 1 1 λ n 1 2 + i λ n 1 3 ) .
Here, the eigenvalues are λ j { 1 , i , 1 , i } ; j = 0 , 1 , , n 1 . Thus, Equations (10)–(13) can be written as:
Q = { λ j 0 + λ j 1 + λ j 2 + λ j 3 = 4 i f λ j = 1 λ j 0 + λ j 1 + λ j 2 + λ j 3 = 0 i f λ j = i , 1 , i ,
W = { λ j 0 + i λ j 1 λ j 2 i λ j 3 = 4 i f λ j = i λ j 0 + i λ j 1 λ j 2 i λ j 3 = 0 i f λ j = 1 , i , 1 ,
E = { λ j 0 λ j 1 + λ j 2 λ j 3 = 4 i f λ j = 1 λ j 0 λ j 1 + λ j 2 λ j 3 = 0 i f λ j = 1 , i , i ,
and
R = { λ j 0 i λ j 1 λ j 2 + i λ j 3 = 4 i f λ j = i λ j 0 i λ j 1 λ j 2 + i λ j 3 = 0 i f λ j = 1 , 1 , i .
Then, the eigenvalues of the FT are preserved in the quantum state, as shown in Equation (18):
{ λ j = 1 | 00 λ j = i | 01 λ j = 1 | 10 λ j = i | 11 .
We know that the eigenvector | u can be written in standard orthogonal basis form, as shown in Equation (19):
| u = j b j | u j ,
where b j is the projection length. Thus, Equation (8) can also be written as:
| ψ 3 = j | ϕ j b j | u j ,
where | ϕ j { | 00 , | 01 , | 10 , | 11 } ; j = 0 , 1 , , n 1 . Next, the phase gates R Z ( π α ) and R Z ( π α / 2 ) are expressed as:
R Z ( π α ) = ( 1 0 0 e π i α )
and
R Z ( π α / 2 ) = ( 1 0 0 e π i α / 2 ) .
We can obtain the quantum state | ψ 4 by
| ψ 4 = j e π i α ϕ j / 2 | ϕ j b j | u j .
The 2-qubit QFT can be expressed as:
Q F T = 1 4 ( 1 1 1 1 1 i 1 i 1 1 1 1 1 i 1 i ) .
After the QFT of | ϕ j , the eigenvalues are restored:
| ψ 5 = 1 2 ( j e π i α ϕ j / 2 λ j 0 b j | u j | 00 + j e π i α ϕ j / 2 λ j 1 b j | u j | 01 + j e π i α ϕ j / 2 λ j 2 b j | u j | 10 + j e π i α ϕ j / 2 λ j 3 b j | u j | 11 ) = 1 2 ( j e π i α ϕ j / 2 b j | u j | 00 + j e π i α ϕ j / 2 F b j | u j | 01 + j e π i α ϕ j / 2 F 2 b j | u j | 10 + j e π i α ϕ j / 2 F 3 b j | u j | 11 ) ,
where ϕ j { 0 , 1 , 2 , 3 } . Then,
| ψ 6 = 1 2 ( j e π i α ϕ j / 2 λ j 0 b j | u j | 00 + j e π i α ϕ j / 2 I F F b j | u j | 01 + j e π i α ϕ j / 2 I F 2 F 2 b j | u j | 10 + j e π i α ϕ j / 2 I F 3 F 3 b j | u j | 11 ) = 1 2 ( j e π i α ϕ j / 2 b j | u j | 00 + j e π i α ϕ j / 2 b j | u j | 01 + j e π i α ϕ j / 2 b j | u j | 10 + j e π i α ϕ j / 2 b j | u j | 11 ) = j e π i α ϕ j / 2 b j | u j 1 2 ( | 0 + | 1 ) 1 2 ( | 0 + | 1 ) ,
where IF denotes the inverse Fourier transform. Finally, through the H gate, we obtain | ψ 7 :
| ψ 7 = | 0 | 0 j e π i α ϕ j / 2 b j | u j .
From Equations (1) and (2), we know that
j e π i ϕ j / 2 b j | u j = F | u ,
j e π i α ϕ j / 2 b j | u j = F α | u .
Then, Equation (27) can be expressed as:
| ψ 7 = | 0 | 0 F α | u .
Therefore, the QFrFT is obtained.
In Figure 2, the upper half is the control register, and the lower half is the target register. As the eigenvalues of QFT are only 4, the scale of the QFT and IQFT in the control register is 2-qubit, while the scale of the QFT and IQFT in the target register can be arbitrary. The QFrFT uses two additional auxiliary qubits (i.e., control registers) in space. In time, when the target register is large, only the quantum gate complexity is considered. The complexity of the QFT is O ( n 2 ) , where n is the number of qubits of the target register [1].

3. Classical Algorithm Corresponding to the QFrFT

There exist diverse definitions for the DFrFT [29]. This section analyses the correspondence between the proposed QFrFT and classical algorithms. The design idea underlying the proposed QFrFT is the fractional power of the eigenvalues. We know that the DFT can be expressed as eigenvalues and eigenvectors, as shown in Equation (31):
F = V D V T ,
where D is the eigenvalue matrix and V = [ v 0 | v 1 | v N 2 | v N 1 ] is the eigenvector. Thus, the DFrFT can be expressed as:
F α = V D α V T .
The eigenvalues of the FT are 1, i, −1, and −i. Therefore, Equations (31) and (32) can be further expressed as:
F = k E 1 ( 1 ) v k v k T + k E 2 ( i ) v k v k T + k E 3 ( 1 ) v k v k T + k E 4 ( i ) v k v k T
and
F α = k E 1 ( 1 ) α v k v k T + k E 2 ( i ) α v k v k T + k E 3 ( 1 ) α v k v k T + k E 4 ( i ) α v k v k T ,
where E 1 , E 2 , E 3 , and E 4 represent the set of subscripts corresponding to the respective eigenvalues and eigenvectors. For the QFrFT in Figure 2, the important calculation steps are
I Q F T ( I | 00 + F | 01 + F 2 | 10 + F 3 | 11 ) = I Q F T ( I F F 2 F 3 ) = 1 4 ( 1 1 1 1 1 i 1 i 1 1 1 1 1 i 1 i ) ( I F F 2 F 3 ) .
Let
{ Y 0 = I + F + F 2 + F 3 Y 1 = I + i F F 2 i F 3 Y 2 = I F + F 2 F 3 Y 3 = I i F F 2 + i F 3 .
At this point, we can obtain
{ Y 0 | u = ( D 0 + D 1 + D 2 + D 3 ) | u Y 1 | u = ( D 0 + i D 1 D 2 i D 3 ) | u Y 2 | u = ( D 0 D 1 + D 2 D 3 ) | u Y 3 | u = ( D 0 i D 1 D 2 + i D 3 ) | u .
Then,
{ Y 0 | u = Q | u Y 1 | u = W | u Y 2 | u = E | u Y 3 | u = R | u .
From the above results in Equations (14)–(17) for Q, W, E, and R, Equation (36) can be expressed as:
{ Y 0 = k E 1 v k v k T Y 1 = k E 2 v k v k T Y 2 = k E 3 v k v k T Y 3 = k E 4 v k v k T .
From the phase gates shown in Figure 2, we can obtain
( 1 0 0 e π i α ) ( 1 0 0 e π i α / 2 ) = ( 1 0 e π i α / 2 e 2 π i α / 2 0 e 3 π i α / 2 ) = ( ( 1 ) α 0 ( i ) α ( 1 ) α 0 ( i ) α ) .
Thus, Equation (34) can be obtained by combining Equations (39) and (40). Thus, the mathematical principle of the QFrFT is explained.
Considering the many definitions of the DFrFT [29], we found that Shih’s FrFT is consistent with the proposed QFrFT. In 1995, Shih proposed the FrFT using the weighted combination of four functions [30]. Its definition can be expressed as:
( F α f ) ( t ) = l = 0 3 A l α f l ( t ) ,
where f l ( t ) = F l [ f ( t ) ] ; l = 0 , 1 , 2 , 3 . The weighting coefficient A l α can be further written as:
A l α = cos ( ( α l ) π 4 ) cos ( 2 ( α l ) π 4 ) exp ( 3 ( α l ) i π 4 ) = 1 2 × [ exp ( ( α l ) π i 4 ) + exp ( ( α l ) π i 4 ) ] × 1 2 × [ exp ( 2 ( α l ) π i 4 ) + exp ( 2 ( α l ) π i 4 ) ] × exp ( 3 ( α l ) i π 4 ) = 1 4 ( 1 + exp ( 2 ( α l ) π i 4 ) + exp ( 4 ( α l ) π i 4 ) + exp ( 6 ( α l ) π i 4 ) ) = 1 4 k = 0 3 exp ( 2 π i 4 ( α l ) k ) = 1 4 k = 0 3 exp ( 2 π i α k 4 ) exp ( 2 π i k l 4 ) .
Therefore, the weighting coefficient A l α can be written in matrix form
( A 0 α A 1 α A 2 α A 3 α ) = 1 4 ( 1 1 1 1 1 i 1 i 1 1 1 1 1 i 1 i ) ( B 0 α B 1 α B 2 α B 3 α ) ,
where B k α = exp ( 2 π i k α 4 ) ; k = 0 , 1 , 2 , 3 . Thus, Shih’s FrFT can be expressed as:
F α [ f ( t ) ] = ( I , F , F 2 , F 3 ) ( A 0 α A 1 α A 2 α A 3 α ) f ( t ) = 1 4 ( I , F , F 2 , F 3 ) ( 1 1 1 1 1 i 1 i 1 1 1 1 1 i 1 i ) ( B 0 α B 1 α B 2 α B 3 α ) f ( t ) .
From Equation (36), Shih’s FrFT can be further expressed as:
F α [ f ( t ) ] = 1 4 ( Y 0 , Y 1 , Y 2 , Y 3 ) ( B 0 α B 1 α B 2 α B 3 α ) f ( t ) ,
where B k α = exp ( 2 π i k α 4 ) ; k = 0 , 1 , 2 , 3 and
( B 0 α B 1 α B 2 α B 3 α ) = ( ( 1 ) α ( i ) α ( 1 ) α ( i ) α ) .
Thus, we obtain Equation (34). Therefore, Shih’s FrFT coincides with the proposed QFrFT.

4. Implementation of Quantum Multi-Fractional Fourier Transform

Shih’s FrFT has many extended definitions [31,32,33,34], which are collectively referred to as the multi-fractional Fourier transform (mFrFT). Next, we focus on the design of quantum algorithms for the mFrFT. The implementation of its quantum algorithm will be of great help in the fields of quantum computing and quantum information security. In 2000, Zhu et al. proposed the definition of the mFrFT and applied it to image encryption [31]. This definition can be expressed as:
F M α [ f ( t ) ] = l = 0 M 1 A l α f l ( t ) ,
where the basic functions are f l ( t ) = F 4 l / M [ f ( t ) ] ; l = 0 , 1 , , M 1 . The weighting coefficient A l α can be expressed as:
A l α = 1 M k = 0 M 1 exp [ 2 π i k ( α l ) M ] = 1 M k = 0 M 1 exp ( 2 π i k l M ) exp ( 2 π i α k M ) = I D F T [ e ( 2 π i α k / M ) ] k = 0 , 1 , 2 , , M 1 .
where IDFT is the inverse discrete Fourier transform. Furthermore, the weighting coefficients A l α can be expressed in matrix form as:
( A 0 α A 1 α A M 1 α ) = 1 M ( w 0 × 0 w 0 × 1 w 0 × ( M 1 ) w 1 × 0 w 1 × 1 w 1 × ( M 1 ) w ( M 1 ) × 0 w ( M 1 ) × 1 w ( M 1 ) × ( M 1 ) ) ( B 0 α B 1 α B M 1 α ) ,
where w = exp ( 2 π i / M ) and B k α = exp ( 2 π i k α M ) , k = 0 , 1 , , M 1 . Therefore, Equation (47) can be expressed as:
F M α [ f ( t ) ] = A 0 α f 0 ( t ) + A 1 α f 1 ( t ) + + A M 1 α f M 1 ( t ) = A 0 α F 0 M [ f ( t ) ] + A 1 α F 4 M [ f ( t ) ] + + A M 1 α F 4 ( M 1 ) M [ f ( t ) ] = ( I , F 4 M , , F 4 ( M 1 ) M ) ( A 0 α A 1 α A M 1 α ) f ( t ) ,
Further, we can obtain
F M α [ f ( t ) ] = 1 M ( I , F 4 M , , F 4 ( M 1 ) M ) ( w 0 × 0 w 0 × 1 w 0 × ( M 1 ) w 1 × 0 w 1 × 1 w 1 × ( M 1 ) w ( M 1 ) × 0 w ( M 1 ) × 1 w ( M 1 ) × ( M 1 ) ) ( B 0 α B 1 α B M 1 α ) f ( t ) .
Here, let
{ Y 0 = w 0 × 0 I + w 1 × 0 F 4 M + + w ( M 1 ) × 0 F 4 ( M 1 ) M Y 1 = w 0 × 1 I + w 1 × 1 F 4 M + + w ( M 1 ) × 1 F 4 ( M 1 ) M Y 2 = w 0 × 2 I + w 1 × 2 F 4 M + + w ( M 1 ) × 2 F 4 ( M 1 ) M   Y M 1 = w 0 × ( M 1 ) I + w 1 × ( M 1 ) F 4 M + + w ( M 1 ) × ( M 1 ) F 4 ( M 1 ) M .
In this way, the mFrFT can also be expressed as Equation (53):
F M α [ f ( t ) ] = 1 M ( Y 0 , Y 1 , , Y M 1 ) ( B 0 α B 1 α B M 1 α ) f ( t ) = 1 M k = 0 M 1 Y k B k α f ( t ) ,
where B k α = exp ( 2 π i k α M ) ; k = 0 , 1 , , M 1 . Next, analyzing the results of Equation (52), from Equation (37), we obtain
{ Y 0 | u = ( w 0 × 0 D 0 + w 1 × 0 D 4 M + + w ( M 1 ) × 0 D 4 ( M 1 ) M ) | u Y 1 | u = ( w 0 × 1 D 0 + w 1 × 1 D 4 M + + w ( M 1 ) × 1 D 4 ( M 1 ) M ) | u Y 2 | u = ( w 0 × 2 D 0 + w 1 × 2 D 4 M + + w ( M 1 ) × 2 D 4 ( M 1 ) M ) | u Y M 1 | u = ( w 0 × ( M 1 ) D 0 + w 1 × ( M 1 ) D 4 M + + w ( M 1 ) × ( M 1 ) D 4 ( M 1 ) M ) | u .
Let
{ Q 0 = w 0 × 0 D 0 + w 1 × 0 D 4 M + + w ( M 1 ) × 0 D 4 ( M 1 ) M Q 1 = w 0 × 1 D 0 + w 1 × 1 D 4 M + + w ( M 1 ) × 1 D 4 ( M 1 ) M Q 2 = w 0 × 2 D 0 + w 1 × 2 D 4 M + + w ( M 1 ) × 2 D 4 ( M 1 ) M Q M 1 = w 0 × ( M 1 ) D 0 + w 1 × ( M 1 ) D 4 M + + w ( M 1 ) × ( M 1 ) D 4 ( M 1 ) M .
Then,
Q 0 = w 0 × 0 λ j 0 + w 1 × 0 λ j 4 M + + w ( M 1 ) × 0 λ j 4 ( M 1 ) M = λ j 0 + λ j 4 M + + λ j 4 ( M 1 ) M ,
where λ j ( 1 , i , 1 , i ) . For proof, we set λ j ( e π i 0 / 2 , e π i 1 / 2 , e π i 2 / 2 , e π i 3 / 2 ) .
When λ j = 1 ,
λ j 0 + λ j 4 M + + λ j 4 ( M 1 ) M = M .
When λ j = i ,
λ j 0 + λ j 4 M + + λ j 4 ( M 1 ) M = 1 + e 2 π i 1 / M + e 2 π i 2 / M + + e 2 π i ( M 1 ) / M = 1 ( e 2 π i / M ) M 1 e 2 π i / M = 1 1 1 e 2 π i / M = 0 .
When λ j = 1 ,
λ j 0 + λ j 4 M + + λ j 4 ( M 1 ) M = 1 + e 4 π i 1 / M + e 4 π i 2 / M + + e 4 π i ( M 1 ) / M = 1 ( e 4 π i / M ) M 1 e 2 π i / M = 1 1 1 e 2 π i / M = 0 .
When λ j = 1 ,
λ j 0 + λ j 4 M + + λ j 4 ( M 1 ) M = 1 + e 6 π i 1 / M + e 6 π i 2 / M + + e 6 π i ( M 1 ) / M = 1 ( e 6 π i / M ) M 1 e 2 π i / M = 1 1 1 e 2 π i / M = 0 .
So, we can obtain
Q 0 = { M , λ j = 1 0 , λ j = i , 1 , i .
Next, we calculate Q 1 as
Q 1 = w 0 × 1 D 0 + w 1 × 1 D 4 M + + w ( M 1 ) × 1 D 4 ( M 1 ) M = k = 0 M 1 e 2 π i k / M λ j 4 k / M ,
when λ j ( e π i 0 / 2 , e π i 1 / 2 , e π i 2 / 2 , e π i 3 / 2 ) , we have Q 1 = 0 . Moreover, our calculation shows that, for Equation (55), Q 2 = Q 3 = = Q M 4 = 0 .
So, we calculate Q M 3 as
Q M 3 = w 0 × ( M 3 ) D 0 + w 1 × ( M 3 ) D 4 M + + w ( M 1 ) × ( M 3 ) D 4 ( M 1 ) M = k = 0 M 1 e 2 π i k ( M 3 ) / M λ j 4 k / M ,
where λ j ( e π i 0 / 2 , e π i 1 / 2 , e π i 2 / 2 , e π i 3 / 2 ) . Then, we obtain
Q 3 = k = 0 M 1 e 2 π i k ( M 3 ) / M λ j 4 k / M = k = 0 M 1 e ( 2 π i k M 6 π i k + 2 π i k h ) / M ,
where h = 0 , 1 , 2 , 3 . For Equation (64), when h = 3 , the result is non-zero and the eigenvalue is e π i 3 / 2 = i ,
Q M 3 = { M λ = i 0 λ = 1 , i , 1 .
For Q M 2 ,
Q M 2 = w 0 × ( M 2 ) D 0 + w 1 × ( M 2 ) D 4 M + + w ( M 1 ) × ( M 2 ) D 4 ( M 1 ) M = k = 0 M 1 e 2 π i k ( M 2 ) / M λ j 4 k / M = k = 0 M 1 e ( 2 π i k M 4 π i k + 2 π i k h ) / M .
For Equation (66), when h = 2 , the result is non-zero and the eigenvalue is e π i 2 / 2 = 1 ,
Q M 2 = { M λ = 1 0 λ = 1 , i , i .
For Q M 1 ,
Q M 1 = w 0 × ( M 1 ) D 0 + w 1 × ( M 1 ) D 4 M + + w ( M 1 ) × ( M 1 ) D 4 ( M 1 ) M = k = 0 M 1 e 2 π i k ( M 1 ) / M λ j 4 k / M = k = 0 M 1 e ( 2 π i k M 2 π i k + 2 π i k h ) / M .
For Equation (68), when h = 1 , the result is non-zero and the eigenvalue is e π i / 2 = i ,
Q M 1 = { M λ = i 0 λ = 1 , 1 , i .
Therefore, for Equation (52), only Y 0 , Y M 3 , Y M 2 , and Y M 1 are non-zero, while the other terms are zero. If the eigenvalue λ j ( e π i 0 / 2 , e π i 1 / 2 , e π i 2 / 2 , e π i 3 / 2 ) , in Equation (55), only Q 0 , Q 1 , Q 2 , and Q 3 are non-zero, while the other results are zero. For Equation (53), there are only four effective weighting terms for the mFrFT.
However, for Equation (53), when B k α takes different values, we obtain different definitions [35]. For example, B k α = exp ( 2 π i k α 4 ) ; k = 0 , 1 , 2 , 3 gives Shih’s FrFT [30]; B k α = exp ( 2 π i k α M ) ; k = 0 , 1 , , M 1 gives the mFrFT of Zhu et al. [31]; B k α = exp [ 2 π i ( m k M + 1 ) α ( n k M + k ) M ] ; k = 0 , 1 , , M 1 gives the multi-parameter fractional Fourier transform of Tao et al. [32]; B k α = exp [ 2 π i α ( r k M + l ) M ] ; k = 0 , 1 , , M 1 gives the modified multi-parameter fractional Fourier transform of Ran et al. [33]; and B k α ¯ = exp [ 2 π i α l ( r l M + l ) M ] ; k = 0 , 1 , , M 1 gives the vector power multi-parameter fractional Fourier transform of Ran et al. [34].
Obviously, all of these extended definitions also have only four effective weighting terms, which has a great impact on the application of these algorithms in signal processing and image encryption.

5. Discussion

5.1. Phase Gates

For the QFrFT shown in Figure 2, we select the phase gates R Z ( π α ) and R Z ( π α / 2 ) . If the phase gates R Z ( π α ) and R Z ( π α / 2 ) are selected, the QFrFT circuit is as shown in Figure 3.
From Figure 3, we can see that the phase gates have changed, and the order of the QFT and the IQFT regarding the 2-qubits in the control register has also changed. This is to ensure that the eigenvalues correspond to the eigenvectors. The correspondence between the eigenvalues of the control register is:
{ λ j = 1 | 00 λ j = i | 01 λ j = 1 | 10 λ j = i | 11 .
From the phase gates R Z ( π α ) and R Z ( π α / 2 ) , we obtain
( 1 0 0 e π i α ) ( 1 0 0 e π i α / 2 ) = ( 1 0 e π i α / 2 e 2 π i α / 2 0 e 3 π i α / 2 ) = ( ( 1 ) α 0 ( i ) α ( 1 ) α 0 ( i ) α ) .
This ensures the correspondence between eigenvalues and eigenvectors. Thus, we have
F = k E 1 ( 1 ) v k v k T + k E 2 ( i ) v k v k T + k E 3 ( 1 ) v k v k T + k E 4 ( i ) v k v k T
and
F α = k E 1 ( 1 ) α v k v k T + k E 2 ( i ) α v k v k T + k E 3 ( 1 ) α v k v k T + k E 4 ( i ) α v k v k T .
The QFrFT and Shih’s FrFT are the same. To ensure that the eigenvalues correspond to the eigenvectors, when the eigenvalues become fractional powers, they also correspond to the eigenvectors. Therefore, the key idea of the algorithm’s design is to retain the eigenvector of the FT, such that the eigenvalue takes a fractional power.
Considering Figure 2, if we only change the phase gates, the result is shown in Figure 4. The relevant correspondence between the eigenvalues and eigenvectors is shown in Equations (33) and (73). In this way, the obtained results do not correspond to the feature vector of the FT and, so, the design scheme shown in Figure 4 is not feasible.

5.2. Flaws of the mFrFT

From the design of the quantum mFrFT, we find that there are many invalid weighting terms in the definition of classical mFrFT. The reason for this problem is due to the eigenvalues of the algorithm. We know that the eigenvalues of the FT are 1, i, −1 and −i, and the time–frequency spatial distribution is shown in Figure 5. For Equation (54), the sum of eigenvalues can be expressed as a change in angle in Figure 5.
For angle θ , we set
w = e i θ ,
where θ = 2 π / M . We can obtain
1 + w + w 2 + + w M 1 = 0 .
Furthermore, we have
1 + w h + w 2 h + + w ( M 1 ) h = { M h = 0   0   h 0 .
Thus, Equation (64) can be expressed as
k = 0 M 1 e ( 2 π i k M 6 π i k + 2 π i k h ) / M = k = 0 M 1 e 2 π i k e 2 π i k M ( h 3 ) .
Here, e 2 π i k = 1 . Therefore, the result of Equation (77) is non-zero only when h = 3 . The other results are also confirmed.
Since 2000, Zhu’s mFrFT [31] and previous results [36], as well as some extended definitions, have been shown to have flaws [32,33,34], which we discussed in terms of the eigenvalues of the QFT. In the extended definition, the number of effective weighting terms depends on the period of the matrix; for example, the period of the FT matrix is 4, so its effective weighting is 4, while the period of the Hadamard transform is 2, so the effective weighting term is defined as 2 accordingly.

6. Conclusions

In this paper, we used the QFT and QPE to design the QFrFT. We retained the eigenvectors of the QFT and forced its eigenvalues to take fractional powers. As diverse definitions for the DFrFT exist, we determined the correspondence between the proposed QFrFT and classical algorithms, which lays the foundation for future application of the QFrFT to solve practical problems. Regarding the design of the quantum mFrFT, we observed a flaw in the definition of classical algorithms—namely, there are only four effective weighting terms in a series of mFrFT definitions—which has a significant impact on the applicability of the algorithms. Finally, we discussed the correspondence between eigenvalues and eigenvectors in the QFrFT circuit and analyzed the definition of the algorithm in depth. Furthermore, the defects of the mFrFT were mathematically explained.

Author Contributions

Methodology, T.Z.; validation, T.Z. and Y.C.; investigation, T.Z.; writing—original draft preparation, T.Z.; writing—review and editing, T.Z.; supervision, T.Z.; project administration, T.Z.; funding acquisition, T.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by the Fundamental Research Funds for the Central Universities (N2123016); and the Scientific Research Projects of Hebei colleges and universities (QN2020511).

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Nielsen, M.A.; Chuang, I.L. Quantum computation and quantum information. Phys. Today 2001, 54, 60. [Google Scholar]
  2. Feynman, R.P. Simulating physics with computers. Int. J. Theor. Phys. 1982, 21, 467–488. [Google Scholar] [CrossRef]
  3. Shor, P.W. Algorithms for Quantum Computation: Discrete Logarithms and Factoring. In Proceedings of the 35th Annual Symposium on Foundation of Computer Science, Washington, DC, USA, 20–22 November 1994; pp. 124–134. [Google Scholar]
  4. Coppersmith, D. An approximate Fourier transform useful in quantum factoring. arXiv 2002. [Google Scholar] [CrossRef]
  5. Flamini, F.; Spagnolo, N.; Sciarrino, F. Photonic quantum information processing: A review. Rep. Prog. Phys. 2018, 82, 016001. [Google Scholar] [PubMed]
  6. Biamonte, J.; Wittek, P.; Pancotti, N.; Rebentrost, P.; Wiebe, N.; Lloyd, S. Quantum machine learning. Nature 2017, 549, 195–202. [Google Scholar]
  7. Beer, K.; Bondarenko, D.; Farrelly, T.; Osborne, T.J.; Salzmann, R.; Scheiermann, D.; Wolf, R. Training deep quantum neural networks. Nat. Commun. 2020, 11, 808. [Google Scholar] [CrossRef]
  8. Venegas-Andraca, S.E.; Bose, S. Quantum computation and image processing: New trends in artificial intelligence. In Proceedings of the IJCAI, Acapulco, Mexico, 9–15 August 2003. [Google Scholar]
  9. Latorre, J.I. Image compression and entanglement. arXiv 2005. [Google Scholar] [CrossRef]
  10. Venegas-Andraca, S.E.; Ball, J. Processing images in entangled quantum systems. Quantum Inf. Process. 2010, 9, 1–11. [Google Scholar] [CrossRef]
  11. Le, P.Q.; Iliyasu, A.M.; Dong, F.; Hirota, K. A flexible representation and invertible transformations for images on quantum computers. In New Advances in Intelligent Signal Processing. Studies in Computational Intelligence; Ruano, A.E., Várkonyi-Kóczy, A.R., Eds.; Springer: Berlin/Heidelberg, Germany, 2011; Volume 372, pp. 179–202. [Google Scholar]
  12. Zhang, Y.; Lu, K.; Gao, Y.; Wang, M. NEQR: A novel enhanced quantum representation of digital images. Quantum Inf. Process. 2013, 12, 2833–2860. [Google Scholar] [CrossRef]
  13. Wang, Z.; Xu, M.; Zhang, Y. Review of quantum image processing. Arch. Comput. Methods Eng. 2022, 29, 737–761. [Google Scholar]
  14. Zhang, Y.; Lu, K.; Xu, K.; Gao, Y.; Wilson, R. Local feature point extraction for quantum images. Quantum Inf. Process. 2015, 14, 1573–1588. [Google Scholar] [CrossRef]
  15. Li, P.; Liu, X.; Xiao, H. Quantum image weighted average filtering in spatial domain. Int. J. Theor. Phys. 2017, 56, 3690–3716. [Google Scholar] [CrossRef]
  16. Liu, K.; Zhang, Y.; Lu, K.; Wang, X. Restoration for noise removal in quantum images. Int. J. Theor. Phys. 2017, 56, 2867–2886. [Google Scholar] [CrossRef]
  17. Yin, H.; Lu, D.; Zhang, R. Quantum Windowed Fourier Transform and its Application to Quantum Signal Processing. Int. J. Theor. Phys. 2021, 60, 3896–3918. [Google Scholar] [CrossRef]
  18. Asaka, R.; Sakai, K.; Yahagi, R. Quantum circuit for the fast Fourier transform. Quantum Inf. Process. 2020, 19, 277. [Google Scholar] [CrossRef]
  19. Ruiz-Perez, L.; Garcia-Escartin, J.C. Quantum arithmetic with the quantum Fourier transform. Quantum Inf. Process. 2017, 16, 152. [Google Scholar] [CrossRef]
  20. Lohmann, A.W. Image rotation, Wigner rotation, and the fractional Fourier transform. JOSA A 1993, 10, 2181–2186. [Google Scholar] [CrossRef]
  21. Bernardo, L.M.; Soares, O.D. Fractional Fourier transforms and optical systems. Opt. Commun. 1994, 110, 517–522. [Google Scholar] [CrossRef]
  22. Hennelly, B.M.; Sheridan, J.T. Image encryption and the fractional Fourier transform. Optik 2003, 114, 251–265. [Google Scholar] [CrossRef]
  23. Sejdić, E.; Djurović, I.; Stanković, L. Fractional Fourier transform as a signal processing tool: An overview of recent developments. Signal Process. 2011, 91, 1351–1369. [Google Scholar]
  24. Gaglione, D.; Clemente, C.; Ilioudis, C.V.; Persico, A.R.; Proudler, I.K.; Soraghan, J.J.; Farina, A. Waveform design for communicating radar systems using fractional Fourier transform. Digit. Signal Process. 2018, 80, 57–69. [Google Scholar] [CrossRef]
  25. Parasa, V.; Perkowski, M. In Quantum pseudo-fractional fourier transform using multiple-valued logic. In Proceedings of the 2012 IEEE 42nd International Symposium on Multiple-Valued Logic, Victoria, BC, Canada, 14–16 May 2012; pp. 311–314. [Google Scholar]
  26. Lv, C.-H.; Fan, H.-Y.; Li, D.-W. From fractional Fourier transformation to quantum mechanical fractional squeezing transformation. Chin. Phys. B 2015, 24, 020301. [Google Scholar] [CrossRef]
  27. Weimann, S.; Perez-Leija, A.; Lebugle, M.; Keil, R.; Tichy, M.; Gräfe, M.; Heilmann, R.; Nolte, S.; Moya-Cessa, H.; Weihs, G. Implementation of quantum and classical discrete fractional Fourier transforms. Nat. Commun. 2016, 7, 11027. [Google Scholar] [CrossRef] [PubMed]
  28. Madrid, Y.; Molina, M.; Torres, R. Quantum Fractional Fourier Transform. In Frontiers in Optics; Optica Publishing Group: Washington, DC, USA, 2018. [Google Scholar]
  29. Gómez-Echavarría, A.; Ugarte, J.P.; Tobón, C. The fractional Fourier transform as a biomedical signal and image processing tool: A review. Biocybern. Biomed. Eng. 2020, 40, 1081–1093. [Google Scholar]
  30. Shih, C.-C. Fractionalization of Fourier transform. Opt. Commun. 1995, 118, 495–498. [Google Scholar] [CrossRef]
  31. Zhu, B.; Liu, S.; Ran, Q. Optical image encryption based on multifractional Fourier transforms. Opt. Lett. 2000, 25, 1159–1161. [Google Scholar] [CrossRef]
  32. Tao, R.; Lang, J.; Wang, Y. Optical image encryption based on the multiple-parameter fractional Fourier transform. Opt. Lett. 2008, 33, 581–583. [Google Scholar] [CrossRef]
  33. Ran, Q.; Zhang, H.; Zhang, J.; Tan, L.; Ma, J. Deficiencies of the cryptography based on multiple-parameter fractional Fourier transform. Opt. Lett. 2009, 34, 1729–1731. [Google Scholar] [CrossRef]
  34. Ran, Q.; Zhao, T.; Yuan, L.; Wang, J.; Xu, L. Vector power multiple-parameter fractional Fourier transform of image encryption algorithm. Opt. Laser Eng. 2014, 62, 80–86. [Google Scholar] [CrossRef]
  35. Zhao, T.; Ran, Q. The weighted fractional Fourier transform and its application in image encryption. Math. Probl. Eng. 2019, 2019, 4789194. [Google Scholar] [CrossRef]
  36. Liu, S.; Jiang, J.; Zhang, Y.; Zhang, J. Generalized fractional Fourier transforms. J. Phys. A Math. Gen. 1997, 30, 973. [Google Scholar] [CrossRef]
Figure 1. Circuit for quantum phase estimation.
Figure 1. Circuit for quantum phase estimation.
Fractalfract 07 00743 g001
Figure 2. Circuit for the QFrFT.
Figure 2. Circuit for the QFrFT.
Fractalfract 07 00743 g002
Figure 3. QFrFT circuit with both phase gates and QFT changed (the colored parts have changed relative to Figure 2).
Figure 3. QFrFT circuit with both phase gates and QFT changed (the colored parts have changed relative to Figure 2).
Fractalfract 07 00743 g003
Figure 4. QFrFT circuit with only the phase gates changed (the colored parts have changed relative to Figure 2).
Figure 4. QFrFT circuit with only the phase gates changed (the colored parts have changed relative to Figure 2).
Fractalfract 07 00743 g004
Figure 5. Time–frequency distribution of eigenvalues.
Figure 5. Time–frequency distribution of eigenvalues.
Fractalfract 07 00743 g005
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zhao, T.; Chi, Y. Circuit of Quantum Fractional Fourier Transform. Fractal Fract. 2023, 7, 743. https://doi.org/10.3390/fractalfract7100743

AMA Style

Zhao T, Chi Y. Circuit of Quantum Fractional Fourier Transform. Fractal and Fractional. 2023; 7(10):743. https://doi.org/10.3390/fractalfract7100743

Chicago/Turabian Style

Zhao, Tieyu, and Yingying Chi. 2023. "Circuit of Quantum Fractional Fourier Transform" Fractal and Fractional 7, no. 10: 743. https://doi.org/10.3390/fractalfract7100743

APA Style

Zhao, T., & Chi, Y. (2023). Circuit of Quantum Fractional Fourier Transform. Fractal and Fractional, 7(10), 743. https://doi.org/10.3390/fractalfract7100743

Article Metrics

Back to TopTop