%function [Gstack_AB_g] = load_data_Gstack_AB(dirname); %function [Gstack_AB_g] = load_data_Gstack_AB(dirname, dl); function [Gstack_AB_g] = load_data_Gstack_AB_fixed(dirname, dl, ls); % load_data_Gstack_fixed.m をコピーして編集。 %[A1_spk1_g, A2_spk1_g, B1_spk1_g, B2_spk1_g] = load_data_raw_AB(strcat(dirname, '/spk1.dat')); [A1_spk1_g, A2_spk1_g, B1_spk1_g, B2_spk1_g] = load_data_raw_AB_fixed(strcat(dirname, '/spk1.dat'), dl, ls); %[A1_spk2_g, A2_spk2_g, B1_spk2_g, B2_spk2_g] = load_data_raw_AB(strcat(dirname, '/spk2.dat')); [A1_spk2_g, A2_spk2_g, B1_spk2_g, B2_spk2_g] = load_data_raw_AB_fixed(strcat(dirname, '/spk2.dat'), dl, ls); w = A1_spk1_g.freq; n = size(w, 1); for k = 1:n tmp = [A2_spk1_g.resp(k), A2_spk2_g.resp(k); B1_spk1_g.resp(k), B1_spk2_g.resp(k)] / [A1_spk1_g.resp(k), A1_spk2_g.resp(k); B2_spk1_g.resp(k), B2_spk2_g.resp(k)]; resp(1,1,k) = tmp(1,1); resp(1,2,k) = tmp(1,2); resp(2,1,k) = tmp(2,1); resp(2,2,k) = tmp(2,2); end Gstack_AB_g = frd(resp, w);