Scribd Logo
Cargar

¡Te damos la bienvenida a Scribd!

  • 96 vistas

    Tutorial 3 de Matlab

    Cargado por

    JAVIER ALONSO MURILLO MURILLO
    Este documento describe cómo calcular el espectro de potencia de diferentes señales (senoidal, cuadrada y triangular) utilizando MATLAB. Explica que la transformada de Fourier descompone las señales en componentes de frecuencia pura y que el espectro de potencia muestra la distribución de potencia de la señal en función de la frecuencia. A continuación, presenta el código MATLAB para calcular y graficar el espectro de potencia de cada tipo de señal.

    Copyright:

    © All Rights Reserved
    Descargue como PDF, TXT o lea en línea desde Scribd

    Tutorial 3 de Matlab

    Cargado por

    JAVIER ALONSO MURILLO MURILLO
    96 vistas5 páginas
    Este documento describe cómo calcular el espectro de potencia de diferentes señales (senoidal, cuadrada y triangular) utilizando MATLAB. Explica que la transformada de Fourier descompone las señales en componentes de frecuencia pura y que el espectro de potencia muestra la distribución de potencia de la señal en función de la frecuencia. A continuación, presenta el código MATLAB para calcular y graficar el espectro de potencia de cada tipo de señal.

    Descripción original:

    Introducción a Matlab

    Derechos de autor

    © © All Rights Reserved

    Formatos disponibles

    PDF, TXT o lea en línea desde Scribd

    ¿Le pareció útil este documento?

    ¿Este contenido es inapropiado?

    Este documento describe cómo calcular el espectro de potencia de diferentes señales (senoidal, cuadrada y triangular) utilizando MATLAB. Explica que la transformada de Fourier descompone las señales en componentes de frecuencia pura y que el espectro de potencia muestra la distribución de potencia de la señal en función de la frecuencia. A continuación, presenta el código MATLAB para calcular y graficar el espectro de potencia de cada tipo de señal.

    Copyright:

    © All Rights Reserved
    Descargue como PDF, TXT o lea en línea desde Scribd
    Descargar como pdf o txt
    96 vistas5 páginas

    Tutorial 3 de Matlab

    Cargado por

    JAVIER ALONSO MURILLO MURILLO
    Este documento describe cómo calcular el espectro de potencia de diferentes señales (senoidal, cuadrada y triangular) utilizando MATLAB. Explica que la transformada de Fourier descompone las señales en componentes de frecuencia pura y que el espectro de potencia muestra la distribución de potencia de la señal en función de la frecuencia. A continuación, presenta el código MATLAB para calcular y graficar el espectro de potencia de cada tipo de señal.

    Copyright:

    © All Rights Reserved
    Descargue como PDF, TXT o lea en línea desde Scribd
    Descargar como pdf o txt
    de 5

    UNIVERSIDAD SANTIAGO DE CALI

    FACULTAD INGENIERIAS
    PROGRAMAS DE INGENIERIA ELECTRONICA Y BIOINGENIERIA
    SEALES Y SISTEMAS
    Profesor: MSc. JAVIER A. MURILLO M.

    VI. ESPECTRO DE POTENCIA USANDO MATLAB

    La seal ms simple o seal de


    frecuencia pura es la funcin sinusoidal.

    Sumando seales de frecuencia pura


    se puede obtener cualquier tipo de
    seal real.

    Los anlisis de Fourier (series y


    transformadas) permiten calcular las
    amplitudes, frecuencias y fases que es
    necesario dar a dichas seales puras
    para obtener cualquier otra seal.

    La siguiente figura muestra la aproximacin a una seal cuadrada de 1 kHz. Se


    han utilizado 7 componentes (C0 ..C7), aunque realmente C0=C2=C4=C6=0 (para
    una seal cuadrada, el valor de los componentes Cn=4/n/, si n impar y Cn=0 si n
    par).

    41
    UNIVERSIDAD SANTIAGO DE CALI
    FACULTAD INGENIERIAS
    PROGRAMAS DE INGENIERIA ELECTRONICA Y BIOINGENIERIA
    SEALES Y SISTEMAS
    Profesor: MSc. JAVIER A. MURILLO M.

    En secciones anteriores vimos como graficar los componentes de frecuencia de la


    seal anterior:

    Veamos ahora la forma de obtener en Matlab otra grafica muy importante como es
    la del espectro de potencias.

    Espectro de potencia para una seal senoidal:


    Cdigo en Matlab
    % Calculo de la FFT de una onda senoidal
    Fs = 150; % frecuencia de sampleo (en Hz)
    t = 0:1/Fs:0.1; % vector tiempo de 1 segundo
    f = 10; % frecuencia de la onda senoidal
    x=sin(2*pi*t*f);%onda senoidal
    nfft=1024;%el numero de puntos de la fft
    Y=fft(x,nfft);% tomar la FFT, y llenando con ceros, de manera que el
    largo de la FFT sea nfft
    Y = Y(1:nfft/2); % la FFT es simtrica, as que se tira la mitad
    my = abs(Y).^2;% tomar la potencia espectral, mdulo alcuadrado de la FFT
    f = (0:nfft/2-1)*Fs/nfft; %construccion del vector de frecuencias
    % Genera los plots, titulos y nombres.
    figure(1);
    plot(t,x);
    title('seal');
    xlabel('Tiempo (s)');
    ylabel('Amplitud ')
    figure(2);
    plot(f,my);
    title('Espectro de potencia');
    xlabel('Frecuencia (Hz)');
    ylabel('Potencia');

    42
    UNIVERSIDAD SANTIAGO DE CALI
    FACULTAD INGENIERIAS
    PROGRAMAS DE INGENIERIA ELECTRONICA Y BIOINGENIERIA
    SEALES Y SISTEMAS
    Profesor: MSc. JAVIER A. MURILLO M.

    Espectro de potencia para una seal cuadrada:


    Cdigo en Matlab

    % Grafica del espectro de potencia para una onda cuadrada


    A=1;
    w=10*pi;
    rho=0.5;
    Fs = 150; % frecuencia de sampleo (en Hz)
    t = 0:1/Fs:1; % vector tiempo de 1 segundo
    f = 10; % frecuencia de la onda cuadrada
    x=square(w*t+rho);%onda cuadrada
    nfft=1024;%el numero de puntos de la fft
    Y=fft(x,nfft);% tomar la FFT, y llenando con ceros, de manera que el
    largo de la FFT sea nfft
    Y = Y(1:nfft/2); % la FFT es simtrica, as que se tira la mitad
    my = abs(Y).^2;% tomar la potencia espectral, mdulo alcuadrado de la FFT
    f = (0:nfft/2-1)*Fs/nfft; %construccion del vector de frecuencias
    % Genera los plots, titulos y nombres.
    figure(1);
    plot(t,x);
    title('seal');
    xlabel('Tiempo (s)');
    ylabel('Amplitud ')
    figure(2);
    plot(f,my);
    title('Espectro de potencia');
    xlabel('Frecuencia (Hz)');
    ylabel('Potencia');

    43
    UNIVERSIDAD SANTIAGO DE CALI
    FACULTAD INGENIERIAS
    PROGRAMAS DE INGENIERIA ELECTRONICA Y BIOINGENIERIA
    SEALES Y SISTEMAS
    Profesor: MSc. JAVIER A. MURILLO M.

    Espectro de potencia para una seal triangular:


    Cdigo en Matlab

    % Grafica del espectro de potencia para una onda triangular


    A=1;
    w=10*pi;
    rho=0.5;
    Wdt=0.2;
    Fs = 150; % frecuencia de sampleo (en Hz)
    t = 0:1/Fs:1; % vector tiempo de 1 segundo
    f = 10; % frecuencia de la onda triangular
    x=A * sawtooth(w * t + Wdt);%onda triangular
    nfft=1024;%el numero de puntos de la fft
    Y=fft(x,nfft);% tomar la FFT, y llenando con ceros, de manera que el
    largo de la FFT sea nfft
    Y = Y(1:nfft/2); % la FFT es simtrica, as que se tira la mitad
    my = abs(Y).^2;% tomar la potencia espectral, mdulo alcuadrado de la FFT
    f = (0:nfft/2-1)*Fs/nfft; %construccion del vector de frecuencias
    % Genera los plots, titulos y nombres.
    figure(1);
    plot(t,x);
    title('seal');
    xlabel('Tiempo (s)');
    ylabel('Amplitud ')
    figure(2);
    plot(f,my);
    title('Espectro de potencia');
    xlabel('Frecuencia (Hz)');
    ylabel('Potencia');

    44
    UNIVERSIDAD SANTIAGO DE CALI
    FACULTAD INGENIERIAS
    PROGRAMAS DE INGENIERIA ELECTRONICA Y BIOINGENIERIA
    SEALES Y SISTEMAS
    Profesor: MSc. JAVIER A. MURILLO M.

    Precisamente la representacin espectral de un proceso es una herramienta


    especialmente adecuada para sacar a la luz componentes peridicas, ya que
    consiste en expresar el proceso como una suma de sucesiones peridicas,
    concretamente trigonomtricas. Esta representacin, por tanto, nos permite
    centrar la atencin en las frecuencias que ms nos interesen.

    Una serie temporal puede entenderse como una sucesin de nmeros. Las
    sucesiones de nmeros reales se pueden escribir como combinaciones lineales de
    senos y cosenos (o exponenciales complejas). Esta representacin se hace
    mediante la Transformada de Fourier Discreta (TFD).

    45

    También podría gustarte