Но как оно на практике, до меня не доходит.

Короче если кому интересно.

clear all;

for i=1:8192;
nf(i) = sin(i/6*pi);
mf(i) = sin(i/100*pi);

ssb(i) = sin(i/6*pi).*sin(i/100*pi) - cos(i/6*pi).*cos(i/100*pi);
end

f = nf.*mf;
%ssb = ssbmod(mf,150,1000);
%ssb =

subplot(4,1,1);
spectrogram(nf,2048,2000);
subplot(4,1,2);
spectrogram(mf,2048,2000);
subplot(4,1,3);
spectrogram(nf.*mf,2048,2000);
subplot(4,1,4);
spectrogram(ssb,2048,2000);

----------------------------------------------------

clear all;
%x = 0:0.5:1;
%f=round(sin(x));
%x(0)= 0;
for i=1:1000;
if(sin(i/300*pi)>0) x(i) = 1;
else x(i) = -1.0;
end
end
fc = ssbmod(x,40,1000);
fr = ssbdemod(fc,40,1000);
subplot(4,1,1);
plot(x);
subplot(4,1,2);
plot(fc);
subplot(4,1,3);
plot(fr);
subplot(4,1,4);
spectrogram(fc,32,31,'yaxis');