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.

Statistics Download as Zip
| Branch: | Revision:

root / multithreshold 1.46 / plotProfile.m @ 29:b51bf546ca3f

History | View | Annotate | Download (3.92 KB)

1 28:02aa9826efe0 rmeddis
function plotProfile(longTone,shortTone,gaps,BFs,TMC,offBFs,IFMCs)
2
3
%% plot profile
4
if nargin<1
5
    load profile
6
end
7
8 29:b51bf546ca3f rmeddis
% comparison data (e.g. participants)
9
% rows are BFs
10 28:02aa9826efe0 rmeddis
11 29:b51bf546ca3f rmeddis
% -------------------------------------------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 28:02aa9826efe0 rmeddis
19 29:b51bf546ca3f rmeddis
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 28:02aa9826efe0 rmeddis
];
25 29:b51bf546ca3f rmeddis
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 28:02aa9826efe0 rmeddis
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 29:b51bf546ca3f rmeddis
semilogx(compareBFs,compareLongTone,'ko:'); hold on
95
semilogx(compareBFs,compareShortTone,'bo:'); hold on
96 28:02aa9826efe0 rmeddis
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 29:b51bf546ca3f rmeddis
    plot(gaps,TMC(:,BFno),'b','lineWidth',3), hold on
103
    ylim([-10 110])
104 28:02aa9826efe0 rmeddis
    xlim([0.01 0.1])
105 29:b51bf546ca3f rmeddis
    grid on
106 28:02aa9826efe0 rmeddis
    if BFno==1
107
        ylabel('masker dB SL')
108
        xlabel('gap')
109 29:b51bf546ca3f rmeddis
%         text(0.02,80,' TMC','backgroundColor','w')
110 28:02aa9826efe0 rmeddis
    end
111 29:b51bf546ca3f rmeddis
    title([num2str(BFs(BFno)) ' Hz'])
112
    set(gca,'XTick',[ 0.1],'xTickLabel', { '0.1'})
113 28:02aa9826efe0 rmeddis
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 29:b51bf546ca3f rmeddis
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