To check out this repository please hg clone the following URL, or open the URL using EasyMercurial or your preferred Mercurial client.
The primary repository for this project is hosted at git://github.com/rmeddis/MAP.git .
This repository is a read-only copy which is updated automatically every hour.
root / multithreshold 1.46 / plotProfile.m @ 31:c54a34161e4a
History | View | Annotate | Download (3.92 KB)
| 1 |
function plotProfile(longTone,shortTone,gaps,BFs,TMC,offBFs,IFMCs) |
|---|---|
| 2 |
|
| 3 |
%% plot profile |
| 4 |
if nargin<1 |
| 5 |
load profile |
| 6 |
end |
| 7 |
|
| 8 |
% comparison data (e.g. participants) |
| 9 |
% rows are BFs |
| 10 |
|
| 11 |
% -------------------------------------------JSan |
| 12 |
compareBFs=[250 500 1000 2000 3000 |
| 13 |
]; |
| 14 |
compareLongTone= [67.46485 56.95655 65.01985 61.46655 73.33265 |
| 15 |
]; |
| 16 |
compareShortTone=[ 72.3185 63.2818 69.0373 65.2853 76 |
| 17 |
]; |
| 18 |
|
| 19 |
compareGaps=[0.02 0.05 0.08]; |
| 20 |
compareTMC= [ |
| 21 |
84 69 77 75 93 |
| 22 |
88 73 81 79 97 |
| 23 |
95 79 85 83 98 |
| 24 |
]; |
| 25 |
|
| 26 |
compareMaskerFreqs=[0.7 0.9 1 1.1 1.3 ]; |
| 27 |
compareIFMCs=[83.1698 77.3165 79.8474 82.9074 82.3294 |
| 28 |
80.9667 73.6653 80.9446 80.7005 79.0022 |
| 29 |
82.0135 71.2284 78.7345 74.3342 84.126 |
| 30 |
79.3348 70.3347 78.5769 79.7274 91.9849 |
| 31 |
79.2308 76.3164 83.6881 86.2105 NaN |
| 32 |
]; |
| 33 |
|
| 34 |
% % -------------------------------------------JE |
| 35 |
% BFs=[250 500 1000 2000 4000 8000]; |
| 36 |
% compareLongTone= [32 30 31 40 54 NaN]; |
| 37 |
% compareShortTone=[ 49 50 47 56 63 NaN]; |
| 38 |
% |
| 39 |
% compareGaps=0.01:0.01:0.09; |
| 40 |
% compareTMC= [ |
| 41 |
% 69 83 82 NaN NaN NaN NaN NaN NaN |
| 42 |
% 61 68 79 88 93 NaN NaN NaN NaN |
| 43 |
% 63 69 79 84 92 NaN NaN NaN NaN |
| 44 |
% 67 71 75 80 82 84 88 93 NaN |
| 45 |
% 82 82 86 86 NaN 83 88 90 75 |
| 46 |
% NaN NaN NaN NaN NaN NaN NaN NaN NaN |
| 47 |
% ]; |
| 48 |
% compareTMC=compareTMC'; |
| 49 |
% |
| 50 |
% compareMaskerFreqs=[0.5 0.7 0.9 1 1.1 1.3 1.6]; |
| 51 |
% compareIFMCs=[ |
| 52 |
% 64 60 58 58 57 60 64 |
| 53 |
% 65 63 58 55 54 59 69 |
| 54 |
% 68 64 60 59 62 73 79 |
| 55 |
% 76 75 71 67 68 71 77 |
| 56 |
% 79 71 68 69 73 75 77 |
| 57 |
% 76 73 75 75 76 80 NaN |
| 58 |
% ]; |
| 59 |
% compareIFMCs=compareIFMCs'; |
| 60 |
|
| 61 |
% -------------------------------------------CMR |
| 62 |
% CMR |
| 63 |
% BFs=[250 500 1000 2000 4000 8000]; |
| 64 |
% compareLongTone=[ 11.4 1.55 -13.5 -6.35 -6.4 7.45]; |
| 65 |
% compareShortTone=[ 23.85 18.9 9.85 10.6 9.55 21.9]; |
| 66 |
% |
| 67 |
% compareGaps=0.01:0.01:0.09; |
| 68 |
% compareTMC= [ |
| 69 |
% 28.5 35.0 49.3 70.1 80.5 85.5 NaN NaN NaN; |
| 70 |
% 31.1 44.3 48.4 59.5 56.4 76.7 70.2 82.4 76.3; |
| 71 |
% 33.4 38.4 48.8 55.8 64.5 78.7 84.2 88.3 90.3; |
| 72 |
% 25.4 37.0 49.2 49.7 58.2 69.6 87.7 95.8 93.0; |
| 73 |
% 18.2 23.5 27.4 41.5 64.3 82.1 86.7 91.2 NaN; |
| 74 |
% 32.5 35.8 43.5 52.1 69.1 78.6 86.6 86.0 NaN; |
| 75 |
% ]; |
| 76 |
% compareTMC=compareTMC'; |
| 77 |
% |
| 78 |
% compareMaskerFreqs=[0.5 0.7 0.9 1 1.1 1.3 1.6]; |
| 79 |
% compareIFMCs=[ |
| 80 |
% 50 42 34 35 34 33 37; |
| 81 |
% 58 51 38 33 28 41 49; |
| 82 |
% 57 41 27 20 28 37 66; |
| 83 |
% 61 49 27 20 34 68 79; |
| 84 |
% 67 45 27 22 46 74 87; |
| 85 |
% 62 62 43 22 47 56 83; |
| 86 |
% ]; |
| 87 |
% compareIFMCs=compareIFMCs'; |
| 88 |
|
| 89 |
% absolute thresholds |
| 90 |
figure(90), clf |
| 91 |
subplot(2,1,2) |
| 92 |
semilogx(BFs,longTone,'ko-','lineWidth',2); hold on |
| 93 |
semilogx(BFs,shortTone,'bo-','lineWidth',2); hold on |
| 94 |
semilogx(compareBFs,compareLongTone,'ko:'); hold on |
| 95 |
semilogx(compareBFs,compareShortTone,'bo:'); hold on |
| 96 |
ylim([0 100]) |
| 97 |
|
| 98 |
% TMC |
| 99 |
for BFno=1:length(BFs) |
| 100 |
subplot(2,6,BFno) |
| 101 |
plot(gaps,TMC(:,BFno)-longTone(BFno),'r','lineWidth',3), hold on |
| 102 |
plot(gaps,TMC(:,BFno),'b','lineWidth',3), hold on |
| 103 |
ylim([-10 110]) |
| 104 |
xlim([0.01 0.1]) |
| 105 |
grid on |
| 106 |
if BFno==1 |
| 107 |
ylabel('masker dB SL')
|
| 108 |
xlabel('gap')
|
| 109 |
% text(0.02,80,' TMC','backgroundColor','w') |
| 110 |
end |
| 111 |
title([num2str(BFs(BFno)) ' Hz']) |
| 112 |
set(gca,'XTick',[ 0.1],'xTickLabel', { '0.1'})
|
| 113 |
end |
| 114 |
|
| 115 |
% IFMCs |
| 116 |
for BFno=1:length(BFs) |
| 117 |
freq=offBFs'*BFs(BFno); |
| 118 |
subplot(2,1,2) |
| 119 |
semilogx(freq,IFMCs(:,BFno),'r','lineWidth',3), hold on |
| 120 |
ylim([0 100]) |
| 121 |
xlim([100 12000]) |
| 122 |
grid on |
| 123 |
end |
| 124 |
xlabel('frequency (Hz)')
|
| 125 |
ylabel('masker dB / probe dB')
|
| 126 |
set(gca,'XTick',BFs) |
| 127 |
|
| 128 |
for BFno=1:length(compareBFs) |
| 129 |
subplot(2,6,BFno) |
| 130 |
plot(compareGaps,compareTMC(:,BFno)-longTone(BFno),'k:') |
| 131 |
plot(compareGaps,compareTMC(:,BFno),'k:') |
| 132 |
ylim([-10 110]) |
| 133 |
xlim([0.01 0.1]) |
| 134 |
grid on |
| 135 |
if BFno==1 |
| 136 |
ylabel('masker dB SL')
|
| 137 |
xlabel('gap')
|
| 138 |
% text(0.02,80,' TMC','backgroundColor','w') |
| 139 |
end |
| 140 |
title([num2str(BFs(BFno)) ' Hz']) |
| 141 |
set(gca,'XTick',[ 0.1],'xTickLabel', { '0.1'})
|
| 142 |
end |
| 143 |
|
| 144 |
% IFMCs |
| 145 |
for BFno=1:length(compareBFs) |
| 146 |
compareFreq=compareMaskerFreqs'*BFs(BFno); |
| 147 |
subplot(2,1,2) |
| 148 |
semilogx(compareFreq,compareIFMCs(:,BFno),'k:') |
| 149 |
ylim([0 100]) |
| 150 |
xlim([100 12000]) |
| 151 |
grid on |
| 152 |
end |