FDM.rar_FDM modulation_audio am matlab_audio filter_fdm_matlab F

% Modulaci鏮 anal鏬ica de informaci鏮 de audio: %archivo original: close all %-------------------------------------------------------------------------- %Archivo1 [y_original,Fs,Nbits,Format] = wavread('here put a wav song'); % sound(y_original,Fs) %extraemos 1 solo canal y a鎙dimos un DC de 0.3: canal_izquierdo_original=y_original(:,1)+0.3; BB = canal_izquierdo_original; % Filtramos para evitar posibles complicaciones en aliassing: Fmax=Fs/2; orden = 40; wn = 5000/ Fmax; % al variar, dio 5000 que nos deja oir sin alteraciones b = fir1(orden,wn,'low'); BB_filt=filter(b,1,BB); informacion = BB; inf_size = size(informacion) % Portadora: fc= 4000; t= 0:1/Fs:(length(canal_izquierdo_original)-1)/Fs; portadora = cos(2*pi*fc*t); % La portadora. port_size= size(portadora) % 1x151 % MODULACI粍 EN AM: MULTIPLICACI粍 DE DOS SE哻LES. onda_modulada = informacion .* portadora'; %transpuesta para que sean de la misma dimensi鏮 Fc_10 = 4000; % Carrier frequency in Hz mod_10 = onda_modulada; % Para los filtros pasa banda: Fmax=Fs/2; orden = 20; % Filtramos espectro INFERIOR desde Fc_10-5000 hasta Fc_10: wn_inf_10 = [Fc_10-2000 Fc_10]/ Fmax; % Filtramos espectro SUPERIOR desde Fc_10 hasta Fc_10 + 5000: wn_sup_10 = [Fc_10 Fc_10+2000]/ Fmax; % Filtramos espectro INFERIOR M糜 SUPERIOR desde Fc_10-5000 hasta Fc_10 + 5000: wn_inf_sup_10 = [Fc_10-2000 Fc_10+2000]/ Fmax; % Coeficientes: [b_inf_10,a_inf_10]= butter(orden,wn_inf_10); [b_sup_10,a_sup_10]= butter(orden,wn_sup_10); [b_inf_sup_10,a_inf_sup_10]= butter(orden,wn_inf_sup_10); % filtramos: mod_inf_filt_10=filter(b_inf_10,a_inf_10,onda_modulada); mod_sup_filt_10=filter(b_sup_10,a_sup_10,onda_modulada); mod_inf_sup_filt_10=filter(b_inf_sup_10,a_inf_sup_10,onda_modulada); % amplificamos: mod_inf_filt_amp_10=2*mod_inf_filt_10; mod_sup_filt_amp_10=2*mod_sup_filt_10; mod_inf_sup_filt_amp_10=2*mod_inf_sup_filt_10; %Espectro de potencias: nFFT=1024; LongVentana= nFFT/4; Solapamiento =LongVentana/2; % Espectro de la BB: [P_welch_BB_con_mod,F_welch_BB_con_mod] = pwelch(onda_modulada,LongVentana,Solapamiento,nFFT,Fs); informacion_uno=P_welch_BB_con_mod; %-------------------------------------------------------------------------- %Archivo 2 [y_original,Fs,Nbits,Format] = wavread('Chispa'); % sound(y_original,Fs) %extraemos 1 solo canal y a鎙dimos un DC de 0.3: canal_izquierdo_original=y_original(:,1)+0.3; BB = canal_izquierdo_original; % Filtramos para evitar posibles complicaciones en aliassing: Fmax=Fs/2; orden = 40; wn = 5000/ Fmax; % al variar, dio 5000 que nos deja oir sin alteraciones b = fir1(orden,wn,'low'); BB_filt=filter(b,1,BB); informacion = BB; inf_size = size(informacion) % Portadora: fc= 8000; t= 0:1/Fs:(length(canal_izquierdo_original)-1)/Fs; portadora = cos(2*pi*fc*t); % La portadora. port_size= size(portadora) % 1x151 % MODULACI粍 EN AM: MULTIPLICACI粍 DE DOS SE哻LES. onda_modulada = informacion .* portadora'; %transpuesta para que sean de la misma dimensi鏮 Fc_10 = 8000; % Carrier frequency in Hz mod_10 = onda_modulada; % Para los filtros pasa banda: Fmax=Fs/2; orden = 20; % Filtramos espectro INFERIOR desde Fc_10-5000 hasta Fc_10: wn_inf_10 = [Fc_10-4000 Fc_10]/ Fmax; % Filtramos espectro SUPERIOR desde Fc_10 hasta Fc_10 + 5000: wn_sup_10 = [Fc_10 Fc_10+4000]/ Fmax; % Filtramos espectro INFERIOR M糜 SUPERIOR desde Fc_10-5000 hasta Fc_10 + 5000: wn_inf_sup_10 = [Fc_10-4000 Fc_10+4000]/ Fmax; % Coeficientes: [b_inf_10,a_inf_10]= butter(orden,wn_inf_10); [b_sup_10,a_sup_10]= butter(orden,wn_sup_10); [b_inf_sup_10,a_inf_sup_10]= butter(orden,wn_inf_sup_10); % filtramos: mod_inf_filt_10=filter(b_inf_10,a_inf_10,onda_modulada); mod_sup_filt_10=filter(b_sup_10,a_sup_10,onda_modulada); mod_inf_sup_filt_10=filter(b_inf_sup_10,a_inf_sup_10,onda_modulada); % amplificamos: mod_inf_filt_amp_10=2*mod_inf_filt_10; mod_sup_filt_amp_10=2*mod_sup_filt_10; mod_inf_sup_filt_amp_10=2*mod_inf_sup_filt_10; %Espectro de potencias: nFFT=1024; LongVentana= nFFT/4; Solapamiento =LongVentana/2; % Espectro de la BB: [P_welch_BB_con_mod,F_welch_BB_con_mod] = pwelch(onda_modulada,LongVentana,Solapamiento,nFFT,Fs); informacion_dos=P_welch_BB_con_mod; %--------------------------------------------------------------------------------------------------------- %Archivo 3 [y_original,Fs,Nbits,Format] = wavread('Pinguinos'); % sound(y_original,Fs) %extraemos 1 solo canal y a鎙dimos un DC de 0.3: canal_izquierdo_original=y_original(:,1)+0.3; BB = canal_izquierdo_original; % Filtramos para evitar posibles complicaciones en aliassing: Fmax=Fs/2; orden = 40; wn = 5000/ Fmax; % al variar, dio 5000 que nos deja oir sin alteraciones b = fir1(orden,wn,'low'); BB_filt=filter(b,1,BB); informacion = BB; inf_size = size(informacion) % Portadora: fc= 12000; t= 0:1/Fs:(length(canal_izquierdo_original)-1)/Fs; portadora = cos(2*pi*fc*t); % La portadora. port_size= size(portadora) % 1x151 % MODULACI粍 EN AM: MULTIPLICACI粍 DE DOS SE哻LES. onda_modulada = informacion .* portadora'; %transpuesta para que sean de la misma dimensi鏮 Fc_10 = 12000; % Carrier frequency in Hz mod_10 = onda_modulada; % Para los filtros pasa banda: Fmax=Fs/2; orden = 20; % Filtramos espectro INFERIOR desde Fc_10-5000 hasta Fc_10: wn_inf_10 = [Fc_10-6000 Fc_10]/ Fmax; % Filtramos espectro SUPERIOR desde Fc_10 hasta Fc_10 + 5000: wn_sup_10 = [Fc_10 Fc_10+6000]/ Fmax; % Filtramos espectro INFERIOR M糜 SUPERIOR desde Fc_10-5000 hasta Fc_10 + 5000: wn_inf_sup_10 = [Fc_10-6000 Fc_10+6000]/ Fmax; % Coeficientes: [b_inf_10,a_inf_10]= butter(orden,wn_inf_10); [b_sup_10,a_sup_10]= butter(orden,wn_sup_10); [b_inf_sup_10,a_inf_sup_10]= butter(orden,wn_inf_sup_10); % filtramos: mod_inf_filt_10=filter(b_inf_10,a_inf_10,onda_modulada); mod_sup_filt_10=filter(b_sup_10,a_sup_10,onda_modulada); mod_inf_sup_filt_10=filter(b_inf_sup_10,a_inf_sup_10,onda_modulada); % amplificamos: mod_inf_filt_amp_10=2*mod_inf_filt_10; mod_sup_filt_amp_10=2*mod_sup_filt_10; mod_inf_sup_filt_amp_10=2*mod_inf_sup_filt_10; %Espectro de potencias: nFFT=1024; LongVentana= nFFT/4; Solapamiento =LongVentana/2; % Espectro de la BB: [P_welch_BB_con_mod,F_welch_BB_con_mod] = pwelch(onda_modulada,LongVentana,Solapamiento,nFFT,Fs); informacion_tres=P_welch_BB_con_mod; informacion_total=informacion_uno+informacion_dos+informacion_tres; plot(F_welch_BB_con_mod,informacion_total,'b') %-------------------------------------------------------------------------- %FILTRANDO SE哻LES Y RECUPERANDO % ---------------------------------------- % % Para los filtros pasa banda: % % % filtramos: % recuperacion_music1=filter(4000,5000,informacion_total); % recuperacion_music2=filter(8000,9000,informacion_total); % recuperacion_music3=filter(12000,12500,informacion_total); % % amplificamos: % recuperacion_music1_amp=2*recuperacion_music1; % recuperacion_music2_amp=2*recuperacion_music2; % recuperacion_music3_amp=2*recuperacion_music3; % % figure % subplot(3,1,1) % plot(recuperacion_music1_amp,'b') % subplot(3,1,2) % plot(recuperacion_music2_amp,'y') % subplot(3,1,3) % plot(recuperacion_music3_amp,'k')