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