Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- %WORKING WITH HILBERT Data
- %Import Hilber
- fileID = 'vowdata.txt';
- formatSpec = '%s %4.1f %4.1f %4.1f %4.1f %4.1f %4.1f%4.1f %4.1f %4.1f %4.1f %4.1f %4.1f %4.1f %4.1f%4.1f';
- [filenames,dur,F0s,F1s,F2s,F3s,F4s,F120,F220,F320,F150,F250,F350,F180,F280,F380] = textread(fileID,formatSpec);
- filenames=char(filenames);
- %vowel (ae="had", ah="hod", aw="hawed", eh="head", er="heard",
- % ei="haid", ih="hid", iy="heed", oa=/o/ as in "boat",
- % oo="hood", uh="hud", uw="who'd")
- %Config
- vowel_codes = char('ae','ah','aw','eh','er','ei','ih','iy','oa','oo','uh','uw');
- vowels = char('had','hod','hawed','head','head','heard','haid','hid','heed','boat','hood','hud','whod');
- gender = 'm';
- figure;
- hold on;
- sd = 1;
- color = 'r'
- for ivow=1:12;
- vowel = vowels(ivow,:)
- vowel_code = vowel_codes(ivow,:)
- ind = find(filenames(:,1) == gender & filenames(:,end-1) == vowel_code(1) & filenames(:,end) == vowel_code(2)); % had
- plot( F150( ind,: ), F250( ind,: ), ['.',color], 'MarkerSize', 2 ); % red dots, size 12
- text( mean(F150( ind,: )), mean(F250( ind,: )), vowel);
- z = ellipse( F150( ind,: ), F250( ind,: ), sd);
- plot( z(:,1), z(:,2) , ['-',color]);
- end
- %WORKING ON OWN Data
- %Importing
- fileID = 'vowel_forments_data.txt';
- formatSpec = '%s\t%.6f\t%.6f\t%.6f';
- [ad_vowels,ad_F1,ad_F2,ad_F3] = textread(fileID,formatSpec);
- ad_vowels=char(ad_vowels);
- %Config
- color = 'b';
- ad_F1 = ad_F1 * 1000;
- ad_F2 = ad_F2 * 1000;
- ad_F3 = ad_F3 * 1000;
- for ivow=1:12;
- vowel = ad_vowels(ivow,:);
- disp(vowel);
- plot( ad_F1(ivow,:), ad_F2(ivow,:), ['.',color], 'MarkerSize', 30 );
- text( ad_F1(ivow,:), ad_F2(ivow,:), vowel);
- end
- set(gca, 'XScale', 'log', 'YScale', 'log');
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement