load spk1.dat; load spk2.dat; fl = 10; % the lowest freq.(Hz) fh = 2000; % the highest freq.(Hz) f = spk1(:,1); % frequency data w = 2*pi*f; n = size(spk1,1); % number of frequencies Gyw_jw = spk1(:,2).*exp(1i*pi/180.0*spk1(:,3)); Gzw_jw = spk1(:,4).*exp(1i*pi/180.0*spk1(:,5)); Gyu_jw = spk2(:,2).*exp(1i*pi/180.0*spk2(:,3)); Gzu_jw = spk2(:,4).*exp(1i*pi/180.0*spk2(:,5)); for k = 1:n resp(1,1,k) = Gzw_jw(k); resp(1,2,k) = Gzu_jw(k); resp(2,1,k) = Gyw_jw(k); resp(2,2,k) = Gyu_jw(k); end G_g = frd(resp, w); figure(1); bode(G_g, 'b-');