Induction motor model in rotor flux frame with three phase sinusoidal excitation

This model simulates three phase induction motor model in rotor frame of reference.It is helpful in simulation of induction motor drives.
Electrical system equations:

Vs = Rs is + (1/wb) (dfs/dt) + wf M(pi/2) fs

where the variables is, vs, and fs are 2-dimensional space vectors and M(pi/2) represents a 90 degree space rotator namely [0 -1;1 0]

wf = wm + wr;
wr = Rr Lm/Lr (iqs / fdr)
1/wb Lr/Rr (dfdr/dt) + fdr = Lm ids
fqr = 0

The flux linkage current relations:
fs = Ls is + Lm ir
fr = Lm is + Lr ir
Ls= Lm + Lls;
Lr = Lm + Lrl

Mechanical system equations:
Te = 2H (dwm/dt) + Bm wm + Tl
where Te = Lm/Lr (iqs fdr)