Code covered by the BSD License  

Highlights from
Contemporary Communications Systems Matlab Files

image thumbnail

Contemporary Communications Systems Matlab Files

by

Omar Ruiz (view profile)

 

Matlab Files in this book

dsb3.m
% dsb1.m
% Matlab demonstration script for DSB-AM modulation. The message signal
% is +1 for 0 < t < t0/3, -2 for t0/3 < t < 2t0/3 and zero otherwise.
echo on
t0=.15;                              	% signal duration
ts=0.0005;                           	% sampling interval
fc=250;                              	% carrier frequency
snr=10;                              	% SNR in dB (logarithmic)
fs=1/ts;                             	% sampling frequency
t_long=20*t0;                        	% extended signal duration 
n0=floor(t_long/ts)+1;               	% length of extended signal vector 
n00=2^(ceil(log2(n0)));
df=fs/n00;                            	% frequency resolution
t_long1=[0:ts:t_long];               	% extended time vector
f=[-fs/2+df:df:fs/2];                	% frequency vector 
% message signal
m=[ones(1,t0/(3*ts)),-2*ones(1,t0/(3*ts)),zeros(1,t0/(3*ts)+1)];
m_long=[m,zeros(1,n0-length(m))];    	% extended message signal 
c_long=cos(2*pi*fc.*t_long1);        	% extended carrier signal 
u_long=m_long.*c_long;               	% extended modulated signal
du_long=u_long.*c_long;
DU=fft(du_long,n00)/fs;
M=fft(m_long,n00)/fs;                 	% spectrum of the extended message signal 
U=fft(u_long,n00)/fs;                  	% spectrum of the extended modulated signal 
pause   % Press any key to see a plots of the magnitude of the message and the
	% modulated signal in the frequency domain.
subplot(2,1,1)
plot(f,abs(fftshift(M)))
xlabel('Frequency')
title('Spectrum of the message signal')
subplot(2,1,2)
plot(f,abs(fftshift(U)))
title('Spectrum of the modulated signal')
xlabel('Frequency') 
pause  % Press a key to see the modulated and modulated X carrier in freq. domain
subplot(2,1,1)
plot(f,abs(fftshift(U)))
title('Modulated signal spectrum')
xlabel('Frequency')
subplot(2,1,2)
plot(f,abs(fftshift(DU))) 
title('Modulated signal X Carrier spectrum')
xlabel('Frequency')

Contact us