%EE7780 Array Signal Processing
%Tarea 1 semana 1
function [st]=fCreateNBSig(fs,Tlength,fAt,fbin)

%Vamos a crear Tlength segundos de señal
Ts = 1/fs; %periodo de muestreo

t = 0:Ts:Tlength-Ts; %vector de tiempos, Tlength segundos Notese: se elimina Ts del 
                %final para que de un número redondo, en el caso de 2 segundos 
                %debe dar la longitud del vector t 88mil muestras

%A(t) que es una señal con variación lenta comparada a la frecuencia donde está montada                
omegaAt = 2*pi*fAt;     %frecuencia angular           
At = 1+ sin(omegaAt*t); %amplitud de la señal de interes
figure
ax(1) = subplot(3,1,1)
plot(t,At,'Parent',ax(1)), grid on
ylabel('A(t) slow mov. signal')

%Tono donde está montada A(t)
omegabin = 2*pi*fbin; %frecuencia angular
expomegat = exp(j*omegabin*t); %exponencial compleja donde está el tono
ax(2) = subplot(3,1,2)
plot(t,real(expomegat),'Parent',ax(2)), grid on
ylabel('exp(j omega t)')

%Montar A(t) con la exponencial para generar s(t)
st =At.*expomegat; %SOI
ax(3) = subplot(3,1,3)
plot(t,real(st),'Parent',ax(3)), grid on
xlabel('time, sec')
linkaxes(ax,'x'); %encadena los ejes x para hacer zoom juntos
linkaxes(ax,'y'); %encadena los ejes y para comparar figuras
axis([0 Tlength -4 4]) %Fijar ejes en la figura
ylabel('s(t), SOI')

disp('>>>>>>>>>>>>>>>>>>>>>>>>>end')
