例えば、100Hzの搬送波を10Hzの正弦波で変調する場合、以下の様に書けると思います。
%%
fs = 1000;
t = 0 :1/fs : 1;
y1 = sin(2*pi*10*t); %ベースバンド信号
y2 = sin(2*pi*100*t); %搬送波（中間周波数、IF）
y = y1.*y2; % mixing
Y = fft(y);
f = 0 : fs/length(Y) : fs*(1-1/length(Y));
plot(f,abs(Y)), xlim([0,fs/2]), grid
%%
この例の場合、100-10=90[Hz]と、100+10=[110]Hzが、ビート周波数に相当します。
Signal Processing Toolboxや、Communications Toolboxをお持ちであれば、以下も参考になると思います。
