сигналу. Для початку було проведено декілька досліджень роботи алгоритму виділення обвідної для типових сигналів: синусоїда, радіоімпульс, амплітудно-модульоване коливання (АМК) і АМК з домішкою шуму в модулююча коливання. p align="justify"> Програма виділення обвідної гармонійного сигналу:
s = zeros (1,256); (1:256) = cos (pi/8 * (1:256)); = hilbert (s); = imag (sa); = real (sa ); = sqrt ((Xi. ^ 2) + (Xr. ^ 2)); = (-128:127)/128; (2,1,1); (f, Xr); (2,1,2 ); (f, abs (A)); YLIM ([0 2]);
В
Рис. 2.1 Гармонійний сигнал
В
Рис. 2.2 Що огинає гармонійного сигналу
Програма виділення обвідної радіоімпульсу:
s = zeros (1,256);
s (64:196) = cos (pi/8 * (64:196)); = hilbert (s); = imag (sa); = real (sa); = sqrt (( Xi. ^ 2) + (Xr. ^ 2)); = (-128:127)/128; (2,1,1); (f, Xr); (2,1,2); (f, abs (A)); YLIM ([0 2]);
В
Рис. 2.3 Радіоімпульс
В
Рис. 2.4 Що огинає радиоимпульса
Програма виділення обвідної АМК:
s = zeros (1,256);
s (1:256) = (1 +0.5 * cos (pi/64 * (1:256))). * cos (pi/4 * (1:256)); = hilbert (s); = imag (sa); = real (sa); = sqrt ((Xi. ^ 2) + (Xr. ^ 2)); = (-128:127)/128; (2,1,1 ); (f, Xr); (2,1,2); (f, abs (A)); YLIM ([0 2]);
В
Рис. 2.5 Амплітудно-модульоване коливання
В
Рис. 2.6 Що огинає АМК
В
Рис. 2.7 Амплітудно-модульоване коливання з іншого несучою частотою
В
Рис. 2.8. Що огинає АМК
Як видно з останніх графіків частота несучого коливання не впливає на вигляд обвідної сигналу.
Програма виділення обвідної АМК з домішкою шуму в модулююча коливання:
s = zeros (1,256), = 1 * rand (1); = +1 * rand (1); (1:256) = ((B * 0.2 * sin (pi/8 * ( 1:256))) + (A * 0.101 * cos (pi/4 * (1:256))) +
+ (1 +0.5 * cos (pi/64 * (1:256)))). * cos (pi/4 * (1:256)); = hilbert (s); = imag (sa); = (-128:127)/128; (2,1,1); (f, Xr); (2,1,2); (f, abs (A)); YLIM ([0 2]);
В
Рис. 2.9 АМК з домішкою шуму в модулююча коливання
В
Рис. 2.10 Що огинає АМК з домішкою шуму в модулююча коливання
.2.2 Програми виділення обвідної одиночних звуків
Проведемо аналіз простих одиночних звуків (глассних) на прикладі звуків [а] [і] [е]. За отриманими огинають цих звуків можна судити про такі перехідних процесах, як атака і згасання звуку. p align="justify"> Програма виділення обвідної звуку [a]:
s = WAVREAD ('a.wav'); = hilbert (s); = imag (sa); = real (sa); = sqrt ((Xi. ^ 2) + (Xr. ^ 2)); = 12000; = 8;
[b, a] = ellip (4, 1, 80, F0/Fs * 2); = length (s) = filter (b, a, abs (A)); = (- q/2: (q/2) -1)./q, (3,1,1), (s), XLIM ([80000 140000]), (3,1,2), (abs (y)) ; XLIM ([80000 140000]); (3,1,3); (f, abs (fftshift (fft (s)))); XLIM ([-0.03 0.03]);
В
Рис. 2.11 Тимчасова реалізація звуку [a]
В
Рис. 2.12 Що огинає звуку [a]
В
Рис. 2.13 Нормований спектр звуку [a]
Програма виділення обвідної звуку [і]:
s = WAVREAD ('i.wav'); = hilbert (s); = imag (sa); = real (sa); = sqrt ((Xi. ^ 2) + (Xr. ^ 2)); = 12000; = 8;
[b, a] = ellip (4, 1, 80, F0/Fs * 2); = length (s) = filter (b, a, abs (A)); = (- q/2: (q/2) -1)./q; (3,1,1); (s); XLIM ([50000 110000]); (3,1,2); (abs (y)) ; XLIM ([50000 110000]); (3,1,3);
В
рис.2.14. Тимчасова реалізація звуку [і]
В
Рис. 2.15 Що огинає звуку [і]
В
Рис. 2.16 Нормований спектр звуку [і]
Програма виділення обвідної звуку [е]:
s = WAVREAD ('e.wav'); = hilbert (s); = imag (sa); = real (sa); = sqrt ((Xi. ^ 2) + (Xr. ^ 2)); = 12000; = 8;
[b, a] = ellip (4, 1, 80, F0/Fs * 2); = length (s) = filter (b, a, abs (A)), = (- q/2: (q/2) -1)./q; (3,1,1); (s); XLIM ([39000 65000]); (3,1,2); (abs (y)) ; XLIM ([39000 65000]); (3,1,3); (f, abs (fftshift (fft (s)))); XLIM ([-0.02 0.02]);
В
Рис. 2.17 Тимчасова реалізація звуку [е]
В
Рис. 2.18 Що огинає звуку [е]
В
Рис. 2.19 Нормований спектр звуку [е]
.2.3 Програма виділення обвідної складних звуків
При аналізі обвідної складних звуків використовув...
