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