emf_B = ppval (time, dpp_flux_B);
emf_C = ppval (time, dpp_flux_C);
figure (1);
plot (time * 1000, flux_link_A, 'r-');
hold on;
plot (time * 1000, flux_link_B, 'b-');
plot (time * 1000, flux_link_C, 'g-');
title ('Flux linkage - under noload');
xlabel ('Time [ms]'), ylabel ('Flux linkage [WbT]')
grid;
figure (2);
plot (time * 1000, emf_A, 'r-');
hold on;
plot (time * 1000, emf_B, 'b-');
plot (time * 1000, emf_C, 'g-');
title ('Back Emf - under noload');
xlabel ('Time [ms]'), ylabel ('Back EMF [V]')
grid;
x = length (flux_link_A);
FLux_RMS = norm (flux_link_A)/sqrt (x)
y = length (emf_A);
EMF_RMS = norm (emf_A)/sqrt (y)