%Parameters
FM = 1000;
AmpFM = 2.5;
Fc = 10000;
AmpFc = 2.5;
Kf = 80;
SampRate = 1000000;
Time = -1/FM:1/SampRate:1/FM;
betaValues = [1 2 5 10 20];
DeltaW = Kf*AmpFM;
%Signals
SineSignal = sin(2*pi*FM*Time);

for i = 1:5
    Wm = DeltaW/betaValues(i);
    beta = DeltaW/Wm;
    OutFM = AmpFM*cos(2*pi*Fc*Time+beta*sin(2*pi*FM*Time));
    OutFM_S = fftshift(fft(OutFM));
    figure(i);
    plot(Time,abs(OutFM_S));
    title(['Constant DeltaW / Beta = ' num2str(betaValues(i))]);

end
clear;