CreateTFfromFreqResponse.m

pkg load control;
pkg load signal;
clear all
x = csvread('buckVfreqResp.csv');
SampleFrequency = 20000;
T = 1/SampleFrequency;
f = x(:,1);
r = x(:,2);
phase_rad = x(:,3).*(pi/180);
h = r.*cos(phase_rad) + r.*(sin(phase_rad)*i);
w = f.*(2*pi*T);
[b,a] = invfreqz(h,w,2,3)
sys = tf(b,a)
bode(sys)