C Sai Bhavana CSlab MATLAB
C Sai Bhavana CSlab MATLAB
WARANGAL
CS lab report - I
A modulated signal is obtained by multiplying message signal and its carrier wave. Two
modulated signals are added together to form FDM signal.
MATLAB code:
clear all
close all
clc
fs=1000;%sampling frequency
t=0:1/fs:1;
a1=2;%amplitude of mesaage signal 1
f1=4;%frequency of message signal 1
m1=a1*cos(2*pi*f1*t);%message signal 1
a2=2;%amplitude of mesaage signal 2
f2=5;%frequency of message signal 2
m2=a2*cos(2*pi*f2*t);%message signal 2
%carrier signal generation
ac1=2;%amplitude of carrier signal 1
fc1=50;%frequency of carrier signal 1
c1=ac1*cos(2*pi*fc1*t);%carrier signal 1
ac2=2;%amplitude of carrier signal 2
fc2=150;%frequency ofcarrier signal 2
c2=ac2*cos(2*pi*fc2*t);%carrier signal 2
%composite signal
x=m1.*c1+m2.*c2;
figure(1);
plot(t,x);
xlabel('Time');
ylabel('Composite signal');
%Absolute frequency on X-axis vs.Magniitude on Y axis
N=length(x);
X=fftshift(fft(x,N));
f=fs*[-N/2:N/2-1]/N;
figure(2);
plot(f,abs(X));
xlabel('Absolute frequency');
ylabel('DFT values');
title('Composite signal spectrum');
%Recovering message signal 1
%Bandpass filter
[num den]=butter(5,[(fc1-f1)*2/fs (fc1+f1)*4/fs]);
bpf_1=filter(num,den,x);
%mixing
z1=2*bpf_1.*c1;
[num den]=butter(5,f1*4/fs);
rec1=filter(num,den,z1);
figure(3);
plot(t,rec1);
xlabel('Time');
ylabel('Recovered signal 1');
%Absolute frequency on X-axis vs.Magniitude on Y axis
N=length(rec1);
r1=fftshift(fft(rec1,N));
f=fs*[-N/2:N/2-1]/N;
figure(4);
plot(f,abs(r1));
xlabel('Absolute frequency');
ylabel('DFT values');
title('Spectrum of Recovered signal 1');