１

Size: px

Start display at page:

Download "１"

ありおきおおはし
5 years ago
Views:

1 KC KC H.13

2 MATLAB GUI (Graphical User Interface)

3 3

4 [1] 2 1 [cm] 4

5 3[KHz]

6 [ms] 50[ms]

7 TV. 2 7

8 MATLAB

9 (10 x) Ar Al Ar = x ( x ) 2 x Al = (10 x) ( x ) 2 x A Al 2 Ar 2 + Al 2 = x 2 2 x + (10 x) (10 x) 2 x + (10 x) 2 = 1 9

10 Technics RP-DJ Hz 1000Hz 2000Hz

11 MATLAB MATLAB % sin_amp Ver 2.0 % (amplitude) ( ) % % Fs=8000; x=1:16000; y=zeros(24000,2); % % freq=500; t=x/fs; % % 0 left 10 right 5 center M=5; AR=M; AL=10-AR; % ( ) %AR %AL % seigen=sin(2.*pi*(x'-1)*(freq/fs)); %Left yleft=(al/sqrt(al*al+ar*ar))*seigen; %Right yright=(ar/sqrt(al*al+ar*ar))*seigen; y=[yleft yright]; 11

12 sound(y,fs); %wavwrite(y, Fs,'d: sin_amp_test1.wav'); return 1. % sin_amp2 Ver 1.1 % (amplitude) ( ) % sin_amp % center 1 % Fs=8000; x=1:16000; y=zeros(24000,2); % % freq=500; t=x/fs; % %0 10 M=2; L=randn(11*M,1); [hh,kk]=sort(l); kk=fix((kk+1)/m)-1; %2 % % seigen=sin(2.*pi*(x'-1)*(freq/fs)); 12

13 yleft1=(5/sqrt(5*5+5*5))*seigen; yright1=yleft1; % %center for n=1:11*m AR=kk(n); AL=10-AR; % AR 0 left 10 right 5 center %Left yleft=[yleft1;zeros(8000,1);(al/sqrt(al*al+ar*ar))*seigen]; %Right yright=[yright1;zeros(8000,1);(ar/sqrt(al*al+ar*ar))*seigen]; y=[yleft yright]; input('?') sound(y,fs); amp=kk(n) end return 13

14 % noise_amp Ver 2.2 % (amplitude) ( ) % % Fs=8000; x=1:32000; y=zeros(32000,2); % % % 0 left 10 right 5 center M=1/10*5; AR=M; AL=1.0-AR; % ( ) % 1/10 %AR %AL % noise=randn(12000,1); %Left yleft=(al/sqrt(al*al+ar*ar))*noise; %Right yright=(ar/sqrt(al*al+ar*ar))*noise; y=[yleft yright]; sound(y,fs); %wavwrite(y, Fs,'c: noise_amp_test1.wav'); 14

15 return 1/10 % noise_amp2 Ver 2.0 % (amplitude) ( ) % noise_amp % center 1 % Fs=8000; x=1:32000; y=zeros(32000,2); % % %0 10 M=2; L=randn(11*M,1); [hh,kk]=sort(l); kk=fix((kk+1)/m)-1; %2 %22 % % noise=randn(12000,1); yleft1=sqrt(0.5)*noise; yright1=yleft1; %center % for n=1:11*m 15

16 AR=1/10*kk(n); % AR 0 left 10 right AL=1-AR; %Left yleft=[yleft1;zeros(8000,1);(al/sqrt(al*al+ar*ar))*noise]; %Right yright=[yright1;zeros(8000,1);(ar/sqrt(al*al+ar*ar))*noise]; y=[yleft yright]; input('?') sound(y,fs); amp=kk(n) end return 16

17 % wav_amp Ver 1.0 % (amplitude) ( ) % % Fs=8000; x=1:24000; yy=zeros(24000,2); % % % 0 left 10 right 5 center M=5; AR=M; AL=10-AR; % ( ) %AR %AR % [y,fs]=wavread('mono_voice.wav'); %wav % mono %Left yleft=(al/sqrt(al*al+ar*ar))*y; %Right yright=(ar/sqrt(al*al+ar*ar))*y; yy=[yleft yright]; sound(yy,fs); %wavwrite(y, Fs,'c: wav_amp_test1.wav'); 17

18 return Good Morning wav wav % wav_amp2 Ver 1.0 % (amplitude) ( ) % wav_amp % center 1 % Fs=8000; x=1:24000; yy=zeros(24000,2); % % %0 (10*M) M=2; L=randn(11*M,1); [hh,kk]=sort(l); kk=fix((kk+1)/m)-1; %2 %22 % % [y,fs]=wavread('mono_voice.wav'); %wav % mono yleft1=(5/sqrt(5*5+5*5))*y; yright1=yleft1; %center % 18

19 for n=1:11*m AR=kk(n); AL=10-AR; % AR 0 left 10 right 5 center %Left yleft=[yleft1;zeros(8000,1);(al/sqrt(al*al+ar*ar))*y]; %Right yright=[yright1;zeros(8000,1);(ar/sqrt(al*al+ar*ar))*y]; yy=[yleft yright]; input('?') sound(yy,fs); amp=kk(n) end return 19

20 R L

21 Hz X= (-3 4)+(-2 4)+(-1 2) / Hz [ ] ( )

22 2 1000Hz [ ] ( ) Hz [ ] ( )

23 4 [ ] ( ) [ ] ( )

24 Hz Hz 24

25 7 2000Hz 8 25

26 x=2 4 x=6 8 =

27 10 27

28 MATLAB 6. GUI (Graphical User Interface) [ms] -1 1[ms] 50[ms] -1 1[ms] 0.2[ms]

29 12 R L 29

30 Technics RP-DJ Hz 1000Hz 2000Hz [ms] 1[ms] -1[ms] 1[ms] 0.2[ms]

31 [ms] 1[ms] Hz [ ] [ms]

32 7 1000Hz [ ] [ms] Hz [ ] [ms]

33 9 [ ] [ms] [ ] [ms]

34 13 500Hz Hz 34

35 Hz 16 35

36 [ms] 0.6[ms] 1.0[ms] 500[Hz] 36

37 [ms] 50[ms] Center S 37

39 1.0[ms] 39

40 5.1 (BOSE 101M) [m] 20 40

41 41

42 [ms] 50[ms] 5[ms] 10 25[ms] 20[ms] 2[ms] 10 10[ms] 10[ms] 2[ms] 0.2[ms] [ms] 1.0[ms] 1.0[ms] good morning 0.85[sec] 1000[Hz] 0.5[sec] 3000[Hz] 0.5[sec] 1.0[ms] 0.1[ms] 20 42

44 21 44

45 Good morning 0.85[sec] [Hz] 0.5[sec] 45

46 [Hz] 0.5[sec] 6 0.3[ms] {(1 2)+(5 1)+(6 0)} [ms] 0 0.1[ms] -0.1[ms] 46

47 [ms] 6 [ ] ( ) CENTER Good morning 0.85[sec] 47

48 [ms] 7 [ ] ( ) CENTER Good morning 0.85[sec] 48

49 [ms] 8 [ ] ( ) CENTER [Hz] 0.5[sec] 49

50 [ms] 9 [ ] ( ) CENTER [Hz] 0.5[sec] 50

51 [ms] 10 [ ] ( ) CENTER [Hz] 0.5[sec] 51

52 [ms] 11 [ ] ( ) CENTER [Hz] 0.5[sec] 52

53 25 53

54 26 54

55 27 55

56 28 56

57 [KHz] 1[KHz] 3[KHz] 1.0[ms] 1.0[ms] 1.0[ms] 100[Hz] 3[KHz] 57

58 GUI (Graphical User Interface) MATLAB GUI 6.1. GUI 29 58

59 figure(1); set(gcf,'menubar','none'); addmenu=uimenu(gcf,'label','test'); doglid=uimenu(addmenu,'label','grid','callback','grid'); vl=0.75; x=1:8001; Fs=8000; t=-3025; fr=400 br=8; tm=1.0; lr=0.5; tou=0; tout=0; ra1=uicontrol('style','radiobutton','string','8000hz','position',[ ],'value',1,... 'callback',['if get(ra2,''value'') == 1,set(ra2,''value'',0);end; if get(ra3,''value'') == 1,set(ra3,''value'',0);end; if get(ra4,''value'') == 1,set(ra4,''value'',0);end; Fs=8000;t=-8000;']); 59

60 ra2=uicontrol('style','radiobutton','string','16000hz','position',[ ],... 'callback',['if get(ra1,''value'') == 1,set(ra1,''value'',0);end; if get(ra3,''value'') == 1,set(ra3,''value'',0);end; if get(ra4,''value'') == 1,set(ra4,''value'',0);end;Fs=16000;t=0;']); br1=uicontrol('style','radiobutton','string','8bit','position',[ ],'value',1,... 'callback',['if get(br2,''value'') == 1,set(br2,''value'',0);end;br=8;']); br2=uicontrol('style','radiobutton','string','16bit','position',[ ],... 'callback',['if get(br1,''value'') == 1,set(br1,''value'',0);end;br=16;']); tm1=uicontrol('style','slider','position',[ ],'min',0.01,'max',10.01,'value',1. 0,'sliderstep',[0.01,0.1],...'callback','tm=get(tm1,''value'');set(time,''string'',num2str(t m));'); f2=uicontrol('style','edit','string',num2str(fr),'backgroundcolor',[1 1 1],'position',[ ],...'callback','fr=eval(get(f2,''string''));'); time=uicontrol('style','edit','string',num2str(tm),'backgroundcolor',[1 1 1],'position',[ ],...'callback','tm=eval(get(time,''string''));'); v1=uicontrol('style','slider','position',[ ],'Min',0,'Max',1,'Value',0.75,... 'callback','vl=get(v1,''value'');'); 60

61 lr1=uicontrol('style','slider','position',[150,125,100,15],'min',0,'max',1,'value',0.5,'slid erstep',[0.01,0.1],...'callback','lr=get(lr1,''value'');set(lr2,''string'',num2str(lr));'); lr2=uicontrol('style','edit','string',num2str(lr),'backgroundcolor',[1 1 1],'position',[ ],...'callback','lr=eval(get(lr1,''string''));'); toutime1=uicontrol('style','slider','position',[150,340,50,15],'min',-10,'max',10,'value',0,'sliderstep',[0.05,0.1],... 'callback','tout=get(toutime1,''value'');set(toutime2,''string'',num2str(tout));tou=tout *(Fs/1000);set(tousam,''string'',num2str(tou));'); toutime2=uicontrol('style','edit','string',num2str(tout),'backgroundcolor',[1,1,1],'posit ion',[ ],... 'callback','tout=eval(get(toutime2,''string''));tou=tout*(fs/1000);set(tousam,''string'', num2str(tou));'); tousam=uicontrol('style','edit','string',num2str(tou),'backgroundcolor',[1,1,1],'positio n', [ ],...'callback','tou=eval(get(tousam,''string''));'); tx1=uicontrol('style','text','string','vol','position',[20,360,20,10]); tx10=uicontrol('style','text','string','sample','position',[210,321,40,15]); so2=uicontrol('string','noise playback','position',[ ],'callback',... 'y=(vl-0.4)*randn(((16000+t)*tm),1);sound(y,fs,br);'); so3=uicontrol('string','sinwave playback','position',[ ],'callback',... 'x=1:fs*tm;y=vl*sin(2*pi*(x-1)/(fs/fr));sound(y,fs,br);'); 61

62 ww=uicontrol('string','save this sound','position',[ ],'callback',... 'wavwrite(y,fs,br,''newsound'');'); stes1=uicontrol('string','sinwave playback','position',[ ],'callback',... 'x=1:fs*tm;y=vl*[((1-lr)*(sin(2*pi*(x-1)/(fs/fr))));(lr*(sin(2*pi*(x-1)/(fs/fr))))]'';sound (y,fs,br);'); sten1=uicontrol('string','noize playback','position',[ ],'callback',... 'no=(vl-0.4)*randn(((16000+t)*tm),1);y=[(1-lr)*no lr*no];sound(y,fs,br);'); stes2=uicontrol('string','sinwave playback','position',[ ],'callback',... 'q=zeros(1,(fs/100));q1=zeros(1,((fs/100)-tou));q2=zeros(1,((fs/100)+tou));xa0=1:fs* tm;xa1=[q1,xa0,q2];xa2=[q,xa0,q];y=vl*[(0.5*(sin(2*pi*(xa1-1)/(fs/fr))));(0.5*(sin(2* pi*(xa2-1)/(fs/fr))))]'';sound(y,fs,br);'); sten2=uicontrol('string','noise playback','position',[ ],'callback',... 'q=zeros(1,(fs/100));q1=zeros(1,((fs/100)-tou));q2=zeros(1,((fs/100)+tou));no=(vl-0.4 )*randn(((16000+t)*tm),1);xa1=[q1 no q2];xa2=[q no q];y=[xa1,xa2]'';sound(y,fs,br);'); pl=uicontrol('string','plot this sound','position',[ ],'callback',... 'figure(2);plot(y);'); 62

63 ( ) time C,l or r/2,l L L/2 C user ms [ms] 63

64 yl=wavread('voice00m.wav'); fr=11025; tl01=0; tl02=0; tl03=0; tl04=0; tl05=0; tl06=0; tl07=0; tl08=0; tl09=0; tl10=0; tr01=0; tr02=0; tr03=0; tr04=0; tr05=0; tr06=0; tr07=0; tr08=0; tr09=0; tr10=0; jl51=0; jr51=0; jl80=0; jr80=0; 64

65 ce=uicontrol('string','c','position',[ ],'callback',... 'y=[yl yl];sound(y,fr,16);'); l01=uicontrol('string','0.1','position',[ ],'callback',... 'ms =-0.1; dl=round(fr/(1/(0.001*abs(ms)))); d2=zeros(1,dl); yl2=[d2 yl'']; yr2=[yl'' d2]; y=[yr2;yl2]''; sound(y,fr,16); tl01=tl01+1;set(el01,''string'',num2str(tl01));'); r01=uicontrol('string','0.1','position',[ ],'callback',... 'ms = 0.1; dl=round(fr/(1/(0.001*abs(ms))));d2=zeros(1,dl);yl2=[d2 yl''];yr2=[yl'' d2]; y=[yl2;yr2]''; sound(y,fr,16); tr01=tr01+1;set(er01,''string'',num2str(tr01));'); el01=uicontrol('style','edit','backgroundcolor',[1 1 1],'position',[ ]); er01=uicontrol('style','edit','backgroundcolor',[1 1 1],'position',[ ]); bl01=uicontrol('string','c','position',[ ],'callback',... 'jl51=jl51+1;set(el51,''string'',num2str(jl51));'); bl30=uicontrol('string','l','position',[ ],'callback',... 'jl80=jl80+1;set(el80,''string'',num2str(jl80));'); % br01 br30 % editbox el51=uicontrol('style','edit','backgroundcolor',[1 1 1],'position',[ ]); 65

66 el80=uicontrol('style','edit','backgroundcolor',[1 1 1],'position',[ ]); % ez01=uicontrol('style','edit','backgroundcolor',[1 1 1],'position',[ ]); tz01=uicontrol('style','text','string','ms','position',[52,100,20,15]); bz01=uicontrol('string','user','position',[ ],'callback',... 'ms=eval(get(ez01,''string''));dl=round(fr/(1/(0.001*abs(ms))));d2=zeros(1,dl);yl2=[d 2 yl''];yr2=[yl'' d2]; if ms < 0,y=[yr2;yl2]''; else y=[yl2;yr2]''; end; sound(y,fr,16);'); v1=uicontrol('style','slider','position',[ ],'Min',0.2,'Max',2,'Value',1,... 'callback','vl=get(v1,''value''); yl = yl*vl;'); el21=uicontrol('style','edit','backgroundcolor',[1 1 1],'position',[ ]); el40=uicontrol('style','edit','backgroundcolor',[1 1 1],'position',[ ]); % er21 er40 % 66

67 2 3000[Hz] 1.0[ms]

68 [1], pp ,

TSP信号を用いた音響系評価の研究

TSP信号を用いた音響系評価の研究 1 TSP 98kc068 2 1. 4 2. TSP 2.1 2.2 TSP 2.2.1 ATSP 2.2.2 OATSP 2.3 N 3. 5 5 6 6 7 8 2.3.1 N 8 2.3.2 TSP 9 2.3.3 2.3.4 m 3.1 3.1.1 3.1.2 3.2 3.3 15 18 20 20 20 21 22 24 3.3.1 24 3.3.2 3.4 3.4.1 3.4.2 4.