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: