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