[y,Fs]=wavread('alex.wav');
[L,M]=size(y);
%Fs = 1000; % Sampling frequency
T = 1/Fs; % Sample time
%L = 1000; % Length of signal
t = (0:L-1)*T; % Time vector
% Sum of a 50 Hz sinusoid and a 120 Hz sinusoid
%x = 0.7*sin(2*pi*200*t);% + sin(2*pi*120*t);
%y = x; %+ 2*randn(size(t)); % Sinusoids plus noise
figure(1)
plot(Fs*t(1:3000),y(1:3000))
title('Signal Corrupted with Zero-Mean Random Noise')
xlabel('time (milliseconds)')
Portadora=sin(2*pi*20000*t);
yf=y(:,1);
yf=yf';
modulada=Portadora.*yf;
figure(1)
plot(Fs*t(1:3000),modulada(1:3000))
title('Signal Corrupted with Zero-Mean Random Noise')
xlabel('time (milliseconds)')
y=modulada;
figure(3)
plot(modulada)
NFFT = 2^nextpow2(L); % Next power of 2 from length of y
Y = fft(y,NFFT)/L;
f = Fs/2*linspace(0,1,NFFT/2+1);
% Plot single-sided amplitude spectrum.
figure(2)
plot(f,2*abs(Y(1:NFFT/2+1)))
title('Single-Sided Amplitude Spectrum of y(t)')
xlabel('Frequency (Hz)')
ylabel('|Y(f)|')
Es un espacio donde publico notas que creo interesantes desde octubre de 2007.
domingo, 31 de julio de 2011
Código de Matlab señal de 200 Hz.
Fs = 1000; % Sampling frequency
T = 1/Fs; % Sample time
L = 1000; % Length of signal
t = (0:L-1)*T; % Time vector
% Sum of a 50 Hz sinusoid and a 120 Hz sinusoid
x = 0.7*sin(2*pi*200*t);% + sin(2*pi*120*t);
y = x; %+ 2*randn(size(t)); % Sinusoids plus noise
figure(1)
plot(Fs*t(1:50),y(1:50))
title('Signal Corrupted with Zero-Mean Random Noise')
xlabel('time (milliseconds)')
NFFT = 2^nextpow2(L); % Next power of 2 from length of y
Y = fft(y,NFFT)/L;
f = Fs/2*linspace(0,1,NFFT/2+1);
% Plot single-sided amplitude spectrum.
figure(2)
plot(f,2*abs(Y(1:NFFT/2+1)))
title('Single-Sided Amplitude Spectrum of y(t)')
xlabel('Frequency (Hz)')
ylabel('|Y(f)|')
T = 1/Fs; % Sample time
L = 1000; % Length of signal
t = (0:L-1)*T; % Time vector
% Sum of a 50 Hz sinusoid and a 120 Hz sinusoid
x = 0.7*sin(2*pi*200*t);% + sin(2*pi*120*t);
y = x; %+ 2*randn(size(t)); % Sinusoids plus noise
figure(1)
plot(Fs*t(1:50),y(1:50))
title('Signal Corrupted with Zero-Mean Random Noise')
xlabel('time (milliseconds)')
NFFT = 2^nextpow2(L); % Next power of 2 from length of y
Y = fft(y,NFFT)/L;
f = Fs/2*linspace(0,1,NFFT/2+1);
% Plot single-sided amplitude spectrum.
figure(2)
plot(f,2*abs(Y(1:NFFT/2+1)))
title('Single-Sided Amplitude Spectrum of y(t)')
xlabel('Frequency (Hz)')
ylabel('|Y(f)|')
miércoles, 27 de julio de 2011
jueves, 21 de julio de 2011
miércoles, 20 de julio de 2011
martes, 19 de julio de 2011
martes, 12 de julio de 2011
lunes, 4 de julio de 2011
Suscribirse a:
Entradas (Atom)