Next Article in Journal
Performance Evaluation of a New Sport Watch in Sleep Tracking: A Comparison against Overnight Polysomnography in Young Adults
Previous Article in Journal
3D Galileo Reference Antenna Pattern for Space Service Volume Applications
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Fast Prediction Method for Scattering Parameters of Rigid-Flex PCBs Based on ANN

1
Shanghai Institute of Technical Physics, Chinese Academy of Sciences, Shanghai 200083, China
2
School of Electronic, Electrical and Communication Engineering, University of Chinese Academy of Sciences, Beijing 100049, China
3
Key Laboratory of Infrared System Detection and Imaging Technologies, Shanghai Institute of Technical Physics, Chinese Academy of Sciences, Shanghai 200083, China
4
Key Laboratory of Ultra High Speed Circuit Design and Electromagnetic Compatibility, Ministry of Education, Xidian University, Xi’an 710071, China
*
Author to whom correspondence should be addressed.
Sensors 2024, 24(7), 2221; https://doi.org/10.3390/s24072221
Submission received: 18 January 2024 / Revised: 20 March 2024 / Accepted: 28 March 2024 / Published: 30 March 2024
(This article belongs to the Section Sensor Networks)

Abstract

:
InGaAs detection systems have been increasingly used in the aerospace field, and due to the high signal-to-noise ratio requirements of short-wave infrared quantitative payloads, there is an urgent need for methods for the rapid and precise evaluation and the optimal design of these systems. The rigid-flex printed circuit board (PCB) is a vital component of InGaAs detectors, as its grid ground plane design parameters impact parasitic capacitance and thus affect weak infrared analog signals. To address the time-intensive and costly nature of design optimization achieved with simulations and experimental measurements, we propose an innovative method based on a neural network to predict the scattering parameters of rigid-flex boards for InGaAs detection links. This is the first study in which the effects of rigid-flex boards on weak infrared detection signals have been considered. We first obtained sufficient samples with software simulation. A backpropagation (BP) neural network prediction model was trained on existing sample sets and then verified on a rigid-flex board used in a crucial aerospace short-wave infrared quantitative mission. The model efficiently and accurately predicted high-speed interconnect scattering parameters under various rigid-flex board grid plane parameter conditions. The prediction error was less than 1% compared with a 3D field solver, indicating an overcoming of the iterative optimization inefficiency and showing improved design quality for InGaAs detection circuits.

1. Introduction

In recent decades, short-wave infrared InGaAs focal plane array detectors [1] have been extensively used in aerospace remote sensing. With the continuous improvement in remote sensing payload indicators, the total number of remote sensing data has increased rapidly [2]. It has thus become common to see the data transfer capacity of high-resolution satellites in the order of GHz, a trend that in turn has led to the increased use of large-area array detectors and a corresponding increase in readout frequency. Quantification, which is essential in quantitative payloads, requires an improvement in signal fidelity. However, in practice, there is often a need to balance the design requirements of the optical system with those of the mechanical structure. Generally, the InGaAs detectors cannot be directly coupled to the driver acquisition circuit board. However, rigid-flex boards, with their compact, lightweight, and highly reliable mechanical structure [3], play a vital role in the connection between short-wave infrared photodetectors [4] and information acquisition circuits. Nevertheless, high-speed signal transmission on rigid-flex boards can significantly distort infrared analog signals [5] before they reach the driver acquisition module, presenting potential signal integrity challenges in a low-noise infrared detection system design. To enable accurate quantization under high signal-to-noise ratio requirements, it is imperative to design each component of the detection link precisely. Therefore, modeling and optimizing the transmission characteristics of rigid-flex boards are necessary.
Grid ground planes are quite commonly used in a modern rigid-flex PCB design to meet the requirements of impedance matching and interlayer interconnections [6,7,8,9,10]. Traces over the grid ground plane are associated with signal integrity (SI) degradation, electromagnetic interference (EMI), and ground bounce in general [11]. Hence, it is vital to investigate how the electrical properties of a transmission line change as a function of the geometry of the grid ground planes. Scattering parameters (S-parameters), as a beneficial means for describing the behavior of linear passive electrical properties, contain much information about interconnect transmission characteristics such as the impedance, crosstalk, reflection, and attenuation of the transmitted signals [12]. Moreover, recent studies on rigid-flex PCBs have focused on how various design parameters in grid grounding can be applied to impedance matching [13,14,15], but how this affects the S-parameters has not yet been investigated. To achieve optimal performance in the signal link transmission, the link design parameters must be adjusted repeatedly for S-parameter optimization, which can improve rigid-flex PCBs’ analysis and optimization capabilities in infrared fine detection systems, resulting in better design effects and lower noise levels.
In the past, S-parameters were frequently obtained in the early design phase with time-consuming and labor-intensive experimental measurements [16,17,18,19] which required precise test equipment. It is now widely recognized that the transmission line referenced to a grid ground plane is a periodically varying structure [20,21]. Therefore, with the advancement of modern electronic design automation (EDA), 3D field solvers [22,23] have been widely adopted to simulate experimental measurements. Substituting expensive and time-consuming physical measurements with simulated results obtained with 3D field solvers significantly reduces the cost and design cycle associated with fabrication and measurements. Macromodeling is a fundamental tool for the characterization of high-frequency interconnects and has been applied to the improvement in the efficiency of 3D field solvers, the modeling of lossy transmission lines, and design optimization [24,25,26,27]. Many popular macromodeling techniques exploit vector fitting (VF) [24,25] to build a rational function approximation for interconnect transfer functions. A novel, low-cost multivariate parametric macro modeling technique based on the kriging method for the residuals and poles of the reflection coefficients and the zeros of the transmission coefficients of microwave filters was proposed in [26]. A technique to derive macro-models of lumped-element networks for distributed microwave circuits by using system identification was presented in [27]. Nevertheless, both traditional 3D field solvers and those based on macromodeling optimization are becoming increasingly costly in terms of computer resources and simulation time, even more significantly so when the number of design choices and parameters increases [28].
In recent decades, powerful techniques of machine learning (ML) models such as artificial neural networks (ANNs) [29,30] and support vector machines (SVMs) [31,32] for the analysis and modeling of complex nonlinear problems have been developed as promising alternatives to the 3D field solver approach. ML models have high simulation speeds and do not need to fully understand the internal structure of modeling objects. Thus, S-parameter modeling technology based on ANNs has attracted increasing attention, and preliminary feasibility studies have been recently carried out [33,34,35,36]. While both SVMs and ANNs demonstrate efficacy in fitting, they have yet to be extensively applied in grid ground plane design parameters. Therefore, a method exhibiting low resource utilization and high convergence speed to achieve a balance between resource consumption and accuracy in predicting S-parameters would be beneficial to optimizing rigid-flex PCBs for InGaAs detection links.
In this study, we established a method for rapidly predicting the S-parameters of rigid-flexed PCBs based on a BP neural network for infrared InGaAs detection systems. We compared the generalized regression neural network (GRNN) model proposed in [36] and a traditional 3D field solver. The proposed method was implemented with a simple BP neural network structure, which can effectively control and prevent overfitting. Furthermore, the BP neural network only requires the storage of the weight matrix between each pair of layers, resulting in minimal memory usage. Compared with the traditional 3D field solver methods and the GRNN model, this approach significantly enhanced computational efficiency, reduced costs, and maintained the prediction error to less than 1%. This article is organized as follows: In Section 2, we briefly review the application and problems of rigid-flex PCBs in infrared InGaAs detection links and introduce BP networks. The details of our S-parameter modeling method for a rigid-flex board are also elucidated. In Section 3, we describe the data acquisition process and provide a detailed analysis of the prediction results. Finally, we provide the discussion and concluding thoughts in Section 4 and Section 5, respectively.

2. Prediction of Rigid-Flex Printed Circuit Board’s Scattering Parameters

2.1. Rigid-Flex Printed Circuit Boards and Problems of Infrared InGaAs Detection Links

While collecting terrestrial object information from an optical system to an InGaAs detector, as a single pixel output is a fixed level value, the ideal readout circuit outputs an analog signal waveform of an amplitude-modulated square wave sequence. The amplitude of the latter reflects the frequency at which ground information can be obtained at that spatial resolution. Rigid-flex printed circuit boards combine rigid and flexible printed circuits. Due to the limitations on the position and size of infrared photodetectors in the structure, they cannot be directly integrated into the driver acquisition circuit board. Therefore, foldable rigid-flex printed circuit boards are used to complete the signal transmission between the infrared photodetectors and the information acquisition module, which results in a lightweight and compact device [37]. However, when the sensor’s output signal passes through the rigid-flex board, before reaching the information acquisition module, the ideal square wave sequence experiences distortion and produces unexpected noise, as shown in Figure 1.
Circuit designers often opt for solid copper as the reference ground plane for rigid boards. Although solid copper provides superior shielding, using cross-hatched copper in flexible areas provides better bending and flexibility to the ground plane, ultimately reducing the probability of circuit damage and fracture. In addition, the cross-hatched copper results in better heat dissipation, solving the problem of uneven heating in solid copper, which causes PCB warpage. Consequently, cross-hatched copper is widely used in rigid-flex printed PCB manufacturing. In conventional cross-hatched reference planes, two types of grids are primarily used, i.e., 45° uniform oblique and vertical crossed grids, as illustrated in Figure 2. However, there is currently a lack of research on grid size, proportion, and form optimization, which often leads to poor shielding of the produced grid reference surface, which in turn can result in poor signal integrity and electromagnetic magnetic compatibility (EMC). The grid design affects the wiring’s parasitic capacitance and its projection area on the grid reference plane, where the latter also affects the former. In high-frequency circuits, it is essential to maintain the parasitic capacitance values as low as possible [38], as even small variations in this parameter can cause a significant frequency shift [39].

2.2. BP Neural Network

A BP neural network is an intelligent algorithmic model that simulates the thought process of the human brain. Learning inherent mappings from a large amount of data provides a powerful tool for solving complex nonlinear problems [40]. A BP neural network is a multi-layer feedforward network trained by an error backpropagation algorithm, and it is one of the most widely used and successful neural network models [41].

2.3. Modeling of Rigid-Flex PCB S-Parameters Based on BP Neural Network

Considering the mechanical structure and impedance matching problem of rigid-flex boards, the line width, thickness, length, and relative dielectric constant of the transmission line should be fixed. The main variation was in the design of the copper grid. At present, the design of the grid size and ratio of flexible boards is based on experience and lacks strict theoretical guidance, which may lead to more serious signal integrity and EMC problems. The above problem can further improve the signal transmission rate of rigid-flex circuit boards. To optimize the grid structure at the beginning of the design and improve the reliability of the device, this study analyzed the S-parameters based on a BP neural network to explore the effect of the grid copper design on signal integrity.

2.3.1. Network Parameter Settings

The BP neural network architecture employed in this study comprised an input layer, an output layer, and two hidden layers. The dimension of the input layer was [1, 6], corresponding to the optimization variables for the grid ground plane. The dimension of the output layer was [1, 201], matching the number of frequency points. The S-parameters in the frequency domain were represented by a complex-valued function depending on frequency. To simplify the output, we decomposed it into real and imaginary parts, on which the network was trained separately to produce fitted curves at the output layer corresponding to the two parts at 201 frequency points. In Figure 3, the relationship between the settings for root mean square error (RMSE), time, and number of neurons in predicting the imaginary part of S21 is depicted. The horizontal coordinates are [number of neurons in the first hidden layer, number of neurons in the second hidden layer]. As shown in the figure, when the number of neurons increases, the RMSE does not decrease significantly, while the time cost increases. Therefore, we determined that utilizing two hidden layers with two neurons each provides an effective trade-off between model efficiency and accuracy. The network structure allows it to capture more complex nonlinear relationships and intricate nuances within the data, enhancing its ability to model the task while maintaining feasible resource requirements accurately.
Furthermore, the fitting effect of this configuration was deemed excellent. The loss function was set as the mean squared error (MSE), and the threshold was set to 0.001. The maximum number of training steps was 1000, and activation functions f1 and f2 were set as Tan-Sigmoid and Purelin functions, respectively. In Figure 4, the specific structure of the network prediction model is illustrated.

2.3.2. Model Training

Several training algorithms are available for BP neural networks. Among the classical ones are the gradient descent algorithm and the Levenberg–Marquardt (LM) algorithm. The former is easy to implement, but its training convergence is slow, and it may even stop training prematurely. The LM algorithm combines the advantages of the gradient descent and Gauss–Newton methods, and its convergence speed is significantly higher than that of the gradient descent method while maintaining high accuracy. Therefore, the LM algorithm was adopted in this study to train the model. The error function is expressed as follows:
E ( w ) = 1 2 i = 1 n ( x i y i ) 2 = 1 2 i = 1 n e i 2 ( w )
where n represents the number of samples, w represents the weight vector, and e i 2 ( w ) denotes the error value of the ith node in the output layer. The weight vector at the (m + 1)th iteration can be written as follows:
w m + 1 = w m + Δ w
where the weight increment is given as:
Δ w = F T ( w ) F ( w ) + μ I 1 F T ( w ) e ( w )
The Jacobian matrix of the weights is calculated as follows:
F ( w ) = e 1 ( w ) w 1 e 1 ( w ) w 2 e 1 ( w ) w n e 2 ( w ) w 1 e 2 ( w ) w 2 e 2 ( w ) w n e N ( w ) w 1 e N ( w ) w 2 e N ( w ) w n
where w m stands for the weight vector at the mth iteration, e ( w ) stands for the error vector, I stands for the identity matrix, and e N ( w ) w n stands for the partial derivative of the error vector over the weight vector.

3. Simulation

3.1. Data Acquisition and Analysis

An automated simulation platform was developed by using PowerSI (Cadence, San Jose, CA, USA) to obtain the S-parameter dataset. In this study, the ANN training was performed by using Matlab R2023b (MathWorks, Natick, MA, USA). The optimization parameters of the reference grid ground plane included line width, line spacing, and angle, as depicted in Figure 5. All simulations were performed on a personal computer equipped with an Intel Core i7-1360P (Intel Corporation, Santa Clara, CA, USA) CPU operating at 2.2 GHz with 32 GB memory.
When designing the reference ground grid copper as a 45° uniform oblique grid, θ1 and θ2 were set to 45° and 135°, respectively. Grid line width w1 was initially set to 6 mil and increased stepwise by 1 mil up to 15 mil, resulting in 10 different parameter settings grouped together. Moreover, w2, s1, and s2 were kept consistent within each group. However, when the reference ground grid copper was designed as a vertical crossed grid, the grid copper angles were altered to 90° and 180°. Spacings s1 and s2 and line widths w1 and w2 were set to be consistent with the 45° uniform oblique grid design. The specific parameter settings are listed in Table 1.
To verify the universality and robustness of the model proposed in this study for predicting the S-parameters of transmission lines in rigid-flex PCBs for aerospace payloads, a single-ended transmission line (A) and a differential transmission line pair (B) in the rigid-flex PCBs were selected respectively for simulation validation.
The main parameter settings of Line A were as follows: the line width was 5.1 mil, the copper thickness was 0.708 mil, the transmission length was 7874 mil, the relative dielectric constant was 3.3, and the dielectric thickness was 2.953 mil. The differential transmission line pair (B) consisted of two Lines A with a 12 mil spacing, as shown in Figure 6.
The settings for all the simulations were consistent and included the following: the simulation frequency range was 1 kHz~2 GHz, and the linear scanning mode step size was 10 MHz; the frequency was divided into 201 points.
Finally, 480 different sets of S-parameter data were obtained with simulations for a single-ended transmission line (A) and a differential transmission line pair (B), respectively, as shown in Figure 7. To better verify the generalization ability of the proposed method and perform an unbiased evaluation, the samples obtained from the last 80 design parameters in Table 1 were selected as fixed test samples and were not used in model training. The first 400 samples in Table 1, each with an input dimension of 6 corresponding to the design parameters and an output dimension of 201 frequency points, were randomly divided into a test set and a validation set according to a 4:1 ratio, and a total of five experiments were conducted.

3.2. Prediction Results and Analysis

3.2.1. Overfitting, Robustness and Generalization Ability Analysis

Overfitting is a common issue in neural network training that leads to poor generalization performance on new data [42]. To investigate the susceptibility to overfitting of the proposed model, we reserved 20% of the samples in the training sets of the single-ended transmission line and differential transmission lines for validation data checking, respectively. To assess the generalization ability and robustness of the proposed model, five experiments were conducted with different training set compositions. As illustrated in Figure 8, despite variations in the training data, the mean MSEs on the training and validation sets converged to similar values and stabilized within each experiment. Notably, model performance improved, while the MSE curve on the training set did not further decrease, indicating that overfitting was effectively avoided. These results demonstrate that the proposed approach exhibits strong generalization capability and robustness, without being susceptible to overfitting issues.

3.2.2. Comparison of Goodness-of-Fit between BP Model and GRNN Model

We selected the same sample set from the test set and predicted the S-parameters by using the BP and GRNN models; the prediction curves are compared in Figure 9, Figure 10, Figure 11 and Figure 12. It can be clearly seen that the S-parameters predicted by using the BP model were closer to the actual values than those predicted by the GRNN model in both the real and imaginary parts.
The amplitude and phase can also be calculated from the real and imaginary parts, which can help us understand the physical meaning of the S-parameters more intuitively. According to the above comparison results, we used the proposed BP model for curve fitting. A test sample was randomly selected, and Figure 13a–d illustrates the differences between the actual and predicted values of the amplitude and phase for S11 and S21, respectively. It can be seen from the four sub-figures that the predicted values of the amplitude and phase fit well with the actual values, which demonstrates the superior prediction performance of the proposed method.
Comparing curves from a single test sample is not sufficient to demonstrate the accuracy and robustness of a method. For nonlinear regression models, the goodness-of-fit is a critical evaluative metric that quantifies how closely the model aligns with the observed data. Greater goodness-of-fit values indicate that the regression function more accurately describes the data.
A common approach to assessing the goodness-of-fit is calculating the adjusted R2 [43], which accounts for model complexity and degrees of freedom to mitigate overfitting. This metric ranges from 0 to 1, where values nearing 1 indicate that the model adequately explains variation in the data, while values approaching 0 denote poor explanatory power. Moreover, the goodness-of-fit analysis was expanded to include the RMSE, providing additional insights into model performance from a different perspective. While the adjusted R2 quantifies the variance explained by the model, the RMSE measures residual differences between predicted and actual values. Evaluating nonlinear regression models by using the adjusted R2 and RMSE provides a complementary quantification of model fit quality.
Therefore, we compared the prediction performance of the proposed BP model and the GRNN model by evaluating the goodness-of-fit metrics in all 80 test samples. In Figure 14 and Figure 15, we show the comparison of averaging the RMSE and the adjusted R2, respectively, of the real and imaginary parts of the 80 identical test sets in one experiment for the single-ended transmission line and the differential line pair. It can be seen that our method outperformed the GRNN model in 80 samples, as evidenced by smaller RMSEs and adjusted R2 values closer to 1.
To further demonstrate the robustness and high fit of the proposed method, five experiments were conducted as described in Section 3.1. The quantitative analysis summarized in Table 2 and Table 3 reinforced the superior goodness-of-fit of the BP model, as the average adjusted R2 was closer to 1 in the five experiments on 80 test sets.

3.2.3. Time Efficiency

The time periods required to obtain 40 and 80 S-parameter curves by using PowerSI software simulation were 200 min and 400 min, respectively. However, by training the two-port network model (the real and imaginary parts of S11, S21, S12, and S22) once with the BP model and predicting the 40 and 80 S-parameter curves, these were reduced to 51.3 s and 51.7 s, respectively. The GRNN model took 68.2 s and 68.7 s, respectively, to do so. It is evident that obtaining S-parameters by using the neural network-trained model significantly reduces the time cost of evaluating the S-parameters, as the prediction time was almost negligible even as the number of S-parameter curves increased, while the time necessary when using traditional simulation software increased linearly. Furthermore, the training and prediction efficiency of the BP model was higher than that of the GRNN model. Therefore, our method can improve the efficiency of iteratively optimized rigid-flex PCBs.

3.2.4. Errors in Eye Height and Eye Width

To further validate the modeling accuracy, time–domain transient simulations were performed for a single-ended transmission line by using the S-parameters predicted from the proposed BP model, the GRNN model [36], and the corresponding test samples, respectively. A pseudorandom binary sequence (PRBS) with a code length of 255 was applied as the excitation source at 1 Gbps with 0.25 ns rise/fall time. The transient simulations were implemented in ADS, as shown in Figure 16, to obtain the simulated eye diagrams. A total of 80 test samples of predicted and actual S-parameters were used in the simulations. Eye height and eye width errors were extracted from the simulated eye diagrams and analyzed. Employing time–domain transient simulations allowed us to quantitatively analyze the modeling accuracy with industry-standard statistical eye diagram parameters.
The results demonstrate excellent concordance between the eye diagrams predicted by using the proposed model and those obtained with simulations. Notably, the errors associated with the eye width and eye height parameters were 0.44% and 0.78%, respectively. These error values represent a significant reduction, being merely 0.68 times and 0.35 times the corresponding errors obtained with the GRNN model. In Figure 17, the trends of eye height and eye width errors for 80 sets of test samples in one experiment are reported. In Table 4, the simulated mean, predicted mean, and error mean of eye height and eye width in the experiment are listed.

4. Discussion

Considering the significance of rigid-flex PCBs in InGaAs detection links for weak infrared detection signals, we investigated the use of BP neural networks to rapidly predict the S-parameters of the grid ground plane of rigid-flex PCBs. We evaluated the results by using robustness, goodness-of-fit, generalization ability, time efficiency, and eye width and eye height errors as performance metrics. The results indicate that the proposed BP neural network can accurately reconstruct the nonlinear behavior of the grid ground plane design for rigid-flex PCBs over the entire frequency range. The constructed neural network model can effectively replace simulations in obtaining S-parameters with good prediction performance.
As presented in Section 3, in-depth research was conducted in order to demonstrate the choice of the model and its accuracy and time efficiency compared with other prediction models. The GRNN method has been used to predict the HBT scattering parameters at different temperatures [36]. However, unlike the approach in [36], our method can predict values within a certain frequency range instead of a single frequency value. From the perspective of the goodness-of-fit, the RMSE of the BP model was smaller (see Table 2), and the adjusted R2 was maximally 0.06 larger than that of the GRNN model (see Table 3). It is evident from the results that using an ANN to predict the S-parameters of the differential transmission line pair was more effective than using it to predict those of the single-ended transmission line. Based on the results of the eye width and eye height errors, the BP model improved them by 0.21% and 1.44%, respectively (see Table 4). In terms of time efficiency, the BP model also showed faster performance.
In addition, we demonstrated the potential of our neural network for designing the grid ground plane by comparing it with a traditional 3D field solver. This confirms the effectiveness of neural networks in electromagnetic modeling and parameter prediction, as previously shown in other studies [29,30,31,32,33,34,35,36]. Given its predictive capability, the model represents a powerful tool for system-level simulation and optimization, eliminating the need for tedious and repetitive simulation processes and saving significant time and computational resources. This result meets the current industrial demand for efficient design optimization methods and provides new ideas for future research in the field of rigid-flex PCBs and InGaAs detector applications. However, despite the results obtained in this study, there is still room for improvement. In future work, we will focus on using the higher-efficiency model to solve more complex problems, thereby broadening its applicability.

5. Conclusions

In the aerospace field, large-array InGaAs focal plane detectors for quantitative payload applications have experienced rapid development. Consequently, there is an increased demand for precise noise control and quantization of high-speed weak signals. Detailed analysis of each detection system component is necessary to achieve this goal. In this study, we designed a BP neural network to predict the scattering parameters of rigid-flex PCBs for InGaAs infrared detection links.
The simulation verification and analysis results demonstrate that the BP network model can achieve high-precision scattering parameter prediction. Moreover, the efficiency of this model was found to be 463 times higher than that of a traditional 3D field solver when 80 S-parameters were obtained. It also exhibited 25% faster performance and 30% and 65% smaller eye errors than the GRNN model.
This method effectively optimizes the grid copper design in early rigid-flex PCB layout stages, substantially reducing the design cycle time. It also improves design quality for aerospace quantitative payloads. The findings of this study represent a foundation for future infrared detection systems requiring high sensitivity and low noise. They can also provide technical support to expand the application prospects of infrared weak-spectrum detection. It is recommended that other influencing factors be considered to improve prediction accuracy in future research. Moreover, extending this modeling approach to other vital electrical parameters would enable more comprehensive PCB evaluation and optimization.

Author Contributions

Conceptualization, X.C. and J.M.; methodology, J.M. and H.Y.; software, J.M. and H.Y.; validation, J.M., H.Y. and X.G.; formal analysis, J.M. and H.Y.; investigation, J.M.; resources, X.C. and L.D.; data curation, H.Y., J.M. and X.G.; writing—original draft preparation, J.M.; writing—review and editing, X.C. and L.D.; supervision, L.D. and X.C.; funding acquisition, L.D. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Key Research and Development Program of China (grant number 2022FYB3904803).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data generated during the study are contained within the article.

Acknowledgments

The authors are grateful to the editors and anonymous reviewers for their helpful comments.

Conflicts of Interest

The authors declare no conflicts of interest.

Abbreviations

The following abbreviations are used in this manuscript:
3DThree-dimensional
ANNArtificial Neural Network
BPBackpropagation
EDAElectronic Design Automation
EMCElectromagnetic Magnetic Compatibility
EMIElectromagnetic Interference
GRNNGeneralized Regression Neural Network
LMLevenberg–Marquardt
MLMachine Learning
MSEMean Squared Error
RMSERoot Mean Squared Error
PCBRigid-flex Printed Circuit Board
PRBSPseudo-random Binary Sequence
S-parametersScattering Parameters
SVMSupport Vector Machine
VFVector Fitting

References

  1. Mu, Y.; Zhang, Z.; Chen, C.; Yuan, D.; Wang, J.; Gao, H.; Chi, Y. The Design of a Low-Noise, High-Speed Readout-Integrated Circuit for Infrared Focal Plane Arrays. Sensors 2023, 23, 8715. [Google Scholar] [CrossRef] [PubMed]
  2. Ma, Y.; Wu, H.; Wang, L.; Huang, B.; Ranjan, R.; Zomaya, A.Y.; Jie, W. Remote Sensing Big Data Computing: Challenges and Opportunities. Future Gener. Comput. Syst. 2015, 51, 47–60. [Google Scholar] [CrossRef]
  3. Wu, T.; Wu, F.; Qiu, C.; Redouté, J.-M.; Yuce, M.R. A Rigid-Flex Wearable Health Monitoring Sensor Patch for IoT-Connected Healthcare Applications. IEEE Internet Things J. 2020, 7, 6932–6945. [Google Scholar] [CrossRef]
  4. Cao, F.; Liu, L.; Li, L. Short-Wave Infrared Photodetector. Mater. Today 2023, 62, 327–349. [Google Scholar] [CrossRef]
  5. Ngamakeur, K.; Yongchareon, S.; Yu, J.; Sheng, Q.Z. Deep CNN-LSTM Network for Indoor Location Estimation Using Analog Signals of Passive Infrared Sensors. IEEE Internet Things J. 2022, 9, 22582–22594. [Google Scholar] [CrossRef]
  6. Pan, G.; Zhu, X.; Gilbert, B.K. Analysis of Transmission Lines of Finite Thickness above a Periodically Perforated Ground Plane at Oblique Orientations. IEEE Trans. Microw. Theory Tech. 1995, 43, 383–393. [Google Scholar] [CrossRef]
  7. Wu, S.; Shi, H.; Herndon, M.; Cornelius, B.; Halligan, M.; Fan, J. Modeling and Analysis of a Trace Referenced to a Meshed Ground Plane. In Proceedings of the 2011 IEEE International Symposium on Electromagnetic Compatibility, Long Beach, CA, USA, 14–19 August 2011; pp. 137–141. [Google Scholar] [CrossRef]
  8. Cicchetti, R.; Faraone, A. A Full-Wave Radiation Model for a Class of Gridded Ground Interconnecting Structures. IEEE Trans. Antennas Propag. 1999, 47, 212–213. [Google Scholar] [CrossRef]
  9. Yu, S.; Kam, D.-G. Improving Signal Integrity by Reducing Mode Conversion in Differential Pairs on a Meshed Ground in Flexible Printed Circuit Boards. J. Electromagn. Eng. Sci. 2022, 22, 229–235. [Google Scholar] [CrossRef]
  10. Xiao, F.; Murano, K.; Kami, Y. Modeling of Differential Line Referenced to a Meshed Ground Plane. In Proceedings of the 2014 International Symposium on Electromagnetic Compatibility, Gothenburg, Sweden, 1–4 September 2014; pp. 735–740. [Google Scholar] [CrossRef]
  11. Xiao, F.; Nakada, Y.; Murano, K.; Kami, Y. Crosstalk Analysis Model for Traces Crossing Split Ground Plane and Its Reduction by Stitching Capacitor. Electron. Commun. Jpn. 2007, 90, 26–34. [Google Scholar] [CrossRef]
  12. Bogatin, E. Signal and Power Integrity—Simplified; Prentice Hall International: Upper Saddle River, NJ, USA, 2009. [Google Scholar]
  13. Cho, K.; Kim, Y.; Lee, D.; Song, J.; Park, J.; Lee, S.; Kim, S.; Park, G.; Son, K.; Kim, J. Signal Integrity Design and Analysis of Differential High-Speed Serial Links in Silicon Interposer with Through-Silicon Via. IEEE Trans. Compon. Packag. Manuf. Technol. 2019, 9, 107–121. [Google Scholar] [CrossRef]
  14. Cho, K.; Kim, Y.; Lee, D.; Kim, H.; Choi, S.; Song, J.; Kim, S.; Park, J.; Lee, S.; Kim, J. Signal Integrity Design and Analysis of Silicon Interposer for GPU-Memory Channels in High-Bandwidth Memory Interface. IEEE Trans. Compon. Packag. Manuf. Technol. 2018, 8, 1658–1671. [Google Scholar] [CrossRef]
  15. Yu, S.; Kam, D.G. Offset Mesh in FPCB for Better Impedance Control. Electron. Lett. 2018, 54, 812–813. [Google Scholar] [CrossRef]
  16. Kuhlmann, K.; Gellersen, F.; Tschauder, M. Comparison of S-Parameter Measurement Methods for Adapters. Adv. Radio Sci. 2023, 20, 1–8. [Google Scholar] [CrossRef]
  17. Stokes, D.; Gellersen, F.; Allal, D.; Skinner, J.F.; Phung, G.N.; Kuhlmann, K. Traceable S-Parameter Measurements up to 90 GHz in 1.35 Mm Coaxial. Meas. Sci. Technol. 2023, 34, 064006. [Google Scholar] [CrossRef]
  18. Pham, T.D.; Allal, D.; Ziadé, F.; Bergeault, E. On-Wafer Coplanar Waveguide Standards for S-Parameter Measurements of Balanced Circuits up to 40 GHz. IEEE Trans. Instrum. Meas. 2019, 68, 2160–2167. [Google Scholar] [CrossRef]
  19. Zhao, W.; Cheng, C.; Li, L. Sensitivity Analysis of S-Parameter Measurements Due to Nonideal SOLR Calibration Standards. IEEE Trans. Instrum. Meas. 2021, 70, 1002910. [Google Scholar] [CrossRef]
  20. Gupta, K.C.; Garg, R.; Bahl, I.; Bhartia, P. Microstrip Lines and Slotlines; Artech House: Norwood, MA, USA, 1996. [Google Scholar]
  21. Mariani, E.A.; Heinzman, C.P.; Agrios, J.P.; Cohn, S.B. Slot Line Characteristics. IEEE Trans. Microw. Theory Tech. 1969, 17, 1091–1096. [Google Scholar] [CrossRef]
  22. Yang, S.L.L. 3D SiP Simulation Technology and Application. In Proceedings of the 2019 20th International Conference on Electronic Packaging Technology, Hong Kong, China, 12–15 August 2019; pp. 1–4. [Google Scholar] [CrossRef]
  23. Hackl, H.; Pommerenke, D.; Ibel, M.; Auinger, B. Extraction of Single Cell Impedance from Battery Module Measurement by Simulation-Based De-Embedding. IEEE Trans. Signal Power Integr. 2022, 1, 112–120. [Google Scholar] [CrossRef]
  24. Grivet-Talocia, S.; Gustavsen, B. Black-box macromodeling and its EMC applications. IEEE Electromagn. Compat. Mag. 2016, 5, 71–78. [Google Scholar] [CrossRef]
  25. Gustavsen, B.; Semlyen, A. Rational approximation of frequency domain responses by vector fitting. IEEE Trans. Power Deliv. 1999, 14, 1052–1061. [Google Scholar] [CrossRef]
  26. Leszczynska, N.; Couckuyt, I.; Dhaene, T.; Mrozowski, M. Low-Cost surrogate models for microwave filters. IEEE Microw. Wirel. Compon. Lett. 2016, 26, 969–971. [Google Scholar] [CrossRef]
  27. Kuznetsov, Y.V.; Baev, A.; Shevgunov, T.; Siart, U.; Yordanov, H.; Russen, P. Generation of network models for planar microwave circuits by system identification methods. In Proceedings of the IEEE 2009 International Conference on Electromagnetics in Advanced Applications, Torino, Italy, 14–18 September 2009; pp. 966–969. [Google Scholar] [CrossRef]
  28. Torun, H.M.; Durgun, A.C.; Aygün, K.; Swaminathan, M. Causal and Passive Parameterization of S-Parameters Using Neural Networks. IEEE Trans. Microw. Theory Tech. 2020, 68, 4290–4304. [Google Scholar] [CrossRef]
  29. Jin, J.; Feng, F.; Zhang, J.; Yan, S.; Na, W.; Zhang, Q. A Novel Deep Neural Network Topology for Parametric Modeling of Passive Microwave Components. IEEE Access 2020, 8, 82273–82285. [Google Scholar] [CrossRef]
  30. Zhang, Q.; Gupta, K.C.; Devabhaktuni, V. Artificial Neural Networks for Rf and Microwave Design-from Theory to Practice. IEEE Trans. Microw. Theory Tech. 2003, 51, 1339–1350. [Google Scholar] [CrossRef]
  31. Trinchero, R.; Larbi, M.; Torun, H.M.; Canavero, F.; Swaminathan, M. Machine Learning and Uncertainty Quantification for Surrogate Models of Integrated Devices with a Large Number of Parameters. IEEE Access 2019, 7, 4056–4066. [Google Scholar] [CrossRef]
  32. Trinchero, R.; Manfredi, P.; Stievano, I.S.; Canavero, F. Machine Learning for the Performance Assessment of High-Speed Links. IEEE Trans. Electromagn. Compat. 2018, 60, 1627–1634. [Google Scholar] [CrossRef]
  33. Garbuglia, F.; Reuschel, T.; Schuster, C.; Deschrijver, D.; Dhaene, T.; Spina, D. Modeling Electrically Long Interconnects Using Physics-Informed Delayed Gaussian Processes. IEEE Trans. Electromagn. Compat. 2023, 65, 1715–1723. [Google Scholar] [CrossRef]
  34. Swaminathan, M.; Torun, H.M.; Yu, H.; Hejase, J.A.; Becker, W.D. Demystifying Machine Learning for Signal and Power Integrity Problems in Packaging. IEEE Trans. Compon. Packag. Manuf. Technol. 2020, 10, 1276–1295. [Google Scholar] [CrossRef]
  35. Schierholz, M.; Erdin, I.; Balachandran, J.; Yang, C.; Schuster, C. Parametric S-Parameters for PCB Based Power Delivery Network Design Using Machine Learning. In Proceedings of the IEEE 26th Workshop Signal Power Integrity, Siegen, Germany, 22–25 May 2022; pp. 1–4. [Google Scholar] [CrossRef]
  36. Lin, Q.; Wang, X.; Wu, H. GRNN-Based Scattering Parameter modeling investigation for HBT at different temperature. IEEE Access 2023, 11, 38002–38009. [Google Scholar] [CrossRef]
  37. Gagnon-Turcotte, G.; Kisomi, A.A.; Ameli, R.; Camaro, C.-O.D.; LeChasseur, Y.; Néron, J.-L.; Bareil, P.B.; Fortier, P.; Bories, C.; De Koninck, Y.; et al. A Wireless Optogenetic Headstage with Multichannel Electrophysiological Recording Capability. Sensors 2015, 15, 22776–22797. [Google Scholar] [CrossRef]
  38. Mulcahy, J.; Shi, S.; Peters, F.H.; Dai, X. Minimisation of Parasitic Capacitance in Lumped-Element Electro-Absorption Modulators for High-Speed Optical Components. Photonics 2023, 10, 885. [Google Scholar] [CrossRef]
  39. Qu, Z.; Zhu, Z.; Liu, Y.; Yu, M.; Ye, T.T. Parasitic Capacitance Modeling and Measurements of Conductive Yarns for E-Textile Devices. Nat. Commun. 2023, 14, 2785. [Google Scholar] [CrossRef]
  40. Wang, Z.; Wu, J.; Wang, H.; Wang, H.; Hao, Y. Optimal Underwater Acoustic Warfare Strategy Based on a Three-Layer GA-BP Neural Network. Sensors 2022, 22, 9701. [Google Scholar] [CrossRef]
  41. Gao, X.; Mou, J.; Banerjee, S.; Zhang, Y. Color-Gray Multi-Image Hybrid Compression–Encryption Scheme Based on BP Neural Network and Knight Tour. IEEE Trans. Cybern. 2023, 53, 5037–5047. [Google Scholar] [CrossRef]
  42. Chu, J.; Li, X.; Zhang, Z.; Zhang, Y.; He, M. A novel method overcomeing overfitting of artificial neural network for accurate prediction: Application on thermophysical property of natural gas. Case Stud. Therm. Eng. 2021, 28, 101406. [Google Scholar] [CrossRef]
  43. Shi, J.; Wang, C.; Xi, X.; Yang, X.; Wang, J.; Ding, X. Retrieving fPAR of maize canopy using artificial neural networks with airborne LiDAR and hyperspectral data. Remote Sens. Lett. 2020, 11, 1002–1011. [Google Scholar] [CrossRef]
Figure 1. Infrared InGaAs detection link.
Figure 1. Infrared InGaAs detection link.
Sensors 24 02221 g001
Figure 2. Two types of grid ground planes. (a) 45° uniform oblique grid. (b) Vertical crossed grid.
Figure 2. Two types of grid ground planes. (a) 45° uniform oblique grid. (b) Vertical crossed grid.
Sensors 24 02221 g002
Figure 3. Number of neurons vs. RMSE and time.
Figure 3. Number of neurons vs. RMSE and time.
Sensors 24 02221 g003
Figure 4. BP neural network model architecture.
Figure 4. BP neural network model architecture.
Sensors 24 02221 g004
Figure 5. Line widths, spacings, angles of the reference ground plane copper.
Figure 5. Line widths, spacings, angles of the reference ground plane copper.
Sensors 24 02221 g005
Figure 6. (a) Topology of Line A and Line pair B. (b) Topology of Line A and Line pair B on a rigid–flex PCBs.
Figure 6. (a) Topology of Line A and Line pair B. (b) Topology of Line A and Line pair B on a rigid–flex PCBs.
Sensors 24 02221 g006
Figure 7. Reflection coefficients for 480 sets of different grid parameters. (a) The S11 of Single-ended transmission line A. (b) The SDD11 of differential transmission line pair B.
Figure 7. Reflection coefficients for 480 sets of different grid parameters. (a) The S11 of Single-ended transmission line A. (b) The SDD11 of differential transmission line pair B.
Sensors 24 02221 g007
Figure 8. Performance of training and validation sets for 5 experiments. (a) Single-ended transmission line A. (b) Differential transmission line pair B.
Figure 8. Performance of training and validation sets for 5 experiments. (a) Single-ended transmission line A. (b) Differential transmission line pair B.
Sensors 24 02221 g008
Figure 9. Comparison curves of real part predictions in 80 test sample sets. (a) Single-ended transmission line A’s S21. (b) Differential transmission line pair B’s SDD21.
Figure 9. Comparison curves of real part predictions in 80 test sample sets. (a) Single-ended transmission line A’s S21. (b) Differential transmission line pair B’s SDD21.
Sensors 24 02221 g009
Figure 10. Comparison curves of imaginary part predictions in 80 test sample sets. (a) Single-ended transmission line A’s S21. (b) Differential transmission line pair B’s SDD21.
Figure 10. Comparison curves of imaginary part predictions in 80 test sample sets. (a) Single-ended transmission line A’s S21. (b) Differential transmission line pair B’s SDD21.
Sensors 24 02221 g010
Figure 11. Comparison curves of real part predictions in 80 test sample sets. (a) Single-ended transmission line A’s S11. (b) Differential transmission line pair B’s SDD11.
Figure 11. Comparison curves of real part predictions in 80 test sample sets. (a) Single-ended transmission line A’s S11. (b) Differential transmission line pair B’s SDD11.
Sensors 24 02221 g011
Figure 12. Comparison curves of imaginary part predictions in 80 test sample sets. (a) Single-ended transmission line A’s S11. (b) Differential transmission line pair B’s SDD11.
Figure 12. Comparison curves of imaginary part predictions in 80 test sample sets. (a) Single-ended transmission line A’s S11. (b) Differential transmission line pair B’s SDD11.
Sensors 24 02221 g012
Figure 13. (a) Comparison results of amplitude prediction for S11. (b) Comparison results of phase prediction for S11. (c) Comparison results of amplitude prediction for S21. (d) Comparison results of phase prediction for S21.
Figure 13. (a) Comparison results of amplitude prediction for S11. (b) Comparison results of phase prediction for S11. (c) Comparison results of amplitude prediction for S21. (d) Comparison results of phase prediction for S21.
Sensors 24 02221 g013
Figure 14. The RMSE of the real and imaginary parts of the same 80 test sample sets. (a) Single-ended transmission line A. (b) Differential transmission line pair B.
Figure 14. The RMSE of the real and imaginary parts of the same 80 test sample sets. (a) Single-ended transmission line A. (b) Differential transmission line pair B.
Sensors 24 02221 g014
Figure 15. The adjusted R2 of the real and imaginary parts of the same 80 test sample sets. (a) Single-ended transmission line A. (b) Differential transmission line pair B.
Figure 15. The adjusted R2 of the real and imaginary parts of the same 80 test sample sets. (a) Single-ended transmission line A. (b) Differential transmission line pair B.
Sensors 24 02221 g015
Figure 16. Transient simulation is built in the ADS.
Figure 16. Transient simulation is built in the ADS.
Sensors 24 02221 g016
Figure 17. Comparison of error trends by the BP and GRNN methods in 80 test sample sets. (a) Eye width. (b) Eye height.
Figure 17. Comparison of error trends by the BP and GRNN methods in 80 test sample sets. (a) Eye width. (b) Eye height.
Sensors 24 02221 g017
Table 1. Settings of mesh parameters.
Table 1. Settings of mesh parameters.
Group Numberw1
(mil)
w2
(mil)
s1
(mil)
s2
(mil)
θ1
(°)
θ2
(°)
16~15661545135
26~15771445135
106~151515645135
116~156151545135
126~157141445135
206~15156645135
216~15661590180
226~15771490180
306~151515690180
316~156151590180
326~157141490180
406~15156690180
416~15681445135
426~158101245135
436~1510121045135
446~151214845135
456~15681490180
466~158101290180
476~1510121090180
486~151214890180
Table 2. Average RMSE over five experiments.
Table 2. Average RMSE over five experiments.
BP GRNN
Real parts of S110.01070.0168
Imaginary parts of S110.01090.0165
Real parts of S210.02100.0214
Imaginary parts of S210.02190.0205
Real parts of SDD110.00420.0086
Imaginary parts of SDD110.00410.0090
Real parts of SDD210.01060.0130
Imaginary parts of SDD210.01050.0138
Table 3. Average adjusted R2 over five experiments.
Table 3. Average adjusted R2 over five experiments.
BP GRNN
Real parts of S110.99610.9899
Imaginary parts of S110.99580.9899
Real parts of S210.99680.9928
Imaginary parts of S210.99740.9932
Real parts of SDD110.99870.9932
Imaginary parts of SDD110.99870.9926
Real parts of SDD210.99960.9971
Imaginary parts of SDD210.99970.9972
Table 4. Comparison of eye width and eye height obtained by different methods.
Table 4. Comparison of eye width and eye height obtained by different methods.
Error of Eye Width
(%)
Error of Eye Height
(%)
BP method0.440.78
GRNN method0.652.22
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

Mei, J.; Yuan, H.; Guo, X.; Chu, X.; Ding, L. Fast Prediction Method for Scattering Parameters of Rigid-Flex PCBs Based on ANN. Sensors 2024, 24, 2221. https://doi.org/10.3390/s24072221

AMA Style

Mei J, Yuan H, Guo X, Chu X, Ding L. Fast Prediction Method for Scattering Parameters of Rigid-Flex PCBs Based on ANN. Sensors. 2024; 24(7):2221. https://doi.org/10.3390/s24072221

Chicago/Turabian Style

Mei, Jingling, Haiyue Yuan, Xinxin Guo, Xiuqin Chu, and Lei Ding. 2024. "Fast Prediction Method for Scattering Parameters of Rigid-Flex PCBs Based on ANN" Sensors 24, no. 7: 2221. https://doi.org/10.3390/s24072221

APA Style

Mei, J., Yuan, H., Guo, X., Chu, X., & Ding, L. (2024). Fast Prediction Method for Scattering Parameters of Rigid-Flex PCBs Based on ANN. Sensors, 24(7), 2221. https://doi.org/10.3390/s24072221

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop