Modeling and Nonlinear Characteristics Analysis of Fluorescent Lamp Driven by a Full-Bridge Inverter

: Addressing the issue of system stability is a crucial step towards the successful integration of memristive devices in power electronics applications. This paper focuses on ﬂuorescent lamps possessing memristive characteristics and investigates the instability phenomena and mechanisms within a full-bridge inverter with a ﬂuorescent lamp load. Based on the memristive characteristics of ﬂuorescent lamps, this paper establishes the averaged model of the system, whose coefﬁcient matrix is nonlinear, periodic, and time-varying. This study identiﬁes the occurrence of low-frequency oscillations within the system and elucidates the fundamental mechanism underlying the emergence of low-frequency oscillations. Furthermore, this paper establishes the stability boundaries of the system across different parameter planes. The research ﬁndings indicate that the low-frequency oscillations within the system are attributed to the occurrence of Hopf bifurcations in a frequency range higher than line frequency but signiﬁcantly lower than switching frequency. Lastly, the PSpice circuit of the system is designed, and simulation results are provided for validation. This study can offer guidance on parameters and control strategies for ensuring the stable operation of a full-bridge inverter with ﬂuorescent lamps. Moreover, it can facilitate the comprehension of instability mechanisms in systems incorporating memristive devices, thereby offering a foundation for the expansion of memristor applications.


Introduction
In 1971, according to the symmetry and completeness principle of circuit theory, Chua LO proposed the concept of the memristor, which describes the relationship between the charge and the magnetic flux and is called the fourth basic circuit element [1].The proposal of the memristor concept increased the number of fundamental circuit elements from three to four, thereby enriching and completing the theoretical framework of circuit theory.Moreover, Chua LO has expanded and refined the concept of the memristor.He classified all memristors into three classes called ideal, Generic, or Extended memristors [2][3][4].However, despite the proposal of the memristor concept, progress in the corresponding research was hindered by the absence of a physical entity.Until 2008, HP LABS successfully manufactured solid-state memristor devices based on nanoscale TiO 2 [5].Since then, the study of memristors and memristive systems has attracted more and more attention, resulting in a wealth of research achievements documented in the open literature.For example, Ventra MD extended the concept of memristors to include meminductors and memcapacitors [6].On this basis, Yin Z Y provided a more detailed explanation of the memristor, meminductor, and memcapacitor [7].Biolek Z studied the mechanism of the hysteretic curve of the memristor [8] and proposed the method for calculating the area in the hysteretic curve of the memristor, which proved that the area of the pinched hysteresis loop of the current-controlled ideal memristor represents the quantity "content" introduced into the theory of nonlinear systems by Millar in 1951 [9,10].In addition, according to the Energies 2023, 16 results on fluorescent lamps in the literature [11], it can be seen that fluorescent lamps show the obvious memristor's characteristics, which arouse people's interest.For example, based on the theoretical analysis of low-pressure discharge in the positive column region and the Francis equation, Dai Guojun proposed a PSpice model of a fluorescent lamp [12].Based on gas discharge theory, Wei Yan established a semi-theoretical model of a fluorescent lamp that fully accords with the definition of a fluid-controlled memristor.The model can predict the electrical characteristics of different types of fluorescent lamps at different operating frequencies and different dimming levels [13].On the basis of the literature [13], D. Lin used a unified mathematical framework to build mathematical models of highpressure and low-pressure gas discharge lamps and pointed out that both high-pressure and low-pressure gas discharge lamps have memristor characteristics [14].Based on [14], S. P. Adhikari briefly analyzed the non-crossing characteristics of the hysteresis loops of fluorescent lamps at the origin [15].Moreover, a considerable number of electrical devices have been experimentally identified as belonging to the category of memristive systems.These devices include; however, are not limited to, thermistors, tungsten filament lamps, current-controlled discharge tubes, and first-generation radio detectors [2,16,17].Additionally, certain commonly encountered circuit systems have exhibited memristive characteristics in their external behavior during energization.For instance, in a significant study conducted in 2012, Corinto and Ascoli demonstrated that a diode bridge post-connected with an LCR filter can manifest memristive external characteristics, leading to its designation as a memristive diode bridge circuit [18].It is foreseeable that in the future, an increasing number of electrical devices exhibiting memristive characteristics will be discovered or manufactured.Memristive electrical devices require operating in an alternating current environment devoid of direct current bias.The inverter circuits are power electronic conversion circuits capable of converting DC to AC power and serve as vital power sources for memristive electrical devices.Due to the transitions between different states of the power switches, the inverter itself operates as a segmented circuit, exhibiting intricate nonlinear phenomena.
Currently, research on complex nonlinear phenomena in inverters has become widespread.In 2002, Robert et al. first extended the study of nonlinear complex behaviors to the realm of DC/AC converters, analyzing boundary-collision bifurcations and chaos phenomena in single-phase H-bridge inverters [19].In 2006, Kousaka et al. conducted an analysis of period orbit stability in single-phase H-bridge inverters using the return map method.By using this map, the stability of the periodic orbit can be analyzed simply using the symbolic sequence [20].In 2009, Wang et al. established a first-order discrete model through iterative techniques and employed frequency-flashing mapping, folding mapping, and bifurcation mapping to describe bifurcation and chaos phenomena in proportional control single-phase H-bridge inverter systems [21,22].In 2012, Wu et al. proposed control methods for fastscale and slow-scale bifurcations in voltage-controlled H-bridge inverter circuits from a frequency-domain perspective [23].In 2013, Liu et al. investigated the symmetric phenomena of fast-scale bifurcations in peak current and valley current-controlled H-bridge inverter circuits.The simulation showed that the symmetrical dynamical phenomenon occurs in the single-phase H-bridge inverter controlled by the peak current or the valley current [24].In 2016, Zhang et al. proposed an improved averaging modeling method and analyzed the fast-scale bifurcation behavior in peak current-controlled Buck-Boost inverters.This method can increase the resolution of the conventional classic averaged model by half the switching frequency [25].However, despite the significant body of research in this area, it is worth noting that the majority of the aforementioned studies have predominantly concentrated on linear loads, and the stability analysis of inverters with memristive loads has not been extensively explored.Existing research findings have demonstrated that memristive loads have a certain impact on the dynamic characteristics of power electronic circuits [26][27][28].Due to the nonlinear characteristics of memristors, the full-bridge inverter with memristive loads exhibits nonlinearity during the periods of power switching on and off, indicating that it is a piecewise nonlinear system.Furthermore, in comparison with Energies 2023, 16, 6622 3 of 20 inverters having linear loads, the memristive nature of memristors leads to an increase in the system's order.This complexity contributes to a range of nonlinear behaviors and challenges, thus amplifying the intricacy and difficulty of system parameter design.Some investigations have resorted to approximating memristive loads as linear loads for the sake of simplicity [29].Such oversimplified approaches inevitably overlook the intricate and unique characteristics inherent in inverters with memristive loads, thereby hindering their ability to provide comprehensive guidance for the precise design and optimization of these specialized systems.
Therefore, the stable operation of inverters with memristive loads becomes a pivotal foundation driving the advancement and application of memristors and is poised to emerge as a crucial issue in the future.Fluorescent lamps, as the most common memristive electrical devices, exhibit significant representativeness.In practical applications, an inverter supplies high-frequency AC power to fluorescent lamps to drive their stable operation [30].However, both fluorescent lamps and inverters are strong nonlinear systems, leading to the full-bridge inverter with the fluorescent lamp being a class of strongly nonlinear systems in which there must be a wealth of complex nonlinear phenomena, such as bifurcation, chaos, and so on, leading to greater noise and instability, not conducive to practical engineering applications.Therefore, a comprehensive investigation into the nonlinear dynamic behavior of the inverter with a fluorescent lamp load not only serves as a valuable reference for better comprehension and enhancement of system design but also contributes to the advancement of memristive device applications in power electronic circuits.
Addressing the aforementioned issues and limitations, this paper establishes a mathematical model for a full-bridge inverter system with a fluorescent lamp, taking into account the memristive characteristics of the fluorescent lamp.Through Matlab/Simulink simulations, we identify the presence of low-frequency oscillations in the system.Employing the harmonic balance method and Floquet theory, it elucidates the fundamental mechanism of these oscillations and outlines the system's stability boundaries across various parameter planes.Furthermore, a PSpice simulation model is constructed, yielding results consistent with theoretical analysis, thus validating the accuracy of the theoretical insights.

Mathematical Modeling and Simulations
In 2007, Wei Yan et al. established the mathematical model of a fluorescent lamp through experiments and a genetic algorithm [14], namely where t is the time, e is the charge on the electron, and e = 1.6 × 10 −19 C, k is the Boltzmann constant, and k = 1.38 × 10 −23 J/K, T is the electron temperature, and the initial value of T is T 0 = 350 K, and a 1 -a 6 are the adjustable model coefficients.
Based on the classification of memristors in Ref. [31], the mathematical model of a fluorescent lamp conforms to the constitutive relation of a current-controlled memristor, whose memristor value is The second and third expressions of (1) are transcendental functions with respect to T, which will bring difficulties to theoretical calculation and analysis.Here, according to the model parameters given in Ref. [15], the approximate fitting is performed on equations a 5 T −0.75 exp(ea 6 /(2kT)) and a 2 exp(−ea 3 /(kT)).That is, a 5 T −0.75 exp( ea 6 2kT ) ≈ 2.145 × 10 15 • T −4.518 a 2 exp( −ea 3 kT ) ≈ 5.569 × 10 −56 •T 16.28   (3) As can be seen from Figure 1, the fitting R-square of a 5 T −0.75 exp(ea 6 /(2kT)) is 0.9944, and the fitting R-square of a 2 exp(−ea 3 /(kT)) is 0.9998, which shows that the fitting accuracy is very high.Thus, it is effective to use the simplified model instead of a 5 T −0.75 exp(ea 6 /(2kT)) and a 2 exp(−ea 3 /(kT)).
a5T −0.75 exp(ea6/(2kT)) and a2exp(−ea3/(kT)).That is, As can be seen from Figure 1, the fitting R-square of a5T −0.75 exp(ea6/(2kT)) is 0.9944, and the fitting R-square of a2exp(−ea3/(kT)) is 0.9998, which shows that the fitting accuracy is very high.Thus, it is effective to use the simplified model instead of a5T −0.75 exp(ea6/(2kT)) and a2exp(−ea3/(kT)).Generally speaking, the operational frequency of fluorescent lamps typically ranges from 20 kHz to 100 kHz.A high-frequency AC power supply enables fluorescent lamps to produce more stable and uniform illumination, thereby reducing visible flickering.In electronic ballasts, the high-frequency AC power is generated by the inverter circuit.Therefore, the inverter is the most basic and critical component of the electronic ballast, and the full-bridge inverter is one of the most common topologies.Therefore, here, the PI current feedback-controlled full-bridge inverter is used to drive a fluorescent lamp as the research object, and its circuit model is shown in Figure 2.This circuit is mainly composed of the DC voltage source E, power switches Q1~Q4, LCC (LCC is a circuit topology consisting of an inductor and two capacitors, hence its name LCC) branch, the fluorescent lamp load, and the current feedback controller.The control section adopts the SPWM (Sinusoidal Pulse Width Modulation) current control mode, where the output current i L is compared with the reference current i ref .
The error signal resulting from this Energies 2023, 16, 6622 5 of 20 comparison is fed into a PI controller to obtain the modulating signal i c .This modulating signal i c is then compared with a triangular wave i t to generate a driving signal, which is used to control the on-off of the power switches.When the driving signal is high, the corresponding power switches are on; otherwise, they are off.Note that Q1 and Q4 have the same drive signals, as do Q2 and Q3.Note that the insulation resistance R c2 of capacitor C 2 is considered (note: R c2 stands for loss in the dielectric, and R c2 is usually large, generally above 1 MΩ [32,33]).
research object, and its circuit model is shown in Figure 2.This circuit is mainly compo of the DC voltage source E, power switches Q1~Q4, LCC (LCC is a circuit topolo consisting of an inductor and two capacitors, hence its name LCC) branch, the fluoresc lamp load, and the current feedback controller.The control section adopts the SPW (Sinusoidal Pulse Width Modulation) current control mode, where the output current compared with the reference current iref.The error signal resulting from this compari is fed into a PI controller to obtain the modulating signal ic.This modulating signal then compared with a triangular wave it to generate a driving signal, which is used control the on-off of the power switches.When the driving signal is high, corresponding power switches are on; otherwise, they are off.Note that Q1 and Q4 h the same drive signals, as do Q2 and Q3.Note that the insulation resistance Rc2 of capac C2 is considered (note: Rc2 stands for loss in the dielectric, and Rc2 is usually large, gener above 1 MΩ [32,33]).The bipolar SPWM generator in Figure 2 consists of a COM(Comparator) triangular carrier it, and a power switch driver.Note that the chosen operating freque f0 for the fluorescent lamp is f0 = 25 kHz, which results in ω0 = 2πf0 = 50,000π, And triangular carrier it and reference current iref are presented as follows: The bipolar SPWM generator in Figure 2 consists of a COM (Comparator), a triangular carrier i t , and a power switch driver.Note that the chosen operating frequency f 0 for the fluorescent lamp is f 0 = 25 kHz, which results in ω 0 = 2πf 0 = 50,000π, And the triangular carrier i t and reference current i ref are presented as follows: where ω sw and I tri are the angular frequency and amplitude of the triangular carrier respectively.
Energies 2023, 16, 6622 6 of 20 The output of a comparator is high; otherwise, it is low.According to KCL (Kirchhoff's Current Law) and KVL (Kirchhoff's Voltage Law) [34], the mathematical model of the system described by the differential equation can be derived as follows: When S = 1, means Q1 and Q4 are turned on, and Q2 and Q3 are turned off.On the contrary, Q2 and Q3 are turned on, and Q1 and Q4 are turned off.According to Formula (5), the MATLAB/Simulink simulation model of the full-bridge inverter system with the fluorescent lamp is established, as shown in Figure 3.As indicated by Equation ( 5) and Figure 3 in the original manuscript, it can be observed that when a fluorescent lamp is connected, compared to an inverter with a linear load, the memristive nature of the lamp leads to an increase in the system's order.Usually, power electronic systems can be regarded as segmented linear switching circuits.However, due to the memristive characteristics of fluorescent lamps, the full-bridge inverter with fluorescent lamps exhibits nonlinearity during the periods of power switching on and off.In other words, the full-bridge inverter with fluorescent lamp is a segmented nonlinear system that is more complicated and more prone to generating oscillations and other nonlinear behaviors.In actual electronic ballasts, the amplitude of the rectified DC voltage is around 300 V; hence, we opted for E = 300 V.The inductance L of the LCC branch is on the order of millihenries (mH), while both C1 and C2 have values in the nanofarad (nF) range.The electronic ballast achieves high-voltage ignition of the fluorescent lamp through the resonance of the LCC branch.The resonance frequency for ignition (fres1) should be higher than the operating frequency of the fluorescent lamp.And to ensure that the LCC branch operates in an inductive state, the operating frequency of the fluorescent lamp needs to be higher than the steady-state resonance frequency (fres2).In actual electronic ballasts, the amplitude of the rectified DC voltage is around 300 V; hence, we opted for E = 300 V.The inductance L of the LCC branch is on the order of millihenries (mH), while both C operates in an inductive state, the operating frequency of the fluorescent lamp needs to be higher than the steady-state resonance frequency (f res2 ).
Choosing C 1 = 10 nF, C 2 = 100 nF, and L = 1.5 mH can meet the requirements.The other circuit parameters are selected as Figure 4 shows the Matlab/Simulink simulation results under R = 5 kΩ.In Figure 4a, U c1 is the time-domain waveform of the voltage across the fluorescent lamp, while Figure 4b is the time-domain waveform of the current through the fluorescent lamp.From Figure 4, it is obvious that the system is in stable operation.However, the system is unstable and in a state of low-frequency oscillation when R = 50 Ω, which is shown in Figure 5. Obviously, the resistance R has an important effect on the occurrence of low-frequency oscillation in the fluorescent lamp system driven by the full bridge inverter.

Analysis of Low-Frequency Oscillation Mechanism of The System
Firstly, the approximate solution of a full-bridge inverter with a fluorescent lamp is obtained by the harmonic balance method.Then, the Floquet theory is used to identify the mechanism responsible for the occurrence of low-frequency oscillation.

Approximate Solution
Since the angular frequency of the triangular carrier it is much larger than that of the modulated wave, the switching frequency is much larger than that of the modulated wave's frequency.Therefore, based on the average modeling method of power electronic systems [35], the average model of the system can be established, and the expression of the duty cycle d   is

Analysis of Low-Frequency Oscillation Mechanism of the System
Firstly, the approximate solution of a full-bridge inverter with a fluorescent lamp is obtained by the harmonic balance method.Then, the Floquet theory is used to identify the mechanism responsible for the occurrence of low-frequency oscillation.

Approximate Solution
Since the angular frequency of the triangular carrier i t is much larger than that of the modulated wave, the switching frequency is much larger than that of the modulated wave's frequency.Therefore, based on the average modeling method of power electronic systems [35], the average model of the system can be established, and the expression of the duty cycle d is By substituting d for the term S in the first formula of Equation ( 7), the average model of the full-bridge inverter with fluorescent lamp controlled by PI current feedback control can be established as where, i L , i c , u C1 and T are respectively the average values of i L , i c , u C1 and T in a switching cycle.
It can be seen that the current mode-controlled full-bridge inverter with fluorescent lamp described in ( 11) is a nonlinear and non-autonomous system.Equation ( 11) involves five variables.In order to obtain the approximate solution of the full-bridge inverter with a fluorescent lamp, it is necessary to reduce the number of variables through equation manipulation.By combining the first, third, and fifth formulas of Equation ( 11), the variables u C2 and i c can be eliminated.The specific process is as follows: From the first formula of (11), we can obtain After differentiating both sides of ( 12) with respect to time t, we can obtain After substituting the third expression of ( 11) into ( 13), the variable u C2 can be eliminated, resulting in where, exp is the exponential function.
By combining the fifth formula of ( 11) with ( 14), the variable i c can be eliminated.we can obtain Equation ( 15) contains only the variables i L and u C1 , as well as their differentials and integrals, which can be abbreviated as where, By combining the second and fourth formulas of Equation ( 11), the variables T can be eliminated.The specific process is as follows: From the second formula of (11), we can obtain Energies 2023, 16, 6622 10 of 20 After differentiating both sides of ( 19) with respect to time t, we can obtain By combining (19), (20), and the fourth formula of (11), we can obtain 547.6( Equation ( 15) contains only the variables i L and u C1 , as well as their differentials, which can be abbreviated as where, For this system, the average value of the output voltage u C1 and the inductor current i L , which are expected to contain as few harmonic components as possible, can be approximated as a series of sine waves.Therefore, suppose Defining that M = [1, cos(ω 0 t), . . . ,cos(Nω 0 t), sin(ω 0 t), . . . ,sin(Nω 0 t)] where, Vectors A and B contain coefficients to be determined, with a total of (4N + 2) coefficients in both vectors combined.N is the chosen harmonic number.The larger the value of N, the more accurate the result.
Therefore, i L , u c1 and their first and second derivatives can be written as By substituting ( 27) into ( 16) and ( 22), it can be obtained Since Equation ( 28) involves (4N + 2) unknowns, it is necessary to establish a set of (4N + 2) algebraic equations to solve for the unknowns.Considering the orthogonality of trigonometric functions, (4N + 2) algebraic equations can be obtained by multiplying both sides of Equation ( 28) by M T and integrating over time t from 0 to T 1 .Where, T 1 = 2π/ω 0 .
Figure 6 shows the comparison between the simulation results obtained from the simulation model in MATLAB/Simulink and the approximate solutions obtained from the harmonic balance method.It can be seen that the two results are consistent, indicating that the harmonic balance method can effectively calculate the approximate solution of a current-mode-controlled full-bridge inverter with a fluorescent lamp.

Theoretical Analysis
Suppose the periodic equilibrium solution of the system is When the periodic state of the system is subjected to a small disturbance

Theoretical Analysis
Suppose the periodic equilibrium solution of the system is When the periodic state of the system is subjected to a small disturbance ∆X(t), there is By substituting ( 32) into (11) and linearizing the average model equation at the periodic equilibrium solution X 0 (t), the following linear time-varying differential equation can be obtained: where, (36) According to (34), H(t) is a periodic matrix, and its period is the same as that of X 0 (t).Therefore, ∆X (t) = H(t)∆X(t) has a fundamental solution matrix, which can be expressed as where P i (t) represents the solution vector of (33), with By solving differential Equation ( 40), it can be obtained Equation ( 41) satisfies the initial condition P(0) = I, and I is the fifth-order identity matrix.Since P(t + T 1 ) = P(t), each column vector of P(t + T 1 ) is also a fundamental solution to Equation (33).Therefore, there is a linear correlation between P(t + T 1 ) and P(t), and their relationship can be expressed as where, Z is referred to as the transfer matrix.Set t = 0, the transfer matrix Z can be written as The eigenvalues of the transfer matrix Z are known as Floquet multipliers.According to Floquet theory, the full-bridge inverter with the fluorescent lamp is stable only if all Floquet multipliers are located within the unit circle; otherwise, the type of bifurcation can be determined by the way in which the Floquet multipliers break through the unit circle.
• If one of the Floquet multipliers of the transition matrix Z passes through the unit circle from the negative real axis and the rest remains in the unit circle, the periodic time-varying system has a period-doubling bifurcation; • If one of the Floquet multipliers of the transition matrix Z passes through the unit circle from the positive real axis and the rest remains in the unit circle, the periodic time-varying system has a pitch-fork bifurcation; • If a pair of conjugate Floquet multipliers of the transition matrix Z exit the unit circle and the rest remain in the unit circle, the periodic time-varying system has a Hopf bifurcation.
Therefore, it is necessary to calculate the eigenvalues of the transfer matrix Z to determine the stability and instability mechanisms of the full-bridge inverter with the fluorescent lamp.The Jacobian matrix of the averaged model of a full-bridge inverter with a resistive load remains constant.However, as indicated by Equation ( 34), H(t) it is a time-varying matrix.This is a distinctive difference between a full-bridge inverter with a memristive load and one with a resistive load.Consequently, numerical methods are required to compute the transfer matrix Z.
Divide the period [0, T 1 ] into N k subintervals, then the length of each subinterval is And the kth subinterval starts at The time variable matrix P(t) is continuous with respect to time t.Therefore, when N k is sufficiently large, the value of the matrix P(t) in the kth subinterval can be replaced by the mean value of P(t) in this subinterval, then where By substituting (46) into (43), the transition matrix Z can be expressed as Energies 2023, 16, 6622 14 of 20 The power series expansion formula of the exponential function is known as where N j represents the number of terms in the exponential series expansion.Clearly, as the value of N j increases, the computational accuracy of the results improves.
The Floquet multipliers can be obtained by calculating |λI − Z| = 0.The circuit parameters are the same as in Section 3.1.Change the resistance R and observe the change in the Floquet multipliers of the system.The calculation results are shown in Table 1 and Figure 7.The Floquet multipliers can be obtained by calculating |λI − Z| = 0.The circuit parameters are the same as in Section 3.1.Change the resistance R and observe the change in the Floquet multipliers of the system.The calculation results are shown in Table 1 and Figure 7.As shown in Table 1 and Figure 7, as the resistance R decreases from 5000 Ω to 50, Ω, a pair of complex conjugate Floquet multipliers penetrate from the inside of the unit circle and move to the outside, indicating that the full-bridge inverter with the fluorescent lamp transitions from a stable to an unstable working state.The system has a pair of conjugate Floquet multipliers that penetrate out of the unit circle, while the rest of the Floquet multipliers remain in the unit circle, which means that Hopf bifurcation occurs in this system.Therefore, the low-frequency oscillation phenomenon in the system is caused by Hopf bifurcation.As shown in Table 1 and Figure 7, as the resistance R decreases from 5000 Ω to 50 Ω, a pair of complex conjugate Floquet multipliers penetrate from the inside of the unit circle and move to the outside, indicating that the full-bridge inverter with the fluorescent lamp transitions from a stable to an unstable working state.The system has a pair of conjugate Floquet multipliers that penetrate out of the unit circle, while the rest of the Floquet multipliers remain in the unit circle, which means that Hopf bifurcation occurs in this system.Therefore, the low-frequency oscillation phenomenon in the system is caused by Hopf bifurcation.

PSpice Simulations
In order to verify the theoretical analysis, a simulation circuit model of a full-bridge inverter with a fluorescent lamp was constructed in PSpice.The circuit parameters are the same as those in Section 3.1, where LF356 is selected for the operational amplifier and LM311 for the comparator.
Figure 8 shows the time-domain waveform of the voltage across the fluorescent lamp and the current through it when R = 5 kΩ, and the system is in stable operation.Figure 9 shows the PSpice simulations for R = 50 Ω.It is clear that the system is operating in an unstable state.For the purpose of comparison, the black dashed lines in Figures 8 and 9 represent the simulation results from Matlab/Simulink in Figures 4 and 5.After comparing Figures 8 and 9 with Figures 4 and 5, respectively, it is found that the PSpice simulations are consistent with those of MATLAB/Simulink simulations, which effectively verifies the correctness of the theoretical analysis.Figure 8 shows the time-domain waveform of the voltage across the fluorescent lamp and the current through it when R = 5 kΩ, and the system is in stable operation.Figure 9 shows the PSpice simulations for R = 50 Ω.It is clear that the system is operating in an unstable state.For the purpose of comparison, the black dashed lines in Figures 8 and 9 represent the simulation results from Matlab/Simulink in Figures 4 and 5.After comparing Figures 8 and 9 with Figures 4 and 5, respectively, it is found that the PSpice simulations are consistent with those of MATLAB/Simulink simulations, which effectively verifies the correctness of the theoretical analysis.In addition, by importing the data into MATLAB and selecting a fundamental frequency of 25 kHz for spectrum analysis, the harmonic spectrum of voltage Uc1 can be obtained.As depicted in Figure 10a, it can be observed that when the resistance R is 5 KΩ, the total harmonic distortion (THD) of the output voltage Uc1 is measured at 1.21% within the frequency range of 0-500 kHz.Remarkably, the harmonics exhibit smaller amplitudes compared to the fundamental frequency, with each harmonic amplitude approximating 0 V.However, when the resistance R is 50 Ω, the THD of the output voltage Uc1 significantly increases to 27.70%.Under this condition, the system exhibits oscillations; the oscillation frequency is mainly around 110 KHz which is between the operating frequency of 25 kHz and the switching frequency of 500 kHz.This oscillation phenomenon can be called ʺlowfrequency oscillation.ʺIn the region of the oscillation frequency, the voltage Uc1 has larger harmonics.Therefore, the occurrence of low-frequency oscillation will make the fluorescent lamp work in harsh environments.In addition, by importing the data into MATLAB and selecting a fundamental frequency of 25 kHz for spectrum analysis, the harmonic spectrum of voltage U c1 can be obtained.As depicted in Figure 10a, it can be observed that when the resistance R is 5 KΩ, the total harmonic distortion (THD) of the output voltage U c1 is measured at 1.21% within the frequency range of 0-500 kHz.Remarkably, the harmonics exhibit smaller amplitudes compared to the fundamental frequency, with each harmonic amplitude approximating 0 V.However, when the resistance R is 50 Ω, the THD of the output voltage U c1 significantly increases to 27.70%.Under this condition, the system exhibits oscillations; the oscillation frequency is mainly around 110 KHz which is between the operating frequency of 25 kHz and the switching frequency of 500 kHz.This oscillation phenomenon can be called "lowfrequency oscillation".In the region of the oscillation frequency, the voltage U c1 has larger harmonics.Therefore, the occurrence of low-frequency oscillation will make the fluorescent lamp work in harsh environments.Figure 11a shows the stability boundary of the system obtained by PSpice simulations and theoretical calculations in the C1-R parameter plane under different parameters a4, indicating that with the decrease of capacitor C1 or resistance R, the fullbridge inverter with the fluorescent lamp is more prone to instability.In addition, when a4 = 0.07266, the stable region is A1 + A2, while when a4 = 5.07266, the stable region decreases to region A1.Therefore, the stability region of the system will decrease with the increase of the parameter a4. Figure 11b shows the stability boundary of the system obtained by PSpice simulations and theoretical calculations in the R0-R parameter plane under different capacitance C1, indicating that the full-bridge inverter with the fluorescent lamp is more prone to instability with the decrease of resistance R0 or R. In Figure 11b, when C1= 15 nF, the stable region is A3 + A4, while when C1= 10 nF, the stable region is reduced to region A3.Thus, the stable region of the system will increase with the increase of the capacitance C1.Based on the stability boundaries of the system, it is possible to optimize the parameter design of the system, thereby reducing the risk of system oscillation or even failure, and enhancing the reliability and stability of the system.Figure 11a shows the stability boundary of the system obtained by PSpice simulations and theoretical calculations in the C 1 -R parameter plane under different parameters a 4 , indicating that with the decrease of capacitor C 1 or resistance R, the full-bridge inverter with the fluorescent lamp is more prone to instability.In addition, when a 4 = 0.07266, the stable region is A1 + A2, while when a 4 = 5.07266, the stable region decreases to region A1.Therefore, the stability region of the system will decrease with the increase of the parameter a 4 .Figure 11b shows the stability boundary of the system obtained by PSpice simulations and theoretical calculations in the R 0 -R parameter plane under different capacitance C 1 , indicating that the full-bridge inverter with the fluorescent lamp is more prone to instability with the decrease of resistance R 0 or R. In Figure 11b, when C 1 = 15 nF, the stable region is A3 + A4, while when C 1 = 10 nF, the stable region is reduced to region A3.Thus, the stable region of the system will increase with the increase of the capacitance C 1 .Based on the stability boundaries of the system, it is possible to optimize the parameter design of the system, thereby reducing the risk of system oscillation or even failure, and enhancing the reliability and stability of the system.

Conclusions
This paper presents the first analysis of the impact of fluorescent lamp memristive characteristics on the stability of a full-bridge inverter.It identifies the presence of lowfrequency oscillations in a full-bridge inverter with a fluorescent lamp.By utilizing harmonic balance methods and Floquet theory, the essential mechanism generating lowfrequency oscillations is unveiled as the occurrence of Hopf bifurcation in the system, along with establishing stability boundaries across different parameter planes.Additionally, PSpice models for both the fluorescent lamp and the entire system are developed, and the obtained simulation results are consistent with the theoretical analysis.Experiments will be conducted in future work.However, in general, under the premise of consistency between theoretical analysis results and PSpice simulation results, the theoretical analysis outcomes possess a higher level of reliability.The conclusion of this paper can offer guidance on parameters and control strategies for the design of a full-

Conclusions
This paper presents the first analysis of the impact of fluorescent lamp memristive characteristics on the stability of a full-bridge inverter.It identifies the presence of low-frequency oscillations in a full-bridge inverter with a fluorescent lamp.By utilizing harmonic balance methods and Floquet theory, the essential mechanism generating low-frequency oscillations is unveiled as the occurrence of Hopf bifurcation in the system, along with establishing stability boundaries across different parameter planes.Additionally, PSpice models for both the fluorescent lamp and the entire system are developed, and the obtained simulation results are consistent with the theoretical analysis.Experiments will be conducted in future work.However, in general, under the premise of consistency between theoretical analysis results and PSpice simulation results, the theoretical analysis outcomes possess a higher level of reliability.The conclusion of this paper can offer guidance on parameters and control strategies for the design of a full-bridge inverter with a fluorescent lamp load.Moreover, it can facilitate the comprehension of instability mechanisms in

Figure 1 .
Figure 1.Comparison between the original formula and the result after fitting an approximation.(a) The fitting effect of a 5 T −0.75 exp(ea 6 /(2kT)); (b) The fitting effect of a 2 exp(−ea 3 /(kT)).

Figure 2 .
Figure 2. Schematic diagram of a current feedback-controlled full-bridge inverter with fluores lamp load.

Figure 2 .
Figure 2. Schematic diagram of a current feedback-controlled full-bridge inverter with fluorescent lamp load.

Figure 3 .
Figure 3. Simulation model for the system in the MATLAB/Simulink.

Figure 3 .
Figure 3. Simulation model for the system in the MATLAB/Simulink.
1 and C 2 have values in the nanofarad (nF) range.The electronic ballast achieves high-voltage ignition of the fluorescent lamp through the resonance of the LCC branch.The resonance frequency for ignition (f res1 ) should be higher than the operating frequency of the fluorescent lamp.And to ensure that the LCC branch Energies 2023, 16, 6622 7 of 20 Energies 2023, 16, x FOR PEER REVIEW 8 of 22 (a) (b)

Figure 6 .
Figure 6.Comparison between harmonic balance method and Simulink simulation.

Figure 6 .
Figure 6.Comparison between harmonic balance method and Simulink simulation.

Figure 7 .
Figure 7. Trajectory of Floquet multipliers λ1,2 and λ3 with R in the complex plane.

Figure 11 .
Figure 11.Under different parameters, the system stability boundary is compared by theoretical calculations and PSpice simulations (a) Different parameters a4 in the C1-R parameter plane; (b) Different parameters in C1 in the R0-R parameter plane.

Figure 11 .
Figure 11.Under different parameters, the system stability boundary is compared by theoretical calculations and PSpice simulations (a) Different parameters a 4 in the C 1 -R parameter plane; (b) Different parameters in C 1 in the R 0 -R parameter plane.

Table 1 .
Floquet multipliers and operating conditions of a full-bridge inverter with fluorescent lamp.

Table 1 .
Floquet multipliers and operating conditions of a full-bridge inverter with fluorescent lamp.