Fm = 1000;
Am = 2.5;
Fc = 10000;
Ac = 2.5;
s = 1000000;
t= -1/Fm:1/s:1/Fm;

betas = [0.01 1 2.4 10 50];
%Sine Signal Input
signal = Ac*cos(2*pi*Fc*t);

%Modulating Signal
modulating = Am*cos(2*pi*Fm*t);


for i=1:5
   fprintf('\n');
   modulated_signal = Am*cos(2*pi*Fc*t+betas(i)*sin(2*pi*Fm*t));
   fast_fourier = fftshift(fft(modulated_signal));
   len = length(modulated_signal);
   freqs = -s/2:s/len:s/2-s/len;
   delta = obw(ab(fast_fourier)) * betas(i);
   kf = delta/ Am;
   fprintf('The delta w for B: %d is equal to %d\n', betas(i), delta);
   fprintf('The Scaling factor for B: %d is equal to %d\n', betas(i), kf);
%    figure(i), plot(t,modulated_signal,'g', t,fast_fourier, 'r'), title(['Angle Modulation at Beta : ', num2str(betas(i))])
   figure (i), plot(freqs,abs(fast_fourier)), title(['FFT generated at Beta : ', num2str(betas(i))])
end

clear;
clc;