function [kAR_hat NsNf_hat] = estimateParam(Ia_shunt,Ia_serie,If_shunt,If_eff_shunt,If_eff_serie)
% TSFS04: Elektriska drivsystem. Script till Laboration 2.
% Datum: 111207
% estimateParam skattar parametrarna kAR och Ns/Nf med
% minsta-kvadratanpassning
% Input: 
%   2 x ankarström (Ia)
%   1 x fälström (If)
%   2 x shunt-ekvivalent ström (If_eff)
%   Samtliga vektorer ska vara KOLUMNVEKTORER!
%   Modell: If_eff = If_shunt + Ns/Nf*If_series - kAR*Ia.^2
% Output:
%   Skattat kAR
%   Skattat Ns/Nf

A = [[Ia_serie; zeros(length(Ia_shunt),1)] ...
    -[Ia_serie.^2; Ia_shunt.^2]];
b = [If_eff_serie; ...
    [If_eff_shunt - If_shunt]];
xsol = A\b;
NsNf_hat = xsol(1);
kAR_hat = xsol(2);

end