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