I believe your using the wrong equations for your transfer function. The equations that you are referencing to in the above article are for levitating a soft magnetic material (such as steel). If you are using a hard magnetic material (such as a permanent magnet) the equations are going to be different. The reasoning is that the mechanism that generates force is different. The soft magnetic material's force is due to the change in inductance vs change in position. Whereas the permanent magnets flux couples with the coil and the force mechanism is the change in magnetic flux vs the change in position.
In addition, a permanent magnet has the same permeability as air, so the inductance that the coil sees does not change with respect to the position of the magnet. Steel has a permeability around 2000 times that of air. Hence the inductance changes when the object gets closer or further from the coil.
The amount of flux that couples with the coil is going to change with respect to the position of the coil you are using. I would highly recommend modeling your system in FEMM and solving for the flux linkage at several magnet positions.
I will solve out for you so you can see the difference
First perform a power balance between mechanical power and electrical power. \begin{equation} P_m=\frac{Fx}{t}\\ P_e=vi\\ P_m=P_e=\frac{Fx}{t}=vi \end{equation} Next apply faradays law.
\begin{equation} v=\frac{d\lambda}{dt}\\ \end{equation}
For a soft magnetic object (steel) total magnetic flux \$\lambda\$ is a function of the inductance. The inductance is a function of the distance from the core. \begin{equation} \lambda = L(x)i \end{equation}
Now put it all together. I believe the chain rule is used here. \begin{equation} P_m=P_e=\frac{Fx}{t}=i\frac{d\lambda}{dt}\\ \frac{Fx}{t}=i^2\frac{d(L(x))}{dt} \end{equation}\begin{equation} P_m=P_e=\frac{Fx}{t}=\frac{d\lambda}{dt}\\ \end{equation} Break down the derivative into partials. \begin{equation} \frac{Fx}{t}=i^2\frac{d(L(x))}{dx}\frac{dx}{dt} \end{equation}\begin{equation} \frac{Fx}{t}=\frac{i^2}{2}\frac{d(L(x))}{dx}\frac{dx}{dt} \end{equation} Finally canceling out the velocity term \$\frac{x}{t}\$ gives you the following relationship. \begin{equation} F=i^2\frac{d(L(x))}{dx} \end{equation}\begin{equation} F=\frac{i^2}{2}\frac{d(L(x))}{dx} \end{equation}
Now look at the units to make sure they work out to be in Newtons \begin{equation} kg*m*s^2=(A^2)(kg*m^2*s^{−2}*A^{−2})(m^{-1}) \end{equation}## Heading ##
For a hard magnet material the total magnetic flux \$\lambda\$ is a function of the magnetic flux that couples with the coil. The inductance does not change with position \$L(x) = 0\$ so we can leave it out, as it contributes zero torque. \$\lambda_m\$ is the B field normal to the coil times the area of the coil (it comes from gauss law)
\begin{equation} \lambda = \lambda_m(x) + L(x)i\\ \lambda = \lambda_m(x) = B*A \end{equation}
Now put it all together \begin{equation} P_m=P_e=\frac{Fx}{t}=i\frac{d\lambda_m}{dt} \end{equation} Break down the derivative into partials. \begin{equation} \frac{Fx}{t}=i\frac{d(\lambda_m(x))}{dx}\frac{dx}{dt} \end{equation} Finally canceling out the velocity term \$\frac{x}{t}\$ gives you the following relationship. \begin{equation} F=i\frac{d(\lambda_m(x))}{dx} \end{equation}
Now look at the units to make sure they work out to be in Newtons \begin{equation} kg*m*s^2=(A)(kg*m^2*s^{−2}*A^{−1})(m^{-1}) \end{equation}
In conclusion in order to model \$\frac{d(\lambda_m(x))}{dx}\$ you likely need to use finite element analysis. I would recommend you make a simple model using FEMM. But the biggest difference is that you will see that the current is squared in one equation and not in the other equation. This kind of assumption would make a huge difference in your model.