alexbrandmeyer@678: % Author: Alex Brandmeyer
alexbrandmeyer@678: %
alexbrandmeyer@678: % This Matlab file is part of an implementation of Lyon's cochlear model:
alexbrandmeyer@678: % "Cascade of Asymmetric Resonators with Fast-Acting Compression"
alexbrandmeyer@678: % to supplement Lyon's upcoming book "Human and Machine Hearing"
alexbrandmeyer@678: %
alexbrandmeyer@678: % Licensed under the Apache License, Version 2.0 (the "License");
alexbrandmeyer@678: % you may not use this file except in compliance with the License.
alexbrandmeyer@678: % You may obtain a copy of the License at
alexbrandmeyer@678: %
alexbrandmeyer@678: %     http://www.apache.org/licenses/LICENSE-2.0
alexbrandmeyer@678: %
alexbrandmeyer@678: % Unless required by applicable law or agreed to in writing, software
alexbrandmeyer@678: % distributed under the License is distributed on an "AS IS" BASIS,
alexbrandmeyer@678: % WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
alexbrandmeyer@678: % See the License for the specific language governing permissions and
alexbrandmeyer@678: % limitations under the License.
alexbrandmeyer@678: 
alexbrandmeyer@678: function CARFAC_Compare_CPP_Test_Data()
alexbrandmeyer@678: data_dir = '../../../carfac/test_data/';
alexbrandmeyer@678: filename = 'long_test_nap1.txt';
alexbrandmeyer@678: matlab_nap1 = dlmread([data_dir filename]);
alexbrandmeyer@678: filename = 'cpp_nap_output_1_long_test.txt';
alexbrandmeyer@678: cpp_nap1 = dlmread([data_dir filename]);
alexbrandmeyer@678: filename = 'long_test_nap2.txt';
alexbrandmeyer@678: matlab_nap2 = dlmread([data_dir filename]);
alexbrandmeyer@678: filename = 'cpp_nap_output_2_long_test.txt';
alexbrandmeyer@678: cpp_nap2 = dlmread([data_dir filename]);
alexbrandmeyer@678: 
alexbrandmeyer@678: factor = 10;
alexbrandmeyer@678: figure(1)
alexbrandmeyer@678: subplot(2,2,1)
alexbrandmeyer@678: image(matlab_nap1' * factor);
alexbrandmeyer@678: title('Long Test: Matlab NAP, Ear 1');
alexbrandmeyer@678: ylabel('Channel');
alexbrandmeyer@678: xlabel('Sample Index');
alexbrandmeyer@678: subplot(2,2,2)
alexbrandmeyer@678: image(matlab_nap2' * factor);
alexbrandmeyer@678: title('Long Test: Matlab NAP, Ear 2')
alexbrandmeyer@678: ylabel('Channel');
alexbrandmeyer@678: xlabel('Sample Index');
alexbrandmeyer@678: subplot(2,2,3)
alexbrandmeyer@678: image(cpp_nap1' * factor);
alexbrandmeyer@678: title('Long Test: C++ NAP, Ear 1')
alexbrandmeyer@678: ylabel('Channel');
alexbrandmeyer@678: xlabel('Sample Index');
alexbrandmeyer@678: subplot(2,2,4)
alexbrandmeyer@678: image(cpp_nap2' * factor);
alexbrandmeyer@678: title('Long Test: C++ NAP, Ear 2')
alexbrandmeyer@678: ylabel('Channel');
alexbrandmeyer@678: xlabel('Sample Index');
alexbrandmeyer@678: 
alexbrandmeyer@678: filename = 'binaural_test_nap1.txt';
alexbrandmeyer@678: matlab_nap1 = dlmread([data_dir filename]);
alexbrandmeyer@678: filename = 'cpp_nap_output_1_binaural_test.txt';
alexbrandmeyer@678: cpp_nap1 = dlmread([data_dir filename]);
alexbrandmeyer@678: filename = 'binaural_test_nap2.txt';
alexbrandmeyer@678: matlab_nap2 = dlmread([data_dir filename]);
alexbrandmeyer@678: filename = 'cpp_nap_output_2_binaural_test.txt';
alexbrandmeyer@678: cpp_nap2 = dlmread([data_dir filename]);
alexbrandmeyer@678: 
alexbrandmeyer@678: factor = 10;
alexbrandmeyer@678: figure(2)
alexbrandmeyer@678: subplot(2,2,1)
alexbrandmeyer@678: image(matlab_nap1' * factor);
alexbrandmeyer@678: title('Binaural Test: Matlab NAP, Ear 1');
alexbrandmeyer@678: ylabel('Channel');
alexbrandmeyer@678: xlabel('Sample Index');
alexbrandmeyer@678: subplot(2,2,2)
alexbrandmeyer@678: image(matlab_nap2' * factor);
alexbrandmeyer@678: title('Binaural Test: Matlab NAP, Ear 2')
alexbrandmeyer@678: ylabel('Channel');
alexbrandmeyer@678: xlabel('Sample Index');
alexbrandmeyer@678: subplot(2,2,3)
alexbrandmeyer@678: image(cpp_nap1' * factor);
alexbrandmeyer@678: title('Binaural Test: C++ NAP, Ear 1')
alexbrandmeyer@678: ylabel('Channel');
alexbrandmeyer@678: xlabel('Sample Index');
alexbrandmeyer@678: subplot(2,2,4)
alexbrandmeyer@678: image(cpp_nap2' * factor);
alexbrandmeyer@678: title('Binaural Test: C++ NAP, Ear 2')
alexbrandmeyer@678: ylabel('Channel');
alexbrandmeyer@678: xlabel('Sample Index');