List of Simulation Experiments For Analog Modulation Techniques
This document contains Matlab programs and simulations for analog modulation techniques. It includes:
1) A DS-B modulation simulation that generates a modulated signal, analyzes its spectrum, and calculates the power of the modulated and message signals.
2) A simulation of USB modulation that generates a modulated signal, calculates the power content, and analyzes the spectrum of the carrier, unmodulated, and modulated signals.
3) Additional simulations generate modulated signals, analyze the power in the demodulated signal as a function of phase, and generate Hilbert transforms of message signals.
List of Simulation Experiments For Analog Modulation Techniques
This document contains Matlab programs and simulations for analog modulation techniques. It includes:
1) A DS-B modulation simulation that generates a modulated signal, analyzes its spectrum, and calculates the power of the modulated and message signals.
2) A simulation of USB modulation that generates a modulated signal, calculates the power content, and analyzes the spectrum of the carrier, unmodulated, and modulated signals.
3) Additional simulations generate modulated signals, analyze the power in the demodulated signal as a function of phase, and generate Hilbert transforms of message signals.
1) Matlab programme:- clc; clear all; close all; ts=0.001; %sampling time df=0.3; %desired frequeny resolution fs=1/ts; %sampling frequency t=0:0.001:2; %time interval x=zeros(1,2001); for i=1:2001 %defining msg signal
if(t(i)>=0.1 && t(i)<1)
x(i)=t(i);
elseif(t(i)>=1 && t(i)<2)
x(i)=-t(i)+2; else x(i)=0.1;
end;
end;
1
u=x.*sin(2*pi*25*t); %dsb modulated signal subplot(221); plot(t,u); xlabel('time');ylabel('modulated signal'); title('modulated signal'); M=fft(u); % fourier transform of modulated signal f=[0:1/2:(length(M)-1)/2]-fs/2; % frequency vector subplot(222); plot(f,fftshift(abs(M))); xlabel('frequency'); title('spectrum of modulated signal'); hp = spectrum.periodogram('hamming'); hpopts = psdopts(hp,u); hpsd1 = psd(hp,u,hpopts);% calculate psd of modulated signal subplot(223); plot(hpsd1) power_modulated = avgpower(hpsd1); % calculate avg power of modulated signal display(power_modulated); hpsd2 = psd(hp,x,hpopts);% calculate psd of msg signal subplot(224); plot(hpsd2) power_msg = avgpower(hpsd1);% power of msg signal display(power_msg);
power_modulated = 1
0.3000 Q2) Matlab program:- clc; clear all; close all; ts=0.001;%sampling time df=0.3;%desired frequeny resolution fs=1/ts;%sampling frequency t=0:0.001:2;%time interval x=zeros(1,2001); for i=1:2001
if(t(i)>=0 && t(i)<1)
x(i)=t(i);
elseif(t(i)>=1 && t(i)<=2)
x(i)=-t(i)+2;
end;
end;
u=x.*sin(2*pi*25*t); %dsb modulated signal subplot(221); plot(t,u); xlabel('time');ylabel('modulated signal'); title('modulated signal'); M=fft(u); % fourier transform of modulated signal f=[0:1/2:(length(M)-1)/2]-fs/2; % frequency vector subplot(222); plot(f,fftshift(abs(M))); xlabel('frequency'); title('spectrum of modulated signal'); hp = spectrum.periodogram('hamming'); hpopts = psdopts(hp,u); hpsd1 = psd(hp,u,hpopts);% calculate psd of modulated signal subplot(223); plot(hpsd1) power_modulated = avgpower(hpsd1); % calculate avg power of modulated signal display(power_modulated); hpsd2 = psd(hp,x,hpopts);% calculate psd of msg signal subplot(224); plot(hpsd2) power_msg = avgpower(hpsd1);% power of msg signal display(power_msg); 1
a)modulated signal
syms t t1=(sinc(10*t)*sinc(10*t))^2;pw1=int(t1,t,-2,1); pm=(1/p)*(pw1); display('Q_b)power content of the modulated signal') pow_u_mod=((ac^2)/2)*pm %b)power content of the modulated signal
Warning: Explicit integral could not be found. Q_b)power content of the modulated signal
pow_u_mod =
int(sin(10*pi*t)^4/(10000*pi^4*t^4), t == -2..1)/8
syms t t1=(sinc(10*t)*sinc(10*t))^2;pw1=int(t1,t,-2,1); pm=(1/p)*(pw1); 3
display('Q_b)power content of the modulated signal') pow_u_mod=((ac^2)/2)*pm %b)power content of the modulated signal
Warning: Explicit integral could not be found. Q_b)power content of the modulated signal
pow_u_mod =
int(sin(10*pi*t)^4/(10000*pi^4*t^4), t == -2..1)/8
c)spectrum of the carrier and modulated
Fs = 1000; tc = 0:1/Fs:1-(1/Fs); % x = 14*cos(200*pi*t-pi/3);% x=2*cos(2*pi*100*t); c=ac*cos(2*pi*fc*tc); xdft = (1/length(c))*fft(c); freq = -500:(Fs/length(c)):500-(Fs/length(c)); figure,plot(freq,abs(fftshift(xdft))); %spectrum of the carrier signal title('spectrum of the carrier signal') le=length(tc); m=zeros(1,le); for i=1:le if (tc(i)>=-2 && tc(i)<=2) m(i)=m(i)+(sinc(10*tc(i))*sinc(10*tc(i))); else
end end u_mod=m.*c; m(i)=m(i)+0; xdft = (1/length(u_mod))*fft(u_mod); freq = -500:(Fs/length(u_mod)):500-(Fs/length(u_mod)); figure,plot(freq,abs(fftshift(xdft))); title('spectrum of the modulated signal') %c)spectrum of the modulated s
4
d)Power spectral density PSD carrier,modulated and message
Hs=spectrum.cov; psd(Hs,c,'Fs',Fs); %% PSD of carrier figure,psd(Hs,u_mod,'Fs',Fs);hold on; %% PSD of modulated signal figure,psd(Hs,m,'Fs',Fs); %% PSD of msg
% peak amplitude is at f=100 Hz
5 1
1
3.5
a)USB Modulated
clc;clear all;close all; t=-50:0.1:50;p=2;ac=1;fc=25; le=length(t); m=zeros(1,le); for i=1:le if (t(i)>=0 && t(i)<0.1) m(i)=m(i)+0.1; else if (t(i)>=0.1 && t(i)<1) m(i)=m(i)+t(i); else if (t(i)>=1 && t(i)<1.9) m(i)=m(i)-t(i)+2; else m(i)=m(i)+0.1; end end end end m; plot(t,m,'r');axis([0 2 0 1.5]);title('Given i/p sig.'); % t1 = (0:1/1023:1); y = hilbert(m); plot(t,real(y));axis([-5 5 -3 3]);grid on; hold on plot(t,imag(y),':');axis([-5 5 -3 3]) % axis([0 0.05 -1.1 2]); legend('Real-Part-i/p','Imag-Part-hilb-trans','location','best'); im=imag(y); c=ac*cos(2*pi*fc*t);c1=ac*sin(2*pi*fc*t); %carrier u_mod=(ac/2)*((m.*c)-(im.*c1)); %a)USB Modulated sig figure,plot(t,u_mod);axis([-5 5 -3 3]);grid on;title('Modulated sig');
3
4
b)power content of the modulated
syms t t1=0.1;pw1=int(t1,t,0,0.1); t2=t^2;pw2=int(t2,t,0.1,1); t3=(-t+2)^2;pw3=int(t3,t,1,1.9); t4=0.1;pw4=int(t4,t,1.9,2); pm=(1/p)*(pw1+pw2+pw3+pw4); display('Q_b)power content of the USB modulated signal') pow_u_mod=((ac^2)/4)*pm %b)power content of the modulated
Q_b)power content of the USB modulated signal pow_u_mod = 343/4000
c)spectrum of the carrier ,unmodulated modulated
Fs = 1000; tc = 0:1/Fs:1-(1/Fs); % x = 14*cos(200*pi*t-pi/3);% x=2*cos(2*pi*100*t); c=ac*cos(2*pi*fc*tc); xdft = (1/length(c))*fft(c); freq = -500:(Fs/length(c)):500-(Fs/length(c)); figure,plot(freq,abs(fftshift(xdft)));grid on; %spectrum of the carri title('spectrum of the carrier signal') le=length(tc); m1=zeros(1,le); for i=1:le if (tc(i)>=0 && tc(i)<0.1) m1(i)=m1(i)+0.1; else if (tc(i)>=0.1 && tc(i)<1) m1(i)=m1(i)+tc(i); else if (tc(i)>=1 && tc(i)<1.9) m1(i)=m1(i)-tc(i)+2; else m1(i)=m1(i)+0.1; end end end end m1; xdft = (1/length(m1))*fft(m1); freq = -500:(Fs/length(m1)):500-(Fs/length(m1)); figure,plot(freq,abs(fftshift(xdft)));grid on; title('spectrum of the unmodulated signal') %spectrum of the unmodulat y1 = hilbert(m1); 5
im=imag(y1); c=ac*cos(2*pi*fc*tc);c1=ac*sin(2*pi*fc*tc); u_mod=(ac/2)*((m1.*c)-(im.*c1)); xdft = (1/length(u_mod))*fft(u_mod); freq = -500:(Fs/length(u_mod)):500-(Fs/length(u_mod)); figure,plot(freq,abs(fftshift(xdft)));grid on; title('spectrum of the modulated signal') %c)spectrum of the modulated s
5
d)Power spectral density PSD carrier,modulated and message
Hs=spectrum.cov; psd(Hs,c,'Fs',Fs); %% PSD of carrier figure,psd(Hs,u_mod,'Fs',Fs); %% PSD of modulated signal figure,psd(Hs,m,'Fs',Fs); %% PSD of msg
% peak amplitude is at f=25hz
6
PROGRAM-3.8
a)Modulated signal
clc;clear all;close all; t=-50:0.1:50;p=2;ac=1;fc=25; le=leng th(t); m=zeros (1,le); for i=1:le if (t(i)>=0 && t(i)<0.1) m(i)=m(i)+0.1; else if (t(i)>=0.1 && t(i)<1) m(i)=m(i)+t(i); else if (t(i)>=1 && t(i)<1.9) m(i)=m(i)-t(i)+2; else m(i)=m(i)+0; end end end end plot(t,m,'r');axis([0 2 0 1.5]);grid on;title('Given i/p sig.');
syms t f t1=0.1;pw1=int(t1,t,0,0.1); t2=t^2;pw2=int(t2,t,0.1,1); t3=(-t+2)^2;pw3=int(t3,t,1,1.9); t4=0.1;pw4=int(t4,t,1.9,2); pm=(1/p)*(pw1+pw2+pw3+pw4); display('Q_b)power content of the demodulated signal') ac=ac^2;pf=cos(f); pow_u_mod=((ac*pf)^2)*pm %power content of the demodulated pow_u_mod(f)=pow_u_mod; p=0; f=0:0.01:pi/2; p=[p,pow_u_mod(f)]; p=[p(2:length(p))]; figure,plot(f,p,'r');grid on;title('PWR in demod sig VS f') %b)PWR in demod sig legend('0<=theta<=pi/2','location','Best')
q_b)power content of the demodulated signal pow_u_mod = (343*cos(f)^2)/1000 1
PROGRAM-3.9
a)Modulated
clc;clear all;close all; t=-50:0.1:50;p=2;ac=1;fc=25; le=length(t); m=zeros(1,le); for i=1:le if (t(i)>=0 && t(i)<0.1) m(i)=m(i)+0.1; else if (t(i)>=0.1 && t(i)<1) m(i)=m(i)+t(i); else if (t(i)>=1 && t(i)<1.9) m(i)=m(i)-t(i)+2; else m(i)=m(i)+0.1; end end end end plot(t,m,'r');axis([0 2 0 1.5]); 2
y = hilbert(m); plot(t,real(y));axis([-5 5 -3 3]); hold on plot(t,imag(y),':');axis([-5 5 -3 3]);grid on;title('Given i/p sig.'); legend('Real-Part-i/p','Imag-Part-hilb-trans','location','best'); im=imag(y); max_amp=max(im); c=ac*cos(2*pi*fc*t);c1=ac*sin(2*pi*fc*t); %carrier u_mod=(ac/2)*((m.*c)-(im.*c1)); %a)USB Modulated sig figure,plot(t,u_mod);axis([-5 5 -3 3]);grid on;title('Modulated sig');
syms t f t1=0.1;pw1=int(t1,t,0,0.1); t2=t^2;pw2=int(t2,t,0.1,1); t3=(-t+2)^2;pw3=int(t3,t,1,1.9); t4=0.1;pw4=int(t4,t,1.9,2); pm=(1/p)*(pw1+pw2+pw3+pw4); %power content of the msg signal display('power content of the demodulated signal') pow_u_mod=((ac^2)/4)*pm ; pf=cos(f); pow_u_mod1=((pf)^2)*pow_u_mod; pow_u_mod1(f)=pow_u_mod1 pwh=int((max_amp)^2,t,0,2); p_hilb=(1/p)*pwh ; %power content of the hilbert trans. display('power content of the hilbert trans. of msg signal') pow_u_hmod=((ac^2)/4)*p_hilb; pf1=sin(f); pow_u_mod2=((pf1)^2)*pow_u_hmod; pow_u_mod2(f)=pow_u_mod2 p_ussb=pow_u_mod1(f)-pow_u_mod2(f);
p_ussb(f)=p_ussb
power content of the demodulated signal
pow_u_mod1(f) = (343*cos(f)^2)/4000 power content of the hilbert trans. of msg signal pow_u_mod2(f) = (4586806831851157*sin(f)^2)/72057594037927936