rmeddis@9: function testPhaseLocking
rmeddis@9: 
rmeddis@9: testFrequencies=[250 500 1000 2000 4000 8000];
rmeddis@9: levels=-20:10:80;
rmeddis@9: figure(14), clf
rmeddis@9: allStrengths=zeros(length(testFrequencies), length(levels));
rmeddis@9: peakVectorStrength=zeros(1,length(testFrequencies));
rmeddis@9: freqCount=0;
rmeddis@9: for targetFrequency=testFrequencies;
rmeddis@9:     %single test
rmeddis@9:     freqCount=freqCount+1;
rmeddis@9:     vectorStrength=testAN(targetFrequency,targetFrequency, levels);
rmeddis@9:     allStrengths(freqCount,:)=vectorStrength';
rmeddis@9:     peakVectorStrength(freqCount)=max(vectorStrength');
rmeddis@9: end
rmeddis@9: %% plot results
rmeddis@9: figure(14)
rmeddis@9: subplot(2,1,2)
rmeddis@9: plot(levels,allStrengths)
rmeddis@9: xlabel('levels')
rmeddis@9: ylabel('vector strength')
rmeddis@9: legend (num2str(testFrequencies'),'location','eastOutside')
rmeddis@9: 
rmeddis@9: subplot(2,1,1)
rmeddis@9: semilogx(testFrequencies,peakVectorStrength)
rmeddis@9: grid on
rmeddis@9: title ('peak vector strength')
rmeddis@9: xlabel('frequency')
rmeddis@9: ylabel('vector strength')
rmeddis@9: 
rmeddis@9: johnson=[250	0.79
rmeddis@9: 500	0.82
rmeddis@9: 1000	0.8
rmeddis@9: 2000	0.7
rmeddis@9: 4000	0.25
rmeddis@9: 5500	0.05
rmeddis@9: ];
rmeddis@9: hold on
rmeddis@9: plot(johnson(:,1),johnson(:,2),'o')
rmeddis@9: legend({'model','Johnson 80'},'location','eastOutside')
rmeddis@9: hold off
rmeddis@9: 
rmeddis@9: