Hola Beto:
Yo también estoy haciendo cosillas en este plan con el Matlab, y he conseguido generar la señal modulada en FM. Básicamente, lo que he hecho es partir de la ecuación general de la modulación en frecuencia:
y(t) = Ac*cos(wc*t + Kf * Integral( xn(t)dt ) )
Teniendo en cuenta que el ancho de banda aproximado que va a ocupar tu señal modulada en FM (de acuerdo con la regla de Carson) es B_FM = 2(fd + fm), siendo fd la máxima desviación en frecuencia y fm el ancho de banda de la señal moduladora x(t), deberás definir tu frecuencia de muestreo de manera que sea por lo menos el doble (y yo te recomiendo que sea apreciablemente mayor que ese mínimo de Nyquist) que la máxima frecuencia a representar.
Asimismo, para implementar la ecuación anterior necesitas desarrollar la integral. No te aconsejo uisar las funciones de integración de Matlab, sino más bien aproximar la integral por una acumulación. Para ello he programado el bloque acumula.m, cuyo código es el siguiente:
function [sum] = acumula(x,t)
N=length(x);
acum=zeros(1,N);
delta_t=t(2)-t(1);
for i=1:N
if i==1
acum(i)=x(i)*delta_t;
else
acum(i)=acum(i-1)+delta_t*x(i);
end
end
sum=acum;
offset=0.5*(max(sum)+min(sum));
sum=sum-offset;
Posteriormente he definido el bloque llamado Modula_FM.m, que es:
function [y] = Modula_FM(x,t,Kf,wp,Amp)
xac=acumula(x,t);
y = Amp*cos(wp*t+2*pi*Kf*xac);
donde Kf es la constante del modulador (en electrónica, se mide en Hz / V), wp la pulsación de la portadora (que no la frecuencia, por lo que deberás introducir 2*pi*fp si define tu frecuencia de portadora en Hz) y Amp la amplitud de la portadora.
Habrás tenido que definir previamente los siguientes parámetros: frecuencia de muestreo y ventana de tiempos. A partir de ellos, defines los ejes de tiempos y frecuencias. Posteriormente creas tu señal moduladora x(t), lo normal es que sea una sinusoide de frecuencia mucho menor que la de la portadora.
Si tienes alguna duda más, aquí me tienes. Por cierto, si además te interesa visualizar el espectro de la señal modulada en FM, he desarrollado otro fichero .m llamado Analiza_Espectro.m, que, dada una señal en el tiempo y un eje de frecuencias, muestra el espectro de dicha señal de una manera similar a como aparece en un analizador de espectros de radiocomunicaciones. El código es:
function [output_args] = Analiza_Espectro(x,f)
X = fftshift(fft(x));
Xm = abs(X);
plot(f,Xm)
Un saludo.
LuisDa.
