## 2 Motor Model and Brain Emotional Limbic System Calculation Model
### SynRM Model
Furthermore, in order to build a more accurate SynRM model, the concept of magnetic saturation parameter \(K_{s}\) is introduced and defined as follows:  
\[K_{s}(I_{m})=\frac{L^{\prime}_{d}(I_{m})}{L_{d}}=\frac{L^{\prime}_{q}(I_{m})}{ L_{q}}, \tag{4}\]  
where \(L^{\prime}_{d}(I_{m})\) and \(L^{\prime}_{q}(I_{m})\) are respectively the inductances of the d-axis and q-axis under magnetic saturation, whose values are affected by \(I_{m}\). According to Formula (4), \(L_{d}/L_{q}\) is a constant value; that is, \(K_{s}\) is a constant value, regardless of whether the click is affected by magnetic saturation during the whole operation range of the motor. Therefore, the relationship between the magnetic saturation parameter \(K_{s}\) and the equivalent excitation current \(I_{m}\) can be found and the value of \(K_{s}\) can be corrected by \(I_{m}\) under different current conditions. The inductance values of the \(d\)-axis and \(q\)-axis of the motor under magnetic saturation can be calculated using the formula (4). The \(d-q\) axis voltage equation of the synchronous reluctance motor under magnetic saturation can be expressed as:  
\[\left\{ \begin{array}{l}u^{\prime}_{d}(K_{s})=R_{s}i_{d}-\omega_{r}\psi^{ \prime}_{d}(K_{s})+\frac{d\psi^{\prime}_{d}(K_{s})}{dt}\\ u^{\prime}_{q}(K_{s})=R_{s}i_{q}-\omega_{r}\psi^{\prime}_{q}(K_{s})+\frac{d \psi^{\prime}_{q}(K_{s})}{dt}\end{array} \right.. \tag{5}\]  
Among them, \(u^{\prime}_{d}(K_{s})\) and \(u^{\prime}_{q}(K_{s})\) are the \(d\)- and \(q\)-axis voltage components under magnetic saturation, \(R_{s}\) is the stator phase winding. \(\psi^{\prime}_{d}(K_{s})\) and \(\psi^{\prime}_{q}(K_{s})\) are the \(d\)- and \(q\)-axis flux components under magnetic saturation. The expressions are:  
\[\left\{ \begin{array}{l}\psi^{\prime}_{d}(K_{s})=K_{s}(I_{m})L_{d}i_{d}\\ \psi^{\prime}_{q}(K_{s})=K_{s}(I_{m})L_{q}i_{q}\end{array} \right.. \tag{6}\]  
The torque equation of the SynRM under magnetic saturation is:  
\[T^{\prime}_{em}(K_{s})=pK_{s}(L_{d}-L_{q})i_{d}i_{q}. \tag{7}\]  
In the formula, \(T^{\prime}_{em}(K_{s})\) is the torque of the motor under magnetic saturation. Based on the variation trend of \(L_{d}\) and \(L_{q}\) with \(i_{d}\) and \(i_{q}\) shown in Figure 2 and the relationship between magnetic saturation parameter and \(d\)-axis and \(q\)-axis inductance in Formula (4), the curve of magnetic saturation parameter \(K_{s}\) with \(I_{m}\) in Figure 4 can be obtained. Figure 3: Flux density diagram of the motor under different currents: (**a**) \(i_{1}\); (**b**) \(i_{2}\); (**c**) \(i_{3}\). By fitting the magnetic saturation parameter data in Figure 4, the expression of magnetic saturation parameter \(K_{\rm s}\) is obtained as follows:  
\[K_{\rm s}(I_{m})=p_{1}I_{m}^{4}+p_{2}I_{m}^{3}+p_{3}I_{m}^{2}+p_{4}I_{m}+p_{5}. \tag{8}\]  
The values of the coefficients in the fitting function are shown in Table 1.  
By fitting the relation between \(K_{\rm s}\) and equivalent magnetization current \(I_{m}\), it can be found that the critical point of motor saturation is \(I_{m}\) = 10 A and \(I_{d}\) = 8 A. When \(I_{m}\) < 10 A (\(I_{d}\) < 8 A), \(K_{\rm s}\) = 1, it is clear that the motor has not entered the magnetic saturation state. When \(I_{m}\) > 10 A (\(I_{d}\) > 8 A) and \(K_{\rm s}\) < 1, the motor is seen to be in magnetic saturation state. With the increase of \(I_{m}\), the motor is more seriously affected by magnetic saturation and the value of \(K_{\rm s}\) is smaller.