Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Decoupling control algorithm based on numerical orthogonal polynomials for a woofer-tweeter adaptive optics system

Open Access Open Access

Abstract

To resolve cooperative control issues of dual wavefront correctors in generalized irregular pupil regions, we propose a decoupling control algorithm based on numerical orthogonal polynomials (NOP). The proposed algorithm firstly deduces NOP from Zernike polynomials in generalized irregular pupil regions. Then, according to wavefront restoration, different spatial frequency aberrations to different wavefront correctors are assigned precisely. Finally, the algorithm calculates and eliminates the cross-coupling between dual wavefront correctors. As observed in numerical simulations and experiments based on a typical woofer-tweeter (W-T) adaptive optics system, NOP decoupling control algorithm restrains the cross-coupling between woofer and tweeter in generalized irregular pupil regions. Moreover, there are obvious advantages over Zernike polynomials decoupling control algorithm in cross-coupling suppression for various scenarios in irregular pupil regions and restoration orders.

© 2021 Optical Society of America under the terms of the OSA Open Access Publishing Agreement

1. Introduction

Adaptive optics (AO) is an essential technology for detecting and correcting wavefront aberrations in real-time [1]. Nowadays, growing number of applications start to use AO systems to correct large amplitude aberrations with lower spatial frequency and small aberrations with higher spatial frequency, such as astronomical observation [2], atmospheric corrections [3], and beam clean-up [4]. Unfortunately, limited by the spatial frequency and actuator stroke of the wavefront corrector [5], it is challenging for AO systems with single wavefront corrector to compensate aberrations which has large amplitude and high spatial frequency simultaneously. Accordingly, a solution to the problem is provided by using dual wavefront correctors [6,7]. For example, the W-T system is a typical double deformable mirrors (DMs) AO system [7], whose woofer is a mirror with large stroke and low spatial frequency and tweeter provides high spatial frequency and small stroke. With appropriate control methods, low order and large amplitude wavefront aberrations can be assigned to the woofer for correcting, while the tweeter is responsible for the residual wavefront aberrations with high orders and small amplitude. However, there is cross-coupling between the dual DMs in the W-T system, which often lead to aberrations correction disorder. Such as, the coupling suppression is tested in the Ref. [8] using Lagrange-multiplier (LM)-based damped least-squares (DLS) control method. Although the coupling between the dual DMs is suppressed by this method, it is not eliminated, which leads to worse stability performance of the AO system. Besides, coupling suppression is not employed in the Ref. [9], which caused a waste of available stroke of the dual DMs.

Therefore, in order to eliminate the cross-coupling, several decoupling control algorithms for dual DMs are proposed and proved, including mode decomposition control algorithm [10,11], zonal decoupling control algorithm [12,13], direct slop-based decoupling control algorithm [1416], and so on. It is worth highlighting that Zernike mode decomposition control algorithm and slope-based decoupling control algorithm demonstrate great effectiveness in suppressing the correlation error of dual DMs within unit circle pupil, among the aforementioned algorithms, according to Ref. [10,15]. However, causing by non-ideality of system device and light intensity distribution, irregular pupil regions often exist. For instance, telescopes and other optical systems are annular systems with central shading [17,18]. Besides, light consistently blinks in real optical systems, which is a result of uneven light intensity distribution [19] caused by light source defects or unstable transmission process [2022]. Where, traditional decoupling control algorithms are no longer applicable. Such as, the orthogonality of Zernike polynomials no longer exists in irregular pupil regions [23]. It prevents the wavefront aberrations from being decomposed orthogonally, resulting in Zernike mode decomposition control algorithm is no longer able to assign accurate wavefront aberrations to the woofer and the tweeter respectively. Moreover, the cross-coupling between the dual DMs is not suppressed validly any more. In addition, slope-based decoupling control algorithm is improved to apply to irregular pupil regions in Ref. [24]. However, according to the experimental results, the aberrations with different spatial frequency distribution require manual selection of the woofer eigen mode, which is a complicated and inflexible process. Thus, it is necessary to develop decoupling control algorithms for generalized irregular pupil regions in W-T systems.

In summary, since Zernike polynomials can characterize aberrations conveniently, it is reasonable to realize the assignment of aberrations with different spatial frequency in this approach. Therefore, a decoupling control algorithm based on NOP is proposed in this paper by modifying Zernike polynomials in generalized irregular pupil regions. Firstly, this algorithm constructs the NOP corresponding to the irregular pupil region. Then, the wavefront aberrations assigned to woofer and tweeter are calculated. Finally, to keep the tweeter from generating components coupled to the woofer, the NOP is applied to derive a coupling constraint matrix to modify the tweeter voltage control signal.

The paper is organized as follows: The principle of the decoupling control algorithm based on the NOP is presented in Section 2. The validity of NOP decoupling algorithm is proved in theory by a numerical simulation model established in section 3 and compared with Zernike polynomials decoupling control algorithm. In Section 4, the feasibility of the proposed algorithm is verified by a real W-T AO system, and compared with Zernike polynomials decoupling control algorithm. Section 5 is the conclusion.

2. Principle

In the W-T AO system, the greatest challenge is to correctly assign various aberrations to the woofer and tweeter according to deformable mirror correction capabilities. Firstly, it is assumed that the light intensity in the pupil region is fully filled up under uniform distribution. Then the residual wavefront slope g(2n×1) detected by an n sub-apertures Hartmann wavefront sensor can be restored using an m-order Zernike polynomials as:

$${\boldsymbol g} = {\boldsymbol {Za}}.$$
where Z(2n×m) is the wavefront restoration matrix corresponding to Zernike polynomials and a(m×1) denotes the coefficient vector for Zernike polynomials. a can be rewritten as:
$${\boldsymbol a} = {{\boldsymbol Z}^ + }{\boldsymbol g}.$$
where Z+ represents the pseudo-inverse matrix of Z. In fact, the non-uniform distribution of light intensity causes a lack of light or saturation of Hartmann sub-apertures, which brings in irregular pupil regions. At this moment, orthogonality of Zernike polynomials loses its effectiveness in irregular pupil regions, which affect the result of wavefront restoration. Therefore, Z should be modified to achieve consistency with the irregular pupil region as:
$${{\boldsymbol Z}_{{\boldsymbol {new}}}}({{{\boldsymbol i}_{\boldsymbol x}},{\boldsymbol j}} )= \left\{ {\begin{array}{{c}} {Z({{{\boldsymbol i}_{\boldsymbol x}},{\boldsymbol j}} ),({{{\boldsymbol I}_{{\boldsymbol {min}}}} \le {{\boldsymbol I}_{\boldsymbol i}} \le {{\boldsymbol I}_{{\boldsymbol {max}}}}} )}\\ {\mathbf{0},({{{\boldsymbol I}_{\boldsymbol i}} \ge {{\boldsymbol I}_{{\boldsymbol {max}}}}\; \; {\boldsymbol {or}}\; \; {{\boldsymbol I}_{\boldsymbol i}} \le {{\boldsymbol I}_{{\boldsymbol {min}}}}} )} \end{array}} \right..$$
where Z(ix,j) is the value in x direction at the ith sub-aperture of the jth Zernike polynomials, Ii is the sum of light intensity at the ith sub-aperture, and the minimum and maximum thresholds of light intensity denote as Imin and Imax respectively. Besides, Z(iy,j) is modified in the same way as Z(ix,j), and then we can get the modified matrix Znew(2n×m) from Z.

According to the modification process and Ref. [23], we derive the NOP applicable to irregular pupil regions by Znew:

$${{\boldsymbol C}_{\boldsymbol Z}} = {\boldsymbol Z}_{{\boldsymbol {new}}}^{\boldsymbol T}{{\boldsymbol Z}_{{\boldsymbol {new}}}}/2{\boldsymbol n},$$
$${\boldsymbol Q} = {\boldsymbol {chol}}({{{\boldsymbol C}_{\boldsymbol Z}}} ),$$
$${\boldsymbol M} = {\boldsymbol {inv}}({{{\boldsymbol Q}^{\boldsymbol T}}} ),$$
$${\boldsymbol F} = {{\boldsymbol Z}_{{\boldsymbol {new}}}}{{\boldsymbol M}^{\boldsymbol T}}.$$
where ZnewT is the transpose matrix of Znew and CZ represents the correlation of the column vectors of the Znew. Equation (5) can be solved uniquely by Cholesky decomposition method to get Q, and the transform matrix M is calculated by Eq. (6). Finally, the new wavefront restoration matrix F(2n×m) corresponding to the irregular region can be obtained by Eq. (7). And the NOP demonstrates similar spatial frequency characteristics with the Zernike polynomials, which means that its spatial frequency gets higher when the polynomial order increases.

Next, using the pseudo-inverse matrix F+ of F, the coefficients of each order corresponding to NOP are calculated noting from Eq. (2):

$${\boldsymbol r} = {{\boldsymbol F}^ + }{\boldsymbol g}.$$

Later, the low order mode coefficients in the aberrations are selected, according to the diagonal matrix IW(m×m), they are allocated to the woofer as the low spatial frequency aberrations. For instance, if IW=diag(1,1,1,0,0,0,…), the woofer is responsible for correcting the first 3 orders modes rw of wavefront aberrations:

$${{\boldsymbol r}_{\boldsymbol w}} = {{\boldsymbol I}_{\boldsymbol W}}{\boldsymbol r}.$$

Then the wavefront aberrations gw to be corrected by woofer will be:

$${{\boldsymbol g}_{\boldsymbol w}} = {\boldsymbol F}{{\boldsymbol r}_{\boldsymbol w}}.$$

At this time, combining the slope response matrix RW(2n×nw) of woofer with nw actuators, we get the cumulative voltage vw of woofer as:

$${{\boldsymbol v}_{\boldsymbol w}} = {\boldsymbol R}_{{\boldsymbol W}\_{\boldsymbol {new}}}^ + {\boldsymbol F}{{\boldsymbol r}_{\boldsymbol w}}.$$
where RW_new(2n×nw) is modified by RW in the same way as Eq. (3), and RW_new+ is the pseudo-inverse matrix of RW_new. Similarly, the modified slope response matrix RT_new(2n×nt) of tweeter with nt actuators can be obtained. Lastly, we can calculate the driving voltage Vw of woofer according to PI controller:
$${{\boldsymbol V}_{\boldsymbol w}}({{\boldsymbol k} + 1} )= {\boldsymbol a} \times {{\boldsymbol V}_{\boldsymbol w}}({\boldsymbol k} )+ {\boldsymbol b} \times {{\boldsymbol v}_{\boldsymbol w}}.$$

Where a and b are the control parameters of PI controller. Then, according to Eq. (10), gt corrected by the tweeter can be calculated as:

$${{\boldsymbol g}_{\boldsymbol t}} = {\boldsymbol g} - {{\boldsymbol g}_{\boldsymbol w}}.$$

Thus, the cumulative voltage vt and the driving voltage Vt of the tweeter are:

$${{\boldsymbol v}_{\boldsymbol t}} = {\boldsymbol R}_{{\boldsymbol T}\_{\boldsymbol {new}}}^ + {{\boldsymbol g}_{\boldsymbol t}},$$
$${{\boldsymbol V}_{\boldsymbol t}}({{\boldsymbol k} + 1} )= {\boldsymbol a} \times {{\boldsymbol V}_{\boldsymbol t}}({\boldsymbol k} )+ {\boldsymbol b} \times {{\boldsymbol v}_{\boldsymbol t}}.$$

At this moment, to restrict the cross-coupling between woofer and tweeter, the following coupling constraint matrix Rset(nt×nt) is constructed by Eq. (16):

$${{\boldsymbol R}_{{\boldsymbol {set}}}} = {\boldsymbol I} - {\boldsymbol R}_{{\boldsymbol T}\_{\boldsymbol {new}}}^ + {\boldsymbol F}{{\boldsymbol I}_{\boldsymbol W}}{{\boldsymbol F}^ + }{{\boldsymbol R}_{{\boldsymbol T}\_{\boldsymbol {new}}}}.$$
where I(nt×nt) is the unit matrix. So, the actual voltage Vt sent to the tweeter without coupling can be calculated using Vt and Rset as Eq. (17):
$${\boldsymbol V}_{\boldsymbol t}^{\boldsymbol \prime} = {{\boldsymbol R}_{{\boldsymbol {set}}}}{{\boldsymbol V}_{\boldsymbol t}}.$$

3. Numeral simulation

3.1 Woofer-tweeter system simulation model

The NOP decoupling control algorithm proposed in this paper mainly focuses on the aberrations distribution and coupling suppression between woofer and tweeter in the AO system for generalized irregular pupil regions. To prove the validity of this algorithm, a woofer-tweeter system is established, using a 19-actuator deformable mirror (DM) as the woofer to correct the low order aberrations and a DM with 127-actuator as the tweeter to correct the high order aberrations. The configuration between the actuators of DM and the sub-apertures of Hartmann wavefront sensor is shown in Fig. 1.

 figure: Fig. 1.

Fig. 1. The configuration between the actuators of DM and the sub-apertures of Hartmann wavefront sensor.

Download Full Size | PDF

3.2 Validation of the NOP decoupling control algorithm

Firstly, three shading ways are employed to verify the effectiveness of NOP decoupling algorithm: (1) no shading, (2) Semi-circular shading, (3) Irregular shading. And the initial aberrations under the three shading ways are shown in Fig. 2(a1), (a2), and (a3). Then, according to the three shading ways, the new wavefront restoration matrix F can be obtained by Eq. (4)–(7), whose column vectors are verified to be orthogonal in all shading ways, as shown in Fig. 3. Therefore, we can assign the aberration to the dual DMs and restore the wavefront aberration with F.

 figure: Fig. 2.

Fig. 2. Correction results with different shading ways.

Download Full Size | PDF

 figure: Fig. 3.

Fig. 3. Column vectors of F orthogonality under different shading ways.

Download Full Size | PDF

For the given initial aberrations, the first 3 order modes are isolated by IW=diag(1,1,1,0,0,0,…) and assigned to the woofer, and the residual aberrations are corrected by the tweeter. In three shading ways, results of the woofer correction are shown in Fig. 2(b1), (b2), and (b3). Results of the tweeter correction are shown in Fig. 2(c1), (c2), and (c3). And the root mean square (RMS) of wavefront aberrations are 0.013λ, 0.014λ, and 0.015λ, respectively, shown in Fig. 2(d1), (d2), and (d3).

In addition, we present the aberrations distribution and correction of woofer and tweeter under three different shading ways on the first 10 order modes. From Fig. 4(a), (b), and (c), the first 3 orders of aberrations are corrected by the woofer and have no coefficients from tweeter.

 figure: Fig. 4.

Fig. 4. Coupling suppression performance of the algorithm under different shading ways: (a)No shading, (b) Semi-circular shading, (c) Irregular shading, (d) Curves of RMS and c.

Download Full Size | PDF

To visualize the decoupling performance, the coupling coefficient c is defined as Eq. (18):

$${\boldsymbol c} = {\boldsymbol {max}}({{\boldsymbol {abs}}({{{\boldsymbol F}^ + }{{\boldsymbol g}_{\boldsymbol t}}} )} ).$$
gt is decomposed on the mode items, then the maximum value of the first 3 order coefficients in absolute value is the coupling coefficient c. When tweeter has no cross-coupling with woofer, c is 0 at this time. Figure 4(d) shows that the NOP decoupling control algorithm maintains c converge to 0 approximately for all these shading ways, while the RMS of residual wavefront converge to 0.014λ approximately.

According to Fig. 4, it can be seen that NOP decoupling control algorithm drives the woofer and tweeter to correct the aberration and suppress the coupling efficiently where the pupil region has been changed.

3.3 Comparison with Zernike polynomials decoupling control algorithm

We set up two simulations to compare the performance of Zernike polynomials decoupling control algorithm and NOP decoupling control algorithm: (1) Different modes are taken to correct the aberrations in the semi-circular shading way, (2) Closed-loop correction is accomplished in three different shading ways using the first 48 order modes. For both simulations, we use the woofer to correct the first 3 orders aberrations, while the tweeter takes charge of the rest.

3.3.1 Different orders of restoration

The initial aberration in the simulation is shown in Fig. 2(a2), and the wavefront correction is completed by using the first 5, first 10 and first 15 order modes, respectively. The correction results are shown in Fig. 5.

 figure: Fig. 5.

Fig. 5. Correction results with different restoring orders: (a) RMS of residual wavefront, (b) RMS of woofer correction, (c) RMS of tweeter correction, (d) Curves of c.

Download Full Size | PDF

As shown in Fig. 5(a), algorithm and Zernike polynomials decoupling control algorithm. However, compared with Zernike polynomials decoupling control algorithm, the RMS of residual wavefront has better convergence features when using the NOP decoupling control algorithm during the closed-loop process, shown in Fig. 5(a). Moreover, from Fig. 5(b) and (c), Zernike polynomials decoupling control algorithm has different the RMS of woofer and tweeter correction when the restoration order is varying. Meanwhile, using the NOP decoupling control algorithm, the RMS of woofer and tweeter correction are consistent, shown in Fig. 5(b) and (c). This indicates that the aberration assigned to woofer and tweeter by NOP decoupling control algorithm is always the same when the restoration order is changed. In addition, it is shown in Fig. 5(d) that the coupling coefficients of the NOP decoupling control algorithm are consistent at 0 approximately, while those of the Zernike polynomials decoupling control algorithm are large and changeful.

3.3.2 Different shading ways

The wavefront correction is completed with three different shading ways, which are no shading, semi-circular shading, and irregular shading, shown in Fig. 2(a1), (a2) and (a3), respectively. The correction results are shown in Fig. 6. Presented on Fig. 6(a), the RMS of residual wavefront converge to 0.012λ and the woofer and tweeter have the same correction for both NOP decoupling control algorithm and Zernike polynomials decoupling control algorithm when the no shading way is used. However, in case of the two other shading ways, which are shown in Fig. 6(b) and (c), the RMS of residual wavefront and dual DMs correction converge to stabilization with fewer iterations during the closed-loop process when using NOP decoupling control algorithm compared with Zernike polynomials decoupling control algorithm. Furthermore, shown in Fig. 6(d), the coupling coefficients of the NOP decoupling control algorithm are stable at 0 approximately, while the values of the coupling coefficients of Zernike polynomials decoupling control algorithm are larger at different shading ways.

 figure: Fig. 6.

Fig. 6. Correction results with different shading ways: (a)No shading, (b) Semi-circular shading, (c) Irregular shading, (d) Curves of c.

Download Full Size | PDF

In conclusion, NOP decoupling control algorithm in the W-T system can correct the aberration well in the irregular pupil region. Moreover, this algorithm has a more stable coupling suppression performance for different restoration orders and different irregular pupil regions than Zernike polynomials decoupling control algorithm. In other words, there is high robustness in NOP decoupling control algorithm.

4. Experiment

4.1 Experimental conditions

To further test the performance of NOP decoupling control algorithm, the W-T AO system experimental platform is built presented on Fig. 7. A light source with a wavelength of 650nm is employed here in Fig. 7(a), and an aberration plate is placed in front of the source. Two 140mm focal length lenses are placed between woofer and tweeter. And the configuration between the actuators of DM and the sub-apertures of Hartmann wavefront sensor is shown in Fig. 7(b). Figure 7(c) shows the structure of wavefront sensor, which contains Hartmann wavefront sensor, near field CCD camera and far field CCD camera. Moreover, the Hartmann wavefront sensor consists of a 12×12 lens array and an eight-bit CCD camera with a target range of 372×372 pixels, where each sub-aperture contains both x-axis and y-axis directions, and the sub-aperture size is 31×31. Both woofer and tweeter are 59-actuator deformable mirrors.

 figure: Fig. 7.

Fig. 7. A W-T AO system: (a) Experimental platform, (b) The configuration between 59-actuators DM and the sub-apertures of Hartmann wavefront sensor, (c) Structure of wavefront sensor.

Download Full Size | PDF

4.2 Validation of the NOP decoupling control algorithm

To verify the effectiveness of the NOP decoupling algorithm, we employ three shading ways shown in Fig. 8, which are no shading, Semi-circular shading, and Irregular shading. The first 3 modes are isolated and assigned to the woofer, and the residual aberrations are corrected by the tweeter. In three shading ways, results of the woofer correction are shown in Fig. 9(a1), (a2), and (a3). Results of the tweeter correction are shown in Fig. 9(b1), (b2), and (b3). Moreover, the RMS of wavefront aberrations are 0.051λ, 0.063λ, and 0.072λ, respectively, shown in Fig. 9(c1), (c2), and (c3).

 figure: Fig. 8.

Fig. 8. Three shading ways on the CCD.

Download Full Size | PDF

 figure: Fig. 9.

Fig. 9. Correction results with different shading ways.

Download Full Size | PDF

Moreover, Fig. 10(a), (b), and (c) shows the projection of the woofer and tweeter correction at the first 10 order modes for the three shading ways, and it can be seen that tweeter has almost no coefficients on the first 3 orders. From Fig. 10(d), the NOP decoupling control algorithm keeps c stable to 0 approximately, while the RMS of the residual wavefront converges to 0.051λ, 0.063λ, and 0.072λ, respectively. Figure 10 shows that NOP decoupling control algorithm restrains the coupling error effectively between woofer and tweeter for the changing pupil region.

 figure: Fig. 10.

Fig. 10. Coupling suppression performance of the algorithm under different shading ways: (a)No shading, (b) Semi-circular shading, (c) Irregular shading, (d) Curves of RMS and c.

Download Full Size | PDF

4.3 Comparison with Zernike polynomials decoupling control algorithm

To compare the performance of Zernike polynomials decoupling control algorithm and NOP decoupling control algorithm, two comparison experiments are carried out. There are restoration matrices of different orders to complete the aberrations closed-loop correction in the first experiment where a semi-circular shading way is selected. In the second experiment, the first 48 order modes are taken, applying three shading ways. For both experiments, we use the woofer to correct the first 3 orders aberrations, while the tweeter deals with the rest.

4.3.1 Different orders of restoration

The wavefront correction is completed using the first 5, first 10 and first 15 order modes, respectively, under the semi-circular shading way shown in Fig. 8(b). And the correction results are shown in Fig. 11. Shown in Fig. 11(a), the convergence characteristics of the RMS of the residual wavefront using the NOP decoupling control algorithm is better than that using the Zernike polynomials decoupling control algorithm. Meanwhile, with the different restoring order comparisons, it is demonstrated in Fig. 11(b) and (c) that the convergence process of the dual DMs correction remains the same compared with the Zernike polynomials decoupling control algorithm. Accordingly, coupling coefficients of the NOP decoupling control algorithm are consistent and almost 0 in Fig. 11(d), while those for the closed-loop process are large by using the Zernike polynomials decoupling control algorithm.

 figure: Fig. 11.

Fig. 11. Correction results with different restoring orders: (a) RMS of the residual wavefront, (b) RMS of woofer correction, (c) RMS of tweeter correction, (d) Curves of c.

Download Full Size | PDF

4.3.2 Different shading ways

The wavefront correction is completed with three shading ways in Fig. 8, which are no shading, semi-circular shading, and irregular shading, respectively. The correction results are shown in Fig. 12. As shown in Fig. 12(a), the RMS of residual wavefront are 0.051λ approximately and the correction of dual DMs remains stable for both the NOP decoupling control algorithm and the Zernike polynomials decoupling control algorithm when the no shading way is used. However, compared with the Zernike polynomials decoupling control algorithm, shown in Fig. 12(b) and (c), the RMS of residual wavefront and DMs correction converge to stabilization with fewer iterations during aberration correction process using the NOP decoupling control algorithm with the two other shading ways, while the coupling coefficients are 0 approximately, shown in Fig. 12(d). In contrast, the RMS of residual wavefront and DMs correction have not converged during iterations and the coupling coefficients are large as well as not converged by using Zernike polynomials decoupling control algorithm, shown in Fig. 12(b), (c), and (d).

 figure: Fig. 12.

Fig. 12. Correction results with different shading ways:(a) No shading, (b) Semi-circular shading, (c) Irregular shading, (d) Curves of c.

Download Full Size | PDF

In conclusion, NOP decoupling control algorithm in the W-T system can complete closed-loop correction well in irregular pupil regions. Moreover, this algorithm has a stable coupling suppression performance for different restoration orders and different irregular pupil regions. Hence, the robustness of NOP decoupling control algorithm is much better than that of Zernike polynomials decoupling control algorithm. The conclusions of experiments are consistent with those of numerical simulations and confirm the validity of NOP decoupling control algorithm for the W-T AO system. And as illustrated in Ref. [23], the Zernike polynomials lead to problems such as mode confusion due to non-orthogonality in irregular pupil regions. That results in the coupling suppression of the Zernike polynomials decoupling control algorithm are weak and unstable in irregular pupil regions.

5. Conclusion

In conclusion, this paper proposes and discusses a decoupling control algorithm based on NOP for solving the cross-coupling of dual wavefront correctors in the generalized irregular pupil region. The numerical simulations and experiments show that this algorithm can drive woofer and tweeter to correct aberrations with different spatial frequency and restrain cross-coupling under different restoring orders and different shading ways. Meanwhile, NOP decoupling control algorithm demonstrates better performance than Zernike polynomials decoupling control algorithm. In this paper, we have discussed various scenarios under three different irregular pupil regions to verify the feasibility of NOP decoupling control algorithm. In fact, because of the orthogonality of NOP in generalized irregular pupil regions, we can draw a conclusion that the NOP decoupling control algorithm is applicable for the AO system with multiple wavefront correctors in any other generalized irregular pupil regions.

Funding

National Natural Science Foundation of China (61805251, 61875203, 62005285); Youth Innovation Promotion Association of the Chinese Academy of Sciences (2017429); Western Youth Scholar A CAS.

Disclosures

The authors declare no conflicts of interest.

Data availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

References

1. X. Ye and F. Shen, “Analysis of aiming performance limitation for an optical system in atmospheric turbulence,” J. Opt. Technol. 85(10), 603–609 (2018). [CrossRef]  

2. D. Wu, C. Yang, H. Li, P. Zhang, X. Zhang, Z. Cao, Q. Mu, and L. Xuan, “Astronomical observation by 2-meter telescope based on liquid crystal adaptive optics with phase diversity,” Opt. Commun. 439, 129–132 (2019). [CrossRef]  

3. L. Shao, Z. Cao, Q. Mu, P. Zhang, L. Yao, S. Wang, L. Hu, and L. Xuan, “Experimental analysis of the fitting error of diffractive liquid crystal wavefront correctors for atmospheric turbulence corrections,” Opt. Commun. 367, 199–205 (2016). [CrossRef]  

4. Y. Guo, Q. Peng, Y. Bo, Z. Chen, Y. Li, L. Zhang, C. Shao, L. Yuan, B. Wang, J. Xu, J. Xu, H. Gao, Y. Xu, B. Lai, C. Su, S. Ma, and T. Cheng, “24.6 kW near diffraction limit quasi-continuous-wave Nd:YAG slab laser based on a stable–unstable hybrid cavity,” Opt. Lett. 45(5), 1136–1139 (2020). [CrossRef]  

5. H. Yang, G. Liu, C. Rao, Y. Zhang, and W. Jiang, “Combinational-deformable-mirror adaptive optics system for compensation of high-order modes of wavefront,” Chin. Opt. Lett. 5(8), 435–437 (2007).

6. R. Zhou, F. Shen, H. Ye, and X. Li, “Double deformable mirrors’ control based on volatage decoupling for adaptive optics system,” Proc. SPIE 9255, 925538 (2015). [CrossRef]  

7. O. Keskin, P. Hampton, R. Conan, C. Bradley, A. Hilton, and C. Blain, “Woofer-Tweeter Adaptive Optics Test Bench,” in First NASA/ESA Conference on Adaptive Hardware and Systems (AHS'06) (IEEE, 2006), pp. 74–80.

8. W. Zou, X. Qi, and S. A. Burns, “Woofer-tweeter adaptive optics scanning laser ophthalmoscopic imaging based on Lagrange-multiplier damped least-squares algorithm,” Biomed. Opt. Express 2(7), 1986–2004 (2011). [CrossRef]  

9. R. Conan, C. Bradley, P. Hampton, O. Keskin, A. Hilton, and C. Blain, “Distributed modal command for a two-deformable-mirror adaptive optics system,” Appl. Opt. 46(20), 4329–4340 (2007). [CrossRef]  

10. W. Liu, L. Dong, P. Yang, X. Lei, H. Yan, and B. Xu, “A Zernike mode decomposition decoupling control algorithm for dual deformable mirrors adaptive optics system,” Opt. Express 21(20), 23885–23895 (2013). [CrossRef]  

11. S. Hu, B. Xu, X. Zhang, J. Hou, J. Wu, and W. Jiang, “Double-deformable-mirror adaptive optics system for phase compensation,” Appl. Opt. 45(12), 2638–2642 (2006). [CrossRef]  

12. W. Liu, L. Dong, P. Yang, and B. Xu, “Zonal decoupling algorithm for dual deformable mirror adaptive optics system,” Chin. Opt. Lett. 14(2), 020101 (2016). [CrossRef]  

13. W. Zou, X. Qi, and S. A. Burns, “Wavefront-aberration sorting and correction for a dual-deformable-mirror adaptive-optics system,” Opt. Lett. 33(22), 2602–2604 (2008). [CrossRef]  

14. C. Li, N. Sredar, H. Queener, K. Ivers, and J. Porter, “Direct Slope Algorithm for Simultaneous Control of a Dual Deformable Mirror Adaptive Optics System,” Inves. Ophth. Vis. Sci. 51(13), 2309 (2010).

15. T. Cheng, W. Liu, B. Pang, and B. Xu, “A slope-based decoupling algorithm to simultaneously control dual deformable mirrors in a woofer–tweeter adaptive optics system,” Chin. Phys. B 27(7), 070704 (2018). [CrossRef]  

16. C. Li, N. Sredar, K. M. Ivers, H. Queener, and J. Porter, “A correction algorithm to simultaneously control dual deformable mirrors in a woofer-tweeter adaptive optics system,” Opt. Express 18(16), 16671–16684 (2010). [CrossRef]  

17. O. Karcı and M. Ekinci, “Design of a high-precision, 0.5 m aperture Cassegrain collimator,” Appl. Opt. 59(27), 8434–8442 (2020). [CrossRef]  

18. A. C. Dávila, A. P. R. Cortés, and S. V. Y. Montiel, “Exact calculation of conic constants and baffles for any two-mirror aplanatic telescope,” Appl. Opt. 59(22), 6737–6743 (2020). [CrossRef]  

19. L. Xue, Y. Pang, W. Liu, L. Liu, H. Pang, A. Cao, L. Shi, Y. Fu, and Q. Deng, “Fabrication of random microlens array for laser beam homogenization with high efficiency,” Micromachines 11(3), 338 (2020). [CrossRef]  

20. M. C. Roggemann and D. J. Lee, “Two-deformable-mirror concept for correcting scintillation effects in laser beam projection through the turbulent atmosphere,” Appl. Opt. 37(21), 4577–4585 (1998). [CrossRef]  

21. W. Xiang, P. Yang, S. Wang, B. Xu, and H. Liu, “Underwater image enhancement based on red channel weighted compensation and gamma correction model,” Opto-Electronic Adv. 1(10), 18002401 (2018). [CrossRef]  

22. Y. Ren, J. Jian, J. Wang, C. Zhu, and W. Xia, “Analysis of the spatial intensity distribution of a CW laser light through tissue-like turbid media,” J. Mod. Optic. 68(2), 116–123 (2021). [CrossRef]  

23. J. Ye, W. Wang, Z. Gao, Z. Liu, S. Wang, P. Benítez, J. C. Miñano, and Q. Yuan, “Modal wavefront estimation from its slopes by numerical orthogonal transformation method over general shaped aperture,” Opt. Express 23(20), 26208–26220 (2015). [CrossRef]  

24. T. Cheng, Z. Xu, K. Yang, S. Wang, X. He, P. Yang, and B. Xu, “A decoupling control algorithm for Woofer–Tweeter adaptive optics system in generalized irregular pupil region,” Opt. Commun. 472, 125856 (2020). [CrossRef]  

Data availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (12)

Fig. 1.
Fig. 1. The configuration between the actuators of DM and the sub-apertures of Hartmann wavefront sensor.
Fig. 2.
Fig. 2. Correction results with different shading ways.
Fig. 3.
Fig. 3. Column vectors of F orthogonality under different shading ways.
Fig. 4.
Fig. 4. Coupling suppression performance of the algorithm under different shading ways: (a)No shading, (b) Semi-circular shading, (c) Irregular shading, (d) Curves of RMS and c.
Fig. 5.
Fig. 5. Correction results with different restoring orders: (a) RMS of residual wavefront, (b) RMS of woofer correction, (c) RMS of tweeter correction, (d) Curves of c.
Fig. 6.
Fig. 6. Correction results with different shading ways: (a)No shading, (b) Semi-circular shading, (c) Irregular shading, (d) Curves of c.
Fig. 7.
Fig. 7. A W-T AO system: (a) Experimental platform, (b) The configuration between 59-actuators DM and the sub-apertures of Hartmann wavefront sensor, (c) Structure of wavefront sensor.
Fig. 8.
Fig. 8. Three shading ways on the CCD.
Fig. 9.
Fig. 9. Correction results with different shading ways.
Fig. 10.
Fig. 10. Coupling suppression performance of the algorithm under different shading ways: (a)No shading, (b) Semi-circular shading, (c) Irregular shading, (d) Curves of RMS and c.
Fig. 11.
Fig. 11. Correction results with different restoring orders: (a) RMS of the residual wavefront, (b) RMS of woofer correction, (c) RMS of tweeter correction, (d) Curves of c.
Fig. 12.
Fig. 12. Correction results with different shading ways:(a) No shading, (b) Semi-circular shading, (c) Irregular shading, (d) Curves of c.

Equations (18)

Equations on this page are rendered with MathJax. Learn more.

g = Z a .
a = Z + g .
Z n e w ( i x , j ) = { Z ( i x , j ) , ( I m i n I i I m a x ) 0 , ( I i I m a x o r I i I m i n ) .
C Z = Z n e w T Z n e w / 2 n ,
Q = c h o l ( C Z ) ,
M = i n v ( Q T ) ,
F = Z n e w M T .
r = F + g .
r w = I W r .
g w = F r w .
v w = R W _ n e w + F r w .
V w ( k + 1 ) = a × V w ( k ) + b × v w .
g t = g g w .
v t = R T _ n e w + g t ,
V t ( k + 1 ) = a × V t ( k ) + b × v t .
R s e t = I R T _ n e w + F I W F + R T _ n e w .
V t = R s e t V t .
c = m a x ( a b s ( F + g t ) ) .
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.