Matlab - Frecuencia Onda Cuadrada

clc

clear all

close all

A = 1; % AMPLITUD DE LA ONDA CUADRADA

f = 17; % FRECUENCIA DE LA ONDA CUADRADA

w = 2*pi*f; % VELOCIDAD ANGULAR DE LA ONDA CUADRADA

D = 0.5; % CICLO UTIL DE LA ONDA CUADRADA

Fs = 150; % FRECUENCIA DE MUESTREO [Hz]

t = 0:1/Fs:1; % VECTOR DE TIEMPO

sq = A*square(w*t+D); % GENERACIÓN DE LA ONDA CUADRADA

figure(1)

plot(t,sq)

nfft = 1024; % NÚMERO DE PUNTOS DE LA FFT

Y = fft(sq,nfft); % TOMAR LA FFT, Y LLENANDO CON CEROS, DE MANERA QUE EL LARGO DE LA FFT SEA nfft

Y = Y(1:nfft/2); % POR SER SIMÉTRICA LA FFT, SE PUEDE PRESCINDIR DE LA MITAD

my = abs(Y).^2; % POTENCIA ESPECTRAL. MÓDULO DE LA FFT AL CUADRADO

Vf = (0:nfft/2-1)*Fs/nfft; % CONSTRUCCIÓN DEL VECTOR DE FRECUENCIAS

figure(2)

plot(Vf,my)

[Y,I] = max(my)

frecuencia = Vf(I)