format short
format compact
echo on
% Read in Peterson-Barney Data:
fid = fopen('pb.Table1');
PB = textscan(fid, '%s %s %d %s %s %f %f %f %f', 10000);
fclose(fid);
%
% 1: m = man, w=woman, c=child
% 2: m=male, f=female
% 3: N=speaker number
% 4: V=vowel name (heed/iy, hid/ih, head/eh, had/ae, hod/aa, hawed/ao, hood/uh, who'd/uw, hud/ah, heard/er)
% 5: IPA
% 6: F0 in Hz
% 7: F1
% 8: F2
% 9: F3
%
pause
% Convert column 1 (man, woman, child) from cell array to string array:
MWC=char(PB{1});
% Make some useful binary vectors:
isman=(MWC=='m');
iswoman=(MWC=='w');
ischild=(MWC=='c');
% man/woman/child as values from 1 to 3:
MWCints = isman + 2*iswoman + 3*ischild;
% Convert column 2 (male, female) to string array:
MF=char(PB{2});
% More useful binary vectors:
ismale=(MF=='m');
isfemale=(MF=='f');
% Convert column 4 (vowel ID) to integers 1:10 in VID array
% Note that options are 'iy','ih','eh', 'ae','aa','ao','uh','uw','ah','er'
VID=zeros(length(isman),1);
VID = VID + strcmp('iy',PB{4});
VID = VID + 2*strcmp('ih',PB{4});
VID = VID + 3*strcmp('eh',PB{4});
VID = VID + 4*strcmp('ae',PB{4});
VID = VID + 5*strcmp('aa',PB{4});
VID = VID + 6*strcmp('ao',PB{4});
VID = VID + 7*strcmp('uh',PB{4});
VID = VID + 8*strcmp('uw',PB{4});
VID = VID + 9*strcmp('ah',PB{4});
VID = VID + 10*strcmp('er',PB{4});
%
pause
% Make matrix whose rows are vowels, and
%      columns are F0, F1, F2, F3:
F0123 = [PB{6} PB{7} PB{8} PB{9}];