In this study, to solve the hidden limitation of conventional weighted fractional Fourier transform (WFRFT), a random modulation order parameter pool is established by applying chaos technology. Further, a WFRFT secure communication method based on the chaotic parameter pool (CPP) is proposed. Based on the effective characteristics of tent mapping and the sequence output range, the parameter pool constructor is established by parameter transformation. Furthermore, for each parameter selection period, the information can be processed by WFRFT using different modulation orders. The modeling and simulation demonstrate that this method can significantly increase the bit error rate and processing time of unauthorized receivers. This indicates that it can greatly increase the scanning difficulty of unauthorized users and improve the concealment and security of the original information transmission.
National Natural Science Foundation of China61501309China Postdoctoral Science Foundation2017T100185Natural Science Foundation of Liaoning Province2017011002-301Liaoning Provincial Colleges and Universities Innovative Talents Support Program1. Introduction
Fourier transform [1] has an inherent defect: it provides signal processing and analysis in the frequency or time domain based on a single frequency or time domain signal, which restricts increases in system capacity and spectral efficiency. Fractional Fourier transform (FRFT) [2, 3] is a new signal processing and analysis method in the transform domain of wireless communication systems. As a special branch of traditional Fourier transform, FRFT not only exhibits the characteristic of time-frequency conversion analysis of the traditional Fourier transform but also develops a unique joint signal processing ability in the time-frequency domain. FRFT can transform signals into each other in the time domain, frequency domain, and any transform domain in the transition region. Therefore, compared to the traditional Fourier transform, FRFT has gradually gained widespread attention and recognition in the field of nonstationary signal processing. In general, FRFT includes Chirp FRFT (CFRFT) [4, 5] and weighted FRFT (WFRFT) [6, 7]. CFRFT has been used extensively in radar signal processing, optical image processing, quantum mechanics, ranging, and other fields. However, owing to its discrete algorithm implementation, the application of CFRFT in the communication field has been severely limited. Compared to CFRFT, the WFRFT theory appeared at a later stage. Its signal is a form of the time-frequency domain signal fusion. The weighted fraction domain is in the middle of the traditional time and frequency domains. According to the characteristics of WFRFT, new physical characteristics are derived based on retaining the respective advantages of the traditional time-frequency domain signal. Owing to its simple implementation, and the low complexity of the discrete algorithm, WFRFT has been widely used in information optics, image encryption, and other engineering fields. Moreover, it has gradually been extended to image encryption, signal fractional domain sampling and reconstruction, and communication signal processing.
Because of the unique physical characteristics of the WFRFT signal, its application in secure communication has also become a topic of interest. In [8], a novel combinatory strategy was presented to embed WFRFT precoding signals into a transform domain communication system. This system could achieve an improved anti-eavesdropping capacity while maintaining its own communication quality with regard to bit error rate (BER) performances. In [9], a parallel-combinatory-spreading-based WFRFT system aimed at guaranteeing physical layer (PHY) security was proposed. In [10], a novel user cooperation scheme based on WFRFT was proposed to enhance the PHY security of wireless transmissions against eavesdropping. In [11], a novel transmission scheme based on constellation rotation and WFRFT was proposed to enhance the PHY security in polarization modulation-based dual-polarized satellite communications. In [12], a WFRFT-based cooperative overlay system, aiming at guaranteeing PHY security, was presented. In the proposed system, WFRFT was first performed on the secret data, such that the transmitted signal was distorted and could only be neutralized by inverse WFRFT with the same parameter. Thereafter, two streams of the transformed sequences bearing different messages were cooperatively and simultaneously transmitted to two legitimate receivers by means of a beamforming-like method.
In WFRFT communications, although demodulation parameters cannot be obtained for nonscanning methods, and the purpose of confidentiality can be achieved, demodulation parameters may be obtained by unauthorized receivers with scanning capability. As a result, the confidentiality performance of unauthorized receivers with a rapid scanning ability is also threatened. Although the modulation order can be extended, from one to many, for example, in multiple parameters WFRFT [13–15], or the state function can be extended from four to many, as in MWFRFT [16, 17], this will also increase the system complexity. Thus, to solve the hidden limitation of conventional WFRFT, chaos technology [18, 19] is introduced, following which a random modulation order parameter pool is established, so that a WFRFT secure communication method based on the chaotic parameter pool (CPP) is proposed.
2. CPP Method
Based on chaos and WFRFT, the innovation of the CPP method is the establishment of the association and control rules. Based on the effective characteristics of tent mapping and the sequence output range, the parameter pool constructor is established by parameter transformation. Furthermore, for each parameter selection period, the information can be processed by WFRFT using different modulation orders. The establishment of the association and control rules can greatly improve the security performance of the system.
Firstly, the generation trajectory of the chaotic signal is very complex and irregular resulting in superior concealment effect and high unpredictability. The generation of chaotic signals depends on the iteration equation, initial value, and the fractal parameters, which is conducive to the replication and regeneration of the signal. The specific form of chaos is reflected by bifurcation diagram, which reflects the traversal of chaotic sequence. The bifurcation diagram results of classical chaotic mappings are shown in Figure 1. It can be seen that the sequence has reached the best ergodicity state, and the range of fractal parameters can be obtained when each chaotic map reaches the state of full mapping. And tent mapping has the widest range of fractal parameters and the largest application space. As tent mapping presents the best chaotic state when α is greater than 0.4, the mapping can reach full mapping. And considering that the range of fractal parameter α is between 0 and 1, so the optimum range of α is between 0.4 and 1. In this paper, the parameter α is set to 0.4997.
The bifurcation diagram of chaos mapping. (a) Chebyshev mapping. (b) Improved logistic mapping. (c) Segmented mapping. (d) Tent mapping.
Thus, chaotic technology based on tent mapping is imported, and the CPP method is established, whose principle is illustrated in Figure 2.
Schematic of the CPP method.
Then, following substantial preliminary research and analysis, we believe that the Lyapunov exponents of tent mapping are relatively high; that is, they are more sensitive to the initial values, and the fractal parameters of the full mappings of tent maps have a wider range and larger application space. For this reason, a parameter pool constructor based on tent mapping is established.
The tent sequence is generated by the following iterating formula:(1)xm=1−xm−11−α,α<xm−1<1,xm−1α,0<xm−1≤α.
The range of the tent mapping initial value is from 0 to 1, so the sequence value xm ranges from 0 to 1.
Then, the parameters are transformed by using the sequence xm of the tent mapping output, as indicated in formula (2), so that the range of the parameters following transformation is 0,4, which accords with the periodic characteristics of the modulation order 0–4. The parameter ym is a floating-point value whose value is controlled by chaotic sequence xm, which is approximately random. So ym is much better than the integer parameter pool controlled by m/M sequence:(2)ym=4xm−1.
The parameter pool is formed by the parameter ym, and m∈1,N. The value range of m is determined by N, which is the number of switching times of the modulation order and also the number of times of information division to be transmitted. The calculation of N for a transmission cycle is carried out according to the following formula:(3)N=Rb⋅T,where Rb is the rate of information to be transmitted and T is the period of selecting parameters sequentially from the parameter pool.
At each time T, the parameters are selected sequentially from the parameter pool ym, and the selected parameters are used as the modulation order. Thereafter, the information to be transmitted is processed by WFRFT using this modulation order. In the next time T, the above operation is repeated, the next parameter is selected as the modulation order from the parameter pool ym sequentially, and the information to be transmitted by WFRFT is processed using the new modulation order. Thus, at the ith T time, the parameter yi is selected from ym and used as the modulation order, while the transmission information sn is processed by the yi-order WFRFT and is expressed as(4)Fyisn=ω0yisn+ω1yiSn+ω2yis−n+ω3yiS−n.
The yi-order WFRFT processing is denoted by Fyi⋅. In formula (4), the four “state functions” sn, Sn, s−n, and S−n are the results of 0, 1, 2, and 3 times of Fourier transformation of sn, respectively. The weighting coefficient ωlyi is defined as(5)ωlyi=14∑k=03exp2πj4l−yik,l=0,1,2,3.
In receiving, the baseband signal rn following the frequency downconversion, power amplifying, and band-pass filtering can be set to formula (6). Here, λ0n represents the mixed noise item caused by the transmission process and receiving pretreatment:(6)rn=Fyisn+λ0n.
The signal processing of WFRFT can be viewed as the process of rotating the input signal sn in the time-frequency plane, thereby realizing the signal energy redistribution in the time-frequency plane. Therefore, only when the receiver rotates the communication signal Fyisn at the same angle in the opposite direction can the receiver achieve signal energy aggregation. When the demodulation order is selected incorrectly, the signal energy will be lost, which will result in a loss of reception performance and reduce the receiver reception performance.
To synchronize the demodulation order of the inverse process with the modulation order of the transmitter, the receiver generates the sequence according to the same tent chaotic mapping and generates the same parameter pool ym according to formulas (1) and (2). Moreover, the parameters are selected sequentially from the parameter pool ym every time T, and at the ith T time, the selected parameter is yi′. Furthermore, the received signal is processed by WFRFT using yi′ as the demodulation order, as follows:(7)r′n=F−yi′rn,=F−yi′Fyisn+λ0n,=F−yi′Fyisn+F−yi′λ0n,=Fyi−yi′sn+λ0′n.
Here, λ0′n represents the mixed noise in the receiving system, which is related to the channel environment and receiver processing technology. In the following tests, the channel noise is Gauss noise.
When the parameters selected in the parameter pool are in the same order, and the switching time can be synchronized, yi′=yi can be obtained. Thereafter, it is introduced into formula (7) and the following formula (8) is obtained. It can be observed that, if the influence of noise is within a certain extent, the original information can be received and restored correctly:(8)r′n=F0sn+λ0′n,=ω00sn+ω10Sn+ω20s−n+ω30S−n+λ0′n,=1⋅sn+0⋅Sn+0⋅s−n+0⋅S−n+λ0′n,=sn+λ0′n.
For unauthorized users without scanning capabilities, the demodulation order yΔ of the WFRFT inverse transform of the receiver and modulation order yi of the transmitter cannot be exactly the same or synchronized, so the equivalent original information is processed by WFRFT with the modulation order yi−yΔ. Thus, the correct information sn cannot be received and restored correctly, which is expressed as(9)r′n=F−yΔrn,=F−yΔFyisn+λ0n,=F−yΔFyisn+F−yΔλ0n,=Fyi−yΔsn+λ0′n,≠sn.
However, for unauthorized receivers with scanning capabilities, when yΔ is equal to yi, the true transmission signal can be received and restored correctly. Therefore, to meet the requirements of the yΔ=yi equation under the condition that the switching rule is unknown, the aim of yi scanning is achieved by means of yΔ scanning in the main period of 0,4. The difference between the switching rule of the inverse transformation and positive transformation is Δy=yi−yΔ. Combined with the relationship between the scanning interval and BER, a smaller Δy results in a smaller bit error rate of the signal after inverse transform; that is, the error between r′n and the real signal is reduced. However, with a smaller Δy, more time is required for scanning in the 0–4 interval.
3. Performance Analysis
The BER performance of unauthorized receivers and that of authorized receivers is mainly affected by the Eb/N0. For unauthorized receivers, the signal processed by CPP increases the noise energy superimposed on the signal and reduces the signal-to-noise ratio of unauthorized receivers. Thus, Eb/N0 of the unauthorized receivers is also reduced due to the mathematical relationship between Eb/N0 and the signal-to-noise ratio. Therefore, the CPP method can increase the BER of unauthorized receivers.
The average power of signal is defined as P0 and the average power of noise is N1. When the transformation order is known, the SNR of the receiver is written as(10)μ=10⋅logP0N1.
Then, by using formula (5), ω02 is calculated as(11)ω02=ω0⋅ω0∗=cos2Δyπ4cos2Δyπ2.
For unauthorized receivers, the equivalent energy of the signal after WFRFT becomes P0′:(12)P0′=ω02cos2−3Δyπ4⋅P0=ω0cos−3Δyπ42⋅P0=cosΔyπ4cosΔyπ2cos−3Δyπ42⋅P0.
Thus, the equivalent SNR of the unauthorized receivers is expressed as(13)μ1=10⋅logP0′1−P0′+N1=10⋅logω0cos−3Δyπ/421−ω0cos−3Δyπ/42+10−μ/10.
Further, consider the mathematical relationship between Eb/N0 and P0/N1, which is shown in formula (14), then the logarithmic Eb/N0results of the unauthorized receivers are calculated as shown in formula (15), where W is signal bandwidth, Rb is the number of bits transmitted per second, and it is also the transmission rate:(14)EbN0=P0N1⋅WRb,(15)ebn0=10⋅logP0′1−P0′+N1⋅WRb=10⋅logω0cos−3Δyπ/421−ω0cos−3Δyπ/42+10−μ/10+10⋅logWRb.
This formula shows that the Eb/N0 decreases when the transformation order is unknown, which also shows that the BER increases for the unauthorized receivers. Thus, it shows that the theoretical antiscanning performance of the CPP method is superior.
In addition, if the modulation order is fixed, the unauthorized receiver can scan at a step less than 0.01 and illegally receive signals. The number of scans is shown in formula (16). It can be seen that when the number of scans reaches 400, the security of the system is threatened. However, the CPP method also solves this problem. Under the same condition of less than 0.01 step, the number of scans is shown in formula (17). The number of scans is increased by N times, and when the value of N is larger, the degree of improvement of the number of scans will be very large. Thus, CPP method greatly improves the scanning complexity of unauthorized users and the security performance of communication system:(16)χ=40.01=400,(17)χ′=N40.01=400N.
4. Test and Analysis
Considering that the CPP method is a transmission technique for authorized users, the CPP method is compared to the scanning and nonscanning methods for unauthorized users. For the scanning method, conditions 1 and 2 are the cases in which different modulation orders are scanned. Condition 1 is to lock a parameter in the parameter pool with the smallest error. Condition 2 is to lock a random parameter in the parameter pool.
4.1. BER Analysis
Firstly, when L = 102300, the BER is tested under different modulation modes, which are illustrated in Figures 3.
BER results under different modulation modes: (a) N = 3 and (b) N = 12.
It can be observed that with the increase in Eb/N0, the BER of the nonscanning method is always very high; hence, it cannot recover data in any environment. Then, the CPP method has the lowest BER and the best performance. And the BER of the scanning method is also very high, and it cannot recover data. In traditional scanning methods, the demodulation parameters are first estimated, then fixed, and then the demodulation parameters are not re-estimated. Thus, the usual scanning method can only lock the demodulation order at one time, not all the demodulation orders. And the scanning method cannot synchronize the demodulation order switching frequency. Therefore, the BER of the scanning method is higher than that of the CPP method. This also demonstrates that the CPP method can resist the illegal acceptance of unauthorized users under the premise of effective transmission performance.
In addition, under the same methods, the BER performance of BPSK modulation is slightly better than that of QPSK modulation, but the performance difference is not obvious. Therefore, the following tests are performed with BPSK modulation.
Then, the influence of the amount of data L on the BER is tested. When Eb/N0 is fixed at −8 dB, the lengths of the parameter pool N are 3 and 12. With the increase in L, the BER results are illustrated in Figures 4(a) and 4(b), respectively. Obviously, with the increase in L, the BER of each method is relatively stable. However, under the same L condition, the BER of the CPP method is the lowest, followed by the scanning method, while that of the nonscanning method is the highest. Furthermore, when the Eb/N0 is fixed at 6 dB, the lengths of the parameter pool N are 3 and 12. With the increase in the amount of data L, the BER results are presented in Figures 4(c) and 4(d), respectively, and the same conclusion can be obtained. It can be observed that L has little effect on the BER of all methods.
BER results with the change of L. (a) Eb/N0 = −8 dB and N = 3. (b) Eb/N0 = −8 dB and N = 12. (c) Eb/N0 = 6 dB and N = 3. (d) Eb/N0 = 6 dB and N = 12.
Because L has very little influence on the BER of all methods, the following tests are carried out under the condition of a fixed L of 102300. Next, we test the influence of the number of the parameter pool length N on the BER. Then, when the Eb/N0 is −8 and 6 dB, with the increase in N, the BER results are illustrated in Figure 5. It can be observed that, with the increase in N, the BER of the CPP method is the smallest and relatively stable, while the BER of the scanning method increases gradually. However, after N>10, the BER reaches a high level and the change is not obvious. The BER of the nonscanning method is always very high and is not affected by the change in N. Comparatively, the BER of the CPP method is the lowest, and the scanning method is also lower when N is small, although it is still substantially higher than that of the CPP method. The nonscanning method has the highest BER and cannot achieve the correct reception conditions.
BER results with the change of N. (a) Eb/N0 = −8 dB. (b) Eb/N0 = 6 dB.
Thereafter, the influence of the transmission environment on the BER is tested. When N is 2, 3, 6, and 12, with the increase in Eb/N0, the BER results are illustrated in Figure 6, respectively. It can be observed that, with the increase in Eb/N0, the BER of the nonscanning method is always very high; hence, it cannot recover data in any environment.
BER results with the change of Eb/N0. (a) N = 2. (b) N = 3. (c) N = 6. (d) N = 12.
When N = 2, the parameter pool of the CPP method is equivalent to that of m/M sequence, and m/M sequence parameter pool is no longer affected by N. No matter how much N is taken, the result is the same as that of N = 2. The m/M sequence is a binary sequence, whose values are 0 and 1, and it is easy to be measured. Therefore, the parameter pool formed by m/M sequence will only consist of 0 and 1. The antiscanning ability of the method based on m/M sequence parameter pool is the same as that of Figure 6(a). It can be seen that the BER always keeps the trend of Figure 6(a), and the antiscanning ability is greatly reduced.
However, Chaos sequence parameter pool is greatly affected by N, and the larger the N, the better the antiscanning ability, which is illustrated in Figures 6(a) to 6(d). When N is large, the BER of the scanning method is also very high, and it cannot recover data. When N is small, the BER of the scanning method decreases, but it is still higher than that of the CPP method. Because the scanning method can only lock one of the N parameters, the larger the N, the smaller the proportion of locked parameters. Therefore, the smaller the proportion of correct demodulation, the larger the BER. In particular, when N = 1, that is, when the modulation order is not switched, the BER of the scanning method is the same as that of the CPP method. Therefore, apart from the special condition of N = 1, the BER of the CPP method is the lowest, and it is substantially lower than that of the other methods. To sum up, the bigger the N, the better the antiscanning ability of the CPP method, and the antiscanning ability of chaotic sequence parameter pool is much higher than that of m/M sequence parameter pool.
4.2. Detection Probability Analysis
For a secure transmission system, the effective detection probability analysis is also the main index. The detection probability results are shown in Figure 7. It can be seen that the detection probability of the CPP method is the best, followed by the scanning method, and that of the nonscanning method is the worst. These demonstrate that the CPP method can resist the illegal acceptance of unauthorized users under the premise of effective transmission performance.
The detection probability results with the change of Eb/N0. (a) N = 2. (b) N = 3. (c) N = 6. (d) N = 12.
4.3. Speed Analysis
In addition, the processing time of various methods is tested. However, based on the previous test and analysis because the nonscanning method cannot recover data or cannot receive correctly, the speed analysis of the nonscanning method is meaningless. Therefore, speed tests are carried out on the CPP and scanning methods, and the processing time of the scanning method is illustrated in Figure 8(a), and the processing time of the CPP method is illustrated in Figure 8(b). The results indicate that the processing time of the two methods increases with an increase of the data amount L. Moreover, the processing time of the scanning method is not affected by the parameter N, and it is very long. The processing time of the CPP method is affected by N. With the increase of N, the processing time increases, that is, the processing speed decreases. The processing time of the CPP method is obviously lower than that of the scanning method, that is, the processing speed of the CPP method is higher than that of the scanning method, demonstrating that the CPP method exhibits superior time performance. This is also consistent with the theoretical speed performance. For the transmitters and receivers of the authorized user, the chaotic type, the initial chaotic value, and the discarding and switching frequency are known. The CPP method can receive the signal without a large number of scans, so the speed is faster. These parameters are unknown to unauthorized users, they still need to scan, so the speed performance of CPP methods is better than the traditional methods.
The processing time results with the change of L. (a) The scanning method. (b) The CPP method.
Finally, because the scanning complexity is also one of the indicators to evaluate the processing speed, and the complexity is mainly affected by the number of scans, so the scanning complexity is analyzed from the angle of the number of scans. The test results are shown in Table 1. Although the complexity of the nonscanning method is very low, it cannot receive data successfully. Comparatively speaking, the scanning complexity of the scanning method is very large, especially when N is large. The scanning complexity of the CPP method is very low, and it can receive data successfully, which shows that its performance is the best.
The scanning complexity results.
Number of scans
Number of N
N = 1
N = 2
N = 3
N = 5
N = 8
N = 10
N = 16
N = 20
Nonscanning method
1
1
1
1
1
1
1
1
Scanning method
400
800
2400
2000
3200
4000
6400
8000
CPP method
1
1
1
1
1
1
1
1
5. Conclusions
The WFRFT secure communication method based on the chaotic parameter pool has been proposed to improve information transmission security. Compared to the scanning and nonscanning methods, the CPP method exhibits superior performance. Although the BER of the scanning method is lower when N is smaller, it is still much higher than that of the CPP method. The nonscanning method has the highest BER and cannot achieve the correct reception conditions. The BER of the CPP method is the lowest, and detection probability is the best, which also indicates that the CPP method can resist the illegal acceptance of unauthorized users under the premise of effective transmission performance. Moreover, the processing speed of the CPP method is obviously faster than that of the other methods, demonstrating that the CPP method exhibits superior time performance.
Data Availability
The data used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This work was supported in part by the National Natural Science Foundation of China (Grant no. 61501309), the China Postdoctoral Science Foundation (Grant no. 2017T100185), the Liaoning Natural Science Foundation of China (Grant no. 2017011002-301), and the Liaoning Provincial Colleges and Universities Innovative Talents Support Program.
LiuY.NieZ.LiuQ. H.DIFFT: a fast and accurate algorithm for fourier transform integrals of discontinuous functions2008181171671810.1109/lmwc.2008.20051622-s2.0-55649111602WangT.HuanH.TaoR.WangY.Anti-eavesdropping FrFT-OFDM system exploiting multipath channel characteristics20171191371137810.1049/iet-com.2016.12512-s2.0-85022221829LuoY.YuC.ChenS.LiJ.RuanH.El-SheimyN.A novel doppler rate estimator based on fractional fourier transform for high-dynamic GNSS signal20197295752959610.1109/access.2019.2903185ChenR.WangY.Universal FRFT-based algorithm for parameter estimation of chirp signals201223449550110.1109/jsee.2012.000632-s2.0-84865625520WangH.JiangY.Real-time parameter estimation for SAR moving target based on WVD slice and FrFT2018541474910.1049/el.2017.1740YuanL.DaX.WuJ.XuR.ZhangZ.LiuH.WFRFT modulation recognition based on HOC and optimal order searching algorithm2018293462470WangZ.MeiL.ShaX.LeungV. C. M.BER analysis of WFRFT precoded OFDM and GFDM waveforms with an integer time offset201867109097911110.1109/tvt.2018.28556962-s2.0-85049963701DaX.LiangY.HuH.Embedding WFRFT signals into TDCS for secure communications20186549385495110.1109/access.2018.28729362-s2.0-85054351015FangX.ShaX.LiY.Secret communication using parallel combinatory spreading WFRFT2015191626510.1109/lcomm.2014.23592002-s2.0-84921416523FangX.ZhangN.ZhangS.ChenD.ShaX.ShenX.On physical layer security: weighted fractional fourier transform based user cooperation20171685498551010.1109/twc.2017.27121582-s2.0-85029596759LuoZ.WangH.ZhouK.LvW.Combined constellation rotation with weighted FRFT for secure transmission in polarization modulation based dual-polarized satellite communications20175270612707310.1109/access.2017.27676382-s2.0-85033369684FangX.ShaX.MeiL.Guaranteeing wireless communication secrecy via a WFRFT-based cooperative system2015129768210.1109/cc.2015.72752612-s2.0-84960117523ZhangZ.DaX.LiuH.A study of the covered characteristics of MAP-WFRFT satellite signals2017294460467FangX.ShaX.LiY.MP-WFRFT and constellation scrambling based physical layer security system2016132138145LiangY.DaX.XuR.NiD. L.ZhaiZ.PanY.Research on constellation-splitting criterion in multiple parameters WFRFT modulations20186343543436410.1109/access.2018.28489182-s2.0-85048855617RanQ.YeungD. S.TsangE. C. C.General multi fractional fourier trans-form method based on the generalized permutation matrix group20055318398LiJ.ShaX.FangX.MeiL.8-Weighted-type fractional Fourier transform based three-branch transmission method201815914715910.1109/cc.2018.84564592-s2.0-85053310740ManfrediP.GinsteD. V.ZutterD. D.CanaveroF. G.Generalized decoupled polynomial chaos for nonlinear circuits with many random parameters201525850550710.1109/lmwc.2015.24407792-s2.0-85027949262XuY.MiliL.SanduA.von SpakovskyM. R.ZhaoJ.Propagating uncertainty in power system dynamic simulations using polynomial chaos201934133834810.1109/tpwrs.2018.28655482-s2.0-85051731492