function Moisture = SingBaseline(featureIn) % Extract magnitude and phase columns freq = featureIn(:, 1); % Frequency column magnitude_column = featureIn(:, 2); % Magnitude column phase_column = featureIn(:, 3); % Phase column % Compute the complex number e from magnitude and phase e = magnitude_column .* exp(1i * phase_column); % Compute f using the formula fe = -(abs(e) .* sin(angle(e) / 2)); alpha=-0.0005144; %moisture coefficient beta=0.0007338; %suger coefficient zeta=0.0000000002159;% frequency coefficient meu=0.004684; %intercept S=100; %S=82.84; GT=380; BS=beta*S; k=-1/alpha; Moisture=k*(BS + (zeta .* freq) + meu - fe); end