Frequency modulation by MATLAB

I hope to write a program in MATLAB to plot the output of the following function in time domain

s(t) = A sum ( J(B) cos((Wc-n.Wm).t )

where

• A is a constant = 1
• J Bessel function of B
• B = 10 (modulation index)
• Wc = 1000000 carrier frequency
• Wm = 1000
• Limits of summation n = -20 to 20

Hint: all values above changeable

Answer 1

[Vieniava]

From the top of my head:

Fc=1e6;  % carrier 1MHz
Fm=1e3;  % single tone - 1kHz
fs=2.2*Fc; % sampling frequency (Nyquist limitation!)
t=0:1/fs:1;
Y=zeros(1,length(t));
for n=-20:20
    Y=Y+besselj(n, B).*cos( 2*pi*(Fc-n*Fm).*t);
end

Should work (right now I can't vertify this code - lack of Matlab)