18D070050_Group15.m

s = tf('s');

num = 1.0143*(1+s*4.9e-5)*(1+s*2.8e-5+(0.00013*s)^2);
den = (1+s*310e-5)*((1+s*13e-5)^2);
Cs = num/den;

bode(Cs)
h = gcr;
setoptions(h,'FreqUnits','Hz')

w = [251.3274,314.1593,376.9911,439.823,502.6548,565.4867,628.3185,753.9822,879.6459,1256.637,1759.292,2827.433,3769.911,4398.23,5026.548,5654.867,6283.185,8293.805,8796.459,10681.42,12566.37,15707.96,19980.53,21362.83,21991.15,22619.47,23876.1,25132.74, 26577.87, 28274.33, 29719.47, 30787.61, 31415.93, 33300.88, 35123.01, 37699.11, 39646.9, 40840.7, 43793.8, 46495.57, 50265.48, 56548.67, 62831.85, 69115.04 ,75398.22];
Cjw = 1.0143*(1+4.9e-5*1j*w).*(1+2.8e-5*1j*w+(1.3e-4*1j*w).^2)./((1+0.0031*1j*w).*(1+1.3e-4*1j*w).^2);
mag_jw = abs(Cjw');
phase_jw = 180/pi*angle(Cjw)';