I am attempting to model a phase winding of a Switched Reluctance Generator, it involves a nonlinear inductance. I has inputs as reference current and rotor position and an output as a control pulse. has anyone ventured in such?

the equations involved include the following: 1. LK(θ,i) = L0(i) + L1(i)∗(cos(Nrθ )+π) 2. L0(i) = (Lmax(i)+Lmin(i))/2 3. L1(i) = (Lmax(i)-Lmin(i))/2 4. Lmax(i) = 0.136 − 0.0045*i + 0.0056*i^2 − 0.0022*i^3 + 0.00035*i^4

i = current, L = inductance, Nr = number of rotor poles