%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % ZZS - 12.cvičení: Modulační syntézy % R.Čmejla, 8.prosince 2014 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %% http://amber.feld.cvut.cz/vyu/zzs/zzs12/cv12.htm %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % Příklad 1: Amplitudová a kruhová modulace fs = 8000; Ts = 1/fs; Tmax = 1.0; t = [0:Ts:Tmax-Ts]; x = 0.5*cos(2*pi*750*t) + 0.5*cos(2*pi*450*t); x = sin(2*pi*600*t).*sin(2*pi*150*t); %% x = (1 + 1*sin(2*pi*150*t)).*sin(2*pi*600*t); %% mi = t/Tmax; % linearni narust hloubky modulace od 0 do 1 x = (1 + mi.*sin(2*pi*150*t)).*sin(2*pi*600*t); %% mi = t/Tmax % linearni narust hloubky modulace od 0 do 1 x = (1 + mi.*sin(2*pi*4*t)).*sin(2*pi*400*t); %% mi=0.5; x = (1 + mi.*sin(2*pi*4*t)).*sin(2*pi*400*t).*exp(-t); % Příklad 2: Frekvenční modulace % Parametry doba = 2; % trvani signalu [s] fs = 8000; % vzorkovani signalu [Hz] fc = 500; % nosna frekvence (Hz) fm = 50; % modulacni frekvence (Hz) % 2.5 5 10 25 50 100 250 % 4.05 8.09 16.2 40.5 80.9 162 405 mi = 1; % index modulace t = 0 : 1/fs : doba-1/fs; % casovy vektor % Příklad 3: Syntetické zvony % ZVONY fc = [110,220,110,110,250,250]; % nosna frekvence [Hz] fm = [220,440,220,220,350,350]; % modulacni frekvence [Hz] % pomer 1:2 Io = [10,5,10,10,5,3]; % modulacni index tau = [2,2,12,0.3,2,1]; % casova konstanta doba = [6,6,3,3,5,5]; % trvani [s] fs = 4.*[11025,11025,11025,11025,11025,11025]; % vzorkovaci frekvence [Hz] % Příklad 4: Syntetický vítr % VITR %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % od tonu k sumu % vetsiho "rozkmitu" (meluziny) dosahneme % 1. zvetsenim modulacniho indexu % 2. zvysenim poctu frekvencnich slozek signalu, napr. na 500 fmin = 261.6; % c1 fmax = 523.3; % c2 fc = linspace((fmin),(fmax),100); fm = 0.8; % modulacni frekvence [Hz] doba = 4; % trvani [s] fs = 8000; % vzorkovaci frekvence [Hz] t = 0:1/fs:doba-1/fs; % casova osa mi = 1; % index modulace % Příklad 5: Jednoduché syntetické FM nástroje % délka fc H =fm/fc Imax Imin % ZVON 15 200 1.4 10 0 % BICÍ 0.2 80 0.688 25 0 % KLARINET 0.5 900 0.667 4 2 % FAGOT 0.5 1100 0.2 1.5 0 % ŽESTĚ 0.6 440 1.0 5 0 % Příklad 6: Výstrahy fc = [1500,1500,1500,1500,1500]; % nosna frekvence [Hz] fm = [1.25,2,5,5,10,25]; % modulacni frekvence [Hz] Io = [400,200,100,50,20]; % modulacni index doba = [2,2,2,2,2]; % trvani [s] fs = [44100,44100,44100,44100,44100]; % vzorkovaci frekvence [Hz] % Příklad 7: Synteticke kapky vody fs = 44100; doba1 = 0.02; % doba trvani prvni casti doba2 = 0.18; % doba trvani druhe casti t1 = 0:1/fs:doba1-1/fs; t2 = 0:1/fs:doba2-1/fs; t = 0:1/fs:doba1+doba2-1/fs; fm = 7; mi = 20; tau = 0.02; pocet_kapek = 5; % Příklad 7: Vibráto fs = 22050; f0 = 440; fm = 5; doba = 5; t = 0:1/fs:doba-1/fs; alfa = [1000 300 .1 .1 .1 3 .5 1 4]; % vibrato s aditivni syntezou x = zeros(1,length(t)); for n = 1:9 x = x + alfa(n)*sin(2*pi*f0*n.*t + 7*sqrt(n)*sin(2*pi*fm.*t)); end % obalka X = [0 .2 .8 .95 1].*t(end); Y = [0 1 .3 .1 0]; O=interp1(X,Y,t);