msk.zip_Creating_MATLAB msk_msk

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+2)*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); if (mod(i,2) == 0) bow(x(i*100:(i+1)*100)) = b_p(i); bow(x((i+1)*100:(i+2)*100)) = b_p(i); else bew(x(i*100:(i+1)*100)) = b_p(i); bew(x((i+1)*100:(i+2)*100)) = b_p(i); end if (mod(n,2)~= 0) bow(x(n*100:(n+1)*100)) = -1; bow(x((n+1)*100:(n+2)*100)) = -1; end end end bw = bw(100:end); bew = bew(100:(n+1)*100); bow = bow(200:(n+2)*100); wot = 2*pi*t*(5/4); Wt = 2*pi*t/(4*1); st = bow.*sin(wot+(bew.*bow).*Wt); subplot(4,1,1) plot(t,bw) grid on ; axis([0 n -2 +2]) subplot(4,1,2) plot(t,bow) grid on ; axis([0 n -2 +2]) subplot(4,1,3) plot(t,bew) grid on ; axis([0 n -2 +2]) subplot(4,1,4) plot(t,st) grid on ; axis([0 n -2 +2]) Fs=5/4; %figure %periodogram(st) S = fft(st,65); PSS = S.* conj(S) / 65; f = 1000*(-16:16)/65; figure% plot(f,PSS(1:33))