www.pudn.com > digitalmodulate.zip > bfsk.m, change:2006-09-21,size:724b


clear; 
clc; 
b = input('Enter the Bit stream \n '); 
%b = [0 1 0 1 1 1 0]; 
n = length(b); 
t = 0:.01:n; 
x = 1:1:(n+1)*100; 
 
for i = 1:n 
    if (b(i) == 0) 
        b_p(i) = -1; 
    else 
        b_p(i) = 1; 
    end 
    for j = i:.1:i+1 
        bw(x(i*100:(i+1)*100)) = b_p(i);         
    end 
end 
 
bw = bw(100:end); 
 
wo = 2*(2*pi*t); 
W  = 1*(2*pi*t); 
sinHt = sin(wo+W); 
sinLt = sin(wo-W); 
 
st = sin(wo+(bw).*W); 
 
 
subplot(4,1,1) 
plot(t,bw) 
grid on ; axis([0 n -2 +2]) 
subplot(4,1,2) 
plot(t,sinHt) 
grid on ; axis([0 n -2 +2]) 
subplot(4,1,3) 
plot(t,sinLt) 
grid on ; axis([0 n -2 +2]) 
subplot(4,1,4) 
plot(t,st) 
grid on ; axis([0 n -2 +2]) 
 
Fs=1; 
figure 
%pburg(st,10) 
periodogram(st)